Oligodendrocyte gene expression is reduced by and influences effects of chronic social stress in mice.

Oligodendrocyte gene expression is downregulated in stress-related neuropsychiatric disorders, including depression. In mice, chronic social stress (CSS) leads to depression-relevant changes in brain and emotional behavior, and the present study shows the involvement of oligodendrocytes in this model. In C57BL/6 (BL/6) mice, RNA-sequencing (RNA-Seq) was conducted with prefrontal cortex, amygdala and hippocampus from CSS and controls; a gene enrichment database for neurons, astrocytes and oligodendrocytes was used to identify cell origin of deregulated genes, and cell deconvolution was applied. To assess the potential causal contribution of reduced oligodendrocyte gene expression to CSS effects, mice heterozygous for the oligodendrocyte gene cyclic nucleotide phosphodiesterase (Cnp1) on a BL/6 background were studied; a 2 genotype (wildtype, Cnp1+/- ) × 2 environment (control, CSS) design was used to investigate effects on emotional behavior and amygdala microglia. In BL/6 mice, in prefrontal cortex and amygdala tissue comprising gray and white matter, CSS downregulated expression of multiple oligodendroycte genes encoding myelin and myelin-axon-integrity proteins, and cell deconvolution identified a lower proportion of oligodendrocytes in amygdala. Quantification of oligodendrocyte proteins in amygdala gray matter did not yield evidence for reduced translation, suggesting that CSS impacts primarily on white matter oligodendrocytes or the myelin transcriptome. In Cnp1 mice, social interaction was reduced by CSS in Cnp1+/- mice specifically; using ionized calcium-binding adaptor molecule 1 (IBA1) expression, microglia activity was increased additively by Cnp1+/- and CSS in amygdala gray and white matter. This study provides back-translational evidence that oligodendrocyte changes are relevant to the pathophysiology and potentially the treatment of stress-related neuropsychiatric disorders.

Is there evidence of optimisation for carbon efficiency in plant proteomes?

Flowering plants, angiosperms, can be divided into two major clades, monocots and dicots, and while differences in amino acid composition in different species from the two clades have been reported, a systematic analysis of amino acid content and distribution remains outstanding. Here, we show that monocot and dicot proteins have developed distinct amino acid content. In Arabidopsis thaliana and poplar, as in the ancestral moss Physcomitrella patens, the average mass per amino acid appears to be independent of protein length, while in the monocots rice, maize and sorghum, shorter proteins tend to be made of lighter amino acids. An examination of the elemental content of these proteomes reveals that the difference between monocot and dicot proteins can be largely attributed to their different carbon signatures. In monocots, the shorter proteins, which comprise the majority of all proteins, are made of amino acids with less carbon, while the nitrogen content is unchanged in both monocots and dicots. We hypothesise that this signature could be the result of carbon use and energy optimisation in fast-growing annual Poaceae (grasses).

Adaptive changes in HIV-1 subtype C proteins during early infection are driven by changes in HLA-associated immune pressure.

It is unresolved whether recently transmitted human immunodeficiency viruses (HIV) have genetic features that specifically favour their transmissibility. To identify potential "transmission signatures", we compared 20 full-length HIV-1 subtype C genomes from primary infections, with 66 sampled from ethnically and geographically matched individuals with chronic infections. Controlling for recombination and phylogenetic relatedness, we identified 39 sites at which amino acid frequency spectra differed significantly between groups. These sites were predominantly located within Env, Pol and Gag (14/39, 9/39 and 6/39 respectively) and were significantly clustered (33/39) within known immunoreactive peptides. Within 6 months of infection, we detected reversion-to-consensus mutations at 14 sites and potential CTL escape mutations at seven. Here we provide evidence that frequent reversion mutations probably allows the virus to recover replicative fitness which, together with immune escape driven by the HLA alleles of the new hosts, differentiate sequences from chronic infections from those sampled shortly after transmission.

Quantitating the multiplicity of infection with human immunodeficiency virus type 1 subtype C reveals a non-poisson distribution of transmitted variants.

Identifying the specific genetic characteristics of successfully transmitted variants may prove central to the development of effective vaccine and microbicide interventions. Although human immunodeficiency virus transmission is associated with a population bottleneck, the extent to which different factors influence the diversity of transmitted viruses is unclear. We estimate here the number of transmitted variants in 69 heterosexual men and women with primary subtype C infections. From 1,505 env sequences obtained using a single genome amplification approach we show that 78% of infections involved single variant transmission and 22% involved multiple variant transmissions (median of 3). We found evidence for mutations selected for cytotoxic-T-lymphocyte or antibody escape and a high prevalence of recombination in individuals infected with multiple variants representing another potential escape pathway in these individuals. In a combined analysis of 171 subtype B and C transmission events, we found that infection with more than one variant does not follow a Poisson distribution, indicating that transmission of individual virions cannot be seen as independent events, each occurring with low probability. While most transmissions resulted from a single infectious unit, multiple variant transmissions represent a significant fraction of transmission events, suggesting that there may be important mechanistic differences between these groups that are not yet understood.

A Bayesian model comparison approach to inferring positive selection.

A popular approach to detecting positive selection is to estimate the parameters of a probabilistic model of codon evolution and perform inference based on its maximum likelihood parameter values. This approach has been evaluated intensively in a number of simulation studies and found to be robust when the available data set is large. However, uncertainties in the estimated parameter values can lead to errors in the inference, especially when the data set is small or there is insufficient divergence between the sequences. We introduce a Bayesian model comparison approach to infer whether the sequence as a whole contains sites at which the rate of nonsynonymous substitution is greater than the rate of synonymous substitution. We incorporated this probabilistic model comparison into a Bayesian approach to site-specific inference of positive selection. Using simulated sequences, we compared this approach to the commonly used empirical Bayes approach and investigated the effect of tree length on the performance of both methods. We found that the Bayesian approach outperforms the empirical Bayes method when the amount of sequence divergence is small and is less prone to false-positive inference when the sequences are saturated, while the results are indistinguishable for intermediate levels of sequence divergence.

Expansin-like molecules: novel functions derived from common domains.

An Arabidopsis thaliana transcript ( AtPNP-A) encoding an immunoreactant plant natriuretic peptide (irPNP) analog was identified and isolated. The encoded protein shows similarity to CjBAp12, a functionally undefined protein from citrus that is induced in response to blight infection. CjBAp12 shows significant sequence identity to domains found in the cell wall loosening expansins but has tested negative for cell wall loosening activity. We have thus undertaken to establish the evolutionary and functional relationships of irPNP-like molecules within the superfamily of expansins, pollen allergens, and distantly related molecules such as endoglucanases. We show that irPNP-like molecules are related to expansins and fall in two groups; one includes CjBAp12 and the other AtPNP-A. Members of both groups share distinct sequence motifs (K[VI]VD and [LM]SxxAFxxI) but do not contain the tryptophan and tyrosine rich C-terminal putative polysaccharide-binding domain typical of expansins or bacterial cellulases and hemicellulases. We argue that both irPNP-like molecules and expansin have evolved from primitive/ancestral glucanase-like molecules that hydrolysed the cell wall. Importantly, we have previously demonstrated that irPNPs act on protoplasts, that is plant cells without cell walls as well as microsomes, indicating that these novel proteins specifically interact with the plasma membrane. It follows that the cell wall cannot be an obligatory substrate for irPNPs. Thus, both irPNP function and domain structure point to these molecules having a systemic role in H2O and solute homeostasis.

The contribution of exon-skipping events on chromosome 22 to protein coding diversity.

Completion of the human genome sequence provides evidence for a gene count with lower bound 30,000-40,000. Significant protein complexity may derive in part from multiple transcript isoforms. Recent EST based studies have revealed that alternate transcription, including alternative splicing, polyadenylation and transcription start sites, occurs within at least 30-40% of human genes. Transcript form surveys have yet to integrate the genomic context, expression, frequency, and contribution to protein diversity of isoform variation. We determine here the degree to which protein coding diversity may be influenced by alternate expression of transcripts by exhaustive manual confirmation of genome sequence annotation, and comparison to available transcript data to accurately associate skipped exon isoforms with genomic sequence. Relative expression levels of transcripts are estimated from EST database representation. The rigorous in silico method accurately identifies exon skipping using verified genome sequence. 545 genes have been studied in this first hand-curated assessment of exon skipping on chromosome 22. Combining manual assessment with software screening of exon boundaries provides a highly accurate and internally consistent indication of skipping frequency. 57 of 62 exon skipping events occur in the protein coding regions of 52 genes. A single gene, (FBXO7) expresses an exon repetition. 59% of highly represented multi-exon genes are likely to express exon-skipped isoforms in ratios that vary from 1:1 to 1:>100. The proportion of all transcripts corresponding to multi-exon genes that exhibit an exon skip is estimated to be 5%.

Prevalence of small inversions in yeast gene order evolution.

Gene order evolution in two eukaryotes was studied by comparing the Saccharomyces cerevisiae genome sequence to extensive new data from whole-genome shotgun and cosmid sequencing of Candida albicans. Gene order is substantially different between these two yeasts, with only 9% of gene pairs that are adjacent in one species being conserved as adjacent in the other. Inversion of small segments of DNA, less than 10 genes long, has been a major cause of rearrangement, which means that even where a pair of genes has been conserved as adjacent, the transcriptional orientations of the two genes relative to one another are often different. We estimate that about 1,100 single-gene inversions have occurred since the divergence between these species. Other genes that are adjacent in one species are in the same neighborhood in the other, but their precise arrangement has been disrupted, probably by multiple successive multigene inversions. We estimate that gene adjacencies have been broken as frequently by local rearrangements as by chromosomal translocations or long-distance transpositions. A bias toward small inversions has been suggested by other studies on animals and plants and may be general among eukaryotes.

Updated map of duplicated regions in the yeast genome.

We have updated the map of duplicated chromosomal segments in the Saccharomyces cerevisiae genome originally published by Wolfe and Shields in 1997 (Nature 387, 708-713). The new analysis is based on the more sensitive Smith Waterman search method instead of BLAST. The parameters used to identify duplicated chromosomal regions were optimized such as to maximize the amount of the genome placed into paired regions, under the assumption that the hypothesis that the entire genome was duplicated in a single event is correct. The core of the new map, with 52 pairs of regions containing three or more duplicated genes, is largely unchanged from our original map. 39 tRNA gene pairs and one snRNA pair have been added. To find additional pairs of genes that may have been formed by whole genome duplication, we searched through the parts of the genome that are not covered by this core map, looking for putative duplicated chromosomal regions containing only two duplicate genes instead of three, or having lower-scoring gene pairs. This approach identified a further 32 candidate paired regions, bringing the total number of protein-coding genes on the duplication map to 905 (16% of the proteome). The updated map suggests that a second copy of the ribosomal DNA array has been deleted from chromosome IV.

Yeast genome evolution in the post-genome era.

The Saccharomyces cerevisiae genome sequence, augmented by new data on gene expression and function, continues to yield new findings about eukaryote genome evolution. Analysis of the duplicate gene pairs formed by whole-genome duplication indicates that selection for increased levels of gene expression was a significant factor determining which genes were retained as duplicates and which were returned to a single-copy state, possibly in addition to selection for novel gene functions. Proteome comparisons between worm and yeast show that genes for core metabolic processes are shared among eukaryotes and unchanging in function, while comparisons between different yeast species identify 'orphan' genes as the most rapidly evolving fraction of the proteome. Natural hybridisation among yeast species is frequent, but its long-term evolutionary significance is unknown.

G+C content variation along and among Saccharomyces cerevisiae chromosomes.

Past analyses of the genome of the yeast Saccharomyces cerevisiae have revealed substantial regional variation in G+C content. Important questions remain, though, as to the origin, nature, significance, and generality of this variation. We conducted an extensive analysis of the yeast genome to try to answer these questions. Our results indicate that open reading frames (ORFs) with similar G+C contents at silent codon positions are significantly clustered on chromosomes. This clustering can be explained by very short range correlations of silent-site G+C contents at neighboring ORFs. ORFs of high silent-site G+C content are disproportionately concentrated on shorter chromosomes, which causes a negative relationship between chromosome length and G+C content. Contrary to previous reports, there is no correlation between gene density and silent-site G+C content in yeast. Chromosome III is atypical in many regards, and possible reasons for this are discussed.

Extent of genomic rearrangement after genome duplication in yeast.

Whole-genome duplication approximately 10(8) years ago was proposed as an explanation for the many duplicated chromosomal regions in Saccharomyces cerevisiae. Here we have used computer simulations and analytic methods to estimate some parameters describing the evolution of the yeast genome after this duplication event. Computer simulation of a model in which 8% of the original genes were retained in duplicate after genome duplication, and 70-100 reciprocal translocations occurred between chromosomes, produced arrangements of duplicated chromosomal regions very similar to the map of real duplications in yeast. An analytical method produced an independent estimate of 84 map disruptions. These results imply that many smaller duplicated chromosomal regions exist in the yeast genome in addition to the 55 originally reported. We also examined the possibility of determining the original order of chromosomal blocks in the ancestral unduplicated genome, but this cannot be done without information from one or more additional species. If the genome sequence of one other species (such as Kluyveromyces lactis) were known it should be possible to identify 150-200 paired regions covering the whole yeast genome and to reconstruct approximately two-thirds of the original order of blocks of genes in yeast. Rates of interchromosome translocation in yeast and mammals appear similar despite their very different rates of homologous recombination per kilobase.

Evolution of gene order and chromosome number in Saccharomyces, Kluyveromyces and related fungi.

The extent to which the order of genes along chromosomes is conserved between Saccharomyces cerevisiae and related species was studied by analysing data from DNA sequence database. As expected, the extent of gene order conservation decreases with increasing evolutionary distance. About 59% of adjacent gene pairs in Kluyveromyces lactis or K. marxianus are also adjacent in S. cerevisiae, and a further 16% of Kluyveromyces neighbours can be explained in terms of the inferred ancestral gene order in Saccharomyces prior to the occurrence of an ancient whole-genome duplication. Only 13% of Candida albicans linkages, and no Schizosaccharomyces pombe linkages, are conserved. Analysis of gene order arrangements, chromosome numbers, and ribosomal RNA sequences suggests that genome duplication occurred before the divergence of the four species in Saccharomyces sensu stricto (all of which have 16 chromosomes), but after this lineage had diverged from Saccharomyces kluyveri and the Kluyveromyces lactislmarxianus species assemblage.