An extensive candidate gene approach to speciation: diversity, divergence and linkage disequilibrium in candidate pigmentation genes across the European crow hybrid zone

Colouration patterns have an important role in adaptation and speciation. The European crow system, in which all-black carrion crows and grey-coated hooded crows meet in a narrow hybrid zone, is a prominent example. The marked phenotypic difference is maintained by assortative mating in the absence of neutral genetic divergence, suggesting the presence of few pigmentation genes of major effect. We made use of the rich phenotypic and genetic resources in mammals and identified a comprehensive panel of 95 candidate pigmentation genes for birds. Based on functional annotation, we chose a subset of the most promising 37 candidates, for which we developed a marker system that demonstrably works across the avian phylogeny. In total, we sequenced 107 amplicons (B3 loci per gene, totalling 60kb) in population samples of crows (n¼23 for each taxon). Tajima’s D, Fu’s FS, DHEW and HKA (Hudson–Kreitman–Aguade) statistics revealed several amplicons that deviated from neutrality; however, none of these showed significantly elevated differentiation between the two taxa. Hence, colour divergence in this system may be mediated by uncharacterized pigmentation genes or regulatory regions outside genes. Alternatively, the observed high population recombination rate (4NerB0.03), with overall linkage disequilibrium dropping rapidly within the order of few 100bp, may compromise the power to detect causal loci with nearby markers. Our results add to the debate as to the utility of candidate gene approaches in relation to genomic features and the genetic architecture of the phenotypic trait in question

Twisted signatures of GC-biased gene conversion embedded in an evolutionary stable karyotype

The genomes of many vertebrates show a characteristic heterogeneous distribution of GC content, the so-called GC isochore structure. The origin of isochores has been explained via the mechanism of GC-biased gene conversion (gBGC). However, although the isochore structure is declining in many mammalian genomes, the heterogeneity in GC content is being reinforced in the avian genome. Despite this discrepancy, which remains unexplained, examinations of individual substitution frequencies in mammals and birds are both consistent with the gBGC model of isochore evolution. On the other hand, a negative correlation between substitution and recombination rate found in the chicken genome is inconsistent with the gBGC model. It should therefore be important to consider along with gBGC other consequences of recombination on the origin and fate of mutations, as well as to account for relationships between recombination rate and other genomic features. We therefore developed an analytical model to describe the substitution patterns found in the chicken genome, and further investigated the relationships between substitution patterns and several genomic features in a rigorous statistical framework. Our analysis indicates that GC content itself, either directly or indirectly via interrelations to other genomic features, has an impact on the substitution pattern. Further, we suggest that this phenomenon is particularly visible in avian genomes due to their unusually low rate of chromosomal evolution. Because of this, interrelations between GC content and other genomic features are being reinforced, and are as such more pronounced in avian genomes as compared with other vertebrate genomes with a less stable karyotype

Conservation of Neutral Substitution Rate and Substitutional Asymmetries in Mammalian Genes

Local variation in neutral substitution rate across mammalian genomes is governed by several factors, including sequence context variables and structural variables. In addition, the interplay of replication and transcription, known to induce a strand bias in mutation rate, gives rise to variation in substitutional strand asymmetries. Here, we address the conservation of variation in mutation rate and substitutional strand asymmetries using primate- and rodent-specific repeat elements located within the introns of protein-coding genes. We find significant but weak conservation of local mutation rates between human and mouse orthologs. Likewise, substitutional strand asymmetries are conserved between human and mouse, where substitution rate asymmetries show a higher degree of conservation than mutation rate. Moreover, we provide evidence that replication and transcription are correlated to the strength of substitutional asymmetries. The effect of transcription is particularly visible for genes with highly conserved gene expression. In comparison with replication and transcription, mutation rate influences the strength of substitutional asymmetries only marginally

Heterozygosity increases microsatellite mutation rate, linking it to demographic history

<p>Abstract</p> <p>Background</p> <p>Biochemical experiments in yeast suggest a possible mechanism that would cause heterozygous sites to mutate faster than equivalent homozygous sites. If such a process operates, it could undermine a key assumption at the core of population genetic theory, namely that mutation rate and population size are indpendent, because population expansion would increase heterozygosity that in turn would increase mutation rate. Here we test this hypothesis using both direct counting of microsatellite mutations in human pedigrees and an analysis of the relationship between microsatellite length and patterns of demographically-induced variation in heterozygosity.</p> <p>Results</p> <p>We find that microsatellite alleles of any given length are more likely to mutate when their homologue is unusually different in length. Furthermore, microsatellite lengths in human populations do not vary randomly, but instead exhibit highly predictable trends with both distance from Africa, a surrogate measure of genome-wide heterozygosity, and modern population size. This predictability remains even after statistically controlling for non-independence due to shared ancestry among populations.</p> <p>Conclusion</p> <p>Our results reveal patterns that are unexpected under classical population genetic theory, where no mechanism exists capable of linking allele length to extrinsic variables such as geography or population size. However, the predictability of microsatellite length is consistent with heterozygote instability and suggest that this has an important impact on microsatellite evolution. Whether similar processes impact on single nucleotide polymorphisms remains unclear.</p

Mitochondrial DNA phylogeography and population history of the grey wolf Canis lupus

The grey wolf (Canis lupus) and coyote (C. latrans) are highly mobile carnivores that disperse over great distances in search of territories and mates. Previous genetic studies have shown little geographical structure in either species. However, population genetic structure is also influenced by past isolation events and population fluctuations during glacial periods. In this study, control region sequence data from a worldwide sample of grey wolves and a more limited sample of coyotes were analysed. The results suggest that fluctuating population sizes during the late Pleistocene have left a genetic signature on levels of variation in both species. Genealogical measures of nucleotide diversity suggest that historical population sizes were much larger in both species and grey wolves were more numerous than coyotes. Currently, about 300 000 wolves and 7 million coyotes exist. In grey wolves, genetic diversity is greater than that predicted from census population size, reflecting recent historical population declines. By contrast, nucleotide diversity in coyotes is smaller than that predicted by census population size, reflecting a recent population expansion following the extirpation of wolves from much of North America. Both species show little partitioning of haplotypes on continental or regional scales. However, a statistical parsimony analysis indicates local genetic structure that suggests recent restricted gene flow.Peer Reviewe

The role of mutation rate variation and genetic diversity in the architecture of human disease

Background We have investigated the role that the mutation rate and the structure of genetic variation at a locus play in determining whether a gene is involved in disease. We predict that the mutation rate and its genetic diversity should be higher in genes associated with disease, unless all genes that could cause disease have already been identified. Results Consistent with our predictions we find that genes associated with Mendelian and complex disease are substantially longer than non-disease genes. However, we find that both Mendelian and complex disease genes are found in regions of the genome with relatively low mutation rates, as inferred from intron divergence between humans and chimpanzees, and they are predicted to have similar rates of non-synonymous mutation as other genes. Finally, we find that disease genes are in regions of significantly elevated genetic diversity, even when variation in the rate of mutation is controlled for. The effect is small nevertheless. Conclusions Our results suggest that gene length contributes to whether a gene is associated with disease. However, the mutation rate and the genetic architecture of the locus appear to play only a minor role in determining whether a gene is associated with disease

Optimization of sequence alignment for simple sequence repeat regions

Abstract Background Microsatellites, or simple sequence repeats (SSRs), are tandemly repeated DNA sequences, including tandem copies of specific sequences no longer than six bases, that are distributed in the genome. SSR has been used as a molecular marker because it is easy to detect and is used in a range of applications, including genetic diversity, genome mapping, and marker assisted selection. It is also very mutable because of slipping in the DNA polymerase during DNA replication. This unique mutation increases the insertion/deletion (INDELs) mutation frequency to a high ratio - more than other types of molecular markers such as single nucleotide polymorphism (SNPs). SNPs are more frequent than INDELs. Therefore, all designed algorithms for sequence alignment fit the vast majority of the genomic sequence without considering microsatellite regions, as unique sequences that require special consideration. The old algorithm is limited in its application because there are many overlaps between different repeat units which result in false evolutionary relationships. Findings To overcome the limitation of the aligning algorithm when dealing with SSR loci, a new algorithm was developed using PERL script with a Tk graphical interface. This program is based on aligning sequences after determining the repeated units first, and the last SSR nucleotides positions. This results in a shifting process according to the inserted repeated unit type. When studying the phylogenic relations before and after applying the new algorithm, many differences in the trees were obtained by increasing the SSR length and complexity. However, less distance between different linage had been observed after applying the new algorithm. Conclusions The new algorithm produces better estimates for aligning SSR loci because it reflects more reliable evolutionary relations between different linages. It reduces overlapping during SSR alignment, which results in a more realistic phylogenic relationship.</p

Simple sequence repeat variation in the Daphnia pulex genome

Background: Simple sequence repeats (SSRs) are highly variable features of all genomes. Their rapid evolution makes them useful for tracing the evolutionary history of populations and investigating patterns of selection and mutation across gnomes. The recently sequenced Daphnia pulex genome provides us with a valuable data set to study the mode and tempo of SSR evolution, without the inherent biases that accompany marker selection. Results: Here we catalogue SSR loci in the Daphnia pulex genome with repeated motif sizes of 1-100 nucleotides with a minimum of 3 perfect repeats. We then used whole genome shotgun reads to determine the average heterozygosity of each SSR type and the relationship that it has to repeat number, motif size, motif sequence, and distribution of SSR loci. We find that SSR heterozygosity is motif specific, and positively correlated with repeat number as well as motif size. For non-repeat unit polymorphisms, we identify a motif-dependent end-nucleotide polymorphism bias that may contribute to the patterns of abundance for specific homopolymers, dimers, and trimers. Our observations confirm the high frequency of multiple unit variation (multistep) at large microsatellite loci, and further show that the occurrence of multiple unit variation is dependent on both repeat number and motif size. Using the Daphnia pulex genetic map, we show a positive correlation between dimer and trimer frequency and recombination. Conclusions: This genome-wide analysis of SSR variation in Daphnia pulex indicates that several aspects of SSR variation are motif dependent and suggests that a combination of unit length variation and end repeat biased base substitution contribute to the unique spectrum of SSR repeat loci

Interphase chromosome positioning in in vitro porcine cells and ex vivo porcine tissues

Copyright @ 2012 The Authors. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and 85 reproduction in any medium, provided the original author and source are credited. The article was made available through the Brunel University Open Access Publishing Fund.BACKGROUND: In interphase nuclei of a wide range of species chromosomes are organised into their own specific locations termed territories. These chromosome territories are non-randomly positioned in nuclei which is believed to be related to a spatial aspect of regulatory control over gene expression. In this study we have adopted the pig as a model in which to study interphase chromosome positioning and follows on from other studies from our group of using pig cells and tissues to study interphase genome re-positioning during differentiation. The pig is an important model organism both economically and as a closely related species to study human disease models. This is why great efforts have been made to accomplish the full genome sequence in the last decade. RESULTS: This study has positioned most of the porcine chromosomes in in vitro cultured adult and embryonic fibroblasts, early passage stromal derived mesenchymal stem cells and lymphocytes. The study is further expanded to position four chromosomes in ex vivo tissue derived from pig kidney, lung and brain. CONCLUSIONS: It was concluded that porcine chromosomes are also non-randomly positioned within interphase nuclei with few major differences in chromosome position in interphase nuclei between different cell and tissue types. There were also no differences between preferred nuclear location of chromosomes in in vitro cultured cells as compared to cells in tissue sections. Using a number of analyses to ascertain by what criteria porcine chromosomes were positioned in interphase nuclei; we found a correlation with DNA content.This study is partly supported by Sygen International PLC

Genetic Variation and De Novo Mutations in the Parthenogenetic Caucasian Rock Lizard Darevskia unisexualis

Unisexual all-female lizards of the genus Darevskia that are well adapted to various habitats are known to reproduce normally by true parthenogenesis. Although they consist of unisexual lineages and lack effective genetic recombination, they are characterized by some level of genetic polymorphism. To reveal the mutational contribution to overall genetic variability, the most straightforward and conclusive way is the direct detection of mutation events in pedigree genotyping. Earlier we selected from genomic library of D. unisexualis two polymorphic microsatellite containg loci Du281 and Du215. In this study, these two loci were analyzed to detect possible de novo mutations in 168 parthenogenetic offspring of 49 D. unisexualis mothers and in 147 offspring of 50 D. armeniaca mothers . No mutant alleles were detected in D. armeniaca offspring at both loci, and in D. unisexualis offspring at the Du215 locus. There were a total of seven mutational events in the germ lines of four of the 49 D. unisexualis mothers at the Du281 locus, yielding the mutation rate of 0.1428 events per germ line tissue. Sequencing of the mutant alleles has shown that most mutations occur via deletion or insertion of single microsatellite repeat being identical in all offspring of the family. This indicates that such mutations emerge at the early stages of embryogenesis. In this study we characterized single highly unstable (GATA)n containing locus in parthenogenetic lizard species D. unisexualis. Besides, we characterized various types of mutant alleles of this locus found in the D. unisexualis offspring of the first generation. Our data has shown that microsatellite mutations at highly unstable loci can make a significant contribution to population variability of parthenogenetic lizards