Transcriptional landscape of human keratinocyte models exposed to 60-GHz millimeter-waves

The use of millimeter waves (MMW) will exponentially grow in the coming years due to their future utilization in 5G/6G networks. The question of possible biological effects at these frequencies has been raised. In this present study, we aimed to investigate gene expression changes under exposure to MMW using the Bulk RNA Barcoding and sequencing (BRB-seq) technology. To address this issue, three exposure scenarios were performed aiming at: i) comparing the cellular response of two primary culture of keratinocytes (HEK and NHEK) and one keratinocyte derivate cell line (HaCaT) exposed to MMW; ii) exploring the incident power density dose-effect on gene expression in HaCaT cell line; and, iii) studying the exposure duration at the new ICNIRP exposure limit for the general population. With the exception of heat effect induced by high power MMW (over 10 mW/cm2), those exposure scenarios have not enabled us to demonstrate important gene expression changes in the different cell populations studied. Very few differentially genes were observed between MMW exposed samples and heat shock control, and most of them were significantly associated with heat shock response that may reflect small differences in the heat generation. Together these results show that acute exposure to MMW has no effects on the transcriptional landscape of human keratinocyte models under athermal conditions.

Multimodal cartography of human lymphopoiesis reveals B and T/NK/ILC lineages are subjected to differential regulation

The developmental cartography of human lymphopoiesis remains incompletely understood. Here, we establish a multimodal map demonstrating that lymphoid specification follows independent direct or stepwise hierarchic routes converging toward the emergence of newly characterized CD117lo multi-lymphoid progenitors (MLPs) that undergo a proliferation arrest before entering the CD127- (NK/ILC/T) or CD127+ (B) lymphoid pathways. While the differentiation of CD127- early lymphoid progenitors is mainly driven by Flt3 signaling, emergence of their CD127+ counterparts is regulated cell-intrinsically and depends exclusively on the divisional history of their upstream precursors, including hematopoietic stem cells. Further, transcriptional mapping of differentiation trajectories reveals that whereas myeloid granulomonocytic lineages follow continuous differentiation pathways, lymphoid trajectories are intrinsically discontinuous and characterized by sequential waves of cell proliferation allowing pre-commitment amplification of lymphoid progenitor pools. Besides identifying new lymphoid specification pathways and regulatory checkpoints, our results demonstrate that NK/ILC/T and B lineages are under fundamentally distinct modes of regulation. (149 words).

Expression of the endocannabinoid system and response to cannabinoid components by the human fetal testis

BACKGROUND

Cannabis consumption by pregnant women continues to increase worldwide, raising concerns about adverse effects on fetal growth and deleterious impacts on the newborn, in connection with evidence of placental transfer of cannabis compound. Cannabis action is mediated by the endocannabinoid system (ECS), which expression is well established in the brain but unknown in the developing testis. The fetal testis, whose endocrine function orchestrates the masculinization of many distant organs, is particularly sensitive to disruption by xenobiotics. In this context, we aimed to determine whether cannabis exposure has the potential to directly impact the human fetal testis.

METHODS

We determined the expression of components of the ECS in the human fetal testis from 6 to 17 developmental weeks and assessed the direct effects of phytocannabinoids Δ9-trans-tetrahydrocannabinol (THC) and cannabidiol (CBD) on the testis morphology and cell functions ex vivo.

RESULTS

We demonstrate the presence in the human fetal testis of two key endocannabinoids, 2-arachidonylglycerol (2-AG) and to a lower level anandamide (AEA), as well as a range of enzymes and receptors for the ECS. Ex vivo exposure of first trimester testes to CBD, THC, or CBD/THC [ratio 1:1] at 10^-7 to 10^-5 M altered testosterone secretion by Leydig cells, AMH secretion by Sertoli cells, and impacted testicular cell proliferation and viability as early as 72 h post-exposure. Transcriptomic analysis on 72 h-exposed fetal testis explants revealed 187 differentially expressed genes (DEGs), including genes involved in steroid synthesis and toxic substance response. Depending on the molecules and testis age, highly deleterious effects of phytocannabinoid exposure were observed on testis tissue after 14 days, including Sertoli and germ cell death.

CONCLUSIONS

Our study is the first to evidence the presence of the ECS in the human fetal testis and to highlight the potential adverse effect of cannabis consumption by pregnant women onto the development of the male gonad.

Single-cell Deconvolution of a Specific Malignant Cell Population as a Poor Prognostic Biomarker in Low-risk Clear Cell Renal Cell Carcinoma Patients

BACKGROUND

Intratumor heterogeneity (ITH) is a key feature in clear cell renal cell carcinomas (ccRCCs) that impacts outcomes such as aggressiveness, response to treatments, or recurrence. In particular, it may explain tumor relapse after surgery in clinically low-risk patients who did not benefit from adjuvant therapy. Recently, single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool to unravel expression ITH (eITH) and might enable better assessment of clinical outcomes in ccRCC.

OBJECTIVE

To explore eITH in ccRCC with a focus on malignant cells (MCs) and assess its relevance to improve prognosis for low-risk patients.

DESIGN, SETTING, AND PARTICIPANTS

We performed scRNA-seq on tumor samples from five untreated ccRCC patients ranging from pT1a to pT3b. Data were complemented with a published dataset composed of pairs of matched normal and ccRCC samples.

INTERVENTION

Radical or partial nephrectomy on untreated ccRCC patients.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Viability and cell type proportions were determined by flow cytometry. Following scRNA-seq, a functional analysis was performed and tumor progression trajectories were inferred. A deconvolution approach was applied on an external cohort, and Kaplan-Meier survival curves were estimated with respect to the prevalence of malignant clusters.

RESULTS AND LIMITATIONS

We analyzed 54 812 cells and identified 35 cell subpopulations. The eITH analysis revealed that each tumor contained various degrees of clonal diversity. The transcriptomic signatures of MCs in one particularly heterogeneous sample were used to design a deconvolution-based strategy that allowed the risk stratification of 310 low-risk ccRCC patients.

CONCLUSIONS

We described eITH in ccRCCs, and used this information to establish significant cell population-based prognostic signatures and better discriminate ccRCC patients. This approach has the potential to improve the stratification of clinically low-risk patients and their therapeutic management.

PATIENT SUMMARY

We sequenced the RNA content of individual cell subpopulations composed of clear cell renal cell carcinomas and identified specific malignant cells the genetic information of which can be used to predict tumor progression.

The anti-androgenic fungicide triticonazole induces region-specific transcriptional changes in the developing rat perineum and phallus

Intrauterine exposure to endocrine disrupting chemicals can interfere with male reproductive development. This can lead to male reproductive disorders such as hypospadias, cryptorchidism and reduced fertility, as well as shorter anogenital distance (AGD) - a biomarker for incomplete androgen-dependent fetal masculinization. However, it remains challenging to predict adverse in vivo outcomes based on in vitro effect patterns for many chemicals. This is a challenge for modern toxicology, which aims to reduce animal testing for chemical safety assessments. To enable the transition towards higher reliance on alternative test methods, we need to better map underlying mechanisms leading to adverse effects. Herein, we have analyzed the transcriptome of the perineum and phallus of male fetal rats and defined the impacts of exposure to an anti-androgenic fungicide, triticonazole. Previously we have shown that developmental exposure to triticonazole can induce short male AGD, but without a marked effect on the transcriptome of the fetal testes. In contrast, we report here significant changes to the transcriptional landscape of the perineum and phallus, including regional differences between these adjacent tissues. This highlights the importance of analyzing the correct tissue when characterizing mechanisms of complex in vivo effect outcomes. Our results provide a rich resource for the spatiotemporal gene networks that are involved in the development of male external genitalia, and that can be disrupted upon exposure to chemicals that prevent normal masculinization of the perineum and phallus. Such data will be critical in the development of novel alternative test methods to determine the endocrine disrupting potential of existing and emerging chemicals.

FlexDotPlot: a universal and modular dot plot visualization tool for complex multifaceted data

Motivation

Dot plots are heatmap-like charts that provide a compact way to simultaneously display two quantitative information by means of dots of different sizes and colors. Despite the popularity of this visualization method, particularly in single-cell RNA-sequencing (scRNA-seq) studies, existing tools used to make dot plots are limited in terms of functionality and usability.

Results

We developed FlexDotPlot, an R package for generating dot plots from multifaceted data, including scRNA-seq data. It provides a universal and easy-to-use solution with a high versatility. An interactive R Shiny application is also available allowing non-R users to easily generate dot plots with several tunable parameters.

Availability and implementation

Source code and detailed manual are available on CRAN (stable version) and at https://github.com/Simon-Leonard/FlexDotPlot (development version). Code to reproduce figures is available at https://github.com/Simon-Leonard/FlexDotPlot_paper. A Shiny app is available as a stand-alone application within the package.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Transcriptome analysis of fetal rat testis following intrauterine exposure to the azole fungicides triticonazole and flusilazole reveals subtle changes despite adverse endocrine effects

Azoles are used in agriculture and medicine to combat fungal infections. We have previously examined the endocrine disrupting properties of the agricultural azole fungicides triticonazole and flusilazole. Triticonazole displayed strong androgen receptor (AR) antagonism in vitro, whereas in utero exposure resulted in anti-androgenic effects in vivo evidenced by shorter anogenital distance (AGD) in fetal male rats. Flusilazole displayed strong AR antagonism, but less potent than triticonazole, and disrupted steroidogenesis in vitro, whereas in utero exposure disrupted fetal male plasma hormone levels. To elaborate on how these azole fungicides can disrupt male reproductive development by different mechanisms, and to investigate whether feminization effects such as short AGD in males can also be detected at the transcript level in fetal testes, we profiled fetal testis transcriptomes after in utero exposure to triticonazole and flusilazole by 3'Digital Gene Expression (3'DGE). The analysis revealed few transcriptional changes after exposure to either compound at gestation day 17 and 21. This suggests that the observed influence of flusilazole on hormone production may be by directly targeting steroidogenic enzyme activity in the testis at the protein level, whereas observations of shorter AGD by triticonazole may primarily be due to disturbed androgen signaling in androgen-sensitive tissues. Expression of Calb2 and Gsta2 was altered by flusilazole but not triticonazole and may pinpoint novel pathways of disrupted testicular steroid synthesis. Our findings have wider implication for how we integrate omics data in chemical testing frameworks, including selection of non-animal test methods and building of Adverse Outcome Pathways for regulatory purposes.

The ReproGenomics Viewer: a multi-omics and cross-species resource compatible with single-cell studies for the reproductive science community

MOTIVATION

Recent advances in transcriptomics have enabled unprecedented insight into gene expression analysis at a single-cell resolution. While it is anticipated that the number of publications based on such technologies will increase in the next decade, there is currently no public resource to centralize and enable scientists to explore single-cell datasets published in the field of reproductive biology.

RESULTS

Here, we present a major update of the ReproGenomics Viewer, a cross-species and cross-technology web-based resource of manually-curated sequencing datasets related to reproduction. The redesign of the ReproGenomics Viewer's architecture is accompanied by significant growth of the database content including several landmark single-cell RNA-sequencing datasets. The implementation of additional tools enables users to visualize and browse the complex, high-dimensional data now being generated in the reproductive field.

AVAILABILITY AND IMPLEMENTATION

The ReproGenomics Viewer resource is freely accessible at http://rgv.genouest.org. The website is implemented in Python, JavaScript and MongoDB, and is compatible with all major browsers. Source codes can be downloaded from https://github.com/fchalmel/RGV.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Characterization of Glyceollins as Novel Aryl Hydrocarbon Receptor Ligands and Their Role in Cell Migration

Recent studies strongly support the use of the aryl hydrocarbon receptor (AhR) as a therapeutic target in breast cancer. Glyceollins, a group of soybean phytoalexins, are known to exert therapeutic effects in chronic human diseases and also in cancer. To investigate the interaction between glyceollin I (GI), glyceollin II (GII) and AhR, a computational docking analysis, luciferase assays, immunofluorescence and transcriptome analyses were performed with different cancer cell lines. The docking experiments predicted that GI and GII can enter into the AhR binding pocket, but their interactions with the amino acids of the binding site differ, in part, from those interacting with 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Both GI and GII were able to weakly and partially activate AhR, with GII being more potent. The results from the transcriptome assays showed that approximately 10% of the genes regulated by TCDD were also modified by both GI and GII, which could have either antagonistic or synergistic effects upon TCDD activation. In addition, we report here, on the basis of phenotype, that GI and GII inhibit the migration of triple-negative (ER-, PgR-, HER2NEU-) MDA-MB-231 breast cancer cells, and that they inhibit the expression of genes which code for important regulators of cell migration and invasion in cancer tissues. In conclusion, GI and GII are AhR ligands that should be further investigated to determine their usefulness in cancer treatments.

Vascular Delivery of Allogeneic MuStem Cells in Dystrophic Dogs Requires Only Short-Term Immunosuppression to Avoid Host Immunity and Generate Clinical/Tissue Benefits

Growing demonstrations of regenerative potential for some stem cells led recently to promising therapeutic proposals for neuromuscular diseases. We have shown that allogeneic MuStem cell transplantation into Golden Retriever muscular dystrophy (GRMD) dogs under continuous immunosuppression (IS) leads to persistent clinical stabilization and muscle repair. However, long-term IS in medical practice is associated with adverse effects raising safety concerns. Here, we investigate whether the IS removal or its restriction to the transplantation period could be considered. Dogs aged 4-5 months old received vascular infusions of allogeneic MuStem cells without IS (GRMDMU/no-IS) or under transient IS (GRMDMU/tr-IS). At 5 months post-infusion, persisting clinical status improvement of the GRMDMU/tr-IS dogs was observed while GRMDMU/no-IS dogs exhibited no benefit. Histologically, only 9-month-old GRMDMU/tr-IS dogs showed an increased muscle regenerative activity. A mixed cell reaction with the host peripheral blood mononucleated cells (PBMCs) and corresponding donor cells revealed undetectable to weak lymphocyte proliferation in GRMDMU/tr-IS dogs compared with a significant proliferation in GRMDMU/no-IS dogs. Importantly, any dog group showed neither cellular nor humoral anti-dystrophin responses. Our results show that transient IS is necessary and sufficient to sustain allogeneic MuStem cell transplantation benefits and prevent host immunity. These findings provide useful critical insight to designing therapeutic strategies.

Human serum and platelet lysate are appropriate xeno-free alternatives for clinical-grade production of human MuStem cell batches

BACKGROUND

Canine MuStem cells have demonstrated regenerative efficacy in a dog model of muscular dystrophy, and the recent characterization of human counterparts (hMuStem) has highlighted the therapeutic potential of this muscle-derived stem cell population. To date, these cells have only been generated in research-grade conditions. However, evaluation of the clinical efficacy of any such therapy will require the production of hMuStem cells in compliance with good manufacturing practices (GMPs). Because the current use of fetal bovine serum (FBS) to isolate and expand hMuStem cells raises several ethical, safety, and supply concerns, we assessed the use of two alternative xeno-free blood derivatives: human serum (HS) and a human platelet lysate (hPL).

METHODS

hMuStem cells were isolated and expanded in vitro in either HS-supplemented or hPL-supplemented media and the proliferation rate, clonogenicity, myogenic commitment potential, and oligopotency compared with that observed in FBS-supplemented medium. Flow cytometry and high-throughput 3'-digital gene expression RNA sequencing were used to characterize the phenotype and global gene expression pattern of hMuStem cells cultured with HS or hPL.

RESULTS

HS-supplemented and hPL-supplemented media both supported the isolation and long-term proliferation of hMuStem cells. Compared with FBS-based medium, both supplements enhanced clonogenicity and allowed for a reduction in growth factor supplementation. Neither supplement altered the cell lineage pattern of hMuStem cells. In vitro differentiation assays revealed a decrease in myogenic commitment and in the fusion ability of hMuStem cells when cultured with hPL. In return, this reduction of myogenic potential in hPL-supplemented cultures was rapidly reversed by substitution of hPL with HS or fibrinogen-depleted hPL. Moreover, culture of hMuStem cells in hPL hydrogel and fibrinogen-depleted hPL demonstrated that myogenic differentiation potential is maintained in heparin-free hPL derivatives.

CONCLUSIONS

Our findings indicate that HS and hPL are efficient and viable alternatives to FBS for the preparation of hMuStem cell batches in compliance with GMPs.

Quantitative proteome profiling of dystrophic dog skeletal muscle reveals a stabilized muscular architecture and protection against oxidative stress after systemic delivery of MuStem cells

Proteomic profiling plays a decisive role in the elucidation of molecular signatures representative of a specific clinical context. MuStem cell based therapy represents a promising approach for clinical applications to cure Duchenne muscular dystrophy (DMD). To expand our previous studies collected in the clinically relevant DMD animal model, we decided to investigate the skeletal muscle proteome 4 months after systemic delivery of allogenic MuStem cells. Quantitative proteomics with isotope-coded protein labeling was used to compile quantitative changes in the protein expression profiles of muscle in transplanted Golden Retriever muscular dystrophy (GRMD) dogs as compared to Golden Retriever muscular dystrophy dogs. A total of 492 proteins were quantified, including 25 that were overrepresented and 46 that were underrepresented after MuStem cell transplantation. Interestingly, this study demonstrates that somatic stem cell therapy impacts on the structural integrity of the muscle fascicle by acting on fibers and its connections with the extracellular matrix. We also show that cell infusion promotes protective mechanisms against oxidative stress and favors the initial phase of muscle repair. This study allows us to identify putative candidates for tissue markers that might be of great value in objectively exploring the clinical benefits resulting from our cell-based therapy for DMD. All MS data have been deposited in the ProteomeXchange with identifier PXD001768 (http://proteomecentral.proteomexchange.org/dataset/PXD001768).

Differential Gene Expression Profiling of Dystrophic Dog Muscle after MuStem Cell Transplantation

BACKGROUND

Several adult stem cell populations exhibit myogenic regenerative potential, thus representing attractive candidates for therapeutic approaches of neuromuscular diseases such as Duchenne Muscular Dystrophy (DMD). We have recently shown that systemic delivery of MuStem cells, skeletal muscle-resident stem cells isolated in healthy dog, generates the remodelling of muscle tissue and gives rise to striking clinical benefits in Golden Retriever Muscular Dystrophy (GRMD) dog. This global effect, which is observed in the clinically relevant DMD animal model, leads us to question here the molecular pathways that are impacted by MuStem cell transplantation. To address this issue, we compare the global gene expression profile between healthy, GRMD and MuStem cell treated GRMD dog muscle, four months after allogenic MuStem cell transplantation.

RESULTS

In the dystrophic context of the GRMD dog, disease-related deregulation is observed in the case of 282 genes related to various processes such as inflammatory response, regeneration, calcium ion binding, extracellular matrix organization, metabolism and apoptosis regulation. Importantly, we reveal the impact of MuStem cell transplantation on several molecular and cellular pathways based on a selection of 31 genes displaying signals specifically modulated by the treatment. Concomitant with a diffuse dystrophin expression, a histological remodelling and a stabilization of GRMD dog clinical status, we show that cell delivery is associated with an up-regulation of genes reflecting a sustained enhancement of muscle regeneration. We also identify a decreased mRNA expression of a set of genes having metabolic functions associated with lipid homeostasis and energy. Interestingly, ubiquitin-mediated protein degradation is highly enhanced in GRMD dog muscle after systemic delivery of MuStem cells.

CONCLUSIONS

Overall, our results provide the first high-throughput characterization of GRMD dog muscle and throw new light on the complex molecular/cellular effects associated with muscle repair and the clinical efficacy of MuStem cell-based therapy.

Integrated RNA- and protein profiling of fermentation and respiration in diploid budding yeast provides insight into nutrient control of cell growth and development

Diploid budding yeast undergoes rapid mitosis when it ferments glucose, and in the presence of a non-fermentable carbon source and the absence of a nitrogen source it triggers sporulation. Rich medium with acetate is a commonly used pre-sporulation medium, but our understanding of the molecular events underlying the acetate-driven transition from mitosis to meiosis is still incomplete. We identified 263 proteins for which mRNA and protein synthesis are linked or uncoupled in fermenting and respiring cells. Using motif predictions, interaction data and RNA profiling we find among them 28 likely targets for Ume6, a subunit of the conserved Rpd3/Sin3 histone deacetylase-complex regulating genes involved in metabolism, stress response and meiosis. Finally, we identify 14 genes for which both RNA and proteins are detected exclusively in respiring cells but not in fermenting cells in our sample set, including CSM4, SPR1, SPS4 and RIM4, which were thought to be meiosis-specific. Our work reveals intertwined transcriptional and post-transcriptional control mechanisms acting when a MATa/α strain responds to nutritional signals, and provides molecular clues how the carbon source primes yeast cells for entering meiosis.

BIOLOGICAL SIGNIFICANCE

Our integrated genomics study provides insight into the interplay between the transcriptome and the proteome in diploid yeast cells undergoing vegetative growth in the presence of glucose (fermentation) or acetate (respiration). Furthermore, it reveals novel target genes involved in these processes for Ume6, the DNA binding subunit of the conserved histone deacetylase Rpd3 and the co-repressor Sin3. We have combined data from an RNA profiling experiment using tiling arrays that cover the entire yeast genome, and a large-scale protein detection analysis based on mass spectrometry in diploid MATa/α cells. This distinguishes our study from most others in the field-which investigate haploid yeast strains-because only diploid cells can undergo meiotic development in the simultaneous absence of a non-fermentable carbon source and nitrogen. Indeed, we report molecular clues how respiration of acetate might prime diploid cells for efficient spore formation, a phenomenon that is well known but poorly understood.

The conserved histone deacetylase Rpd3 and its DNA binding subunit Ume6 control dynamic transcript architecture during mitotic growth and meiotic development

It was recently reported that the sizes of many mRNAs change when budding yeast cells exit mitosis and enter the meiotic differentiation pathway. These differences were attributed to length variations of their untranslated regions. The function of UTRs in protein translation is well established. However, the mechanism controlling the expression of distinct transcript isoforms during mitotic growth and meiotic development is unknown. In this study, we order developmentally regulated transcript isoforms according to their expression at specific stages during meiosis and gametogenesis, as compared to vegetative growth and starvation. We employ regulatory motif prediction, in vivo protein-DNA binding assays, genetic analyses and monitoring of epigenetic amino acid modification patterns to identify a novel role for Rpd3 and Ume6, two components of a histone deacetylase complex already known to repress early meiosis-specific genes in dividing cells, in mitotic repression of meiosis-specific transcript isoforms. Our findings classify developmental stage-specific early, middle and late meiotic transcript isoforms, and they point to a novel HDAC-dependent control mechanism for flexible transcript architecture during cell growth and differentiation. Since Rpd3 is highly conserved and ubiquitously expressed in many tissues, our results are likely relevant for development and disease in higher eukaryotes.

Developmental stage dependent metabolic regulation during meiotic differentiation in budding yeast

Background

The meiotic developmental pathway in yeast enables both differentiation of vegetative cells into haploid spores that ensure long-term survival, and recombination of the parental DNA to create genetic diversity. Despite the importance of proper metabolic regulation for the supply of building blocks and energy, little is known about the reprogramming of central metabolic pathways in meiotically differentiating cells during passage through successive developmental stages.

Results

Metabolic regulation during meiotic differentiation in budding yeast was analyzed by integrating information on genome-wide transcriptional activity, 26 enzymatic activities in the central metabolism, the dynamics of 67 metabolites, and a metabolic flux analysis at mid-stage meiosis. Analyses of mutants arresting sporulation at defined stages demonstrated that metabolic reprogramming is tightly controlled by the progression through the developmental pathway. The correlation between transcript levels and enzymatic activities in the central metabolism varies significantly in a developmental stage-dependent manner. The complete loss of phosphofructokinase activity at mid-stage meiosis enables a unique setup of the glycolytic pathway which facilitates carbon flux repartitioning into synthesis of spore wall precursors during the co-assimilation of glycogen and acetate. The need for correct homeostasis of purine nucleotides during the meiotic differentiation was demonstrated by the sporulation defect of the AMP deaminase mutant amd1, which exhibited hyper-accumulation of ATP accompanied by depletion of guanosine nucleotides.

Conclusions

Our systems-level analysis shows that reprogramming of the central metabolism during the meiotic differentiation is controlled at different hierarchical levels to meet the metabolic and energetic needs at successive developmental stages.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-014-0060-x) contains supplementary material, which is available to authorized users.

High-resolution profiling of novel transcribed regions during rat spermatogenesis

Mammalian spermatogenesis is a complex and highly orchestrated combination of processes in which male germline proliferation and differentiation result in the production of mature spermatozoa. If recent genome-wide studies have contributed to the in-depth analysis of the male germline protein-encoding transcriptome, little effort has yet been devoted to the systematic identification of novel unannotated transcribed regions expressed during mammalian spermatogenesis. We report high-resolution expression profiling of male germ cells in rat, using next-generation sequencing technology and highly enriched testicular cell populations. Among 20 424 high-confidence transcripts reconstructed, we defined a stringent set of 1419 long multi-exonic unannotated transcripts expressed in the testis (testis-expressed unannotated transcripts [TUTs]). TUTs were divided into 7 groups with different expression patterns. Most TUTs share many of the characteristics of vertebrate long noncoding RNAs (lncRNAs). We also markedly reinforced the finding that TUTs and known lncRNAs accumulate during the meiotic and postmeiotic stages of spermatogenesis in mammals and that X-linked meiotic TUTs do not escape the silencing effects of meiotic sex chromosome inactivation. Importantly, we discovered that TUTs and known lncRNAs with a peak expression during meiosis define a distinct class of noncoding transcripts that exhibit exons twice as long as those of other transcripts. Our study provides new insights in transcriptional profiling of the male germline and represents a high-quality resource for novel loci expressed during spermatogenesis that significantly contributes to rat genome annotation.

An atlas of chromatoid body components

The genome of male germ cells is actively transcribed during spermatogenesis to produce phase-specific protein-coding mRNAs and a considerable amount of different noncoding RNAs. Ribonucleoprotein (RNP) granule-mediated RNA regulation provides a powerful means to secure the quality and correct expression of the requisite transcripts. Haploid spermatids are characterized by a unique, unusually large cytoplasmic granule, the chromatoid body (CB), which emerges during the switch between the meiotic and post-meiotic phases of spermatogenesis. To better understand the role of the CB in male germ cell differentiation, we isolated CBs from mouse testes and revealed its full RNA and protein composition. We showed that the CB is mainly composed of RNA-binding proteins and other proteins involved RNA regulation. The CB was loaded with RNA, including pachytene piRNAs, a diverse set of mRNAs, and a number of uncharacterized long noncoding transcripts. The CB was demonstrated to accumulate nascent RNA during all the steps of round spermatid differentiation. Our results revealed the CB as a large germ cell-specific RNP platform that is involved in the control of the highly complex transcriptome of haploid male germ cells.

Genome-wide identification of Sox8-, and Sox9-dependent genes during early post-natal testis development in the mouse

The SOX8 and SOX9 transcription factors are involved in, among others, sex differentiation, male gonad development and adult maintenance of spermatogenesis. Sox8(-/-) mice lacking Sox9 in Sertoli cells fail to form testis cords and cannot establish spermatogenesis. Although genetic and histological data show an important role for these transcription factors in regulating spermatogenesis, it is not clear which genes depend upon them at a genome-wide level. To identify transcripts that respond to the absence of Sox8 in all cells and Sox9 in Sertoli cells we measured mRNA concentrations in testicular samples from mice at 0, 6 and 18 days post-partum. In total, 621 and 629 transcripts were found at decreased or increased levels, respectively, at different time points in the mutant as compared to the control samples. These mRNAs were categorized as preferentially expressed in Sertoli cells or germ cells using data obtained with male and female gonad samples and enriched testicular cell populations. Five candidate genes were validated at the protein level. Furthermore, we identified putative direct SOX8 and SOX9 target genes by integrating predicted SOX-binding sites present in potential regulatory regions upstream of the transcription start site. Finally, we used protein network data to gain insight into the effects on regulatory interactions that occur when Sox8 and Sox9 are absent in developing Sertoli cells. The integration of testicular samples with enriched Sertoli cells, germ cells and female gonads enabled us to broadly distinguish transcripts directly affected in Sertoli cells from others that respond to secondary events in testicular cell types. Thus, combined RNA profiling signals, motif predictions and network data identified putative SOX8/SOX9 target genes in Sertoli cells and yielded insight into regulatory interactions that depend upon these transcription factors. In addition, our results will facilitate the interpretation of genome-wide in vivo SOX8 and SOX9 DNA binding data.

Profiling of androgen response in rainbow trout pubertal testis: relevance to male gonad development and spermatogenesis

The capacity of testicular somatic cells to promote and sustain germ cell differentiation is largely regulated by sexual steroids and notably androgens. In fish species the importance of androgens is emphasized by their ability to induce sex reversal of the developing fries and to trigger spermatogenesis. Here we studied the influence of androgens on testicular gene expression in trout testis using microarrays. Following treatment of immature males with physiological doses of testosterone or 11-ketotestosterone, 418 genes that exhibit changes in expression were identified. Interestingly, the activity of testosterone appeared stronger than that of 11-ketotestosterone. Expression profiles of responsive genes throughout testis development and in isolated germ cells confirmed androgens to mainly affect gene expression in somatic cells. Furthermore, specific clusters of genes that exhibit regulation coincidently with changes in the natural circulating levels of androgens during the reproductive cycle were highlighted, reinforcing the physiological significance of these data. Among somatic genes, a phylogenetic footprinting study identified putative androgen response elements within the proximal promoter regions of 42 potential direct androgen target genes. Finally, androgens were also found to alter the germ line towards meiotic expression profiles, supporting the hypothesis of a role for the somatic responsive genes in driving germ cell fate. This study significantly increases our understanding of molecular pathways regulated by androgens in vertebrates. The highly cyclic testicular development in trout together with functions associated with regulated genes reveal potential mechanisms for androgen actions in tubule formation, steroid production, germ cell development and sperm secretion.

Transcription of two long noncoding RNAs mediates mating-type control of gametogenesis in budding yeast

The cell-fate decision leading to gametogenesis is essential for sexual reproduction. In S. cerevisiae, only diploid MATa/α but not haploid MATa or MATα cells undergo gametogenesis, known as sporulation. We find that transcription of two long noncoding RNAs (lncRNAs) mediates mating-type control of sporulation. In MATa or MATα haploids, expression of IME1, the central inducer of gametogenesis, is inhibited in cis by transcription of the lncRNA IRT1, located in the IME1 promoter. IRT1 transcription recruits the Set2 histone methyltransferase and the Set3 histone deacetylase complex to establish repressive chromatin at the IME1 promoter. Inhibiting expression of IRT1 and an antisense transcript that antagonizes the expression of the meiotic regulator IME4 allows cells expressing the haploid mating type to sporulate with kinetics that are indistinguishable from that of MATa/α diploids. Conversely, expression of the two lncRNAs abolishes sporulation in MATa/α diploids. Thus, transcription of two lncRNAs governs mating-type control of gametogenesis in yeast.

Global human tissue profiling and protein network analysis reveals distinct levels of transcriptional germline-specificity and identifies target genes for male infertility

BACKGROUND

Mammalian spermatogenesis is a process that involves a complex expression program in both somatic and germ cells present in the male gonad. A number of studies have attempted to define the transcriptome of male meiosis and gametogenesis in rodents and primates. Few human transcripts, however, have been associated with testicular somatic cells and germ cells at different post-natal developmental stages and little is known about their level of germline-specificity compared with non-testicular tissues.

METHODS

We quantified human transcripts using GeneChips and a total of 47 biopsies from prepubertal children diagnosed with undescended testis, infertile adult patients whose spermatogenesis is arrested at consecutive stages and fertile control individuals. These results were integrated with data from enriched normal germ cells, non-testicular expression data, phenotype information, predicted regulatory DNA-binding motifs and interactome data.

RESULTS

Among 3580 genes for which we found differential transcript concentrations in somatic and germ cells present in human testis, 933 were undetectable in 45 embryonic and adult non-testicular tissues, including many that were corroborated at protein level by published gene annotation data and histological high-throughput protein immunodetection assays. Using motif enrichment analyses, we identified regulatory promoter elements likely involved in germline development. Finally, we constructed a regulatory disease network for human fertility by integrating expression signals, interactome information, phenotypes and functional annotation data.

CONCLUSIONS

Our results provide broad insight into the post-natal human testicular transcriptome at the level of cell populations and in a global somatic tissular context. Furthermore, they yield clues for genetic causes of male infertility and will facilitate the identification of novel cancer/testis genes as targets for cancer immunotherapies.

Direct iterative protein profiling (DIPP) - an innovative method for large-scale protein detection applied to budding yeast mitosis

The budding yeast Saccharomyces cerevisiae is a major model organism for important biological processes such as mitotic growth and meiotic development, it can be a human pathogen, and it is widely used in the food-, and biotechnology industries. Consequently, the genomes of numerous strains have been sequenced and a very large amount of RNA profiling data is available. Moreover, it has recently become possible to quantitatively analyze the entire yeast proteome; however, efficient and cost-effective high-throughput protein profiling remains a challenge. We report here a new approach to direct and label-free large-scale yeast protein identification using a tandem buffer system for protein extraction, two-step protein prefractionation and enzymatic digestion, and detection of peptides by iterative mass spectrometry. Our profiling study of diploid cells undergoing rapid mitotic growth identified 86% of the known proteins and its output was found to be widely concordant with genome-wide mRNA concentrations and DNA variations between yeast strains. This paves the way for comprehensive and straightforward yeast proteome profiling across a wide variety of experimental conditions.

Profiling spermatogenic failure in adult testes bearing Sox9-deficient Sertoli cells identifies genes involved in feminization, inflammation and stress

BACKGROUND

Sox9 (Sry box containing gene 9) is a DNA-binding transcription factor involved in chondrocyte development and sex determination. The protein's absence in testicular Sertoli nurse cells has been shown to disrupt testicular function in adults but little is known at the genome-wide level about molecular events concomitant with testicular break-down.

METHODS

To determine the genome-wide effect on mRNA concentrations triggered by the absence of Sox9 in Sertoli cells we analysed adult testicular tissue from wild-type versus mutant mice with high-density oligonucleotide microarrays and integrated the output of this experiment with regulatory motif predictions and protein-protein network data.

RESULTS

We report the genome-wide mRNA signature of adult testes lacking Sox9 in Sertoli cells before and after the onset of late spermatogenic failure as compared to fertile controls. The GeneChip data integrated with evolutionarily conserved Sox9 DNA binding motifs and regulatory network data identified genes involved in feminization, stress response and inflammation.

CONCLUSIONS

Our results extend previous observations that genes required for female gonadogenesis are up-regulated in the absence of Sox9 in fetal Sertoli cells to the adult stage. Importantly, we identify gene networks involved in immunological processes and stress response which is reminiscent of a phenomenon occurring in a sub-group of infertile men. This suggests mice lacking Sox9 in their Sertoli cells to be a potentially useful model for adult human testicular failure.

GermOnline 4.0 is a genomics gateway for germline development, meiosis and the mitotic cell cycle

GermOnline 4.0 is a cross-species database portal focusing on high-throughput expression data relevant for germline development, the meiotic cell cycle and mitosis in healthy versus malignant cells. It is thus a source of information for life scientists as well as clinicians who are interested in gene expression and regulatory networks. The GermOnline gateway provides unlimited access to information produced with high-density oligonucleotide microarrays (3'-UTR GeneChips), genome-wide protein-DNA binding assays and protein-protein interaction studies in the context of Ensembl genome annotation. Samples used to produce high-throughput expression data and to carry out genome-wide in vivo DNA binding assays are annotated via the MIAME-compliant Multiomics Information Management and Annotation System (MIMAS 3.0). Furthermore, the Saccharomyces Genomics Viewer (SGV) was developed and integrated into the gateway. SGV is a visualization tool that outputs genome annotation and DNA-strand specific expression data produced with high-density oligonucleotide tiling microarrays (Sc_tlg GeneChips) which cover the complete budding yeast genome on both DNA strands. It facilitates the interpretation of expression levels and transcript structures determined for various cell types cultured under different growth and differentiation conditions. Database URL: www.germonline.org/

Proteome analysis and genome-wide regulatory motif prediction identify novel potentially sex-hormone regulated proteins in rat efferent ducts

The efferent ducts are a series of tubules that conduct sperm from the rete testis to the epididymis. They absorb most fluid and proteins originating from the rete testis during concentration of spermatozoa prior to their entry into the epididymis. Proteome analysis of micro-dissected efferent duct samples from adult rats was combined with genome-wide computational prediction of conserved hormone response elements to identify factors likely regulated by oestrogens and androgens. We identified 165 proteins and found subsets of the promoters controlling their corresponding genes to contain androgen- and oestrogen response elements (ARE/EREs) at similar frequencies. Moreover, EREs were significantly enriched among the loci identified compared with their genome-wide occurrence. The expression and localization of Anxa6, Ckb, Krt19, Park7, Pdzk1 and Tpt1 in the efferent ducts and other related hormone controlled tissues was further validated at the RNA or protein level. This study identifies many novel proteins predicted to play roles in sperm maturation and male fertility and provides significant computational evidence that the efferent ducts express genes transcriptionally controlled by sex hormones.

Fhl5/Act, a CREM-binding transcriptional activator required for normal sperm maturation and morphology, is not essential for testicular gene expression

BACKGROUND

The LIM domain protein Fhl5 was previously found to interact with CREM, a DNA binding transcriptional regulator necessary for spermiogenesis in mammals. Co-transfection experiments using heterologous promoter constructs indicated a role for Fhl5 in transcriptional up-regulation of CREM-dependent testicular genes. Male mice lacking Fhl5 were reported to be fertile but displayed partially abnormal sperm maturation and morphology.

METHODS

To identify Fhl5 testicular target genes we carried out two whole-genome expression profiling experiments using high-density oligonucleotide microarrays and total testis samples from Fhl5 wild-type versus homozygous mutant mice first in different and then in isogenic strain backgrounds.

RESULTS

Weak signal differences were detected in non-isogenic samples but no statistically significant expression changes were observed when isogenic Fhl5 mutant and wild-type samples were compared.

CONCLUSION

The outcome of these experiments suggests that testicular expression profiling is extremely sensitive to the genetic background and that Fhl5 is not essential for testicular gene expression to a level detected by microarray-based measurements. This might be due to redundant function of the related and similarly expressed protein Fhl4.

The homeobox gene CHX10/VSX2 regulates RdCVF promoter activity in the inner retina

Rod-derived Cone Viability Factor (RdCVF) is a trophic factor with therapeutic potential for the treatment of retinitis pigmentosa, a retinal disease that commonly results in blindness. RdCVF is encoded by Nucleoredoxin-like 1 (Nxnl1), a gene homologous with the family of thioredoxins that participate in the defense against oxidative stress. RdCVF expression is lost after rod degeneration in the first phase of retinitis pigmentosa, and this loss has been implicated in the more clinically significant secondary cone degeneration that often occurs. Here, we describe a study of the Nxnl1 promoter using an approach that combines promoter and transcriptomic analysis. By transfection of selected candidate transcription factors, chosen based upon their expression pattern, we identified the homeodomain proteins CHX10/VSX2, VSX1 and PAX4, as well as the zinc finger protein SP3, as factors that can stimulate both the mouse and human Nxnl1 promoter. In addition, CHX10/VSX2 binds to the Nxnl1 promoter in vivo. Since CHX10/VSX2 is expressed predominantly in the inner retina, this finding motivated us to demonstrate that RdCVF is expressed in the inner as well as the outer retina. Interestingly, the loss of rods in the rd1 mouse, a model of retinitis pigmentosa, is associated with decreased expression of RdCVF by inner retinal cells as well as by rods. Based upon these results, we propose an alternative therapeutic strategy aimed at recapitulating RdCVF expression in the inner retina, where cell loss is not significant, to prevent secondary cone death and central vision loss in patients suffering from retinitis pigmentosa.

Mu-opioid receptor activation induces transcriptional plasticity in the central extended amygdala

Addiction develops from the gradual adaptation of the brain to chronic drug exposure, and involves genetic reprogramming of neuronal function. The central extended amygdala (EAc) is a network formed by the central amygdala and the bed nucleus of the stria terminalis. This key site controls drug craving and seeking behaviors, and has not been investigated at the gene regulation level. We used Affymetrix microarrays to analyze transcriptional activity in the murine EAc, with a focus on mu-opioid receptor-associated events because these receptors mediate drug reward and dependence. We identified 132 genes whose expression is regulated by a chronic escalating morphine regimen in the EAc from wild-type but not mu-opioid receptor knockout mice. These modifications are mostly EAc-specific. Gene ontology analysis reveals an overrepresentation of neurogenesis, cell growth and signaling protein categories. A separate quantitative PCR analysis of genes in the last of these groups confirms the dysregulation of both orphan (Gpr88) and known (DrD1A, Adora2A, Cnr1, Grm5, Gpr6) G protein-coupled receptors, scaffolding (PSD95, Homer1) and signaling (Sgk, Cap1) proteins, and neuropeptides (CCK, galanin). These transcriptional modifications do not occur following a single morphine injection, and hence result from long-term adaptation to excessive mu receptor activation. Proteins encoded by these genes are classically associated with spine modules function in other brain areas, and therefore our data suggest a remodeling of EAc circuits at sites where glutamatergic and monoaminergic afferences interact. Together, mu receptor-dependent genes identified in this study potentially contribute to drug-induced neural plasticity, and provide a unique molecular repertoire towards understanding drug craving and relapse.

Transcriptome analysis identifies genes with enriched expression in the mouse central extended amygdala

The central extended amygdala (EAc) is an ensemble of highly interconnected limbic structures of the anterior brain, and forms a cellular continuum including the bed nucleus of the stria terminalis (BNST), the central nucleus of the amygdala (CeA) and the nucleus accumbens shell (AcbSh). This neural network is a key site for interactions between brain reward and stress systems, and has been implicated in several aspects of drug abuse. In order to increase our understanding of EAc function at the molecular level, we undertook a genome-wide screen (Affymetrix) to identify genes whose expression is enriched in the mouse EAc. We focused on the less-well known BNST-CeA areas of the EAc, and identified 121 genes that exhibit more than twofold higher expression level in the EAc compared with whole brain. Among these, 43 genes have never been described to be expressed in the EAc. We mapped these genes throughout the brain, using non-radioactive in situ hybridization, and identified eight genes with a unique and distinct rostro-caudal expression pattern along AcbSh, BNST and CeA. Q-PCR analysis performed in brain and peripheral organ tissues indicated that, with the exception of one (Spata13), all these genes are predominantly expressed in brain. These genes encode signaling proteins (Adora2, GPR88, Arpp21 and Rem2), a transcription factor (Limh6) or proteins of unknown function (Rik130, Spata13 and Wfs1). The identification of genes with enriched expression expands our knowledge of EAc at a molecular level, and provides useful information to toward genetic manipulations within the EAc.

Gene expression is altered in the lateral hypothalamus upon activation of the mu opioid receptor

The lateral hypothalamus (LH) is a brain structure that controls hedonic properties of both natural rewards and drugs of abuse. Mu opioid receptors are known to mediate drug reward, but whether overstimulation of these receptors impacts on LH function has not been studied. Here we have used a genome-wide microarray approach to identify LH responses to chronic mu opioid receptor activation at the transcriptional level. We have subjected wild-type and mu opioid receptor knockout mice to an escalating morphine regimen, which produces severe physical dependence in wild-type but not mutant animals. We have analyzed gene profiles in LH samples using the 430A.2 Affymetrix array and identified a set of 25 genes whose expression is altered by morphine in wild-type mice only. The regulation was confirmed for a subset of these genes using real-time quantitative PCR on samples from independent treatments. Altered expression of aquaporin 4, apolipoprotein D, and prostaglandin synthase is indicative of modified LH physiology. The regulation of two signaling genes (the serum glucocorticoid kinase and the regulator of G protein signaling 4) suggests that neurotransmission is altered in LH circuitry. Finally, the downregulation of apelin may indicate a potential role for this neuropeptide in opioid signaling and hedonic homeostasis. Altogether, our study shows that chronic mu opioid receptor stimulation induces gene expression plasticity in the LH and provides a unique collection of mu opioid receptor-dependent genes that potentially contribute to alter reward processes in addictive diseases.

Retinoic acid induces TGFbeta-dependent autocrine fibroblast growth

To evaluate the role of murine TFIID subunit TAF4 in activation of cellular genes by all-trans retinoic acid (T-RA), we have characterized the T-RA response of taf4(lox/-) and taf4(-/-) embryonic fibroblasts. T-RA regulates almost 1000 genes in taf4(lox/-) cells, but less than 300 in taf4(-/-) cells showing that TAF4 is required for T-RA regulation of most, but not all cellular genes. We further show that T-RA-treated taf4(lox/-) cells exhibit transforming growth factor (TGF)beta-dependent autocrine growth and identify a set of genes regulated by loss of TAF4 and by T-RA corresponding to key mediators of the TGFbeta signalling pathway. T-RA rapidly and potently induces expression of connective tissue growth factor (CTGF) via a conserved DR2 type response element in its proximal promoter leading to serum-free autocrine growth. These results highlight the role of TAF4 as a cofactor in the cellular response to T-RA and identify the genetic programme of a novel cross talk between the T-RA and TGFbeta pathways that leads to deregulated cell growth.

Rod-derived Cone Viability Factor-2 is a novel bifunctional-thioredoxin-like protein with therapeutic potential

Background

Cone degeneration is the hallmark of the inherited retinal disease retinitis pigmentosa. We have previously identified a trophic factor "Rod-derived Cone Viability Factor (RdCVF) that is secreted by rods and promote cone viability in a mouse model of the disease.

Results

Here we report the bioinformatic identification and the experimental analysis of RdCVF2, a second trophic factor belonging to the Rod-derived Cone Viability Factor family. The mouse RdCVF gene is known to be bifunctional, encoding both a long thioredoxin-like isoform (RdCVF-L) and a short isoform with trophic cone photoreceptor viability activity (RdCVF-S). RdCVF2 shares many similarities with RdCVF in terms of gene structure, expression in a rod-dependent manner and protein 3D structure. Furthermore, like RdCVF, the RdCVF2 short isoform exhibits cone rescue activity that is independent of its putative thiol-oxydoreductase activity.

Conclusion

Taken together, these findings define a new family of bifunctional genes which are: expressed in vertebrate retina, encode trophic cone viability factors, and have major therapeutic potential for human retinal neurodegenerative diseases such as retinitis pigmentosa.

PromAn: an integrated knowledge-based web server dedicated to promoter analysis

PromAn is a modular web-based tool dedicated to promoter analysis that integrates distinct complementary databases, methods and programs. PromAn provides automatic analysis of a genomic region with minimal prior knowledge of the genomic sequence. Prediction programs and experimental databases are combined to locate the transcription start site (TSS) and the promoter region within a large genomic input sequence. Transcription factor binding sites (TFBSs) can be predicted using several public databases and user-defined motifs. Also, a phylogenetic footprinting strategy, combining multiple alignment of large genomic sequences and assignment of various scores reflecting the evolutionary selection pressure, allows for evaluation and ranking of TFBS predictions. PromAn results can be displayed in an interactive graphical user interface, PromAnGUI. It integrates all of this information to highlight active promoter regions, to identify among the huge number of TFBS predictions those which are the most likely to be potentially functional and to facilitate user refined analysis. Such an integrative approach is essential in the face of a growing number of tools dedicated to promoter analysis in order to propose hypotheses to direct further experimental validations. PromAn is publicly available at http://bips.u-strasbg.fr/PromAn.

Polyglutamine expansion causes neurodegeneration by altering the neuronal differentiation program

Huntington's disease (HD) and spinocerebellar ataxia type 7 (SCA7) belong to a group of inherited neurodegenerative diseases caused by polyglutamine (polyQ) expansion in corresponding proteins. Transcriptional alteration is a unifying feature of polyQ disorders; however, the relationship between polyQ-induced gene expression deregulation and degenerative processes remains unclear. R6/2 and R7E mouse models of HD and SCA7, respectively, present a comparable retinal degeneration characterized by progressive reduction of electroretinograph activity and important morphological changes of rod photoreceptors. The retina, which is a simple central nervous system tissue, allows correlating functional, morphological and molecular defects. Taking advantage of comparing polyQ-induced degeneration in two retina models, we combined gene expression profiling and molecular biology techniques to decipher the molecular pathways underlying polyQ expansion toxicity. We show that R7E and R6/2 retinal phenotype strongly correlates with loss of expression of a large cohort of genes specifically involved in phototransduction function and morphogenesis of differentiated rod photoreceptors. Accordingly, three key transcription factors (Nrl, Crx and Nr2e3) controlling rod differentiation genes, hence expression of photoreceptor specific traits, are down-regulated. Interestingly, other transcription factors known to cause inhibitory effects on photoreceptor differentiation when mis-expressed, such as Stat3, are aberrantly re-activated. Thus, our results suggest that independently from the protein context, polyQ expansion overrides the control of neuronal differentiation and maintenance, thereby causing dysfunction and degeneration.

GOAnno: GO annotation based on multiple alignment

GOAnno is a web tool that automatically annotates proteins according to the Gene Ontology (GO) using evolutionary information available in hierarchized multiple alignments. GO terms present in the aligned functional subfamily can be cross-validated and propagated to obtain highly reliable predicted GO annotation based on the GOAnno algorithm.

AVAILABILITY

The web tool and a reduced version for local installation are freely available at http://igbmc.u-strasbg.fr/GOAnno/GOAnno.html

SUPPLEMENTARY INFORMATION

The website supplies a detailed explanation and illustration of the algorithm at http://igbmc.u-strasbg.fr/GOAnno/GOAnnoHelp.html.