Gene expression differences among three Neurospora species reveal genes required for sexual reproduction in Neurospora crassa

Many fungi form complex three-dimensional fruiting bodies, within which the meiotic machinery for sexual spore production has been considered to be largely conserved over evolutionary time. Indeed, much of what we know about meiosis in plant and animal taxa has been deeply informed by studies of meiosis in Saccharomyces and Neurospora. Nevertheless, the genetic basis of fruiting body development and its regulation in relation to meiosis in fungi is barely known, even within the best studied multicellular fungal model Neurospora crassa. We characterized morphological development and genome-wide transcriptomics in the closely related species Neurospora crassa, Neurospora tetrasperma, and Neurospora discreta, across eight stages of sexual development. Despite diverse life histories within the genus, all three species produce vase-shaped perithecia. Transcriptome sequencing provided gene expression levels of orthologous genes among all three species. Expression of key meiosis genes and sporulation genes corresponded to known phenotypic and developmental differences among these Neurospora species during sexual development. We assembled a list of genes putatively relevant to the recent evolution of fruiting body development by sorting genes whose relative expression across developmental stages increased more in N. crassa relative to the other species. Then, in N. crassa, we characterized the phenotypes of fruiting bodies arising from crosses of homozygous knockout strains of the top genes. Eight N. crassa genes were found to be critical for the successful formation of perithecia. The absence of these genes in these crosses resulted in either no perithecium formation or in arrested development at an early stage. Our results provide insight into the genetic basis of Neurospora sexual reproduction, which is also of great importance with regard to other multicellular ascomycetes, including perithecium-forming pathogens, such as Claviceps purpurea, Ophiostoma ulmi, and Glomerella graminicola.

Phylogenetic informativeness reconciles ray-finned fish molecular divergence times

Background

Discordance among individual molecular age estimates, or between molecular age estimates and the fossil record, is observed in many clades across the Tree of Life. This discordance is attributed to a variety of variables including calibration age uncertainty, calibration placement, nucleotide substitution rate heterogeneity, or the specified molecular clock model. However, the impact of changes in phylogenetic informativeness of individual genes over time on phylogenetic inferences is rarely analyzed. Using nuclear and mitochondrial sequence data for ray-finned fishes (Actinopterygii) as an example, we extend the utility of phylogenetic informativeness profiles to predict the time intervals when nucleotide substitution saturation results in discordance among molecular ages estimated.

Results

We demonstrate that even with identical calibration regimes and molecular clock methods, mitochondrial based molecular age estimates are systematically older than those estimated from nuclear sequences. This discordance is most severe for highly nested nodes corresponding to more recent (i.e., Jurassic-Recent) divergences. By removing data deemed saturated, we reconcile the competing age estimates and highlight that the older mtDNA based ages were driven by nucleotide saturation.

Conclusions

Homoplasious site patterns in a DNA sequence alignment can systematically bias molecular divergence time estimates. Our study demonstrates that PI profiles can provide a non-arbitrary criterion for data exclusion to mitigate the influence of homoplasy on time calibrated branch length estimates. Analyses of actinopterygian molecular clocks demonstrate that scrutiny of the time scale on which sequence data is informative is a fundamental, but generally overlooked, step in molecular divergence time estimation.

Detecting rare gene transfer events in bacterial populations

Horizontal gene transfer (HGT) enables bacteria to access, share, and recombine genetic variation, resulting in genetic diversity that cannot be obtained through mutational processes alone. In most cases, the observation of evolutionary successful HGT events relies on the outcome of initially rare events that lead to novel functions in the new host, and that exhibit a positive effect on host fitness. Conversely, the large majority of HGT events occurring in bacterial populations will go undetected due to lack of replication success of transformants. Moreover, other HGT events that would be highly beneficial to new hosts can fail to ensue due to lack of physical proximity to the donor organism, lack of a suitable gene transfer mechanism, genetic compatibility, and stochasticity in tempo-spatial occurrence. Experimental attempts to detect HGT events in bacterial populations have typically focused on the transformed cells or their immediate offspring. However, rare HGT events occurring in large and structured populations are unlikely to reach relative population sizes that will allow their immediate identification; the exception being the unusually strong positive selection conferred by antibiotics. Most HGT events are not expected to alter the likelihood of host survival to such an extreme extent, and will confer only minor changes in host fitness. Due to the large population sizes of bacteria and the time scales involved, the process and outcome of HGT are often not amenable to experimental investigation. Population genetic modeling of the growth dynamics of bacteria with differing HGT rates and resulting fitness changes is therefore necessary to guide sampling design and predict realistic time frames for detection of HGT, as it occurs in laboratory or natural settings. Here we review the key population genetic parameters, consider their complexity and highlight knowledge gaps for further research.

Potential cost-effectiveness of schistosomiasis treatment for reducing HIV transmission in Africa--the case of Zimbabwean women

Background

Epidemiological data from Zimbabwe suggests that genital infection with Schistosoma haematobium may increase the risk of HIV infection in young women. Therefore, the treatment of Schistosoma haematobium with praziquantel could be a potential strategy for reducing HIV infection. Here we assess the potential cost-effectiveness of praziquantel as a novel intervention strategy against HIV infection.

Methods

We developed a mathematical model of female genital schistosomiasis (FGS) and HIV infections in Zimbabwe that we fitted to cross-sectional data of FGS and HIV prevalence of 1999. We validated our epidemic projections using antenatal clinic data on HIV prevalence. We simulated annual praziquantel administration to school-age children. We then used these model predictions to perform a cost-effectiveness analysis of annual administration of praziquantel as a potential measure to reduce the burden of HIV in sub-Saharan Africa.

Findings

We showed that for a variation of efficacy between 30–70% of mass praziquantel administration for reducing the enhanced risk of HIV transmission per sexual act due to FGS, annual administration of praziquantel to school-age children in Zimbabwe could result in net savings of US$16–101 million compared with no mass treatment of schistosomiasis over a ten-year period. For a variation in efficacy between 30–70% of mass praziquantel administration for reducing the acquisition of FGS, annual administration of praziquantel to school-age children could result in net savings of US$36−92 million over a ten-year period.

Conclusions

In addition to reducing schistosomiasis burden, mass praziquantel administration may be a highly cost-effective way of reducing HIV infections in sub-Saharan Africa. Program costs per case of HIV averted are similar to, and under some conditions much better than, other interventions that are currently implemented in Africa to reduce HIV transmission. As a cost-saving strategy, mass praziquantel administration should be prioritized over other less cost-effective public health interventions.

Evidence from epidemiological and clinical studies supports the hypothesis that genital infection with Schistosoma haematobium increases the risk of becoming infected with HIV among women in sub-Saharan Africa. Praziquantel is an oral, nontoxic, inexpensive medication recommended for treatment of schistosomiasis, which might be able to prevent the development of genital schistosomiasis. We constructed a mathematical model of female genital schistosomiasis and HIV infections, which we calibrated using epidemiological data from Zimbabwe. We used this model to investigate the potential cost-effectiveness of mass drug administration with praziquantel as an intervention strategy for reducing HIV transmission in sub-Saharan Africa. We showed that mass drug administration with praziquantel may be a timely, innovative, and cost-saving intervention strategy for HIV prevention in sub-Saharan Africa. As a cost-saving strategy, mass drug administration with praziquantel should be prioritized over other less cost-effective public health interventions. Our findings indicate the possible benefit of scaling up schistosomiasis control efforts in sub-Saharan Africa, and especially in areas were Schistosoma haematobium and HIV are highly prevalent.

Evaluating phylogenetic informativeness as a predictor of phylogenetic signal for metazoan, fungal, and mammalian phylogenomic data sets

Phylogenetic research is often stymied by selection of a marker that leads to poor phylogenetic resolution despite considerable cost and effort. Profiles of phylogenetic informativeness provide a quantitative measure for prioritizing gene sampling to resolve branching order in a particular epoch. To evaluate the utility of these profiles, we analyzed phylogenomic data sets from metazoans, fungi, and mammals, thus encompassing diverse time scales and taxonomic groups. We also evaluated the utility of profiles created based on simulated data sets. We found that genes selected via their informativeness dramatically outperformed haphazard sampling of markers. Furthermore, our analyses demonstrate that the original phylogenetic informativeness method can be extended to trees with more than four taxa. Thus, although the method currently predicts phylogenetic signal without specifically accounting for the misleading effects of stochastic noise, it is robust to the effects of homoplasy. The phylogenetic informativeness rankings obtained will allow other researchers to select advantageous genes for future studies within these clades, maximizing return on effort and investment. Genes identified might also yield efficient experimental designs for phylogenetic inference for many sister clades and outgroup taxa that are closely related to the diverse groups of organisms analyzed.

Sex-linked transcriptional divergence in the hermaphrodite fungus Neurospora tetrasperma

In the filamentous ascomycete Neurospora tetrasperma, a large (approx. 7 Mbp) region of suppressed recombination surrounds the mating-type (mat) locus. While the remainder of the genome is largely homoallelic, this region of recombinational suppression, extending over 1500 genes, is associated with sequence divergence. Here, we used microarrays to examine how the molecular phenotype of gene expression level is linked to this divergent region, and thus to the mating type. Culturing N. tetrasperma on agar media that induce sexual/female or vegetative/male tissue, we found 196 genes significantly differentially expressed between mat A and mat a mating types. Our data show that the genes exhibiting mat-linked expression are enriched in the region genetically linked to mating type, and sequence and expression divergence are positively correlated. Our results indicate that the phenotype of mat A strains is optimized for traits promoting sexual/female development and the phenotype of mat a strains for vegetative/male development. This discovery of differentially expressed genes associated with mating type provides a link between genotypic and phenotypic divergence in this taxon and illustrates a fungal analogue to sexual dimorphism found among animals and plants.

Epidemiological mechanisms of genetic resistance to kuru

Transmissible spongiform encephalopathies (TSEs), such as kuru, are invariably fatal neurodegenerative conditions caused by a malformation of the prion protein. Heterozygosity of codon 129 of the prion protein gene has been associated with increased host resistance to TSEs, although the mechanism by which this resistance is achieved has not been determined. To evaluate the epidemiological mechanism of human resistance to kuru, we developed a model that combines the dynamics of kuru transmission and the population genetics of human resistance. We fitted our model to kuru data from the epidemic that occurred in Papua New Guinea over the last hundred years. To elucidate the epidemiological mechanism of human resistance, we estimated the incubation period and transmission rate of kuru for codon 129 heterozygotes and homozygotes using kuru incidence data and human genotype frequency data from 1957 to 2004. Our results indicate that human resistance arises from a combination of both a longer incubation period and reduced susceptibility to infection. This work provides evidence for balancing selection acting on a human population and the mechanistic basis for the heterozygote resistance to kuru.

AuthorReward: increasing community curation in biological knowledge wikis through automated authorship quantification

Summary: Community curation—harnessing community intelligence in knowledge curation, bears great promise in dealing with the flood of biological knowledge. To exploit the full potential of the scientific community for knowledge curation, multiple biological wikis (bio-wikis) have been built to date. However, none of them have achieved a substantial impact on knowledge curation. One of the major limitations in bio-wikis is insufficient community participation, which is intrinsically because of lack of explicit authorship and thus no credit for community curation. To increase community curation in bio-wikis, here we develop AuthorReward, an extension to MediaWiki, to reward community-curated efforts in knowledge curation. AuthorReward quantifies researchers’ contributions by properly factoring both edit quantity and quality and yields automated explicit authorship according to their quantitative contributions. AuthorReward provides bio-wikis with an authorship metric, helpful to increase community participation in bio-wikis and to achieve community curation of massive biological knowledge.

Availability:http://cbb.big.ac.cn/software.

Contact:zhangzhang@big.ac.cn

Supplementary information:Supplementary data are available at Bioinformatics online.

Optimal targeting of seasonal influenza vaccination toward younger ages is robust to parameter uncertainty

Identification of the optimal vaccine allocation for the control of influenza requires consideration of uncertainty arising from numerous unpredictable factors, including viral evolution and diversity within the human population's immunity as well as variation in vaccine efficacy. The best policy must account for diverse potential outcomes based on these uncertainties. Here we used a mathematical model parametrized with survey-based contact data, demographic, and epidemiological data from seasonal influenza in the United States to determine the optimal vaccine allocation for five outcome measures: infections, hospitalizations, deaths, years of life loss, and contingent valuation. We incorporated uncertainty of epidemiological parameters and derive probability distributions of optimal age- and risk-specific allocation of vaccine. Our analysis demonstrated that previous recommendations of targeting schoolchildren (ages 5-17 years) and young adults (18-44 years) are generally robust in the face of uncertainty. However, when the outcome measure is to minimize deaths, years of life loss, or contingent valuation, uncertainty analysis identified scenarios under which it is optimal to target people at high risk for complications, even when vaccine are in abundance.

Abundant gene-by-environment interactions in gene expression reaction norms to copper within Saccharomyces cerevisiae

Genetic variation for plastic phenotypes potentially contributes phenotypic variation to populations that can be selected during adaptation to novel ecological contexts. However, the basis and extent of plastic variation that manifests in diverse environments remains elusive. Here, we characterize copper reaction norms for mRNA abundance among five Saccharomyces cerevisiae strains to 1) describe population variation across the full range of ecologically relevant copper concentrations, from starvation to toxicity, and 2) to test the hypothesis that plastic networks exhibit increased population variation for gene expression. We find that although the vast majority of the variation is small in magnitude (considerably <2-fold), not just some, but most genes demonstrate variable expression across environments, across genetic backgrounds, or both. Plastically expressed genes included both genes regulated directly by copper-binding transcription factors Mac1 and Ace1 and genes indirectly responding to the downstream metabolic consequences of the copper gradient, particularly genes involved in copper, iron, and sulfur homeostasis. Copper-regulated gene networks exhibited more similar behavior within the population in environments where those networks have a large impact on fitness. Nevertheless, expression variation in genes like Cup1, important to surviving copper stress, was linked with variation in mitotic fitness and in the breadth of differential expression across the genome. By revealing a broader and deeper range of population variation, our results provide further evidence for the interconnectedness of genome-wide mRNA levels, their dependence on environmental context and genetic background, and the abundance of variation in gene expression that can contribute to future evolution.

Yeast response to LA virus indicates coadapted global gene expression during mycoviral infection

Viruses that infect fungi have a ubiquitous distribution and play an important role in structuring fungal communities. Most of these viruses have an unusual life history in that they are propagated exclusively via asexual reproduction or fission of fungal cells. This asexual mode of transmission intimately ties viral reproductive success to that of its fungal host and should select for viruses that have minimal deleterious impact on the fitness of their hosts. Accordingly, viral infections of fungi frequently do not measurably impact fungal growth, and in some instances, increase the fitness of the fungal host. Here we determine the impact of the loss of coinfection by LA virus and the virus-like particle M1 upon global gene expression of the fungal host Saccharomyces cerevisiae and provide evidence supporting the idea that coevolution has selected for viral infection minimally impacting host gene expression.

Controlling antimicrobial resistance through targeted, vaccine-induced replacement of strains

Vaccination has proven effective in controlling many infectious diseases. However, differential effectiveness with regard to pathogen genotype is a frequent reason for failures in vaccine development. Often, insufficient immune response is induced to prevent infection by the diversity of existing serotypes present in pathogenic populations of bacteria. These vaccines that target a too narrow spectrum of serotypes do not offer sufficient prevention of infections, and can also lead to undesirable strain replacements. Here, we examine a novel idea to specifically exploit the narrow spectrum coverage of some vaccines to combat specific, emerging multi- and pan-resistant strains of pathogens. Application of a narrow-spectrum vaccine could serve to prevent infections by some strains that are hard to treat, rather than offer the vaccinated individual protection against infections by the pathogenic species as such. We suggest that vaccines targeted to resistant serotypes have the potential to become important public health tools, and would represent a new approach toward reducing the burden of particular multi-resistant strains occurring in hospitals. Vaccines targeting drug-resistant serotypes would also be the first clinical intervention with the potential to drive the evolution of pathogenic populations toward drug-sensitivity. We illustrate the feasibility of this approach by modeling a hypothetical vaccine that targets a subset of methicillin-resistant Staphylococcus aureus (MRSA) genotypes, in combination with drug treatment targeted at drug-sensitive genotypes. We find that a combined intervention strategy can limit nosocomial outbreaks, even when vaccine efficacy is imperfect. The broader utility of vaccine-based resistance control strategies should be further explored taking into account population structure, and the resistance and transmission patterns of the pathogen considered.

Analyses of expressed sequence tags in Neurospora reveal rapid evolution of genes associated with the early stages of sexual reproduction in fungi

Background

The broadly accepted pattern of rapid evolution of reproductive genes is primarily based on studies of animal systems, although several examples of rapidly evolving genes involved in reproduction are found in diverse additional taxa. In fungi, genes involved in mate recognition have been found to evolve rapidly. However, the examples are too few to draw conclusions on a genome scale.

Results

In this study, we performed microarray hybridizations between RNA from sexual and vegetative tissues of two strains of the heterothallic (self-sterile) filamentous ascomycete Neurospora intermedia, to identify a set of sex-associated genes in this species. We aligned Expressed Sequence Tags (ESTs) from sexual and vegetative tissue of N. intermedia to orthologs from three closely related species: N. crassa, N. discreta and N. tetrasperma. The resulting four-species alignments provided a dataset for molecular evolutionary analyses. Our results confirm a general pattern of rapid evolution of fungal sex-associated genes, compared to control genes with constitutive expression or a high relative expression during vegetative growth. Among the rapidly evolving sex-associated genes, we identified candidates that could be of importance for mating or fruiting-body development. Analyses of five of these candidate genes from additional species of heterothallic Neurospora revealed that three of them evolve under positive selection.

Conclusions

Taken together, our study represents a novel finding of a genome-wide pattern of rapid evolution of sex-associated genes in the fungal kingdom, and provides a list of candidate genes important for reproductive isolation in Neurospora.

Potential for rabies control through dog vaccination in wildlife-abundant communities of Tanzania

Canine vaccination has been successful in controlling rabies in diverse settings worldwide. However, concerns remain that coverage levels which have previously been sufficient might be insufficient in systems where transmission occurs both between and within populations of domestic dogs and other carnivores. To evaluate the effectiveness of vaccination targeted at domestic dogs when wildlife also contributes to transmission, we applied a next-generation matrix model based on contract tracing data from the Ngorongoro and Serengeti Districts in northwest Tanzania. We calculated corresponding values of R(0), and determined, for policy purposes, the probabilities that various annual vaccination targets would control the disease, taking into account the empirical uncertainty in our field data. We found that transition rate estimates and corresponding probabilities of vaccination-based control indicate that rabies transmission in this region is driven by transmission within domestic dogs. Different patterns of rabies transmission between the two districts exist, with wildlife playing a more important part in Ngorongoro and leading to higher recommended coverage levels in that district. Nonetheless, our findings indicate that an annual dog vaccination campaign achieving the WHO-recommended target of 70% will control rabies in both districts with a high level of certainty. Our results support the feasibility of controlling rabies in Tanzania through dog vaccination.

Reevaluation of epidemiological data demonstrates that it is consistent with cross-immunity among human papillomavirus types

BACKGROUND

The degree of cross-immunity between human papillomavirus (HPV) types is fundamental both to the epidemiological dynamics of HPV and to the impact of HPV vaccination. Epidemiological data on HPV infections has been repeatedly interpreted as inconsistent with cross-immunity.

METHODS

We reevaluate the epidemiological data using a model to determine the odds ratios of multiple to single infections expected in the presence or absence of cross-immunity. We simulate a virtual longitudinal survey to determine the effect cross-immunity has on the prevalence of multiple infections. We calibrate our model to epidemiological data and estimate the extent of type replacement following vaccination against specific HPV types.

RESULTS

We find that cross-immunity can produce odds ratios of infection comparable with epidemiological observations. We show that the sample sizes underlying existing surveys have been insufficient to identify even intense cross-immunity. We also find that the removal of HPV type 16, type 18, and types 6 and 11 would increase the prevalence of nontargeted types by 50%, 29%, and 183%, respectively.

CONCLUSIONS

Cross-immunity between HPV types is consistent with epidemiological data, contrary to previous interpretations. Cross-immunity may cause significant type replacement following vaccination, and therefore should be considered in future vaccine studies and epidemiological models.

RBE controls microRNA164 expression to effect floral organogenesis

The establishment and maintenance of organ boundaries are vital for animal and plant development. In the Arabidopsis flower, three microRNA164 genes (MIR164a, b and c) regulate the expression of CUP-SHAPED COTYLEDON1 (CUC1) and CUC2, which encode key transcriptional regulators involved in organ boundary specification. These three miR164 genes are expressed in distinct spatial and temporal domains that are crucial for their function. Here, we show that the C2H2 zinc finger transcriptional repressor encoded by RABBIT EARS (RBE) regulates the expression of all three miR164 genes. Furthermore, we demonstrate that RBE directly interacts with the promoter of MIR164c and negatively regulates its expression. We also show that the role of RBE in sepal and petal development is mediated in part through the concomitant regulation of the CUC1 and CUC2 gene products. These results indicate that one role of RBE is to fine-tune miR164 expression to regulate the CUC1 and CUC2 effector genes, which, in turn, regulate developmental events required for sepal and petal organogenesis.

Transcriptome analyses during fruiting body formation in Fusarium graminearum and Fusarium verticillioides reflect species life history and ecology

Fusarium graminearum and F. verticillioides are devastating cereal pathogens with very different life history and ecological characteristics. F. graminearum is homothallic, and sexual spores are an important component of its life cycle, responsible for disease initiation. F. verticilloides is heterothallic, and produces only modest numbers of fruiting bodies, which are not a significant source of inoculum. To identify corresponding differences in the transcriptional program underlying fruiting body development in the two species, comparative expression was performed, analyzing six developmental stages. To accompany the transcriptional analysis, detailed morphological characterization of F. verticillioides development was performed and compared to a previous morphological analysis of F. graminearum. Morphological development was similar between the two species, except for the observation of possible trichogynes in F. verticillioides ascogonia, which have not been previously reported for any Fusarium species. Expression of over 9000 orthologous genes were measured for the two species. Functional assignments of highly expressed orthologous genes at each time-point revealed the majority of highly expressed genes fell into the "unclassified proteins" category, reflecting the lack of characterization of genes for sexual development in both species. Simultaneous examination of morphological development and stage-specific gene expression suggests that degeneration of the paraphyses during sexual development is an apoptotic process. Expression of mating type genes in the two species differed, possibly reflecting the divergent roles they play in sexual development. Overall, the differences in gene expression reflect the greater role of fruiting bodies in the life cycle and ecology of F. graminearum versus F. verticillioides.

The influence of altruism on influenza vaccination decisions

Game theory is based on the assumption that individuals act according to self-interest and make decisions that maximize their personal payoffs. To test this fundamental assumption, we conducted a survey study in the context of influenza vaccination decisions. Contrary to the assumption of self-interest, we found that altruism plays an important role in vaccination decisions. Nevertheless, altruistic motivation has not yet been considered in epidemiological models, in predictions of vaccination decisions or in the design of vaccination policies. To determine the impact of altruism on the adherence to optimal vaccination policies and on resulting disease burden, we incorporated altruism into a game-theoretic epidemiological model of influenza vaccination. We found that altruism significantly shifted vaccination decisions away from individual self-interest and towards the community optimum, greatly reducing the total cost, morbidity and mortality for the community. Therefore, promoting altruism could be a potential strategy to improve public health outcomes.

Differential impact of nutrition on developmental and metabolic gene expression during fruiting body development in Neurospora crassa

Fungal fruiting body size and form are influenced by the ecology of the species, including diverse environmental stimuli. Accordingly, nutritional resources available to the fungus during development can be vital to successful production of fruiting bodies. To investigate the effect of nutrition, perithecial development of Neurospora crassa was induced on two different media, a chemically sparsely nutritive Synthetic Crossing Medium (SCM) and a natural Carrot Agar (CA). Protoperithecia were collected before crossing, and perithecia were collected at 2, 24, 48, 72, 96, 120, and at full maturity 144 h after crossing. No differences in fruiting body morphology were observed between the two media at any time point. A circuit of microarray hybridizations comparing cDNA from all neighboring stages was performed. For a majority of differentially expressed genes, expression was higher in SCM than in CA, and expression of core metabolic genes was particularly affected. Effects of nutrition were highest in magnitude before crossing, lowering in magnitude during early perithecial development. Interestingly, metabolic effects of the media were also large in magnitude during late perithecial development, at which stage the lower expression in CA presumably reflected the continued intake of diverse complex initial compounds, diminishing the need for expression of anabolic pathways. However, for genes with key regulatory roles in sexual development, including pheromone precursor ccg-4 and poi2, expression patterns were similar between treatments. When possible, a common nutritional environment is ideal for comparing transcriptional profiles between different fungi. Nevertheless, the observed consistency of the developmental program across media, despite considerable metabolic differentiation is reassuring. This result facilitates comparative studies that will require different nutritional resources for sexual development in different fungi.

Codon Deviation Coefficient: a novel measure for estimating codon usage bias and its statistical significance

Background

Genetic mutation, selective pressure for translational efficiency and accuracy, level of gene expression, and protein function through natural selection are all believed to lead to codon usage bias (CUB). Therefore, informative measurement of CUB is of fundamental importance to making inferences regarding gene function and genome evolution. However, extant measures of CUB have not fully accounted for the quantitative effect of background nucleotide composition and have not statistically evaluated the significance of CUB in sequence analysis.

Results

Here we propose a novel measure--Codon Deviation Coefficient (CDC)--that provides an informative measurement of CUB and its statistical significance without requiring any prior knowledge. Unlike previous measures, CDC estimates CUB by accounting for background nucleotide compositions tailored to codon positions and adopts the bootstrapping to assess the statistical significance of CUB for any given sequence. We evaluate CDC by examining its effectiveness on simulated sequences and empirical data and show that CDC outperforms extant measures by achieving a more informative estimation of CUB and its statistical significance.

Conclusions

As validated by both simulated and empirical data, CDC provides a highly informative quantification of CUB and its statistical significance, useful for determining comparative magnitudes and patterns of biased codon usage for genes or genomes with diverse sequence compositions.

Assessing the probability of detection of horizontal gene transfer events in bacterial populations

Experimental approaches to identify horizontal gene transfer (HGT) events of non-mobile DNA in bacteria have typically relied on detection of the initial transformants or their immediate offspring. However, rare HGT events occurring in large and structured populations are unlikely to be detected in a short time frame. Population genetic modeling of the growth dynamics of bacterial genotypes is therefore necessary to account for natural selection and genetic drift during the time lag and to predict realistic time frames for detection with a given sampling design. Here we draw on statistical approaches to population genetic theory to construct a cohesive probabilistic framework for investigation of HGT of exogenous DNA into bacteria. In particular, the stochastic timing of rare HGT events is accounted for. Integrating over all possible event timings, we provide an equation for the probability of detection, given that HGT actually occurred. Furthermore, we identify the key variables determining the probability of detecting HGT events in four different case scenarios that are representative of bacterial populations in various environments. Our theoretical analysis provides insight into the temporal aspects of dissemination of genetic material, such as antibiotic resistance genes or transgenes present in genetically modified organisms. Due to the long time scales involved and the exponential growth of bacteria with differing fitness, quantitative analyses incorporating bacterial generation time, and levels of selection, such as the one presented here, will be a necessary component of any future experimental design and analysis of HGT as it occurs in natural settings.

PhyDesign: an online application for profiling phylogenetic informativeness

Background

The rapid increase in number of sequenced genomes for species across of the tree of life is revealing a diverse suite of orthologous genes that could potentially be employed to inform molecular phylogenetic studies that encompass broader taxonomic sampling. Optimal usage of this diversity of loci requires user-friendly tools to facilitate widespread cost-effective locus prioritization for phylogenetic sampling. The Townsend (2007) phylogenetic informativeness provides a unique empirical metric for guiding marker selection. However, no software or automated methodology to evaluate sequence alignments and estimate the phylogenetic informativeness metric has been available.

Results

Here, we present PhyDesign, a platform-independent online application that implements the Townsend (2007) phylogenetic informativeness analysis, providing a quantitative prediction of the utility of loci to solve specific phylogenetic questions. An easy-to-use interface facilitates uploading of alignments and ultrametric trees to calculate and depict profiles of informativeness over specified time ranges, and provides rankings of locus prioritization for epochs of interest.

Conclusions

By providing these profiles, PhyDesign facilitates locus prioritization increasing the efficiency of sequencing for phylogenetic purposes compared to traditional studies with more laborious and low capacity screening methods, as well as increasing the accuracy of phylogenetic studies. Together with a manual and sample files, the application is freely accessible at http://phydesign.townsend.yale.edu.