Delineating the genetic heterogeneity of ALS using targeted high-throughput sequencing

BACKGROUND

Over 100 genes have been implicated in the aetiology of amyotrophic lateral sclerosis (ALS). A detailed understanding of their independent and cumulative contributions to disease burden may help guide various clinical and research efforts.

METHODS

Using targeted high-throughput sequencing, we characterised the variation of 10 Mendelian and 23 low penetrance/tentative ALS genes within a population-based cohort of 444 Irish ALS cases (50 fALS, 394 sALS) and 311 age-matched and geographically matched controls.

RESULTS

Known or potential high-penetrance ALS variants were identified within 17.1% of patients (38% of fALS, 14.5% of sALS). 12.8% carried variants of Mendelian disease genes (C9orf72 8.78%; SETX 2.48%; ALS2 1.58%; FUS 0.45%; TARDBP 0.45%; OPTN 0.23%; VCP 0.23%. ANG, SOD1, VAPB 0%), 4.7% carried variants of low penetrance/tentative ALS genes and 9.7% (30% of fALS, 7.1% of sALS) carried previously described ALS variants (C9orf72 8.78%; FUS 0.45%; TARDBP 0.45%). 1.6% of patients carried multiple known/potential disease variants, including all identified carriers of an established ALS variant (p<0.01); TARDBP:c.859G>A(p.[G287S]) (n=2/2 sALS). Comparison of our results with those from studies of other European populations revealed significant differences in the spectrum of disease variation (p=1.7×10(-4)).

CONCLUSIONS

Up to 17% of Irish ALS cases may carry high-penetrance variants within the investigated genes. However, the precise nature of genetic susceptibility differs significantly from that reported within other European populations. Certain variants may not cause disease in isolation and concomitant analysis of disease genes may prove highly important.

combi.pl: a computer program to combine data sets with inconsistent microsatellite marker allele size information

Combining two data sets with allele information from overlapping microsatellite markers is often desirable, particularly in population genetic studies where a substantial body of published data exists. When genotyping is performed in different laboratories, allele size calling may not be presumed to be consistent. Our approach solves this problem by assigning allele sizes across studies using maximum-likelihood theory. Using data overlaps in samples and markers, allele shifts between two studies are calculated for each overlapping marker and a single file containing allele frequencies of consistent alleles is produced. The program (combi.pl) is written in PERL and available at http://data40.uni-tz.gwdg.de/~htaeube.

Multiparametric MRI study of ALS stratified for the C9orf72 genotype

OBJECTIVE

To describe the patterns of cortical and subcortical changes in amyotrophic lateral sclerosis (ALS) stratified for the C9orf72 genotype.

METHODS

A prospective, single-center, single-protocol, gray and white matter magnetic resonance case-control imaging study was undertaken with 30 C9orf72-negative patients with ALS, 9 patients with ALS carrying the C9orf72 hexanucleotide repeat expansion, and 44 healthy controls. Tract-based spatial statistics of multiple white matter diffusion parameters, cortical thickness measurements, and voxel-based morphometry analyses were carried out. All patients underwent comprehensive genetic and neuropsychological profiling.

RESULTS

A congruent pattern of cortical and subcortical involvement was identified in those with the C9orf72 genotype, affecting fusiform, thalamic, supramarginal, and orbitofrontal regions and the Broca area. White matter abnormalities in the C9orf72-negative group were relatively confined to corticospinal and cerebellar pathways with limited extramotor expansion. The body of the corpus callosum and superior motor tracts were affected in both ALS genotypes.

CONCLUSIONS

Extensive cortical and subcortical frontotemporal involvement was identified in association with the C9orf72 genotype, compared to the relatively limited extramotor pathology in patients with C9orf72-negative ALS. The distinctive, genotype-specific pathoanatomical patterns are consistent with the neuropsychological profile of the 2 ALS cohorts. Our findings suggest that previously described extramotor changes in ALS could be largely driven by those with the C9orf72 genotype.

Using reference databases of genetic variation to evaluate the potential pathogenicity of candidate disease variants

The potential pathogenicity of genetic variants identified in disease-based resequencing studies is often overlooked where variants have previously been reported in dbSNP, the 1000 genomes project, or the National Heart, Lung and Blood Institute Exome Sequencing Project (ESP). In this work, we estimate that collectively, these databases capture ∼52% of mutations (dbSNP 50.4%; 1000 genomes 4.8%; and ESP 10.2%) reported as disease causing within phenotype-based locus-specific databases (LSDBs). To investigate whether these mutations may simply represent benign population variants, we evaluated whether the carrier frequencies associated with mutations implicated in amyotrophic lateral sclerosis were higher than what could be accounted for by high-penetrance disease models. In doing so, we have questioned the veracity of 51 mutations, but also demonstrated that each of the three databases included credible disease variants. Our results demonstrate the benefits of using databases such as dbSNP, the 1000 genomes project, and the ESP to evaluate the pathogenicity of putative disease variants, and suggest that many disease mutations reported across LSDBs may not actually be pathogenic. However, they also demonstrate that even in the context of rare Mendelian disorders, the potential pathogenicity of variants reported by these databases should not be overlooked without proper evaluation.

Temporal genetic variation of the red fox, Vulpes vulpes, across western Europe and the British Isles

Quaternary climatic fluctuations have had profound effects on the phylogeographic structure of many species. Classically, species were thought to have become isolated in peninsular refugia, but there is limited evidence that large, non-polar species survived outside traditional refugial areas. We examined the phylogeographic structure of the red fox (Vulpes vulpes), a species that shows high ecological adaptability in the western Palaearctic region. We compared mitochondrial DNA sequences (cytochrome b and control region) from 399 modern and 31 ancient individuals from across Europe. Our objective was to test whether red foxes colonised the British Isles from mainland Europe in the late Pleistocene, or whether there is evidence that they persisted in the region through the Last Glacial Maximum. We found red foxes to show a high degree of phylogeographic structuring across Europe and, consistent with palaeontological and ancient DNA evidence, confirmed via phylogenetic indicators that red foxes were persistent in areas outside peninsular refugia during the last ice age. Bayesian analyses and tests of neutrality indicated population expansion. We conclude that there is evidence that red foxes from the British Isles derived from central European populations that became isolated after the closure of the landbridge with Europe.

Novel SNP Discovery in African Buffalo, Syncerus caffer, using high-throughput Sequencing

The African buffalo, Syncerus caffer, is one of the most abundant and ecologically important species of megafauna in the savannah ecosystem. It is an important prey species, as well as a host for a vast array of nematodes, pathogens and infectious diseases, such as bovine tuberculosis and corridor disease. Large-scale SNP discovery in this species would greatly facilitate further research into the area of host genetics and disease susceptibility, as well as provide a wealth of sequence information for other conservation and genomics studies. We sequenced pools of Cape buffalo DNA from a total of 9 animals, on an ABI SOLiD4 sequencer. The resulting short reads were mapped to the UMD3.1 Bos taurus genome assembly using both BWA and Bowtie software packages. A mean depth of 2.7× coverage over the mapped regions was obtained. Btau4 gene annotation was added to all SNPs identified within gene regions. Bowtie and BWA identified a maximum of 2,222,665 and 276,847 SNPs within the buffalo respectively, depending on analysis method. A panel of 173 SNPs was validated by fluorescent genotyping in 87 individuals. 27 SNPs failed to amplify, and of the remaining 146 SNPs, 43-54% of the Bowtie SNPs and 57-58% of the BWA SNPs were confirmed as polymorphic. dN/dS ratios found no evidence of positive selection, and although there were genes that appeared to be under negative selection, these were more likely to be slowly evolving house-keeping genes.

Runs of homozygosity and population history in cattle

BACKGROUND

Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes that are present in an individual due to parents transmitting identical haplotypes to their offspring. The extent and frequency of ROHs may inform on the ancestry of an individual and its population. Here we use high density (n = 777,962) bi-allelic SNPs in a range of cattle breed samples to correlate ROH with the pedigree-based inbreeding coefficients and to validate subsequent analyses using 54,001 SNP genotypes. This study provides a first testing of the inference drawn from ROH through comparison with estimates of inbreeding from calculations based on the detailed pedigree data available for several breeds.

RESULTS

All animals genotyped on the HD panel displayed at least one ROH that was between 1-5 Mb in length with certain regions of the genome more likely to be involved in a ROH than others. Strong correlations (r = 0.75, p < 0.0001) existed between the pedigree-based inbreeding coefficient and a statistic based on sum of ROH of length > 0.5 KB and suggests that in the absence of an animal's pedigree data, the extent of a genome under ROH may be used to infer aspects of recent population history even from relatively few samples.

CONCLUSIONS

Our findings suggest that ROH are frequent across all breeds but differing patterns of ROH length and burden illustrate variations in breed origins and recent management.

Avian resistance to Campylobacter jejuni colonization is associated with an intestinal immunogene expression signature identified by mRNA sequencing

Campylobacter jejuni is the most common cause of human bacterial gastroenteritis and is associated with several post-infectious manifestations, including onset of the autoimmune neuropathy Guillain-Barré syndrome, causing significant morbidity and mortality. Poorly-cooked chicken meat is the most frequent source of infection as C. jejuni colonizes the avian intestine in a commensal relationship. However, not all chickens are equally colonized and resistance seems to be genetically determined. We hypothesize that differences in immune response may contribute to variation in colonization levels between susceptible and resistant birds. Using high-throughput sequencing in an avian infection model, we investigate gene expression associated with resistance or susceptibility to colonization of the gastrointestinal tract with C. jejuni and find that gut related immune mechanisms are critical for regulating colonization. Amongst a single population of 300 4-week old chickens, there was clear segregation in levels of C. jejuni colonization 48 hours post-exposure. RNAseq analysis of caecal tissue from 14 C. jejuni-susceptible and 14 C. jejuni-resistant birds generated over 363 million short mRNA sequences which were investigated to identify 219 differentially expressed genes. Significantly higher expression of genes involved in the innate immune response, cytokine signaling, B cell and T cell activation and immunoglobulin production, as well as the renin-angiotensin system was observed in resistant birds, suggesting an early active immune response to C. jejuni. Lower expression of these genes in colonized birds suggests suppression or inhibition of a clearing immune response thus facilitating commensal colonization and generating vectors for zoonotic transmission. This study describes biological processes regulating C. jejuni colonization of the avian intestine and gives insight into the differential immune mechanisms incited in response to commensal bacteria in general within vertebrate populations. The results reported here illustrate how an exaggerated immune response may be elicited in a subset of the population, which alters host-microbe interactions and inhibits the commensal state, therefore having wider relevance with regard to inflammatory and autoimmune disease.

Genome-wide associations for milk production and somatic cell score in Holstein-Friesian cattle in Ireland

BACKGROUND

Contemporary dairy breeding goals have broadened to include, along with milk production traits, a number of non-production-related traits in an effort to improve the overall functionality of the dairy cow. Increased indirect selection for resistance to mastitis, one of the most important production-related diseases in the dairy sector, via selection for reduced somatic cell count has been part of these broadened goals. A number of genome-wide association studies have identified genetic variants associated with milk production traits and mastitis resistance, however the majority of these studies have been based on animals which were predominantly kept in confinement and fed a concentrate-based diet (i.e. high-input production systems). This genome-wide association study aims to detect associations using genotypic and phenotypic data from Irish Holstein-Friesian cattle fed predominantly grazed grass in a pasture-based production system (low-input).

RESULTS

Significant associations were detected for milk yield, fat yield, protein yield, fat percentage, protein percentage and somatic cell score using separate single-locus, frequentist and multi-locus, Bayesian approaches. These associations were detected using two separate populations of Holstein-Friesian sires and cows. In total, 1,529 and 37 associations were detected in the sires using a single SNP regression and a Bayesian method, respectively. There were 103 associations in common between the sires and cows across all the traits. As well as detecting associations within known QTL regions, a number of novel associations were detected; the most notable of these was a region of chromosome 13 associated with milk yield in the population of Holstein-Friesian sires.

CONCLUSIONS

A total of 276 of novel SNPs were detected in the sires using a single SNP regression approach. Although obvious candidate genes may not be initially forthcoming, this study provides a preliminary framework upon which to identify the causal mechanisms underlying the various milk production traits and somatic cell score. Consequently this will deepen our understanding of how these traits are expressed.

A genome wide association scan of bovine tuberculosis susceptibility in Holstein-Friesian dairy cattle

Background

Bovine tuberculosis is a significant veterinary and financial problem in many parts of the world. Although many factors influence infection and progression of the disease, there is a host genetic component and dissection of this may enlighten on the wider biology of host response to tuberculosis. However, a binary phenotype of presence/absence of infection presents a noisy signal for genomewide association study.

Methodology/Principal Findings

We calculated a composite phenotype of genetic merit for TB susceptibility based on disease incidence in daughters of elite sires used for artificial insemination in the Irish dairy herd. This robust measure was compared with 44,426 SNP genotypes in the most informative 307 subjects in a genome wide association analysis. Three SNPs in a 65 kb genomic region on BTA 22 were associated (i.e. p<10⁻⁵, peaking at position 59588069, p = 4.02×10⁻⁶) with tuberculosis susceptibility.

Conclusions/Significance

A genomic region on BTA 22 was suggestively associated with tuberculosis susceptibility; it contains the taurine transporter gene SLC6A6, or TauT, which is known to function in the immune system but has not previously been investigated for its role in tuberculosis infection.

The peopling of Europe and the cautionary tale of Y chromosome lineage R-M269

Recently, the debate on the origins of the major European Y chromosome haplogroup R1b1b2-M269 has reignited, and opinion has moved away from Palaeolithic origins to the notion of a younger Neolithic spread of these chromosomes from the Near East. Here, we address this debate by investigating frequency patterns and diversity in the largest collection of R1b1b2-M269 chromosomes yet assembled. Our analysis reveals no geographical trends in diversity, in contradiction to expectation under the Neolithic hypothesis, and suggests an alternative explanation for the apparent cline in diversity recently described. We further investigate the young, STR-based time to the most recent common ancestor estimates proposed so far for R-M269-related lineages and find evidence for an appreciable effect of microsatellite choice on age estimates. As a consequence, the existing data and tools are insufficient to make credible estimates for the age of this haplogroup, and conclusions about the timing of its origin and dispersal should be viewed with a large degree of caution.

Ancient hybridization and an Irish origin for the modern polar bear matriline

BACKGROUND

Polar bears (Ursus maritimus) are among those species most susceptible to the rapidly changing arctic climate, and their survival is of global concern. Despite this, little is known about polar bear species history. Future conservation strategies would significantly benefit from an understanding of basic evolutionary information, such as the timing and conditions of their initial divergence from brown bears (U. arctos) or their response to previous environmental change.

RESULTS

We used a spatially explicit phylogeographic model to estimate the dynamics of 242 brown bear and polar bear matrilines sampled throughout the last 120,000 years and across their present and past geographic ranges. Our results show that the present distribution of these matrilines was shaped by a combination of regional stability and rapid, long-distance dispersal from ice-age refugia. In addition, hybridization between polar bears and brown bears may have occurred multiple times throughout the Late Pleistocene.

CONCLUSIONS

The reconstructed matrilineal history of brown and polar bears has two striking features. First, it is punctuated by dramatic and discrete climate-driven dispersal events. Second, opportunistic mating between these two species as their ranges overlapped has left a strong genetic imprint. In particular, a likely genetic exchange with extinct Irish brown bears forms the origin of the modern polar bear matriline. This suggests that interspecific hybridization not only may be more common than previously considered but may be a mechanism by which species deal with marginal habitats during periods of environmental deterioration.

High concordance of bovine single nucleotide polymorphism genotypes generated using two independent genotyping strategies

Single nucleotide polymorphisms (SNPs) represent the most common form of DNA sequence variation in mammalian livestock genomes. While the past decade has witnessed major advances in SNP genotyping technologies, genotyping errors caused, in part, by the biochemistry underlying the genotyping platform used, can occur. These errors can distort project results and conclusions and can result in incorrect decisions in animal management and breeding programs; hence, SNP genotype calls must be accurate and reliable. In this study, 263 Bos spp. samples were genotyped commercially for a total of 16 SNPs. Of the total possible 4,208 SNP genotypes, 4,179 SNP genotypes were generated, yielding a genotype call rate of 99.31% (standard deviation ± 0.93%). Between 110 and 263 samples were subsequently re-genotyped by us for all 16 markers using a custom-designed SNP genotyping platform, and of the possible 3,819 genotypes a total of 3,768 genotypes were generated (98.70% genotype call rate, SD ± 1.89%). A total of 3,744 duplicate genotypes were generated for both genotyping platforms, and comparison of the genotype calls for both methods revealed 3,741 concordant SNP genotype call rates (99.92% SNP genotype concordance rate). These data indicate that both genotyping methods used can provide livestock geneticists with reliable, reproducible SNP genotypic data for in-depth statistical analysis.

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.

Angiogenin levels and ANG genotypes: dysregulation in amyotrophic lateral sclerosis

Objective

To determine whether 5 single nucleotide polymorphisms (SNPs) associate with ALS in 3 different populations. We also assessed the contribution of genotype to angiogenin levels in plasma and CSF.

Methods

Allelic association statistics were calculated for polymorphisms in the ANG gene in 859 patients and 1047 controls from Sweden, Ireland and Poland. Plasma, serum and CSF angiogenin levels were quantified and stratified according to genotypes across the ANG gene. The contribution of SNP genotypes to variance in circulating angiogenin levels was estimated in patients and controls.

Results

All SNPs showed association with ALS in the Irish group. The SNP rs17114699 replicated in the Swedish cohort. No SNP associated in the Polish cohort. Age- and sex-corrected circulating angiogenin levels were significantly lower in patients than in controls (p<0.001). An allele dose-dependent regulation of angiogenin levels was observed in controls. This regulation was attenuated in the ALS cohort. A significant positive correlation between CSF plasma angiogenin levels was present in controls and abolished in ALS.

Conclusions

ANG variants associate with ALS in the Irish and Swedish populations, but not in the Polish. There is evidence of dysregulation of angiogenin expression in plasma and CSF in sporadic ALS. Angiogenin expression is likely to be important in the pathogenesis of ALS.

Cattle demographic history modelled from autosomal sequence variation

The phylogeography of cattle genetic variants has been extensively described and has informed the history of domestication. However, there remains a dearth of demographic models inferred from such data. Here, we describe sequence diversity at 37 000 bp sampled from 17 genes in cattle from Africa, Europe and India. Clearly distinct population histories are suggested between Bos indicus and Bos taurus, with the former displaying higher diversity statistics. We compare the unfolded site frequency spectra in each to those simulated using a diffusion approximation method and build a best-fitting model of past demography. This implies an earlier, possibly glaciation-induced population bottleneck in B. taurus ancestry with a later, possibly domestication-associated demographic constriction in B. indicus. Strikingly, the modelled indicine history also requires a majority secondary admixture from the South Asian aurochs, indicating a complex, more diffuse domestication process. This perhaps involved multiple domestications and/or introgression from wild oxen to domestic herds; the latter is plausible from archaeological evidence of contemporaneous wild and domestic remains across different regions of South Asia.

The differential evolutionary dynamics of avian cytokine and TLR gene classes

The potential for investigating immune gene diversity has been greatly enhanced by recent advances in sequencing power. In this study, variation at two categories of avian immune genes with differing functional roles, pathogen detection and mediation of immune mechanisms, was examined using high-throughput sequencing. TLRs identify and alert the immune system by detecting molecular motifs that are conserved among pathogenic microorganisms, whereas cytokines act as mediators of resulting inflammation and immunity. Nine genes from each class were resequenced in a panel of domestic chickens and wild jungle fowl (JF). Tests on population-wide genetic variation between the gene classes indicated that allele frequency spectra at each group were distinctive. TLRs showed evidence pointing toward directional selection, whereas cytokines had signals more suggestive of frequency-dependent selection. This difference persisted between the distributions considering only coding sites, suggesting functional relevance. The unique patterns of variation at each gene class may be constrained by their different functional roles in the immune response. TLRs identify a relatively limited number of exogeneous pathogenic-related patterns and would be required to adapt quickly in response to evolving novel microbes encountered in new environmental niches. In contrast, cytokines interact with many molecules in mediating the power of immune mechanisms, and accordingly respond to the selective stimuli of many infectious diseases. Analyses also indicated that a general pattern of high variability has been enhanced by widespread genetic exchange between chicken and red JF, and possibly between chicken and gray JF at TLR1LA and TLR2A.

A complete mitochondrial genome sequence from a mesolithic wild aurochs (Bos primigenius)

BACKGROUND

The derivation of domestic cattle from the extinct wild aurochs (Bos primigenius) has been well-documented by archaeological and genetic studies. Genetic studies point towards the Neolithic Near East as the centre of origin for Bos taurus, with some lines of evidence suggesting possible, albeit rare, genetic contributions from locally domesticated wild aurochsen across Eurasia. Inferences from these investigations have been based largely on the analysis of partial mitochondrial DNA sequences generated from modern animals, with limited sequence data from ancient aurochsen samples. Recent developments in DNA sequencing technologies, however, are affording new opportunities for the examination of genetic material retrieved from extinct species, providing new insight into their evolutionary history. Here we present DNA sequence analysis of the first complete mitochondrial genome (16,338 base pairs) from an archaeologically-verified and exceptionally-well preserved aurochs bone sample.

METHODOLOGY

DNA extracts were generated from an aurochs humerus bone sample recovered from a cave site located in Derbyshire, England and radiocarbon-dated to 6,738+/-68 calibrated years before present. These extracts were prepared for both Sanger and next generation DNA sequencing technologies (Illumina Genome Analyzer). In total, 289.9 megabases (22.48%) of the post-filtered DNA sequences generated using the Illumina Genome Analyzer from this sample mapped with confidence to the bovine genome. A consensus B. primigenius mitochondrial genome sequence was constructed and was analysed alongside all available complete bovine mitochondrial genome sequences.

CONCLUSIONS

For all nucleotide positions where both Sanger and Illumina Genome Analyzer sequencing methods gave high-confidence calls, no discrepancies were observed. Sequence analysis reveals evidence of heteroplasmy in this sample and places this mitochondrial genome sequence securely within a previously identified aurochsen haplogroup (haplogroup P), thus providing novel insights into pre-domestic patterns of variation. The high proportion of authentic, endogenous aurochs DNA preserved in this sample bodes well for future efforts to determine the complete genome sequence of a wild ancestor of domestic cattle.

Evidence of the adaptive evolution of immune genes in chicken

The basis for understanding the characteristics of gene functional categories in chicken has been enhanced by the ongoing sequencing of the zebra finch genome, the second bird species to be extensively sequenced. This sequence provides an avian context for examining how variation in chicken has evolved since its divergence from its common ancestor with zebra finch as well as well as a calibrating point for studying intraspecific diversity within chicken. Immune genes have been subject to many selective processes during their evolutionary history: this gene class was investigated here in a set of orthologous chicken and zebra finch genes with functions assigned from the human ortholog. Tests demonstrated that nonsynonymous sites at immune genes were highly conserved both in chicken and on the avian lineage. McDonald-Kreitman tests provided evidence of adaptive evolution and a higher rate of selection on fixation of nonsynonymous substitutions at immune genes compared to that at non-immune genes. Further analyses showed that GC content was much higher in chicken than in zebra finch genes, and was significantly elevated in both species' immune genes. Pathogen challenges are likely to have driven the selective forces that have shaped variation at chicken immune genes, and continue to restrict diversity in this functional class.

Maternal and paternal genealogy of Eurasian taurine cattle (Bos taurus)

Maternally inherited mitochondrial DNA (mtDNA) has been used extensively to determine origin and diversity of taurine cattle (Bos taurus) but global surveys of paternally inherited Y-chromosome diversity are lacking. Here, we provide mtDNA information on previously uncharacterised Eurasian breeds and present the most comprehensive Y-chromosomal microsatellite data on domestic cattle to date. The mitochondrial haplogroup T3 was the most frequent, whereas T4 was detected only in the Yakutian cattle from Siberia. The mtDNA data indicates that the Ukrainian and Central Asian regions are zones where hybrids between taurine and zebu (B. indicus) cattle have existed. This zebu influence appears to have subsequently spread into southern and southeastern European breeds. The most common Y-chromosomal microsatellite haplotype, termed here as H11, showed an elevated frequency in the Eurasian sample set compared with that detected in Near Eastern and Anatolian breeds. The taurine Y-chromosomal microsatellite haplotypes were found to be structured in a network according to the Y-haplogroups Y1 and Y2. These data do not support the recent hypothesis on the origin of Y1 from the local European hybridization of cattle with male aurochsen. Compared with mtDNA, the intensive culling of breeding males and male-mediated crossbreeding of locally raised native breeds has accelerated loss of Y-chromosomal variation in domestic cattle, and affected the contribution of genetic drift to diversity. In conclusion, to maintain diversity, breeds showing rare Y-haplotypes should be prioritised in the conservation of cattle genetic resources.

Evidence of balanced diversity at the chicken interleukin 4 receptor alpha chain locus

BACKGROUND

The comparative analysis of genome sequences emerging for several avian species with the fully sequenced chicken genome enables the genome-wide investigation of selective processes in functionally important chicken genes. In particular, because of pathogenic challenges it is expected that genes involved in the chicken immune system are subject to particularly strong adaptive pressure. Signatures of selection detected by inter-species comparison may then be investigated at the population level in global chicken populations to highlight potentially relevant functional polymorphisms.

RESULTS

Comparative evolutionary analysis of chicken (Gallus gallus) and zebra finch (Taeniopygia guttata) genes identified interleukin 4 receptor alpha-chain (IL-4Ralpha), a key cytokine receptor as a candidate with a significant excess of substitutions at nonsynonymous sites, suggestive of adaptive evolution. Resequencing and detailed population genetic analysis of this gene in diverse village chickens from Asia and Africa, commercial broilers, and in outgroup species red jungle fowl (JF), grey JF, Ceylon JF, green JF, grey francolin and bamboo partridge, suggested elevated and balanced diversity across all populations at this gene, acting to preserve different high-frequency alleles at two nonsynonymous sites.

CONCLUSION

Haplotype networks indicate that red JF is the primary contributor of diversity at chicken IL-4Ralpha: the signature of variation observed here may be due to the effects of domestication, admixture and introgression, which produce high diversity. However, this gene is a key cytokine-binding receptor in the immune system, so balancing selection related to the host response to pathogens cannot be excluded.

The genome sequence of taurine cattle: a window to ruminant biology and evolution

To understand the biology and evolution of ruminants, the cattle genome was sequenced to about sevenfold coverage. The cattle genome contains a minimum of 22,000 genes, with a core set of 14,345 orthologs shared among seven mammalian species of which 1217 are absent or undetected in noneutherian (marsupial or monotreme) genomes. Cattle-specific evolutionary breakpoint regions in chromosomes have a higher density of segmental duplications, enrichment of repetitive elements, and species-specific variations in genes associated with lactation and immune responsiveness. Genes involved in metabolism are generally highly conserved, although five metabolic genes are deleted or extensively diverged from their human orthologs. The cattle genome sequence thus provides a resource for understanding mammalian evolution and accelerating livestock genetic improvement for milk and meat production.

Cytochrome b sequences of ancient cattle and wild ox support phylogenetic complexity in the ancient and modern bovine populations

Mitochondrial DNA has been the traditional marker for the study of animal domestication, as its high mutation rate allows for the accumulation of molecular diversity within the time frame of domestic history. Additionally, it is exclusively maternally inherited and haplotypes become part of the domestic gene pool via actual capture of a female animal rather than by interbreeding with wild populations. Initial studies of British aurochs identified a haplogroup, designated P, which was found to be highly divergent from all known domestic haplotypes over the most variable portion of the D-loop. Additional analysis of a large and geographically representative sample of aurochs from northern and central Europe found an additional, separate aurochs haplotype, E. Until recently, the European aurochs appeared to have no matrilinear descendants among the publicly available modern cattle control regions sequenced; if aurochs mtDNA was incorporated into the domestic population, aurochs either formed a very small proportion of modern diversity or had been subsequently lost. However, a haplogroup P sequence has recently been found in a modern sample, along with a new divergent haplogroup called Q. Here we confirm the outlying status of the novel Q and E haplogroups and the modern P haplogroup sequence as a descendent of European aurochs, by retrieval and analysis of cytochrome b sequence data from twenty ancient wild and domesticated cattle archaeological samples.

Genome-wide survey of SNP variation uncovers the genetic structure of cattle breeds

The imprints of domestication and breed development on the genomes of livestock likely differ from those of companion animals. A deep draft sequence assembly of shotgun reads from a single Hereford female and comparative sequences sampled from six additional breeds were used to develop probes to interrogate 37,470 single-nucleotide polymorphisms (SNPs) in 497 cattle from 19 geographically and biologically diverse breeds. These data show that cattle have undergone a rapid recent decrease in effective population size from a very large ancestral population, possibly due to bottlenecks associated with domestication, selection, and breed formation. Domestication and artificial selection appear to have left detectable signatures of selection within the cattle genome, yet the current levels of diversity within breeds are at least as great as exists within humans.