Mixture models of nucleotide sequence evolution that account for heterogeneity in the substitution process across sites and across lineages

Molecular phylogenetic studies of homologous sequences of nucleotides often assume that the underlying evolutionary process was globally stationary, reversible, and homogeneous (SRH), and that a model of evolution with one or more site-specific and time-reversible rate matrices (e.g., the GTR rate matrix) is enough to accurately model the evolution of data over the whole tree. However, an increasing body of data suggests that evolution under these conditions is an exception, rather than the norm. To address this issue, several non-SRH models of molecular evolution have been proposed, but they either ignore heterogeneity in the substitution process across sites (HAS) or assume it can be modeled accurately using the distribution. As an alternative to these models of evolution, we introduce a family of mixture models that approximate HAS without the assumption of an underlying predefined statistical distribution. This family of mixture models is combined with non-SRH models of evolution that account for heterogeneity in the substitution process across lineages (HAL). We also present two algorithms for searching model space and identifying an optimal model of evolution that is less likely to over- or underparameterize the data. The performance of the two new algorithms was evaluated using alignments of nucleotides with 10 000 sites simulated under complex non-SRH conditions on a 25-tipped tree. The algorithms were found to be very successful, identifying the correct HAL model with a 75% success rate (the average success rate for assigning rate matrices to the tree's 48 edges was 99.25%) and, for the correct HAL model, identifying the correct HAS model with a 98% success rate. Finally, parameter estimates obtained under the correct HAL-HAS model were found to be accurate and precise. The merits of our new algorithms were illustrated with an analysis of 42 337 second codon sites extracted from a concatenation of 106 alignments of orthologous genes encoded by the nuclear genomes of Saccharomyces cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, S. castellii, S. kluyveri, S. bayanus, and Candida albicans. Our results show that second codon sites in the ancestral genome of these species contained 49.1% invariable sites, 39.6% variable sites belonging to one rate category (V1), and 11.3% variable sites belonging to a second rate category (V2). The ancestral nucleotide content was found to differ markedly across these three sets of sites, and the evolutionary processes operating at the variable sites were found to be non-SRH and best modeled by a combination of eight edge-specific rate matrices (four for V1 and four for V2). The number of substitutions per site at the variable sites also differed markedly, with sites belonging to V1 evolving slower than those belonging to V2 along the lineages separating the seven species of Saccharomyces. Finally, sites belonging to V1 appeared to have ceased evolving along the lineages separating S. cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, and S. bayanus, implying that they might have become so selectively constrained that they could be considered invariable sites in these species.

Molecular and pathological insights into Chlamydia pecorum-associated sporadic bovine encephalomyelitis (SBE) in Western Australia

Background Despite its global recognition as a ruminant pathogen, cases of Chlamydia pecorum infection in Australian livestock are poorly documented. In this report, a C. pecorum specific Multi Locus Sequence Analysis scheme was used to characterise the C. pecorum strains implicated in two cases of sporadic bovine encephalomyelitis confirmed by necropsy, histopathology and immunohistochemistry. This report provides the first molecular evidence for the presence of mixed infections of C. pecorum strains in Australian cattle. Case presentation Affected animals were two markedly depressed, dehydrated and blind calves, 12 and 16 weeks old. The calves were euthanized and necropsied. In one calf, a severe fibrinous polyserositis was noted with excess joint fluid in all joints whereas in the other, no significant lesions were seen. No gross abnormalities were noted in the brain of either calf. Histopathological lesions seen in both calves included: multifocal, severe, subacute meningoencephalitis with vasculitis, fibrinocellular thrombosis and malacia; diffuse, mild, acute interstitial pneumonia; and diffuse, subacute epicarditis, severe in the calf with gross serositis. Immunohistochemical labelling of chlamydial antigen in brain, spleen and lung from the two affected calves and brain from two archived cases, localised the antigen to the cytoplasm of endothelium, mesothelium and macrophages. C. pecorum specific qPCR, showed dissemination of the pathogen to multiple organs. Phylogenetic comparisons with other C. pecorum bovine strains from Australia, Europe and the USA revealed the presence of two genetically distinct sequence types (ST). The predominant ST detected in the brain, heart, lung and liver of both calves was identical to the C. pecorum ST previously described in cases of SBE. A second ST detected in an ileal tissue sample from one of the calves, clustered with previously typed faecal bovine isolates. Conclusion This report provides the first data to suggest that identical C. pecorum STs may be associated with SBE in geographically separated countries and that these may be distinct from those found in the gastrointestinal tract. This report provides a platform for further investigations into SBE and for understanding the genetic relationships that exist between C. pecorum strains detected in association with other infectious diseases in livestock.

Evaluation of the relationship between Chlamydia pecorum sequence types and disease using a species-specific multi-locus sequence typing scheme (MLST)

Chlamydia pecorum is globally associated with several ovine diseases including keratoconjunctivitis and polyarthritis. The exact relationship between the variety of C. pecorum strains reported and the diseases described in sheep remains unclear, challenging efforts to accurately diagnose and manage infected flocks. In the present study, we applied C. pecorum multi-locus sequence typing (MLST) to C. pecorum positive samples collected from sympatric flocks of Australian sheep presenting with conjunctivitis, conjunctivitis with polyarthritis, or polyarthritis only and with no clinical disease (NCD) in order to elucidate the exact relationships between the infecting strains and the range of diseases. Using Bayesian phylogenetic and cluster analyses on 62 C. pecorum positive ocular, vaginal and rectal swab samples from sheep presenting with a range of diseases and in a comparison to C. pecorum sequence types (STs) from other hosts, one ST (ST 23) was recognised as a globally distributed strain associated with ovine and bovine diseases such as polyarthritis and encephalomyelitis. A second ST (ST 69) presently only described in Australian animals, was detected in association with ovine as well as koala chlamydial infections. The majority of vaginal and rectal C. pecorum STs from animals with NCD and/or anatomical sites with no clinical signs of disease in diseased animals, clustered together in a separate group, by both analyses. Furthermore, 8/13 detected STs were novel. This study provides a platform for strain selection for further research into the pathogenic potential of C. pecorum in animals and highlights targets for potential strain-specific diagnostic test development.

Structure and vascular tissue expression of duplicated TERMINAL EAR1-like paralogues in poplar

TERMINAL EAR1-like (TEL) genes encode putative RNA-binding proteins only found in land plants. Previous studies suggested that they may regulate tissue and organ initiation in Poaceae. Two TEL genes were identified in both Populus trichocarpa and the hybrid aspen Populus tremula × P. alba, named, respectively, PoptrTEL1-2 and PtaTEL1-2. The analysis of the organisation around the PoptrTEL genes in the P. trichocarpa genome and the estimation of the synonymous substitution rate for PtaTEL1-2 genes indicate that the paralogous link between these two Populus TEL genes probably results from the Salicoid large-scale gene-duplication event. Phylogenetic analyses confirmed their orthology link with the other TEL genes. The expression pattern of both PtaTEL genes appeared to be restricted to the mother cells of the plant body: leaf founder cells, leaf primordia, axillary buds and root differentiating tissues, as well as to mother cells of vascular tissues. Most interestingly, PtaTEL1-2 transcripts were found in differentiating cells of secondary xylem and phloem, but probably not in the cambium itself. Taken together, these results indicate specific expression of the TEL genes in differentiating cells controlling tissue and organ development in Populus (and other Angiosperm species).

Comparative genomic analysis of human Chlamydia pneumoniae isolates from respiratory, brain and cardiac tissues

Chlamydia pneumoniae is an obligate intracellular bacterium implicated in a wide range of human diseases including atherosclerosis and Alzheimer's disease. Efforts to understand the relationships between C. pneumoniae detected in these diseases have been hindered by the availability of sequence data for non-respiratory strains. In this study, we sequenced the whole genomes for C. pneumoniae isolates from atherosclerosis and Alzheimer's disease, and compared these to previously published C. pneumoniae genomes. Phylogenetic analyses of these new C. pneumoniae strains indicate two sub-groups within human C. pneumoniae, and suggest that both recombination and mutation events have driven the evolution of human C. pneumoniae. Further fine-detailed analyses of these new C. pneumoniae sequences show several genetically variable loci. This suggests that similar strains of C. pneumoniae are found in the brain, lungs and cardiovascular system and that only minor genetic differences may contribute to the adaptation of particular strains in human disease.

Nature and extent of genetic diversity of dengue viruses determined by 454 pyrosequencing

Dengue virus (DENV) populations are characteristically highly diverse. Regular lineage extinction and replacement is an important dynamic DENV feature, and most DENV lineage turnover events are associated with increased incidence of disease. The role of genetic diversity in DENV lineage extinctions is not understood. We investigated the nature and extent of genetic diversity in the envelope (E) gene of DENV serotype 1 representing different lineages histories. A region of the DENV genome spanning the E gene was amplified and sequenced by Roche/454 pyrosequencing. The pyrosequencing results identified distinct sub-populations (haplotypes) for each DENV-1 E gene. A phylogenetic tree was constructed with the consensus DENV-1 E gene nucleotide sequences, and the sequences of each constructed haplotype showed that the haplotypes segregated with the Sanger consensus sequence of the population from which they were drawn. Haplotypes determined through pyrosequencing identified a recombinant DENV genome that could not be identified through Sanger sequencing. Nucleotide level sequence diversities of DENV-1 populations determined from SNP analysis were very low, estimated from 0.009-0.01. There were also no stop codon, frameshift or non-frameshift mutations observed in the E genes of any lineage. No significant correlations between the accumulation of deleterious mutations or increasing genetic diversity and lineage extinction were observed (p>0.5). Although our hypothesis that accumulation of deleterious mutations over time led to the extinction and replacement of DENV lineages was ultimately not supported by the data, our data does highlight the significant technical issues that must be resolved in the way in which population diversity is measured for DENV and other viruses. The results provide an insight into the within-population genetic structure and diversity of DENV-1 populations.

The complexities of knowing what it is you are trapping

The effectiveness of any trapping system is highly dependent on the ability to accurately identify the specimens collected. For many fruit fly species, accurate identification (= diagnostics) using morphological or molecular techniques is relatively straightforward and poses few technical challenges. However, nearly all genera of pest tephritids also contain groups of species where single, stand-alone tools are not sufficient for accurate identification: such groups include the Bactrocera dorsalis complex, the Anastrepha fraterculus complex and the Ceratitis FAR complex. Misidentification of high-impact species from such groups can have dramatic consequences and negate the benefits of an otherwise effective trapping program. To help prevent such problems, this chapter defines what is meant by a species complex and describes in detail how the correct identification of species within a complex requires the use of an integrative taxonomic approach. Integrative taxonomy uses multiple, independent lines of evidence to delimit species boundaries, and the underpinnings of this approach from both the theoretical speciation literature and the systematics/taxonomy literature are described. The strength of the integrative approach lies in the explicit testing of hypotheses and the use of multiple, independent species delimitation tools. A case is made for a core set of species delimitation tools (pre- and post-zygotic compatibility tests, multi-locus phylogenetic analysis, chemoecological studies, and morphometric and geometric morphometric analyses) to be adopted as standards by tephritologists aiming to resolve economically important species complexes. In discussing the integrative approach, emphasis is placed on the subtle but important differences between integrative and iterative taxonomy. The chapter finishes with a case study that illustrates how iterative taxonomy applied to the B. dorsalis species complex led to incorrect taxonomic conclusions, which has had major implications for quarantine, trade, and horticultural pest management. In contrast, an integrative approach to the problem has resolved species limits in this taxonomically difficult group, meaning that robust diagnostics are now available.

The dynamics of individual teacher autonomy: A case study of teachers' work in a Queensland secondary school

The emphasis on collegiality and collaboration in the literature on teachers' work and school reform has tended to underplay the significance of teacher autonomy. This thesis explores the dynamics of teachers' understandings and experiences of individual teacher autonomy (as contrasted with collective autonomy) in an independent school in Queensland which promoted itself as a 'teachers' school' with a strong commitment to individual teacher autonomy. The research was a case study which drew on methodological signposts from critical, feminist and traditional ethnography. Intensive fieldwork in the school over five months incorporated the ethnographic techniques of observation, interviews and document analysis. Teachers at Thornton College understood their experience of individual autonomy at three interrelated levels--in terms of their work in the classroom, their working life in the school, and their voice in the decision-making processes of the school. They felt that they experienced a great deal of individual autonomy at each of these three levels. These understandings and experiences of autonomy were encumbered or enabled by a range of internal and external stakeholder groups. There were also a number of structural influences (community perceptions, market forces, school size, time and bureaucracy) emerging from the economic, social and political structures in Australian society which influenced the experience of autonomy by teachers. The experience of individual teacher autonomy was constantly shifting, but there were some emergent patterns. Consensus on educational goals and vision, and strong expressions of trust and respect between teachers and stakeholders in the school, characterised the contexts in which teachers felt they experienced high levels of autonomy in their work. The demand for accountability and desire for relatedness motivated stakeholders and structural forces to influence teacher autonomy. Some significant gaps emerged between the rhetoric of a commitment to individual teacher autonomy and decision-making practices in the school, that gave ultimate power to the co-principals. Despite the rhetoric and promotion of non-hierarchical structures and collaborative decision-making processes, many teachers perceived that their experience of individual autonomy remained subject to the exercise of 'partial democracy' by school leaders.

Playing alone, playing with others: Differences in player experience and indicators of wellbeing

Video game play is becoming an increasingly social experience, yet we have little understanding of how social and solitary modes of play differ in terms of the player experience or interact with player wellbeing. An online survey (n = 446) collected data on players' current mode of play, their game play experience, social capital gained from game play and wellbeing. The results indicate that social and solitary players differ in terms of degree of autonomy, presence and relatedness experienced, while the different types of social play are associated with differences in relatedness and social capital experienced. Different predictors of wellbeing were also present across solitary and social player samples. People who play games on their own experience greater wellbeing when experiencing in-game autonomy. Social players experience greater wellbeing when playing with strangers, and when experiencing in-game bridging social capital. All players experienced increased wellbeing with age and decreased wellbeing with greater amounts of play.

Indicators of wellbeing in recreational video game players

Video game play is a popular entertainment choice, yet we have a limited understanding of the potential wellbeing benefits associated with recreational play. An online survey (final sample, n = 297) addresses this by investigating how the player experience related to wellbeing. The impact of amount of play, game genre, mode of play (social or solitary play) and the psychological experience of play (flow and need satisfaction) on a multi-dimensional measure of wellbeing (emotional, psychological and social) was examined via hierarchical regression. Age, gender, the play of casual games compared to shooters, and in-game experiences of flow, autonomy and relatedness were associated with increases in dimensions of wellbeing.

Genomic inflation factors under polygenic inheritance

Population structure, including population stratification and cryptic relatedness, can cause spurious associations in genome-wide association studies (GWAS). Usually, the scaled median or mean test statistic for association calculated from multiple single-nucleotide-polymorphisms across the genome is used to assess such effects, and 'genomic control' can be applied subsequently to adjust test statistics at individual loci by a genomic inflation factor. Published GWAS have clearly shown that there are many loci underlying genetic variation for a wide range of complex diseases and traits, implying that a substantial proportion of the genome should show inflation of the test statistic. Here, we show by theory, simulation and analysis of data that in the absence of population structure and other technical artefacts, but in the presence of polygenic inheritance, substantial genomic inflation is expected. Its magnitude depends on sample size, heritability, linkage disequilibrium structure and the number of causal variants. Our predictions are consistent with empirical observations on height in independent samples of ~4000 and ~133,000 individuals.

Evidence from molecules and morphology expands Podonomopsis Brundin (Diptera: Chironomidae: Podonominae) to include ‘genus Chile'

The informal taxon ‘genus Chile’ of Brundin, based solely on pupal exuviae of a podonomine Chironomidae, has remained inadequately known for half a century. New collections reveal life associations, and provide molecular data to hypothesise a precise phylogenetic placement in the austral Podonominae. A densely sampled molecular phylogeny based on two nuclear and one mitochondrial DNA markers shows ‘genus Chile’ to be the sister group to Podonomopsis Brundin, 1966. Within Podonomopsis a clade of South American species is sister to all Australian species. We discuss how to rank such a sister group taxon and treat ‘genus Chile’ as a new subgenus Araucanopsis, subg. nov. with the new species, Podonomopsis (Araucanopsis) avelasse, sp. nov. from Chile and Argentina as genotype of the monotypic subgenus. We describe P. (A.) avelasse in all stages and provide an expanded diagnosis and description of Podonomopsis to include Araucanopsis. A dated biogeographic hypothesis (chronogram) infers the most recent common ancestor (tmcra) of expanded Podonomopsis at 95 million years ago (Mya) (68–122 Mya 95% highest posterior density), ‘core’ Podonomopsis at 83 Mya (58–108) and Australian Podonomopsis at 65 Mya (44–87). All dates are before the South America–Australia geological separation through Antarctica, supporting previous conclusions that the taxon distribution is ‘Gondwanan’ in origin. Podonomopsis, even as expanded here, remains unknown from New Zealand or elsewhere on extant Zealandia.

Characterization of sequence and structural features of the Candida krusei Enolase

The incidence of human infections by the fungal pathogen Candida species has been increasing in recent years. Enolase is an essential protein in fungal metabolism. Sequence data is available for human and a number of medically important fungal species. An understanding of the structural and functional features of fungal enolases may provide the structural basis for their use as a target for the development of new anti-fungal drugs. We have obtained the sequence of the enolase of Candida krusei (C. krusei), as it is a significant medically important fungal pathogen. We have then used multiple sequence alignments with various enolase isoforms in order to identify C. krusei specific amino acid residues. The phylogenetic tree of enolases shows that the C. krusei enolase assembles on the tree with the fungal genes. Importantly, C. krusei lacks four amino acids in the active site compared to human enolase, as revealed by multiple sequence alignments. These differences in the substrate binding site may be exploited for the design of new anti-fungal drugs to selectively block this enzyme. The lack of the important amino acids in the active site also indicates that C. krusei enolase might have evolved as a member of a mechanistically diverse enolase superfamily catalying somewhat different reactions.

Directional dominance on stature and cognition in diverse human populations

Homozygosity has long been associated with rare, often devastating, Mendelian disorders1, and Darwin was one of the first to recognize that inbreeding reduces evolutionary fitness2. However, the effect of the more distant parental relatedness that is common in modern human populations is less well understood. Genomic data now allow us to investigate the effects of homozygosity on traits of public health importance by observing contiguous homozygous segments (runs of homozygosity), which are inferred to be homozygous along their complete length. Given the low levels of genome-wide homozygosity prevalent in most human populations, information is required on very large numbers of people to provide sufficient power3, 4. Here we use runs of homozygosity to study 16 health-related quantitative traits in 354,224 individuals from 102 cohorts, and find statistically significant associations between summed runs of homozygosity and four complex traits: height, forced expiratory lung volume in one second, general cognitive ability and educational attainment (P < 1 × 10−300, 2.1 × 10−6, 2.5 × 10−10 and 1.8 × 10−10, respectively). In each case, increased homozygosity was associated with decreased trait value, equivalent to the offspring of first cousins being 1.2 cm shorter and having 10 months’ less education. Similar effect sizes were found across four continental groups and populations with different degrees of genome-wide homozygosity, providing evidence that homozygosity, rather than confounding, directly contributes to phenotypic variance. Contrary to earlier reports in substantially smaller samples5, 6, no evidence was seen of an influence of genome-wide homozygosity on blood pressure and low density lipoprotein cholesterol, or ten other cardio-metabolic traits. Since directional dominance is predicted for traits under directional evolutionary selection7, this study provides evidence that increased stature and cognitive function have been positively selected in human evolution, whereas many important risk factors for late-onset complex diseases may not have been.

Cooperative game play with avatars and agents: Differences in brain activity and the experience of play

The current study sought to identify the impact of whether teammates in a cooperative videogame were controlled by other humans (avatars) or by the game (agents). The impact on player experience was explored through both subjective questionnaire measures and brain wave activity measurement (electroencephalography). Play with human teammates was associated with a greater sense of relatedness, but less competence and flow than play with other computer-controlled teammates. In terms of brain activity, play with human teammates was associated with greater activity in the alpha, theta and beta power bands than play with computer-controlled teammates. Overall, the results suggest that play with human teammates involves greater cognitive activity in terms of 'mentalising' than play with computer-controlled teammates. Additionally, the associations between subjective measures of player experience and brain activity are described. Limitations of the current study are identified and key directions for future research are discussed.