Istvan Ladunga

PANTHER: a browsable database of gene products organized by biological function, using curated protein family and subfamily classification

The PANTHER database was designed for high-throughput analysis of protein sequences. One of the key features is a simplified ontology of protein function, which allows browsing of the database by biological functions. Biologist curators have associated the ontology terms with groups of protein sequences rather than individual sequences. Statistical models (Hidden Markov Models, or HMMs) are built from each of these groups. The advantage of this approach is that new sequences can be automatically classified as they become available. To ensure accurate functional classification, HMMs are constructed not only for families, but also for functionally distinct subfamilies. Multiple sequence alignments and phylogenetic trees, including curator-assigned information, are available for each family. The current version of the PANTHER database includes training sequences from all organisms in the GenBank non-redundant protein database, and the HMMs have been used to classify gene products across the entire genomes of human, and Drosophila melanogaster. The ontology terms and protein families and subfamilies, as well as Drosophila gene c;assifications, can be browsed and searched for free. Due to outstanding contractual obligations, access to human gene classifications and to protein family trees and multiple sequence alignments will temporarily require a nominal registration fee. PANTHER is publicly available on the web at http://panther.celera.com.

The genome of the polar eukaryotic microalga Coccomyxa subellipsoidea reveals traits of cold adaptation

BackgroundLittle is known about the mechanisms of adaptation of life to the extreme environmental conditions encountered in polar regions. Here we present the genome sequence of a unicellular green alga from the division chlorophyta, Coccomyxa subellipsoidea C-169, which we will hereafter refer to as C-169. This is the first eukaryotic microorganism from a polar environment to have its genome sequenced.ResultsThe 48.8 Mb genome contained in 20 chromosomes exhibits significant synteny conservation with the chromosomes of its relatives Chlorella variabilis and Chlamydomonas reinhardtii. The order of the genes is highly reshuffled within synteny blocks, suggesting that intra-chromosomal rearrangements were more prevalent than inter-chromosomal rearrangements. Remarkably, Zepp retrotransposons occur in clusters of nested elements with strictly one cluster per chromosome probably residing at the centromere. Several protein families overrepresented in C. subellipsoidae include proteins involved in lipid metabolism, transporters, cellulose synthases and short alcohol dehydrogenases. Conversely, C-169 lacks proteins that exist in all other sequenced chlorophytes, including components of the glycosyl phosphatidyl inositol anchoring system, pyruvate phosphate dikinase and the photosystem 1 reaction center subunit N (PsaN).ConclusionsWe suggest that some of these gene losses and gains could have contributed to adaptation to low temperatures. Comparison of these genomic features with the adaptive strategies of psychrophilic microbes suggests that prokaryotes and eukaryotes followed comparable evolutionary routes to adapt to cold environments.

The Highly Similar Arabidopsis Homologs of Trithorax ATX1 and ATX2 Encode Proteins with Divergent Biochemical Functions

Gene duplication followed by functional specialization is a potent force in the evolution of biological diversity. A comparative study of two highly conserved duplicated genes, ARABIDOPSIS TRITHORAX-LIKE PROTEIN1 (ATX1) and ATX2, revealed features of both partial redundancy and of functional divergence. Although structurally similar, their regulatory sequences have diverged, resulting in distinct temporal and spatial patterns of expression of the ATX1 and ATX2 genes. We found that ATX2 methylates only a limited fraction of nucleosomes and that ATX1 and ATX2 influence the expression of largely nonoverlapping gene sets. Even when coregulating shared targets, ATX1 and ATX2 may employ different mechanisms. Most remarkable is the divergence of their biochemical activities: both proteins methylate K4 of histone H3, but while ATX1 trimethylates it, ATX2 dimethylates it. ATX2 and ATX1 provide an example of separated K4 di from K4 trimethyltransferase activity.

Practical Guidelines for the Comprehensive Analysis of ChIP-seq Data

Mapping the chromosomal locations of transcription factors, nucleosomes, histone modifications, chromatin remodeling enzymes, chaperones, and polymerases is one of the key tasks of modern biology, as evidenced by the Encyclopedia of DNA Elements (ENCODE) Project. To this end, chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) is the standard methodology. Mapping such protein-DNA interactions in vivo using ChIP-seq presents multiple challenges not only in sample preparation and sequencing but also for computational analysis. Here, we present step-by-step guidelines for the computational analysis of ChIP-seq data. We address all the major steps in the analysis of ChIP-seq data: sequencing depth selection, quality checking, mapping, data normalization, assessment of reproducibility, peak calling, differential binding analysis, controlling the false discovery rate, peak annotation, visualization, and motif analysis. At each step in our guidelines we discuss some of the software tools most frequently used. We also highlight the challenges and problems associated with each step in ChIP-seq data analysis. We present a concise workflow for the analysis of ChIP-seq data in Figure 1 that complements and expands on the recommendations of the ENCODE and modENCODE projects. Each step in the workflow is described in detail in the following sections.

Transcriptome-Wide Changes in Chlamydomonas reinhardtii Gene Expression Regulated by Carbon Dioxide and the CO2-Concentrating Mechanism Regulator CIA5/CCM1

Using powerful tools to investigate the regulation of global gene expression in the model microalga Chlamydomonas reinhardtii, we observed an impact of CO2 and CIA5, a key transcription regulator, on expression of almost 25% of all the genes. We also discovered an array of gene clusters with distinctive expression patterns that provide insight into the regulatory interaction between CIA5 and CO2. We used RNA sequencing to query the Chlamydomonas reinhardtii transcriptome for regulation by CO2 and by the transcription regulator CIA5 (CCM1). Both CO2 and CIA5 are known to play roles in acclimation to low CO2 and in induction of an essential CO2-concentrating mechanism (CCM), but less is known about their interaction and impact on the whole transcriptome. Our comparison of the transcriptome of a wild type versus a cia5 mutant strain under three different CO2 conditions, high CO2 (5%), low CO2 (0.03 to 0.05%), and very low CO2 (<0.02%), provided an entry into global changes in the gene expression patterns occurring in response to the interaction between CO2 and CIA5. We observed a massive impact of CIA5 and CO2 on the transcriptome, affecting almost 25% of all Chlamydomonas genes, and we discovered an array of gene clusters with distinctive expression patterns that provide insight into the regulatory interaction between CIA5 and CO2. Several individual clusters respond primarily to either CIA5 or CO2, providing access to genes regulated by one factor but decoupled from the other. Three distinct clusters clearly associated with CCM-related genes may represent a rich source of candidates for new CCM components, including a small cluster of genes encoding putative inorganic carbon transporters.

The sigma-1 receptor enhances brain plasticity and functional recovery after experimental stroke

Stroke leads to brain damage with subsequent slow and incomplete recovery of lost brain functions. Enriched housing of stroke-injured rats provides multi-modal sensorimotor stimulation, which improves recovery, although the specific mechanisms involved have not been identified. In rats housed in an enriched environment for two weeks after permanent middle cerebral artery occlusion, we found increased sigma-1 receptor expression in peri-infarct areas. Treatment of rats subjected to permanent or transient middle cerebral artery occlusion with 1-(3,4-dimethoxyphenethyl)-4-(3-phenylpropyl)piperazine dihydrochloride, an agonist of the sigma-1 receptor, starting two days after injury, enhanced the recovery of lost sensorimotor function without decreasing infarct size. The sigma-1 receptor was found in the galactocerebroside enriched membrane microdomains of reactive astrocytes and in neurons. Sigma-1 receptor activation increased the levels of the synaptic protein neurabin and neurexin in membrane rafts in the peri-infarct area, while sigma-1 receptor silencing prevented sigma-1 receptor-mediated neurite outgrowth in primary cortical neuronal cultures. In astrocytic cultures, oxygen and glucose deprivation induced sigma-1 receptor expression and actin dependent membrane raft formation, the latter blocked by sigma-1 receptor small interfering RNA silencing and pharmacological inhibition. We conclude that sigma-1 receptor activation stimulates recovery after stroke by enhancing cellular transport of biomolecules required for brain repair, thereby stimulating brain plasticity. Pharmacological targeting of the sigma-1 receptor provides new opportunities for stroke treatment beyond the therapeutic window of neuroprotection.

Activation of the Carbon Concentrating Mechanism by CO2 Deprivation Coincides with Massive Transcriptional Restructuring in Chlamydomonas reinhardtii

The CO2-concentrating mechanism (CCM) is key to allowing robust growth of eukaryotic algae. CCM activation by exposure of cells to low CO2 accompanies induction of genes involved in inorganic carbon uptake and concomitant suppression of many metabolism-related genes. Transcriptome and promoter motif analyses provide insights into CCM-associated regulatory networks. A CO2-concentrating mechanism (CCM) is essential for the growth of most eukaryotic algae under ambient (392 ppm) and very low (<100 ppm) CO2 concentrations. In this study, we used replicated deep mRNA sequencing and regulatory network reconstruction to capture a remarkable scope of changes in gene expression that occurs when Chlamydomonas reinhardtii cells are shifted from high to very low levels of CO2 (≤100 ppm). CCM induction 30 to 180 min post-CO2 deprivation coincides with statistically significant changes in the expression of an astonishing 38% (5884) of the 15,501 nonoverlapping C. reinhardtii genes. Of these genes, 1088 genes were induced and 3828 genes were downregulated by a log2 factor of 2. The latter indicate a global reduction in photosynthesis, protein synthesis, and energy-related biochemical pathways. The magnitude of transcriptional rearrangement and its major patterns are robust as analyzed by three different statistical methods. De novo DNA motif discovery revealed new putative binding sites for Myeloid oncogene family transcription factors potentially involved in activating low CO2–induced genes. The (CA)n repeat (9 ≤ n ≤ 25) is present in 29% of upregulated genes but almost absent from promoters of downregulated genes. These discoveries open many avenues for new research.

RNAi screening reveals requirement for host cell secretory pathway in infection by diverse families of negative-strand RNA viruses

Negative-strand (NS) RNA viruses comprise many pathogens that cause serious diseases in humans and animals. Despite their clinical importance, little is known about the host factors required for their infection. Using vesicular stomatitis virus (VSV), a prototypic NS RNA virus in the family Rhabdoviridae, we conducted a human genome-wide siRNA screen and identified 72 host genes required for viral infection. Many of these identified genes were also required for infection by two other NS RNA viruses, the lymphocytic choriomeningitis virus of the Arenaviridae family and human parainfluenza virus type 3 of the Paramyxoviridae family. Genes affecting different stages of VSV infection, such as entry/uncoating, gene expression, and assembly/release, were identified. Depletion of the proteins of the coatomer complex I or its upstream effectors ARF1 or GBF1 led to detection of reduced levels of VSV RNA. Coatomer complex I was also required for infection of lymphocytic choriomeningitis virus and human parainfluenza virus type 3. These results highlight the evolutionarily conserved requirements for gene expression of diverse families of NS RNA viruses and demonstrate the involvement of host cell secretory pathway in the process.

More complete gene silencing by fewer siRNAs: transparent optimized design and biophysical signature

Highly accurate knockdown functional analyses based on RNA interference (RNAi) require the possible most complete hydrolysis of the targeted mRNA while avoiding the degradation of untargeted genes (off-target effects). This in turn requires significant improvements to target selection for two reasons. First, the average silencing activity of randomly selected siRNAs is as low as 62%. Second, applying more than five different siRNAs may lead to saturation of the RNA-induced silencing complex (RISC) and to the degradation of untargeted genes. Therefore, selecting a small number of highly active siRNAs is critical for maximizing knockdown and minimizing off-target effects. To satisfy these needs, a publicly available and transparent machine learning tool is presented that ranks all possible siRNAs for each targeted gene. Support vector machines (SVMs) with polynomial kernels and constrained optimization models select and utilize the most predictive effective combinations from 572 sequence, thermodynamic, accessibility and self-hairpin features over 2200 published siRNAs. This tool reaches an accuracy of 92.3% in cross-validation experiments. We fully present the underlying biophysical signature that involves free energy, accessibility and dinucleotide characteristics. We show that while complete silencing is possible at certain structured target sites, accessibility information improves the prediction of the 90% active siRNA target sites. Fast siRNA activity predictions can be performed on our web server at .

Transcriptional profiling of the sperm storage organs of Drosophila melanogaster.

The occurrence of female sperm storage across taxa indicates the importance of this complex and dynamic process. Organs responsible for sperm storage (SSOs) and proteins expressed therein, are important in fundamental aspects of reproduction and could play a major role in evolutionary processes such as post‐mating sexual selection. Given the essential role of SSOs, it is surprising that the process of sperm storage is so poorly understood. This study investigated the transcriptome of female Drosophila melanogaster SSOs (seminal receptacle and spermathecae). Spermathecae were enriched for proteases and metabolic enzymes while the seminal receptacle was enriched for genes involved in localization, signaling and ion transport. Differences in functional gene categories indicate that these organs play unique roles in sperm storage.