Hofbauer cell function in the term placenta associates with adult cardiovascular and depressive outcomes.

Pathological placental inflammation increases the risk for several adult disorders, but these mediators are also expressed under homeostatic conditions, where their contribution to adult health outcomes is unknown. Here we define an inflammation-related expression signature, primarily expressed in Hofbauer cells of the term placenta and use expression quantitative trait loci to create a polygenic score (PGS) predictive of its expression. Using this PGS in the UK Biobank we conduct a phenome-wide association study, followed by Mendelian randomization and identify protective, sex-dependent effects of the placental module on cardiovascular and depressive outcomes. Genes differentially regulated by intra-amniotic infection and preterm birth are over-represented within the module. We also identify aspirin as a putative modulator of this inflammation-related signature. Our data support a model where disruption of placental Hofbauer cell function, due to preterm birth or prenatal infection, contributes to the increased risk of depression and cardiovascular disease observed in these individuals.

Sex and cell-specific gene expression in corticolimbic brain regions associated with psychiatric disorders revealed by bulk and single-nuclei RNA sequencing.

BACKGROUND: There are sex-specific differences in the prevalence, symptomology and course of psychiatric disorders. However, preclinical models have primarily used males, such that the molecular mechanisms underlying sex-specific differences in psychiatric disorders are not well established. METHODS: In this study, we compared transcriptome-wide gene expression profiles in male and female rats within the corticolimbic system, including the cingulate cortex, nucleus accumbens medial shell (NAcS), ventral dentate gyrus and the basolateral amygdala (n = 22-24 per group/region). FINDINGS: We found over 3000 differentially expressed genes (DEGs) in the NAcS between males and females. Of these DEGs in the NAcS, 303 showed sex-dependent conservation DEGs in humans and were significantly enriched for gene ontology terms related to blood vessel morphogenesis and regulation of cell migration. Single nuclei RNA sequencing in the NAcS of male and female rats identified widespread sex-dependent expression, with genes upregulated in females showing a notable enrichment for synaptic function. Female upregulated genes in astrocytes, Drd3+MSNs and oligodendrocyte were also enriched in several psychiatric genome-wide association studies (GWAS). INTERPRETATION: Our data provide comprehensive evidence of sex- and cell-specific molecular profiles in the NAcS. Importantly these differences associate with anxiety, bipolar disorder, schizophrenia, and cross-disorder, suggesting an intrinsic molecular basis for sex-based differences in psychiatric disorders that strongly implicates the NAcS. FUNDING: This work was supported by funding from the Hope for Depression Research Foundation (MJM).

Integrative multi-omics landscape of fluoxetine action across 27 brain regions reveals global increase in energy metabolism and region-specific chromatin remodelling.

Depression and anxiety are major global health burdens. Although SSRIs targeting the serotonergic system are prescribed over 200 million times annually, they have variable therapeutic efficacy and side effects, and mechanisms of action remain incompletely understood. Here, we comprehensively characterise the molecular landscape of gene regulatory changes associated with fluoxetine, a widely-used SSRI. We performed multimodal analysis of SSRI response in 27 mammalian brain regions using 310 bulk RNA-seq and H3K27ac ChIP-seq datasets, followed by in-depth characterisation of two hippocampal regions using single-cell RNA-seq (20 datasets). Remarkably, fluoxetine induced profound region-specific shifts in gene expression and chromatin state, including in the nucleus accumbens shell, locus coeruleus and septal areas, as well as in more well-studied regions such as the raphe and hippocampal dentate gyrus. Expression changes were strongly enriched at GWAS loci for depression and antidepressant drug response, stressing the relevance to human phenotypes. We observed differential expression at dozens of signalling receptors and pathways, many of which are previously unknown. Single-cell analysis revealed stark differences in fluoxetine response between the dorsal and ventral hippocampal dentate gyri, particularly in oligodendrocytes, mossy cells and inhibitory neurons. Across diverse brain regions, integrative omics analysis consistently suggested increased energy metabolism via oxidative phosphorylation and mitochondrial changes, which we corroborated in vitro; this may thus constitute a shared mechanism of action of fluoxetine. Similarly, we observed pervasive chromatin remodelling signatures across the brain. Our study reveals unexpected regional and cell type-specific heterogeneity in SSRI action, highlights under-studied brain regions that may play a major role in antidepressant response, and provides a rich resource of candidate cell types, genes, gene regulatory elements and pathways for mechanistic analysis and identifying new therapeutic targets for depression and anxiety.

Genomic modules and intramodular network concordance in susceptible and resilient male mice across models of stress.

The multifactorial etiology of stress-related disorders necessitates a constant interrogation of the molecular convergences in preclinical models of stress that use disparate paradigms as stressors spanning from environmental challenges to genetic predisposition to hormonal signaling. Using RNA-sequencing, we investigated the genomic signatures in the ventral hippocampus common to mouse models of stress. Chronic oral corticosterone (CORT) induced increased anxiety- and depression-like behavior in wild-type male mice and male mice heterozygous for the gene coding for brain-derived neurotrophic factor Val66Met, a variant associated with genetic susceptibility to stress. In a separate set of male mice, chronic social defeat stress (CSDS) led to a susceptible or a resilient population, whose proportion was dependent on housing conditions, namely standard housing or enriched environment. Rank-rank-hypergeometric overlap (RRHO), a threshold-free approach that ranks genes by their p value and effect size direction, was used to identify genes from a continuous gradient of significancy that were concordant across groups. In mice treated with CORT and in standard-housed susceptible mice, differentially expressed genes (DEGs) were concordant for gene networks involved in neurotransmission, cytoskeleton function, and vascularization. Weighted gene co-expression analysis generated 54 gene hub modules and revealed two modules in which both CORT and CSDS-induced enrichment in DEGs, whose function was concordant with the RRHO predictions, and correlated with behavioral resilience or susceptibility. These data showed transcriptional concordance across models in which the stress coping depends upon hormonal, environmental, or genetic factors revealing common genomic drivers that embody the multifaceted nature of stress-related disorders.

Environmental enrichment increases transcriptional and epigenetic differentiation between mouse dorsal and ventral dentate gyrus.

Early life experience influences stress reactivity and mental health through effects on cognitive-emotional functions that are, in part, linked to gene expression in the dorsal and ventral hippocampus. The hippocampal dentate gyrus (DG) is a major site for experience-dependent plasticity associated with sustained transcriptional alterations, potentially mediated by epigenetic modifications. Here, we report comprehensive DNA methylome, hydroxymethylome and transcriptome data sets from mouse dorsal and ventral DG. We find genome-wide transcriptional and methylation differences between dorsal and ventral DG, including at key developmental transcriptional factors. Peripubertal environmental enrichment increases hippocampal volume and enhances dorsal DG-specific differences in gene expression. Enrichment also enhances dorsal-ventral differences in DNA methylation, including at binding sites of the transcription factor NeuroD1, a regulator of adult neurogenesis. These results indicate a dorsal-ventral asymmetry in transcription and methylation that parallels well-known functional and anatomical differences, and that may be enhanced by environmental enrichment.

Regulating the transcriptomes that mediate the conversion of fibroblasts to various nervous system neural cell types.

Our understanding of the mechanism of cell fate transition during the direct reprogramming of fibroblasts into various central nervous system (CNS) neural cell types has been limited by the lack of a comprehensive analysis on generated cells, independently and in comparison with other CNS neural cell types. Here, we applied an integrative approach on 18 independent high throughput expression data sets to gain insight into the regulation of the transcriptome during the conversion of fibroblasts into induced neural stem cells, induced neurons (iNs), induced astrocytes, and induced oligodendrocyte progenitor cells (iOPCs). We found common down-regulated genes to be mostly related to fibroblast-specific functions, and suggest their potential as markers for screening of the silencing of the fibroblast-specific program. For example, Tagln was significantly down-regulated across all considered data sets. In addition, we identified specific profiles of up-regulated genes for each CNS neural cell types, which could be potential markers for maturation and efficiency screenings. Furthermore, we identified the main TFs involved in the regulation of the gene expression program during direct reprogramming. For example, in the generation of iNs from fibroblasts, the Rest TF was the main regulator of this reprogramming. In summary, our computational approach for meta-analyzing independent expression data sets provides significant details regarding the molecular mechanisms underlying the regulation of the gene expression program, and also suggests potentially useful candidate genes for screening down-regulation of fibroblast gene expression profile, maturation, and efficiency, as well as candidate TFs for increasing the efficiency of the reprogramming process.

Genetics of mating in members of the Chaetomiaceae as revealed by experimental and genomic characterization of reproduction in Myceliophthora heterothallica.

Members of the Chaetomiaceae are among the most studied fungi in industry and among the most reported in investigations of biomass degradation in both natural and laboratory settings. The family is recognized for production of carbohydrate-active enzymes and antibiotics. Thermophilic species are of special interest for their abilities to produce thermally stable enzymes and to be grown under conditions that are unsuitable for potential contaminant microorganisms. Such interests led to the recent acquisition of genome sequences from several members of the family, including thermophilic species, several of which are reported here for the first time. To date, however, thermophilic fungi in industry have served primarily as parts reservoirs and there has been no good genetic model for species in the family Chaetomiaceae or for thermophiles in general. We report here on the reproductive biology of the thermophile Myceliophthora heterothallica, which is heterothallic, unlike most described species in the family. We confirmed heterothallism genetically by following the segregation of mating type idiomorphs and other markers. We have expanded the number of known sexually-compatible individuals from the original isolates from Indiana and Germany to include several isolates from New Mexico. An interesting aspect of development in M. heterothallica is that ascocarp formation is optimal at approximately 30 °C, whereas vegetative growth is optimal at 45 °C. Genome sequences obtained from several strains, including isolates of each mating type, revealed mating-type regions whose genes are organized similarly to those of other members of the Sordariales, except for the presence of a truncated version of the mat A-1 (MAT1-1-1) gene in mating-type a (MAT1-2) strains. In M. heterothallica and other Chaetomiaceae, mating-type A (MAT1-1) strains have the full-length version of mat A-1 that is typical of mating-type A strains of diverse Ascomycota, whereas a strains have only the truncated version. This truncated mat A-1 has an intact open reading frame and a derived start codon that is not present in mat A-1 from A strains. The predicted protein contains a region that is conserved across diverse mat A-1 genes, but it lacks the major alpha1 domain, which characterizes proteins in this family and is known to be required for fertility in A strains from other Ascomycota. Finally, we have used genes from M. heterothallica to probe for mating genes in other homothallic and heterothallic members of the Chaetomiaceae. The majority of homothallic species examined have a typical mat A-1,2,3 (MAT1-1-1,2,3) region in addition to an unlinked mat a-1 (MAT1-2-1) gene, reflecting one type of homothallism commonly observed in diverse Ascomycota.

SnowyOwl: accurate prediction of fungal genes by using RNA-Seq and homology information to select among ab initio models.

BACKGROUND: Locating the protein-coding genes in novel genomes is essential to understanding and exploiting the genomic information but it is still difficult to accurately predict all the genes. The recent availability of detailed information about transcript structure from high-throughput sequencing of messenger RNA (RNA-Seq) delineates many expressed genes and promises increased accuracy in gene prediction. Computational gene predictors have been intensively developed for and tested in well-studied animal genomes. Hundreds of fungal genomes are now or will soon be sequenced. The differences of fungal genomes from animal genomes and the phylogenetic sparsity of well-studied fungi call for gene-prediction tools tailored to them. RESULTS: SnowyOwl is a new gene prediction pipeline that uses RNA-Seq data to train and provide hints for the generation of Hidden Markov Model (HMM)-based gene predictions and to evaluate the resulting models. The pipeline has been developed and streamlined by comparing its predictions to manually curated gene models in three fungal genomes and validated against the high-quality gene annotation of Neurospora crassa; SnowyOwl predicted N. crassa genes with 83% sensitivity and 65% specificity. SnowyOwl gains sensitivity by repeatedly running the HMM gene predictor Augustus with varied input parameters and selectivity by choosing the models with best homology to known proteins and best agreement with the RNA-Seq data. CONCLUSIONS: SnowyOwl efficiently uses RNA-Seq data to produce accurate gene models in both well-studied and novel fungal genomes. The source code for the SnowyOwl pipeline (in Python) and a web interface (in PHP) is freely available from http://sourceforge.net/projects/snowyowl/.

Snapshot of the eukaryotic gene expression in muskoxen rumen--a metatranscriptomic approach.

BACKGROUND: Herbivores rely on digestive tract lignocellulolytic microorganisms, including bacteria, fungi and protozoa, to derive energy and carbon from plant cell wall polysaccharides. Culture independent metagenomic studies have been used to reveal the genetic content of the bacterial species within gut microbiomes. However, the nature of the genes encoded by eukaryotic protozoa and fungi within these environments has not been explored using metagenomic or metatranscriptomic approaches. METHODOLOGY/PRINCIPAL FINDINGS: In this study, a metatranscriptomic approach was used to investigate the functional diversity of the eukaryotic microorganisms within the rumen of muskoxen (Ovibos moschatus), with a focus on plant cell wall degrading enzymes. Polyadenylated RNA (mRNA) was sequenced on the Illumina Genome Analyzer II system and 2.8 gigabases of sequences were obtained and 59129 contigs assembled. Plant cell wall degrading enzyme modules including glycoside hydrolases, carbohydrate esterases and polysaccharide lyases were identified from over 2500 contigs. These included a number of glycoside hydrolase family 6 (GH6), GH48 and swollenin modules, which have rarely been described in previous gut metagenomic studies. CONCLUSIONS/SIGNIFICANCE: The muskoxen rumen metatranscriptome demonstrates a much higher percentage of cellulase enzyme discovery and an 8.7x higher rate of total carbohydrate active enzyme discovery per gigabase of sequence than previous rumen metagenomes. This study provides a snapshot of eukaryotic gene expression in the muskoxen rumen, and identifies a number of candidate genes coding for potentially valuable lignocellulolytic enzymes.

Combining proteomics of root and shoot mitochondria and transcript analysis to define constitutive and variable components in plant mitochondria.

Mitochondria undertake respiration in plant cells, but through metabolic plasticity utilize differ proportions of substrates and deliver different proportions of products to cellular metabolic and biosynthetic pathways. In Arabidopsis the mitochondrial proteome from shoots and cell culture have been reported, but there has been little information on mitochondria in roots. We compare the root mitochondrial proteome with mitochondria isolated from photosynthetic shoots to define the role of protein abundance in these differences. The major differences observed were in the abundance and/or activities of enzymes in the TCA cycle and the mitochondrial enzymes involved in photorespiration. Metabolic pathways linked to TCA cycle and photorespiration were also altered, namely cysteine, formate and one-carbon metabolism, as well as amino acid metabolism focused on 2-oxoglutarate generation. Comparisons to microarray analysis of these same tissues showed a positive correlation between mRNA and mitochondrial protein abundance, but still ample evidence for the role of post-transcriptional processes in defining mitochondrial composition. Broader comparisons of transcript abundances for mitochondrial components across Arabidopsis tissues provided additional evidence for specialization of plant mitochondria, and clustering of these data in functional groups showed the constitutive vs variably expressed components of plant mitochondria.

Divalent metal ions in plant mitochondria and their role in interactions with proteins and oxidative stress-induced damage to respiratory function.

Understanding the metal ion content of plant mitochondria and metal ion interactions with the proteome are vital for insights into both normal respiratory function and the process of protein damage during oxidative stress. We have analyzed the metal content of isolated Arabidopsis (Arabidopsis thaliana) mitochondria, revealing a 26:8:6:1 molar ratio for iron:zinc:copper:manganese and trace amounts of cobalt and molybdenum. We show that selective changes occur in mitochondrial copper and iron content following in vivo and in vitro oxidative stresses. Immobilized metal affinity chromatography charged with Cu(2+), Zn(2+), and Co(2+) was used to identify over 100 mitochondrial proteins with metal-binding properties. There were strong correlations between the sets of immobilized metal affinity chromatography-interacting proteins, proteins predicted to contain metal-binding motifs, and protein sets known to be oxidized or degraded during abiotic stress. Mitochondrial respiratory chain pathways and matrix enzymes varied widely in their susceptibility to metal-induced loss of function, showing the selectivity of the process. A detailed study of oxidized residues and predicted metal interaction sites in the tricarboxylic acid cycle enzyme aconitase identified selective oxidation of residues in the active site and showed an approach for broader screening of functionally significant oxidation events in the mitochondrial proteome.

Phage-type RNA polymerase RPOTmp performs gene-specific transcription in mitochondria of Arabidopsis thaliana.

Transcription of mitochondrial genes in animals, fungi, and plants relies on the activity of T3/T7 phage-type RNA polymerases. Two such enzymes, RPOTm and RPOTmp, are present in the mitochondria of eudicotyledonous plants; RPOTmp is additionally found in plastids. We have characterized the transcriptional role of the dual-targeted RNA polymerase in mitochondria of Arabidopsis thaliana. Examination of mitochondrial transcripts in rpoTmp mutants revealed major differences in transcript abundances between wild-type and rpoTmp plants. Decreased levels of specific transcripts were correlated with reduced abundances of the respiratory chain complexes I and IV. Altered transcript levels in rpoTmp were found to result from gene-specific transcriptional changes, establishing that RPOTmp functions in distinct transcriptional processes within mitochondria. Decreased transcription of specific genes in rpoTmp was not associated with changes in promoter utilization; therefore, RPOTmp function is not promoter specific but gene specific. This implies that additional gene-specific elements direct the transcription of a subset of mitochondrial genes by RPOTmp.

Insights into female sperm storage from the spermathecal fluid proteome of the honeybee Apis mellifera.

BACKGROUND: Female animals are often able to store sperm inside their body--in some species even for several decades. The molecular basis of how females keep non-own cells alive is largely unknown, but since sperm cells are reported to be transcriptionally silenced and, therefore, limited in their ability to maintain their own function, it is likely that females actively participate in sperm maintenance. Because female contributions are likely to be of central importance for sperm survival, molecular insights into the process offer opportunities to observe mechanisms through which females manipulate sperm. RESULTS: We used the honeybee, Apis mellifera, in which queens are highly polyandrous and able to maintain sperm viable for several years. We identified over a hundred proteins representing the major constituents of the spermathecal fluid, which females contribute to sperm in storage. We found that the gel profile of proteins from spermathecal fluid is very similar to the secretions of the spermathecal gland and concluded that the spermathecal glands are the main contributors to the spermathecal fluid proteome. A detailed analysis of the spermathecal fluid proteins indicate that they fall into a range of different functional groups, most notably enzymes of energy metabolism and antioxidant defense. A metabolic network analysis comparing the proteins detected in seminal fluid and spermathecal fluid showed a more integrated network is present in the spermathecal fluid that could facilitate long-term storage of sperm. CONCLUSIONS: We present a large-scale identification of proteins in the spermathecal fluid of honeybee queens and provide insights into the molecular regulation of female sperm storage.

Approaches to defining dual-targeted proteins in Arabidopsis.

A variety of approaches were used to predict dual-targeted proteins in Arabidopsis thaliana. These predictions were experimentally tested using GFP fusions. Twelve new dual-targeted proteins were identified: five that were dual-targeted to mitochondria and plastids, six that were dual-targeted to mitochondria and peroxisomes, and one that was dual-targeted to mitochondria and the nucleus. Two methods to predict dual-targeted proteins had a high success rate: (1) combining the AraPerox database with a variety of subcellular prediction programs to identify mitochondrial- and peroxisomal-targeted proteins, and (2) using a variety of prediction programs on a biochemical pathway or process known to contain at least one dual-targeted protein. Several technical parameters need to be taken into account before assigning subcellular localization using GFP fusion proteins. The position of GFP with respect to the tagged polypeptide, the tissue or cells used to detect subcellular localization, and the portion of a candidate protein fused to GFP are all relevant to the expression and targeting of a fusion protein. Testing all gene models for a chromosomal locus is required if more than one model exists.

Novel proteins, putative membrane transporters, and an integrated metabolic network are revealed by quantitative proteomic analysis of Arabidopsis cell culture peroxisomes.

Peroxisomes play key roles in energy metabolism, cell signaling, and plant development. A better understanding of these important functions will be achieved with a more complete definition of the peroxisome proteome. The isolation of peroxisomes and their separation from mitochondria and other major membrane systems have been significant challenges in the Arabidopsis (Arabidopsis thaliana) model system. In this study, we present new data on the Arabidopsis peroxisome proteome obtained using two new technical advances that have not previously been applied to studies of plant peroxisomes. First, we followed density gradient centrifugation with free-flow electrophoresis to improve the separation of peroxisomes from mitochondria. Second, we used quantitative proteomics to identify proteins enriched in the peroxisome fractions relative to mitochondrial fractions. We provide evidence for peroxisomal localization of 89 proteins, 36 of which have not previously been identified in other analyses of Arabidopsis peroxisomes. Chimeric green fluorescent protein constructs of 35 proteins have been used to confirm their localization in peroxisomes or to identify endoplasmic reticulum contaminants. The distribution of many of these peroxisomal proteins between soluble, membrane-associated, and integral membrane locations has also been determined. This core peroxisomal proteome from nonphotosynthetic cultured cells contains a proportion of proteins that cannot be predicted to be peroxisomal due to the lack of recognizable peroxisomal targeting sequence 1 (PTS1) or PTS2 signals. Proteins identified are likely to be components in peroxisome biogenesis, beta-oxidation for fatty acid degradation and hormone biosynthesis, photorespiration, and metabolite transport. A considerable number of the proteins found in peroxisomes have no known function, and potential roles of these proteins in peroxisomal metabolism are discussed. This is aided by a metabolic network analysis that reveals a tight integration of functions and highlights specific metabolite nodes that most probably represent entry and exit metabolites that could require transport across the peroxisomal membrane.

Large-scale structural modeling of protein complexes at low resolution.

Structural aspects of protein-protein interactions provided by large-scale, genome-wide studies are essential for the description of life processes at the molecular level. A methodology is developed that applies the protein docking approach (GRAMM), based on the knowledge of experimentally determined protein-protein structures (DOCKGROUND resource) and properties of intermolecular energy landscapes, to genome-wide systems of protein interactions. The full sequence-to-structure-of-complex modeling pipeline is implemented in the Genome Wide Docking Database (GWIDD) resource. Protein interaction data are imported to GWIDD from external datasets of experimentally determined interaction networks. Essential information is extracted and unified to form the GWIDD database. Structures of individual interacting proteins in the database are retrieved (if available) or modeled, and protein complex structures are predicted by the docking program. All protein sequence, structure, and docking information is conveniently accessible through a Web interface.

Heterogeneity of the mitochondrial proteome for photosynthetic and non-photosynthetic Arabidopsis metabolism.

Heterogeneity of the mitochondrial proteome in plants underlies fundamental differences in the roles of these organelles in different tissues. We quantitatively compared the mitochondrial proteome isolated from a non-photosynthetic cell culture model with more specialized mitochondria isolated from photosynthetic shoots. Differences in intact mitochondrial respiratory rates with various substrates and activities of specific enzymes provided a backdrop of the functional variation between these mitochondrial populations. Proteomics comparisons provided a deep insight into the different steady-state abundances of specific mitochondrial proteins. Combined these data showed the elevated level of the photorespiratory apparatus and its complex interplay with glycolate, cysteine, formate, and one-carbon metabolism as well as the decrease of selected parts of the tricarboxylic acid cycle, alterations in amino acid metabolism focused on 2-oxoglutarate generation, and degradation of branched chain amino acids. Comparisons with microarray analysis of these tissue types showed a positive, mild correlation between mRNA and mitochondrial protein abundance, a tighter correlation for specific biochemical pathways, but over 78% concordance in direction between changes in protein and transcript abundance in the two tissues. Overall these results indicated that the majority of the variation in the plant mitochondrial proteome occurred in the matrix, highlighted the constitutive nature of the respiratory apparatus, and showed the differences in substrate choice and/or availability during photosynthetic and non-photosynthetic metabolism.

On the expansion of the pentatricopeptide repeat gene family in plants.

Pentatricopeptide repeat (PPR) proteins form a huge family in plants (450 members in Arabidopsis and 477 in rice) defined by tandem repetitions of characteristic sequence motifs. Some of these proteins have been shown to play a role in posttranscriptional processes within organelles, and they are thought to be sequence-specific RNA-binding proteins. The origins of this family are obscure as they are lacking from almost all prokaryotes, and the spectacular expansion of the family in land plants is equally enigmatic. In this study, we investigate the growth of the family in plants by undertaking a genome-wide identification and comparison of the PPR genes of 3 organisms: the flowering plants Arabidopsis thaliana and Oryza sativa and the moss Physcomitrella patens. A large majority of the PPR genes in each of the flowering plants are intron less. In contrast, most of the 103 PPR genes in Physcomitrella are intron rich. A phylogenetic comparison of the PPR genes in all 3 species shows similarities between the intron-rich PPR genes in Physcomitrella and the few intron-rich PPR genes in higher plants. Intron-poor PPR genes in all 3 species also display a bias toward a position of their introns at their 5' ends. These results provide compelling evidence that one or more waves of retrotransposition were responsible for the expansion of the PPR gene family in flowering plants. The differing numbers of PPR proteins are highly correlated with differences in organellar RNA editing between the 3 species.

Large-scale characteristics of the energy landscape in protein-protein interactions.

Characterization of intermolecular energy landscapes in protein-protein interactions is important for understanding the mechanisms of these interactions as well as for designing better protein docking methods. A simplified representation of the landscape was used for a systematic study of its large-scale characteristics in a large nonredundant dataset of protein complexes. The focus of the study is on the basic features of the low-resolution energy basins and their distribution on the landscape. The results clearly show that, in general, the number of such basins is small, these basins are well formed, correlated with actual binding modes, and the pattern of basins distribution depends on the type of the complex. For docking studies, the results suggest that adequate starting points for the structural refinement are detectable by low-resolution procedures and the number of such starting points is relatively small.

Genome-wide analysis of mRNA decay rates and their determinants in Arabidopsis thaliana.

To gain a global view of mRNA decay in Arabidopsis thaliana, suspension cell cultures were treated with a transcriptional inhibitor, and microarrays were used to measure transcript abundance over time. The deduced mRNA half-lives varied widely, from minutes to >24 h. Three features of the transcript displayed a correlation with decay rates: (1) genes possessing at least one intron produce mRNA transcripts significantly more stable than those of intronless genes, and this was not related to overall length, sequence composition, or number of introns; (2) various sequence elements in the 3' untranslated region are enriched among short- and long-lived transcripts, and their multiple occurrence suggests combinatorial control of transcript decay; and (3) transcripts that are microRNA targets generally have short half-lives. The decay rate of transcripts correlated with subcellular localization and function of the encoded proteins. Analysis of transcript decay rates for genes encoding orthologous proteins between Arabidopsis, yeast, and humans indicated that yeast and humans had a higher percentage of transcripts with shorter half-lives and that the relative stability of transcripts from genes encoding proteins involved in cell cycle, transcription, translation, and energy metabolism is conserved. Comparison of decay rates with changes in transcript abundance under a variety of abiotic stresses reveal that a set of transcription factors are downregulated with similar kinetics to decay rates, suggesting that inhibition of their transcription is an important early response to abiotic stress.