Metaheuristics for Search Problems in Genomics - New Algorithms and Applications

In this thesis we made the first steps towards the systematic application of a methodology for automatically building formal models of complex biological systems. Such a methodology could be useful also to design artificial systems possessing desirable properties such as robustness and evolvability. The approach we follow in this thesis is to manipulate formal models by means of adaptive search methods called metaheuristics. In the first part of the thesis we develop state-of-the-art hybrid metaheuristic algorithms to tackle two important problems in genomics, namely, the Haplotype Inference by parsimony and the Founder Sequence Reconstruction Problem. We compare our algorithms with other effective techniques in the literature, we show strength and limitations of our approaches to various problem formulations and, finally, we propose further enhancements that could possibly improve the performance of our algorithms and widen their applicability. In the second part, we concentrate on Boolean network (BN) models of gene regulatory networks (GRNs). We detail our automatic design methodology and apply it to four use cases which correspond to different design criteria and address some limitations of GRN modeling by BNs. Finally, we tackle the Density Classification Problem with the aim of showing the learning capabilities of BNs. Experimental evaluation of this methodology shows its efficacy in producing network that meet our design criteria. Our results, coherently to what has been found in other works, also suggest that networks manipulated by a search process exhibit a mixture of characteristics typical of different dynamical regimes.

Mechanismen der Antigenpräsentation bei der Induktion der CD8-T-Zellantwort gegen das murine Cytomegalovirus

Die Kontrolle der Cytomegalovirus(CMV)-Infektion durch CD8 T-Zellen ist abhängig von der effizienten MHC-Klasse-I-Präsentation viraler Peptide auf der Zelloberfläche. Um die Erkennung infizierter Zellen zu unterdrücken, interferieren während der Early (E)-Phase der murinen CMV (mCMV)-Infektion virale Immunevasine mit dem intrazellulären Transport von Peptid-MHC-I (pMHC-I) Komplexen. Den Immunevasinen gelingt es allerdings nicht, ein Priming mCMV-spezifischer CD8 T-Zellen zu verhindern. Daher wurde angenommen, dass die Initiation der antiviralen CD8 T-Zellantwort primär auf der Cross-Präsentation viraler Peptide auf nicht-infizierten, professionellen Antigen-präsentierenden Zellen (profAPC) beruht und damit unabhängig von viralen Immunevasionsmechanismen ist.rnIm Rahmen der vorliegenden Arbeit wurde mittels BAC-Mutagenese eine mCMV-Rekombinante generiert, um die direkte Präsentation viraler Peptide durch die zusätzliche Expression des zentralen Immunevasins m152 bereits in der Immediate Early (IE)-Phase verstärkt zu unterdrücken. Wie erwartet reduzierte die verstärkte m152-Expression sowohl in der IE- als auch in der E-Phase die pMHC-I-Präsentation in vitro. Dies führte überraschenderweise nach Infektion immunkompetenter BALB/c-Mäuse (Haplotyp H-2d) zu einer verminderten CD8 T-Zellantwort und damit zur Verschlechterung der Kontrolle der Infektion im drainierenden Lymphknoten. Diese Beobachtungen weisen erstmals auf einen wichtigen Beitrag der direkten Antigenpräsentation bei der Initiation der mCMV-spezifischen CD8 T-Zellantwort im immunkompetenten Wirt hin. Zusätzlich konnte auch nach mCMV-Infektion von Cross-Präsentations-defizienten Mäusen (Haplotyp H-2b) eine antivirale CD8 T-Zellantwort initiiert werden. Diese Beobachtung bestätigt, dass durch direkte Antigenpräsentation auf infizierten profAPC trotz viraler Immunevasionsmechanismen eine CD8 T-Zellantwort induziert werden kann. Allerdings wurde weder die antivirale CD8 T-Zellantwort noch die Kontrolle der Infektion im Haplotyp H-2b durch die verstärkte m152-Expression moduliert.rnIn einem weiteren Teil der Arbeit konnte im klinisch relevanten Modellsystem der mCMV-Infektion von Knochenmarktransplantations (KMT)-Rezipienten (Haplotyp H-2d) gezeigt werden, dass die verstärkte m152-Expression die Rekrutierung IE1-spezifischer CD8 T-Zellen in die infizierte Lunge unterdrückt. Dies konnte sowohl früh nach Infektion, als auch während der viralen Latenz nachgewiesen werden. Zusätzlich war die Rekrutierung IE1-spezifischer CD8 T-Zellen in die Lunge deutlich vermindert in Ld--Rezipienten von Ld+-hämatopoetischen Zellen, die das IE1-präsentierende MHC-I-Molekül Ld nicht auf den nicht-hämatopoetischen Gewebszellen exprimieren. Diese Beobachtungen zeigen, dass die Rekrutierung antiviraler CD8 T-Zellen in ein peripheres Organ von der direkten Antigenpräsentation auf nicht-hämatopoetischen, infizierten Gewebszellen bestimmt wird.rnIn der vorliegenden Arbeit konnte somit erstmals gezeigt werden, dass trotz viraler Immunevasionsmechanismen nach mCMV-Infektion des immunkompetenten Wirtes und des KMT-Rezipienten die antivirale CD8 T-Zellantwort von der direkten Antigenpräsentation bestimmt wird.

Molecular genetic investigation of the Neolithic population history in the western Carpathian Basin

Der Fokus dieser Dissertation ist die populationsgenetische Analyse der neolithischen Bevölkerungswechsel in den 6.-5. Jahrtausende vor Christus, die im westlichen Karpatenbecken stattfanden. Die Zielsetzung der Studie war, mittels der Analyse von mitochondrialer und Y-chromosomaler aDNA, den Genpool der sechs neolithischen und kupferzeitlichen Populationen zu untersuchen und die daraus resultierenden Ergebnisse mit anderen prähistorischen und modernen genetischen Daten zu vergleichen.rnInsgesamt wurden 323 Individuen aus 32 ungarischen, kroatischen und slowakischen Fundplätzen beprobt und bearbeitet in den archäogenetischen Laboren der Johannes Gutenberg-Universität in Mainz. Die DNA Ergebnisse wurden mit verschiedenen populationsgenetischen Methoden ausgewertet. Vergleichsdaten von prähistorischen und modernen eurasiatischen Populationen wurden dazu gesammelt.rnDie HVS-I Region der mitochondrialen DNA konnten bei 256 Individuen reproduziert und authentifiziert werden (mit einer Erfolgsrate von 85.9%). Die Typisierung der HVS-II Region war in 80 Fällen erfolgreich. Testend alle gut erhaltene Proben, die Y-chromosomale Haplogruppe konnte in 33 männlichen Individuen typisiert werden.rnDie neolithischen, mitochondrialen Haplogruppen deuten auf eine hohe Variabilität des maternalen Genpools hin. Sowohl die mitochondrialen als auch die Y-chromosomalen Daten lassen Rückschlüsse auf eine nah-östliche bzw. südwestasiatische Herkunft der frühen Bauern zu. Die Starčevo- und linearbandkermaischen-Populationen in westlichem Karpatenbecken (letztere abgekürzt als LBKT) und die linearbandkermaischen-Population in Mitteleuropa (LBK) haben so starke genetische Ähnlichkeit, dass die Verbreitung der LBK nach Mitteleuropa mit vorangegangenen Wanderungsereignissen zu erklären ist. Die Transdanubische aDNA Daten zeigen hohe Affinität zu den publizierten prähistorischen aDNA Datensätzen von Mitteleuropa aus den 6.-4. Jahrtausende vor Chr. Die maternal-genetische Variabilität der Starčevo-Population konnte auch innerhalb der nachfolgenden Populationen Transdanubiens festgestellt werden. Nur kleinere Infiltrationen und Immigrationsereignissen konnten während der Vinča-, LBKT-, Sopot- und Balaton-Lasinja-Kultur in Transdanubien identifiziert werden. Zwischen den transdanubischen Regionen konnten mögliche genetische Unterschiede nur in der LBKT und in der Lengyel-Periode beobachtet werden, als sich die nördlichen Gruppen von den südlichen Populationen trennten. rnDie festgestellte Heterogenität der mtDNA in Zusammenhang mit der Y-chromosomalen Homogenität in den Starčevo- und LBK-Populationen, weisen auf patrilokale Residenzregeln und patrilineare Abstammungsregeln in den ersten Bauergemeinschaften hin. rnObwohl die hier präsentierten Daten einen großen Fortschritt in der Forschung von aDNA und Neolithikum des Karpatenbeckens und Mitteleuropas bedeuten, werfen sie auch mehrere Fragen auf, deren Beantwortung durch zukünftige Genomforschungen erbracht werden könnte.

Phylogeny and phylogeography of the family Hyalidae (Crustacea: Amphipoda) along the northeast Atlantic coasts

The family Hyalidae comprises more than one hundred species, distributed worldwide. They are common and abundant in the littoral and shallow sublittoral habitats and they play an important role in the coastal food chain. Most studies about this family were dealing with taxonomy and ecology, while very little is known about phylogenetic relationship among genera and species. In the present study we aim to achieve the first approach of the phylogenetic patterns of this family in NE Atlantic Ocean and Mediterranean Sea, and to perform the first insight into the phylogeography Apohyale prevostii along both the North Atlantic coasts. In order to do that, eight species belonging to the genera Apohyale, Hyale, Serejohyale and Protohyale were investigated using the mitochondrial COI-5P barcode region. Specimens were collected along European and Moroccan Atlantic rocky shores, including Iceland, the British Isles, Macaronesia and in the Mediterranean Sea. Sequences of A. prevostii, from the NW Atlantic Ocean, available in BOLD and GenBank, were retrieved. As expected, phylogenetic analyses showed highly-divergent clades, clearly discriminating among different species clusters, confirming their morphology-based identifications. Although, within A. perieri, A. media, A. stebbingi, P. (Protohyale) schmidtii and S. spinidactylus, high genetic diversity was found, revealing putative cryptic species. The clade of A. prevostii and A. stebbingi appears well supported and divided from the other two congeneric species, and P. (Protohyale) schmidtii shows a basal divergence. The north-western Atlantic coasts were recently colonized by A. prevostii after the last glacial maximum from the European populations showing also a common haplotype in every population analysed. The use of the COI-5P as DNA barcode provided a good tool to underline the necessity of a revision of this emblematic family, as well as to discern taxonomically the possible new species flagged with this molecular device.

Genome-wide association study identifies new HLA class II haplotypes strongly protective against narcolepsy

Narcolepsy is a rare sleep disorder with the strongest human leukocyte antigen (HLA) association ever reported. Since the associated HLA-DRB1*1501-DQB1*0602 haplotype is common in the general population (15-25%), it has been suggested that it is almost necessary but not sufficient for developing narcolepsy. To further define the genetic basis of narcolepsy risk, we performed a genome-wide association study (GWAS) in 562 European individuals with narcolepsy (cases) and 702 ethnically matched controls, with independent replication in 370 cases and 495 controls, all heterozygous for DRB1*1501-DQB1*0602. We found association with a protective variant near HLA-DQA2 (rs2858884; P < 3 x 10(-8)). Further analysis revealed that rs2858884 is strongly linked to DRB1*03-DQB1*02 (P < 4 x 10(-43)) and DRB1*1301-DQB1*0603 (P < 3 x 10(-7)). Cases almost never carried a trans DRB1*1301-DQB1*0603 haplotype (odds ratio = 0.02; P < 6 x 10(-14)). This unexpected protective HLA haplotype suggests a virtually causal involvement of the HLA region in narcolepsy susceptibility.

Hybridization between distant lineages increases adaptive variation during a biological invasion: stickleback in Switzerland

The three-spined stickleback is a widespread Holarctic species complex that radiated from the sea into freshwaters after the retreat of the Pleistocene ice sheets. In Switzerland, sticklebacks were absent with the exception of the far northwest, but different introduced populations have expanded to occupy a wide range of habitats since the late 19th century. A well-studied adaptive phenotypic trait in sticklebacks is the number of lateral plates. With few exceptions, freshwater and marine populations in Europe are fixed for either the low plated phenotype or the fully plated phenotype, respectively. Switzerland, in contrast, harbours in close proximity the full range of phenotypic variation known from across the continent. We addressed the phylogeographic origins of Swiss sticklebacks using mitochondrial partial cytochrome b and control region sequences. We found only five different haplotypes but these originated from three distinct European regions, fixed for different plate phenotypes. These lineages occur largely in isolation at opposite ends of Switzerland, but co-occur in a large central part. Across the country, we found a strong correlation between a microsatellite linked to the high plate ectodysplasin allele and the mitochondrial haplotype from a region where the fully plated phenotype is fixed. Phylogenomic and population genomic analysis of 481 polymorphic amplified fragment length polymorphism loci indicate genetic admixture in the central part of the country. The same part of the country also carries elevated within-population phenotypic variation. We conclude that during the recent invasive range expansion of sticklebacks in Switzerland, adaptive and neutral between-population genetic variation was converted into within-population variation, raising the possibility that hybridization between colonizing lineages contributed to the ecological success of sticklebacks in Switzerland.

Estimating the net contribution of interleukin-28B variation to spontaneous hepatitis C virus clearance

The identification of associations between interleukin-28B (IL-28B) variants and the spontaneous clearance of hepatitis C virus (HCV) raises the issues of causality and the net contribution of host genetics to the trait. To estimate more precisely the net effect of IL-28B genetic variation on HCV clearance, we optimized genotyping and compared the host contributions in multiple- and single-source cohorts to control for viral and demographic effects. The analysis included individuals with chronic or spontaneously cleared HCV infections from a multiple-source cohort (n = 389) and a single-source cohort (n = 71). We performed detailed genotyping in the coding region of IL-28B and searched for copy number variations to identify the genetic variant or haplotype carrying the strongest association with viral clearance. This analysis was used to compare the effects of IL-28B variation in the two cohorts. Haplotypes characterized by carriage of the major alleles at IL-28B single-nucleotide polymorphisms (SNPs) were highly overrepresented in individuals with spontaneous clearance versus those with chronic HCV infections (66.1% versus 38.6%, P = 6 × 10(-9) ). The odds ratios for clearance were 2.1 [95% confidence interval (CI) = 1.6-3.0] and 3.9 (95% CI = 1.5-10.2) in the multiple- and single-source cohorts, respectively. Protective haplotypes were in perfect linkage (r(2) = 1.0) with a nonsynonymous coding variant (rs8103142). Copy number variants were not detected. CONCLUSION: We identified IL-28B haplotypes highly predictive of spontaneous HCV clearance. The high linkage disequilibrium between IL-28B SNPs indicates that association studies need to be complemented by functional experiments to identify single causal variants. The point estimate for the genetic effect was higher in the single-source cohort, which was used to effectively control for viral diversity, sex, and coinfections and, therefore, offered a precise estimate of the net host genetic contribution.

Common CFTR haplotypes and susceptibility to chronic pancreatitis and congenital bilateral absence of the vas deferens

CFTR mutations enhance susceptibility for idiopathic chronic pancreatitis (ICP) and congenital bilateral absence of the vas deferens (CBAVD); however, it is unknown why CFTR heterozygotes are at increased disease risk. We recently showed that common CFTR variants are associated with aberrantly spliced transcripts. Here, we genotyped for common CFTR variants and tested for associations in two ICP (ICP-A: 126 patients, 319 controls; ICP-B: 666 patients, 1,181 controls) and a CBAVD population (305 patients, 319 controls). Haplotype H10 (TG11-T7-470V) conferred protection (ICP-A: OR 0.19, P<0.0001; ICP-B: OR 0.78, P = 0.06; CBAVD OR 0.08, P<0.001), whereas haplotype H3 (TG10-T7-470M) increased disease risk (ICP-A: OR 8.34, P = 0.003; ICP-B: OR 1.88, P = 0.007; CBAVD: OR 5.67, P = 0.01). The risk of heterozygous CFTR mutations carriers for ICP (OR 2.44, P<0.001) and CBAVD (OR 14.73, P<0.001) was fully abrogated by the H10/H10 genotype. Similarly, ICP risk of heterozygous p.Asn34Ser SPINK1 mutation carriers (OR 10.34, P<0.001) was compensated by H10/H10. Thus, common CFTR haplotypes modulate ICP and CBAVD susceptibility alone and in heterozygous CFTR and p.Asn34Ser mutation carriers. Determination of these haplotypes helps to stratify carriers into high- and low-risk subjects, providing helpful information for genetic counseling.

An unusual splice defect in the mitofusin 2 gene (MFN2) is associated with degenerative axonopathy in Tyrolean Grey cattle

Tyrolean Grey cattle represent a local breed with a population size of approximately 5000 registered cows. In 2003, a previously unknown neurological disorder was recognized in Tyrolean Grey cattle. The clinical signs of the disorder are similar to those of bovine progressive degenerative myeloencephalopathy (weaver syndrome) in Brown Swiss cattle but occur much earlier in life. The neuropathological investigation of an affected calf showed axonal degeneration in the central nervous system (CNS) and femoral nerve. The pedigrees of the affected calves suggested a monogenic autosomal recessive inheritance. We localized the responsible mutation to a 1.9 Mb interval on chromosome 16 by genome-wide association and haplotype mapping. The MFN2 gene located in this interval encodes mitofusin 2, a mitochondrial membrane protein. A heritable human axonal neuropathy, Charcot-Marie-Tooth disease-2A2 (CMT2A2), is caused by MFN2 mutations. Therefore, we considered MFN2 a positional and functional candidate gene and performed mutation analysis in affected and control Tyrolean Grey cattle. We did not find any non-synonymous variants. However, we identified a perfectly associated silent SNP in the coding region of exon 20 of the MFN2 gene. This SNP is located within a putative exonic splice enhancer (ESE) and the variant allele leads to partial retention of the entire intron 19 and a premature stop codon in the aberrant MFN2 transcript. Thus we have identified a highly unusual splicing defect, where an exonic single base exchange leads to the retention of the preceding intron. This splicing defect represents a potential explanation for the observed degenerative axonopathy. Marker assisted selection can now be used to eliminate degenerative axonopathy from Tyrolean Grey cattle.

Dual origins of dairy cattle farming - evidence from a comprehensive survey of European Y-chromosomal variation

BACKGROUND: Diversity patterns of livestock species are informative to the history of agriculture and indicate uniqueness of breeds as relevant for conservation. So far, most studies on cattle have focused on mitochondrial and autosomal DNA variation. Previous studies of Y-chromosomal variation, with limited breed panels, identified two Bos taurus (taurine) haplogroups (Y1 and Y2; both composed of several haplotypes) and one Bos indicus (indicine/zebu) haplogroup (Y3), as well as a strong phylogeographic structuring of paternal lineages. METHODOLOGY AND PRINCIPAL FINDINGS: Haplogroup data were collected for 2087 animals from 138 breeds. For 111 breeds, these were resolved further by genotyping microsatellites INRA189 (10 alleles) and BM861 (2 alleles). European cattle carry exclusively taurine haplotypes, with the zebu Y-chromosomes having appreciable frequencies in Southwest Asian populations. Y1 is predominant in northern and north-western Europe, but is also observed in several Iberian breeds, as well as in Southwest Asia. A single Y1 haplotype is predominant in north-central Europe and a single Y2 haplotype in central Europe. In contrast, we found both Y1 and Y2 haplotypes in Britain, the Nordic region and Russia, with the highest Y-chromosomal diversity seen in the Iberian Peninsula. CONCLUSIONS: We propose that the homogeneous Y1 and Y2 regions reflect founder effects associated with the development and expansion of two groups of dairy cattle, the pied or red breeds from the North Sea and Baltic coasts and the spotted, yellow or brown breeds from Switzerland, respectively. The present Y1-Y2 contrast in central Europe coincides with historic, linguistic, religious and cultural boundaries.

Molecular analysis of carbohydrate N-acetylgalactosamine 4-O sulfotransferase 8 (CHST8) as a candidate gene for bovine spongiform encephalopathy susceptibility

Endogenous prion proteins (PrP) play the central role in the pathogenesis of transmissible spongiform encephalopathies. The carbohydrate N-acetylgalactosamine 4-O sulfotransferase 8 (CHST8) promotes the conversion of the cellular PrP(C) into the pathogenic PrP(d). Six sequence variants within the CHST8 gene were identified by comparative sequencing and genotyped for a sample of 623 animals comprising bovine spongiform encephalopathy (BSE)-affected and healthy control cows representing German Fleckvieh (German Simmental), German Holstein (Holstein-Friesian) and Brown Swiss. Significant differences in the allele, genotype and haplotype frequencies between BSE-affected and healthy cows indicate an association of sequence variant g.37254017G>T with the development of the disease in Brown Swiss cattle.

Genetic variations in bile acid homeostasis are not overrepresented in alcoholic cirrhosis compared to patients with heavy alcohol abuse and absent liver disease

Increased serum bile salt levels have been associated to a single-nucleotide polymorphism in the bile salt export pump (BSEP; ABCB11) in several acquired cholestatic liver diseases but there is little evidence in alcoholic liver disease (ALD). Furthermore, a crosstalk between vitamin D and bile acid synthesis has recently been discovered. Whether this crosstalk has an influence on the course of ALD is unclear to date. Our aim was to analyse the role of genetic polymorphisms in BSEP and the vitamin D receptor gene (NR1I1) on the emergence of cirrhosis in patients with ALD. Therefore, 511 alcoholic patients (131 with cirrhosis and 380 without cirrhosis) underwent ABCB11 genotyping (rs2287622). Of these, 321 (131 with cirrhosis and 190 without cirrhosis) were also tested for NR1I1 polymorphisms (bat-haplotype: BsmI rs1544410, ApaI rs7975232 and TaqI rs731236). Frequencies of ABCB11 and NR1I1 genotypes and haplotypes were compared between alcoholic patients with and without cirrhosis and correlated to serum bile salt, bilirubin and aspartate aminotransferase levels in those with cirrhosis. Frequencies of ABCB11 and NR1I1 genotypes and haplotypes did not differ between the two subgroups and no significant association between genotypes/haplotypes and liver function tests could be determined for neither polymorphism. We conclude that ABCB11 and NR1I1 polymorphisms are obviously not associated with development of cirrhosis in patients with ALD.

β-Arrestin2 influences the response to methadone in opioid-dependent patients

β-Arrestin2 (ARRB2) is a component of the G-protein-coupled receptor complex and is involved in μ-opioid and dopamine D(2) receptor signaling, two central processes in methadone signal transduction. We analyzed 238 patients in methadone maintenance treatment (MMT) and identified a haplotype block (rs34230287, rs3786047, rs1045280 and rs2036657) spanning almost the entire ARRB2 locus. Although none of these single nucleotide polymorphisms (SNPs) leads to a change in amino-acid sequence, we found that for all the SNPs analyzed, with exception of rs34230287, homozygosity for the variant allele confers a nonresponding phenotype (n=73; rs1045280C and rs2036657G: OR=3.1, 95% CI=1.5-6.3, P=0.004; rs3786047A: OR=2.5, 95% CI=1.2-5.1, P=0.02) also illustrated by a 12-fold shorter period of negative urine screening (P=0.01). The ARRB2 genotype may thus contribute to the interindividual variability in the response to MMT and help to predict response to treatment.

The Impact of PPARγ Genetic Variants on IBD Susceptibility and IBD Disease Course

PPARγ is a nuclear receptor that regulates numerous pathways including cytokine expression and immune responses and plays an important role in controlling colon inflammation. We aimed at determining the occurring PPARγ SNPs, at predicting the haplotypes, and at determining the frequency outcome in inflammatory bowel disease (IBD) patients in comparison with healthy controls. We determined genetic variants in the coding exons and flanking intronic sequences of the NR1C3 gene in 284 IBD patients and 194 controls and predicted NR1C3 haplotypes via bioinformatic analysis. We investigated whether certain NR1C3 variants are associated with susceptibility to IBD or its disease course. None of the detected 22 NR1C3 variants were associated with IBD. Two variants with allelic frequencies over 1% were included in haplotype/diplotype analyses. None of the NR3C1 haplotypes showed association with IBD development or disease course. We conclude that NR1C3 haplotypes are not related to IBD susceptibility or IBD disease activity.

Globally dispersed Y chromosomal haplotypes in wild and domestic sheep

To date, investigations of genetic diversity and the origins of domestication in sheep have utilised autosomal microsatellites and variation in the mitochondrial genome. We present the first analysis of both domestic and wild sheep using genetic markers residing on the ovine Y chromosome. Analysis of a single nucleotide polymorphism (oY1) in the SRY promoter region revealed that allele A-oY1 was present in all wild bighorn sheep (Ovis canadensis), two subspecies of thinhorn sheep (Ovis dalli), European Mouflon (Ovis musimon) and the Barbary (Ammontragis lervia). A-oY1 also had the highest frequency (71.4%) within 458 domestic sheep drawn from 65 breeds sampled from Africa, Asia, Australia, the Caribbean, Europe, the Middle East and Central Asia. Sequence analysis of a second locus, microsatellite SRYM18, revealed a compound repeat array displaying fixed differences, which identified bighorn and thinhorn sheep as distinct from the European Mouflon and domestic animals. Combined genotypic data identified 11 male-specific haplotypes that represented at least two separate lineages. Investigation of the geographical distribution of each haplotype revealed that one (H6) was both very common and widespread in the global sample of domestic breeds. The remaining haplotypes each displayed more restricted and informative distributions. For example, H5 was likely founded following the domestication of European breeds and was used to trace the recent transportation of animals to both the Caribbean and Australia. A high rate of Y chromosomal dispersal appears to have taken place during the development of domestic sheep as only 12.9% of the total observed variation was partitioned between major geographical regions.