Gavin A. Huttley

QIIME allows analysis of high-throughput community sequencing data

Supplementary Figure 1 Overview of the analysis pipeline. Supplementary Table 1 Details of conventionally raised and conventionalized mouse samples. Supplementary Discussion Expanded discussion of QIIME analyses presented in the main text; Sequencing of 16S rRNA gene amplicons; QIIME analysis notes; Expanded Figure 1 legend; Links to raw data and processed output from the runs with and without denoising.

Genetic Restriction of HIV-1 Infection and Progression to AIDS by a Deletion Allele of the CKR5 Structural Gene

The chemokine receptor 5 (CKR5) protein serves as a secondary receptor on CD4+ T lymphocytes for certain strains of human immunodeficiency virus-type 1 (HIV-1). The CKR5 structural gene was mapped to human chromosome 3p21, and a 32-base pair deletion allele (CKR5Δ32) was identified that is present at a frequency of ∼0.10 in the Caucasian population of the United States. An examination of 1955 patients included among six well-characterized acquired immunodeficiency syndrome (AIDS) cohort studies revealed that 17 deletion homozygotes occurred exclusively among 612 exposed HIV-1 antibody-negative individuals (2.8 percent) and not at all in 1343 HIV-1-infected individuals. The frequency of CKR5 deletion heterozygotes was significantly elevated in groups of individuals that had survived HIV-1 infection for more than 10 years, and, in some risk groups, twice as frequent as their occurrence in rapid progressors to AIDS. Survival analysis clearly shows that disease progression is slower in CKR5 deletion heterozygotes than in individuals homozygous for the normal CKR5 gene. The CKR5Δ32 deletion may act as a recessive restriction gene against HIV-1 infection and may exert a dominant phenotype of delaying progression to AIDS among infected individuals.

Dating the Origin of the CCR5-Δ32 AIDS-Resistance Allele by the Coalescence of Haplotypes

The CCR5-Delta32 deletion obliterates the CCR5 chemokine and the human immunodeficiency virus (HIV)-1 coreceptor on lymphoid cells, leading to strong resistance against HIV-1 infection and AIDS. A genotype survey of 4,166 individuals revealed a cline of CCR5-Delta32 allele frequencies of 0%-14% across Eurasia, whereas the variant is absent among native African, American Indian, and East Asian ethnic groups. Haplotype analysis of 192 Caucasian chromosomes revealed strong linkage disequilibrium between CCR5 and two microsatellite loci. By use of coalescence theory to interpret modern haplotype genealogy, we estimate the origin of the CCR5-Delta32-containing ancestral haplotype to be approximately 700 years ago, with an estimated range of 275-1,875 years. The geographic cline of CCR5-Delta32 frequencies and its recent emergence are consistent with a historic strong selective event (e.g. , an epidemic of a pathogen that, like HIV-1, utilizes CCR5), driving its frequency upward in ancestral Caucasian populations.

PyEvolve: a toolkit for statistical modelling of molecular evolution.

BackgroundExamining the distribution of variation has proven an extremely profitable technique in the effort to identify sequences of biological significance. Most approaches in the field, however, evaluate only the conserved portions of sequences – ignoring the biological significance of sequence differences. A suite of sophisticated likelihood based statistical models from the field of molecular evolution provides the basis for extracting the information from the full distribution of sequence variation. The number of different problems to which phylogeny-based maximum likelihood calculations can be applied is extensive. Available software packages that can perform likelihood calculations suffer from a lack of flexibility and scalability, or employ error-prone approaches to model parameterisation.ResultsHere we describe the implementation of PyEvolve, a toolkit for the application of existing, and development of new, statistical methods for molecular evolution. We present the object architecture and design schema of PyEvolve, which includes an adaptable multi-level parallelisation schema. The approach for defining new methods is illustrated by implementing a novel dinucleotide model of substitution that includes a parameter for mutation of methylated CpGs, which required 8 lines of standard Python code to define. Benchmarking was performed using either a dinucleotide or codon substitution model applied to an alignment of BRCA1 sequences from 20 mammals, or a 10 species subset. Up to five-fold parallel performance gains over serial were recorded. Compared to leading alternative software, PyEvolve exhibited significantly better real world performance for parameter rich models with a large data set, reducing the time required for optimisation from ~10 days to ~6 hours.ConclusionPyEvolve provides flexible functionality that can be used either for statistical modelling of molecular evolution, or the development of new methods in the field. The toolkit can be used interactively or by writing and executing scripts. The toolkit uses efficient processes for specifying the parameterisation of statistical models, and implements numerous optimisations that make highly parameter rich likelihood functions solvable within hours on multi-cpu hardware. PyEvolve can be readily adapted in response to changing computational demands and hardware configurations to maximise performance. PyEvolve is released under the GPL and can be downloaded from http://cbis.anu.edu.au/software.

Expansion of circulating T cells resembling follicular helper T cells is a fixed phenotype that identifies a subset of severe systemic lupus erythematosus.

OBJECTIVE In the sanroque mouse model of lupus, pathologic germinal centers (GCs) arise due to increased numbers of follicular helper T (Tfh) cells, resulting in high-affinity anti-double-stranded DNA antibodies that cause end-organ inflammation, such as glomerulonephritis. The purpose of this study was to examine the hypothesis that this pathway could account for a subset of patients with systemic lupus erythematosus (SLE). METHODS An expansion of Tfh cells is a causal, and therefore consistent, component of the sanroque mouse phenotype. We validated the enumeration of circulating T cells resembling Tfh cells as a biomarker of this expansion in sanroque mice, and we performed a comprehensive comparison of the surface phenotype of circulating and tonsillar Tfh cells in humans. This circulating biomarker was enumerated in SLE patients (n = 46), Sjögrens syndrome patients (n = 17), and healthy controls (n = 48) and was correlated with disease activity and end-organ involvement. RESULTS In sanroque mice, circulating Tfh cells increased in proportion to their GC counterparts, making circulating Tfh cells a feasible human biomarker of this novel mechanism of breakdown in GC tolerance. In a subset of SLE patients (14 of 46), but in none of the controls, the levels of circulating Tfh cells (defined as circulating CXCR5+CD4+ cells with high expression of Tfh-associated molecules, such as inducible T cell costimulator or programmed death 1) were increased. This cellular phenotype did not vary with time, disease activity, or treatment, but it did correlate with the diversity and titers of autoantibodies and with the severity of end-organ involvement. CONCLUSION These findings in SLE patients are consistent with the autoimmune mechanism in sanroque mice and identify Tfh effector molecules as possible therapeutic targets in a recognizable subset of patients with SLE.

Loss of ACTN3 gene function alters mouse muscle metabolism and shows evidence of positive selection in humans

More than a billion humans worldwide are predicted to be completely deficient in the fast skeletal muscle fiber protein α-actinin-3 owing to homozygosity for a premature stop codon polymorphism, R577X, in the ACTN3 gene. The R577X polymorphism is associated with elite athlete status and human muscle performance, suggesting that α-actinin-3 deficiency influences the function of fast muscle fibers. Here we show that loss of α-actinin-3 expression in a knockout mouse model results in a shift in muscle metabolism toward the more efficient aerobic pathway and an increase in intrinsic endurance performance. In addition, we demonstrate that the genomic region surrounding the 577X null allele shows low levels of genetic variation and recombination in individuals of European and East Asian descent, consistent with strong, recent positive selection. We propose that the 577X allele has been positively selected in some human populations owing to its effect on skeletal muscle metabolism.

PyCogent: a toolkit for making sense from sequence

We have implemented in Python the COmparative GENomic Toolkit, a fully integrated and thoroughly tested framework for novel probabilistic analyses of biological sequences, devising workflows, and generating publication quality graphics. PyCogent includes connectors to remote databases, built-in generalized probabilistic techniques for working with biological sequences, and controllers for third-party applications. The toolkit takes advantage of parallel architectures and runs on a range of hardware and operating systems, and is available under the general public license from http://sourceforge.net/projects/pycogent.

Adaptive evolution of the tumour suppressor BRCA1 in humans and chimpanzees

Mutations in BRCA1 (ref. 1) confer an increased risk of female breast cancer. In a genome-wide scan of linkage disequilibrium (LD), a high level of LD was detected among microsatellite markers flanking BRCA1 (ref. 3), raising the prospect that positive natural selection may have acted on this gene. We have used the predictions of evolutionary genetic theory to investigate this further. Using phylogeny-based maximum likelihood analysis of the BRCA1 sequences from primates and other mammals, we found that the ratios of replacement to silent nucleotide substitutions on the human and chimpanzee lineages were not different from one another (P=0.8), were different from those of other primate lineages (P=0.004) and were greater than 1 (P=0.04). This is consistent with the historic occurrence of positive darwinian selection pressure on the BRCA1 protein in the human and chimpanzee lineages. Analysis of genetic variation in a sample of female Australians of Northern European origin showed evidence for Hardy-Weinberg (HW) disequilibrium at polymorphic sites in BRCA1, consistent with the possibility that natural selection is affecting genotype frequencies in modern Europeans. The clustering of between-species variation in the region of the gene encoding the RAD51-interaction domain of BRCA1 suggests the maintenance of genomic integrity as a possible target of selection.

Cryptococcus gattii in North American Pacific Northwest: Whole-Population Genome Analysis Provides Insights into Species Evolution and Dispersal

ABSTRACT The emergence of distinct populations of Cryptococcus gattii in the temperate North American Pacific Northwest (PNW) was surprising, as this species was previously thought to be confined to tropical and semitropical regions. Beyond a new habitat niche, the dominant emergent population displayed increased virulence and caused primary pulmonary disease, as opposed to the predominantly neurologic disease seen previously elsewhere. Whole-genome sequencing was performed on 118 C. gattii isolates, including the PNW subtypes and the global diversity of molecular type VGII, to better ascertain the natural source and genomic adaptations leading to the emergence of infection in the PNW. Overall, the VGII population was highly diverse, demonstrating large numbers of mutational and recombinational events; however, the three dominant subtypes from the PNW were of low diversity and were completely clonal. Although strains of VGII were found on at least five continents, all genetic subpopulations were represented or were most closely related to strains from South America. The phylogenetic data are consistent with multiple dispersal events from South America to North America and elsewhere. Numerous gene content differences were identified between the emergent clones and other VGII lineages, including genes potentially related to habitat adaptation, virulence, and pathology. Evidence was also found for possible gene introgression from Cryptococcus neoformans var. grubii that is rarely seen in global C. gattii but that was present in all PNW populations. These findings provide greater understanding of C. gattii evolution in North America and support extensive evolution in, and dispersal from, South America. IMPORTANCE Cryptococcus gattii emerged in the temperate North American Pacific Northwest (PNW) in the late 1990s. Beyond a new environmental niche, these emergent populations displayed increased virulence and resulted in a different pattern of clinical disease. In particular, severe pulmonary infections predominated in contrast to presentation with neurologic disease as seen previously elsewhere. We employed population-level whole-genome sequencing and analysis to explore the genetic relationships and gene content of the PNW C. gattii populations. We provide evidence that the PNW strains originated from South America and identified numerous genes potentially related to habitat adaptation, virulence expression, and clinical presentation. Characterization of these genetic features may lead to improved diagnostics and therapies for such fungal infections. The data indicate that there were multiple recent introductions of C. gattii into the PNW. Public health vigilance is warranted for emergence in regions where C. gattii is not thought to be endemic. Cryptococcus gattii emerged in the temperate North American Pacific Northwest (PNW) in the late 1990s. Beyond a new environmental niche, these emergent populations displayed increased virulence and resulted in a different pattern of clinical disease. In particular, severe pulmonary infections predominated in contrast to presentation with neurologic disease as seen previously elsewhere. We employed population-level whole-genome sequencing and analysis to explore the genetic relationships and gene content of the PNW C. gattii populations. We provide evidence that the PNW strains originated from South America and identified numerous genes potentially related to habitat adaptation, virulence expression, and clinical presentation. Characterization of these genetic features may lead to improved diagnostics and therapies for such fungal infections. The data indicate that there were multiple recent introductions of C. gattii into the PNW. Public health vigilance is warranted for emergence in regions where C. gattii is not thought to be endemic.

Modelling and bioinformatics studies of the human Kappa-class glutathione transferase predict a novel third glutathione transferase family with similarity to prokaryotic 2-hydroxychromene-2-carboxylate isomerases

The Kappa class of GSTs (glutathione transferases) comprises soluble enzymes originally isolated from the mitochondrial matrix of rats. We have characterized a Kappa class cDNA from human breast. The cDNA is derived from a single gene comprising eight exons and seven introns located on chromosome 7q34-35. Recombinant hGSTK1-1 was expressed in Escherichia coli as a homodimer (subunit molecular mass approximately 25.5 kDa). Significant glutathione-conjugating activity was found only with the model substrate CDNB (1-chloro-2,4-ditnitrobenzene). Hyperbolic kinetics were obtained for GSH (parameters: K(m)app, 3.3+/-0.95 mM; V(max)app, 21.4+/-1.8 micromol/min per mg of enzyme), while sigmoidal kinetics were obtained for CDNB (parameters: S0.5app, 1.5+/-1.0 mM; V(max)app, 40.3+/-0.3 micromol/min per mg of enzyme; Hill coefficient, 1.3), reflecting low affinities for both substrates. Sequence analyses, homology modelling and secondary structure predictions show that hGSTK1 has (a) most similarity to bacterial HCCA (2-hydroxychromene-2-carboxylate) isomerases and (b) a predicted C-terminal domain structure that is almost identical to that of bacterial disulphide-bond-forming DsbA oxidoreductase (root mean square deviation 0.5-0.6 A). The structures of hGSTK1 and HCCA isomerase are predicted to possess a thioredoxin fold with a polyhelical domain (alpha(x)) embedded between the beta-strands (betaalphabetaalpha(x)betabetaalpha, where the underlined elements represent the N and C motifs of the thioredoxin fold), as occurs in the bacterial disulphide-bond-forming oxidoreductases. This is in contrast with the cytosolic GSTs, where the helical domain occurs exclusively at the C-terminus (betaalphabetaalphabetabetaalphaalpha(x)). Although hGSTK1-1 catalyses some typical GST reactions, we propose that it is structurally distinct from other classes of cytosolic GSTs. The present study suggests that the Kappa class may have arisen in prokaryotes well before the divergence of the cytosolic GSTs.