Publications

Generic chemical structures in patents — an evaluation of the Sheffield University research work

Fragile X syndrome: genetic localisation by linkage mapping of two microsatellite repeats FRAXAC1 and FRAXAC2 which immediately flank the fragile site.

We report the genetic localisation of the fragile site at Xq27.3 associated with fragile X syndrome. The position of the fragile site within the multipoint linkage map was determined using two polymorphic microsatellite AC repeat markers FRAXAC1 and FRAXAC2. These markers were physically located within 10 kilobases and on either side of the p(CCG)n repeat responsible for the fragile site. FRAXAC1 has five alleles with heterozygosity of 44% and is in strong linkage disequilibrium with FRAXAC2 which has eight alleles and a heterozygosity of 71%. No recombination was observed either between these markers in 40 normal CEPH pedigrees or with the fragile X in affected pedigrees. These markers provide the means for accurate diagnosis of the fragile X genotype in families by rapid polymerase chain reaction analysis and were used to position the fragile X within the multipoint map of the X chromosome to a position 3.7 cM distal to DXS297 and 1.2 cM proximal to DXS296.

Evolutionary Dynamics of Asexual Hypermutators Adapting to a Novel Environment

Abstract How microbes adapt to a novel environment is a central question in evolutionary biology. Although adaptive evolution must be fueled by beneficial mutations, whether higher mutation rates facilitate the rate of adaptive evolution remains unclear. To address this question, we cultured Escherichia coli hypermutating populations, in which a defective methyl-directed mismatch repair pathway causes a 140-fold increase in single-nucleotide mutation rates. In parallel with wild-type E. coli, populations were cultured in tubes containing Luria-Bertani broth, a complex medium known to promote the evolution of subpopulation structure. After 900 days of evolution, in three transfer schemes with different population-size bottlenecks, hypermutators always exhibited similar levels of improved fitness as controls. Fluctuation tests revealed that the mutation rates of hypermutator lines converged evolutionarily on those of wild-type populations, which may have contributed to the absence of fitness differences. Further genome-sequence analysis revealed that, although hypermutator populations have higher rates of genomic evolution, this largely reflects strong genetic linkage. Despite these linkage effects, the evolved population exhibits parallelism in fixed mutations, including those potentially related to biofilm formation, transcription regulation, and mutation-rate evolution. Together, these results are generally inconsistent with a hypothesized positive relationship between the mutation rate and the adaptive speed of evolution, and provide insight into how clonal adaptation occurs in novel environments.

Evolutionary Genomics of Sister Species Differing in Effective Population Sizes and Recombination Rates

Abstract Studies of closely related species with known ecological differences provide exceptional opportunities for understanding the genetic mechanisms of evolution. Here, we compared population-genomics data between D. pulex and D. pulicaria , two reproductively compatible sister species experiencing ecological speciation, the first largely confined to intermittent ponds and the second to permanent lakes in the same geographic region. D. pulicaria has lower genome-wide nucleotide diversity, a smaller effective population size, higher incidence of private alleles, and substantially more linkage-disequilibrium than D. pulex . Functional enrichment analysis revealed that positively selected genes in D. pulicaria are enriched in potentially aging-related categories such as cellular homeostasis, which may explain the extended lifespan in D. pulicaria . We also found that opsin-related genes, which may mediate photoperiodic responses, are under different selection pressures in these two species. Additionally, genes involved in mitochondrial functions, ribosomes, and responses to environmental stimuli are found to be under positive selection in both species. Our results provide insights into the physiological traits that differ within this regionally sympatric sister-species pair that occupies unique microhabitats.

The Linkage-Disequilibrium and Recombinational Landscape in <i>Daphnia pulex</i>

Abstract By revealing the influence of recombinational activity beyond what can be achieved with controlled crosses, measures of linkage disequilibrium (LD) in natural populations provide a powerful means of defining the recombinational landscape within which genes evolve. In one of the most comprehensive studies of this sort ever performed, involving whole-genome analyses on nearly 1,000 individuals of the cyclically parthenogenetic microcrustacean Daphnia pulex, the data suggest a relatively uniform pattern of recombination across the genome. Patterns of LD are quite consistent among populations; average rates of recombination are quite similar for all chromosomes; and although some chromosomal regions have elevated recombination rates, the degree of inflation is not large, and the overall spatial pattern of recombination is close to the random expectation. Contrary to expectations for models in which crossing-over is the primary mechanism of recombination, and consistent with data for other species, the distance-dependent pattern of LD indicates excessively high levels at both short and long distances and unexpectedly low levels of decay at long distances, suggesting significant roles for factors such as nonindependent mutation, population subdivision, and recombination mechanisms unassociated with crossing over. These observations raise issues regarding the classical LD equilibrium model widely applied in population genetics to infer recombination rates across various length scales on chromosomes.

PATTERNS OF GENETIC ARCHITECTURE FOR LIFE-HISTORY TRAITS AND MOLECULAR MARKERS IN A SUBDIVIDED SPECIES

Understanding the utility and limitations of molecular markers for predicting the evolutionary potential of natural populations is important for both evolutionary and conservation genetics. To address this issue, the distribution of genetic variation for quantitative traits and molecular markers is estimated within and among 14 permanent lake populations of Daphnia pulicaria representing two regional groups from Oregon. Estimates of population subdivision for molecular and quantitative traits are concordant, with QST generally exceeding GST. There is no evidence that microsatellites loci are less informative about subdivision for quantitative traits than are allozyme loci. Character-specific comparison of QST and GST support divergent selection pressures among populations for the majority of life-history traits in both coast and mountain regions. The level of within-population variation for molecular markers is uninformative as to the genetic variation maintained for quantitative traits. In D. pulicaria, regional differences in the frequency of sex may contribute to variation in the maintenance of expressed within-population quantitative-genetic variation without substantially impacting diversity at the genic level. These data are compared to an identical dataset for 17 populations of the temporary-pond species, D. pulex.

Periodic Variation of Mutation Rates in Bacterial Genomes Associated with Replication Timing

ABSTRACT The causes and consequences of spatiotemporal variation in mutation rates remain to be explored in nearly all organisms. Here we examine relationships between local mutation rates and replication timing in three bacterial species whose genomes have multiple chromosomes: Vibrio fischeri , Vibrio cholerae , and Burkholderia cenocepacia . Following five mutation accumulation experiments with these bacteria conducted in the near absence of natural selection, the genomes of clones from each lineage were sequenced and analyzed to identify variation in mutation rates and spectra. In lineages lacking mismatch repair, base substitution mutation rates vary in a mirrored wave-like pattern on opposing replichores of the large chromosomes of V. fischeri and V. cholerae , where concurrently replicated regions experience similar base substitution mutation rates. The base substitution mutation rates on the small chromosome are less variable in both species but occur at similar rates to those in the concurrently replicated regions of the large chromosome. Neither nucleotide composition nor frequency of nucleotide motifs differed among regions experiencing high and low base substitution rates, which along with the inferred ~800-kb wave period suggests that the source of the periodicity is not sequence specific but rather a systematic process related to the cell cycle. These results support the notion that base substitution mutation rates are likely to vary systematically across many bacterial genomes, which exposes certain genes to elevated deleterious mutational load. IMPORTANCE That mutation rates vary within bacterial genomes is well known, but the detailed study of these biases has been made possible only recently with contemporary sequencing methods. We applied these methods to understand how bacterial genomes with multiple chromosomes, like those of Vibrio and Burkholderia , might experience heterogeneous mutation rates because of their unusual replication and the greater genetic diversity found on smaller chromosomes. This study captured thousands of mutations and revealed wave-like rate variation that is synchronized with replication timing and not explained by sequence context. The scale of this rate variation over hundreds of kilobases of DNA strongly suggests that a temporally regulated cellular process may generate wave-like variation in mutation risk. These findings add to our understanding of how mutation risk is distributed across bacterial and likely also eukaryotic genomes, owing to their highly conserved replication and repair machinery.

Deleterious mutation accumulation in organelle genomes

Reply to Massey: Drift does influence mutation-rate evolution

Although Massey argues that proteome size (P), not the effective population size (Ne), is the primary determinant of mutation-rate evolution (1), nothing in his comment or in his prior work (2) contradicts our conclusions (3).

Reference genomes

Here are the full assemblies (in fasta format) of the reference genomes used to map variants in <em>P. tetraurelia</em>, <em>P. biaurelia</em>, <em>P. sexaurelia</em> and <em>P. caudatum</em>. These reference genomes were obtained from Mcgrath et al 2014 publications provided below. <br> These files were used to perform analyses in the manuscript - "<strong>A population-genetic lens into the process of gene loss following whole-genome duplication"</strong>.

Guy’s Cancer Cohort: Guy’s Cancer Centre’s Real-World Evidence Programme 5 years later

LTR retroelements in the genome of Daphnia pulex

Abstract Background Long terminal repeat (LTR) retroelements represent a successful group of transposable elements (TEs) that have played an important role in shaping the structure of many eukaryotic genomes. Here, we present a genome-wide analysis of LTR retroelements in Daphnia pulex , a cyclical parthenogen and the first crustacean for which the whole genomic sequence is available. In addition, we analyze transcriptional data and perform transposon display assays of lab-reared lineages and natural isolates to identify potential influences on TE mobility and differences in LTR retroelements loads among individuals reproducing with and without sex. Results We conducted a comprehensive de novo search for LTR retroelements and identified 333 intact LTR retroelements representing 142 families in the D. pulex genome. While nearly half of the identified LTR retroelements belong to the gypsy group, we also found copia (95), BEL/ Pao (66) and DIRS (19) retroelements. Phylogenetic analysis of reverse transcriptase sequences showed that LTR retroelements in the D. pulex genome form many lineages distinct from known families, suggesting that the majority are novel. Our investigation of transcriptional activity of LTR retroelements using tiling array data obtained from three different experimental conditions found that 71 LTR retroelements are actively transcribed. Transposon display assays of mutation-accumulation lines showed evidence for putative somatic insertions for two DIRS retroelement families. Losses of presumably heterozygous insertions were observed in lineages in which selfing occurred, but never in asexuals, highlighting the potential impact of reproductive mode on TE abundance and distribution over time. The same two families were also assayed across natural isolates (both cyclical parthenogens and obligate asexuals) and there were more retroelements in populations capable of reproducing sexually for one of the two families assayed. Conclusions Given the importance of LTR retroelements activity in the evolution of other genomes, this comprehensive survey provides insight into the potential impact of LTR retroelements on the genome of D. pulex , a cyclically parthenogenetic microcrustacean that has served as an ecological model for over a century.

Genome-Wide Estimation of Linkage Disequilibrium from Population-Level High-Throughput Sequencing Data

Abstract Rapidly improving sequencing technologies provide unprecedented opportunities for analyzing genome-wide patterns of polymorphisms. In particular, they have great potential for linkage-disequilibrium analyses on both global and local genetic scales, which will substantially improve our ability to derive evolutionary inferences. However, there are some difficulties with analyzing high-throughput sequencing data, including high error rates associated with base reads and complications from the random sampling of sequenced chromosomes in diploid organisms. To overcome these difficulties, we developed a maximum-likelihood estimator of linkage disequilibrium for use with error-prone sampling data. Computer simulations indicate that the estimator is nearly unbiased with a sampling variance at high coverage asymptotically approaching the value expected when all relevant information is accurately estimated. The estimator does not require phasing of haplotypes and enables the estimation of linkage disequilibrium even when all individual reads cover just single polymorphic sites.

Resource allocation to cell envelopes and the scaling of bacterial growth rate

Abstract Although various empirical studies have reported a positive correlation between the specific growth rate and cell size across bacteria, it is currently unclear what causes this relationship. We conjecture that such scaling occurs because smaller cells have a larger surface-to-volume ratio and thus have to allocate a greater fraction of the total resources to the production of the cell envelope, leaving fewer resources for other biosynthetic processes. To test this theory, we developed a coarse-grained model of bacterial physiology composed of the proteome that converts nutrients into biomass, with the cell envelope acting as a resource sink. Assuming resources are partitioned to maximize the growth rate, the model predicts that the growth rate and ribosomal mass fraction scale negatively, while the mass fraction of envelope-producing enzymes scales positively with surface-to-volume. These relationships are compatible with growth measurements and quantitative proteomics data reported in the literature.

Salmonella tel-el-kebir and terrapins

Differential retention and divergent resolution of duplicate genes following whole-genome duplication

The Paramecium aurelia complex is a group of 15 species that share at least three past whole-genome duplications (WGDs). The macronuclear genome sequences of P. biaurelia and P. sexaurelia are presented and compared to the published sequence of P. tetraurelia . Levels of duplicate-gene retention from the recent WGD differ by &gt;10% across species, with P. sexaurelia losing significantly more genes than P. biaurelia or P. tetraurelia . In addition, historically high rates of gene conversion have homogenized WGD paralogs, probably extending the paralogs’ lifetimes. The probability of duplicate retention is positively correlated with GC content and expression level; ribosomal proteins, transcription factors, and intracellular signaling proteins are overrepresented among maintained duplicates. Finally, multiple sources of evidence indicate that P. sexaurelia diverged from the two other lineages immediately following, or perhaps concurrent with, the recent WGD, with approximately half of gene losses between P. tetraurelia and P. sexaurelia representing divergent gene resolutions (i.e., silencing of alternative paralogs), as expected for random duplicate loss between these species. Additionally, though P. biaurelia and P. tetraurelia diverged from each other much later, there are still more than 100 cases of divergent resolution between these two species. Taken together, these results indicate that divergent resolution of duplicate genes between lineages acts to reinforce reproductive isolation between species in the Paramecium aurelia complex.

Low Base-Substitution Mutation Rate but High Rate of Slippage Mutations in the Sequence Repeat-Rich Genome of <i>Dictyostelium discoideum</i>

Abstract We describe the rate and spectrum of spontaneous mutations for the social amoeba Dictyostelium discoideum, a key model organism in molecular, cellular, evolutionary and developmental biology. Whole-genome sequencing of 37 mutation accumulation lines of D. discoideum after an average of 1,500 cell divisions yields a base-substitution mutation rate of 2.47 × 10−11 per site per generation, substantially lower than that of most eukaryotic and prokaryotic organisms, and of the same order of magnitude as in the ciliates Paramecium tetraurelia and Tetrahymena thermophila. Known for its high genomic AT content and abundance of simple sequence repeats, we observe that base-substitution mutations in D. discoideum are highly A/T biased. This bias likely contributes both to the high genomic AT content and to the formation of simple sequence repeats in the AT-rich genome of Dictyostelium discoideum. In contrast to the situation in other surveyed unicellular eukaryotes, indel rates far exceed the base-substitution mutation rate in this organism with a high proportion of 3n indels, particularly in regions without simple sequence repeats. Like ciliates, D. discoideum has a large effective population size, reducing the power of random genetic drift, magnifying the effect of selection on replication fidelity, in principle allowing D. discoideum to evolve an extremely low base-substitution mutation rate.

Simple sequence repeat variation in the Daphnia pulex genome

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.