Intergenic and genic sequence lengths have opposite relationships with respect to gene expression.

Eukaryotic genomes are mostly composed of noncoding DNA whose role is still poorly understood. Studies in several organisms have shown correlations between the length of the intergenic and genic sequences of a gene and the expression of its corresponding mRNA transcript. Some studies have found a positive relationship between intergenic sequence length and expression diversity between tissues, and concluded that genes under greater regulatory control require more regulatory information in their intergenic sequences. Other reports found a negative relationship between expression level and gene length and the interpretation was that there is selection pressure for highly expressed genes to remain small. However, a correlation between gene sequence length and expression diversity, opposite to that observed for intergenic sequences, has also been reported, and to date there is no testable explanation for this observation. To shed light on these varied and sometimes conflicting results, we performed a thorough study of the relationships between sequence length and gene expression using cell-type (tissue) specific microarray data in Arabidopsis thaliana. We measured median gene expression across tissues (expression level), expression variability between tissues (expression pattern uniformity), and expression variability between replicates (expression noise). We found that intergenic (upstream and downstream) and genic (coding and noncoding) sequences have generally opposite relationships with respect to expression, whether it is tissue variability, median, or expression noise. To explain these results we propose a model, in which the lengths of the intergenic and genic sequences have opposite effects on the ability of the transcribed region of the gene to be epigenetically regulated for differential expression. These findings could shed light on the role and influence of noncoding sequences on gene expression.

Use of 16S ribosomal RNA gene analyses to characterize the bacterial signature associated with poor oral health in West Virginia.

BACKGROUND: West Virginia has the worst oral health in the United States, but the reasons for this are unclear. This pilot study explored the etiology of this disparity using culture-independent analyses to identify bacterial species associated with oral disease. METHODS: Bacteria in subgingival plaque samples from twelve participants in two independent West Virginia dental-related studies were characterized using 16S rRNA gene sequencing and Human Oral Microbe Identification Microarray (HOMIM) analysis. Unifrac analysis was used to characterize phylogenetic differences between bacterial communities obtained from plaque of participants with low or high oral disease, which was further evaluated using clustering and Principal Coordinate Analysis. RESULTS: Statistically different bacterial signatures (P<0.001) were identified in subgingival plaque of individuals with low or high oral disease in West Virginia based on 16S rRNA gene sequencing. Low disease contained a high frequency of Veillonella and Streptococcus, with a moderate number of Capnocytophaga. High disease exhibited substantially increased bacterial diversity and included a large proportion of Clostridiales cluster bacteria (Selenomonas, Eubacterium, Dialister). Phylogenetic trees constructed using 16S rRNA gene sequencing revealed that Clostridiales were repeated colonizers in plaque associated with high oral disease, providing evidence that the oral environment is somehow influencing the bacterial signature linked to disease. CONCLUSIONS: Culture-independent analyses identified an atypical bacterial signature associated with high oral disease in West Virginians and provided evidence that the oral environment influenced this signature. Both findings provide insight into the etiology of the oral disparity in West Virginia.

Gene selection using iterative feature elimination random forests for survival outcomes.

Although many feature selection methods for classification have been developed, there is a need to identify genes in high-dimensional data with censored survival outcomes. Traditional methods for gene selection in classification problems have several drawbacks. First, the majority of the gene selection approaches for classification are single-gene based. Second, many of the gene selection procedures are not embedded within the algorithm itself. The technique of random forests has been found to perform well in high-dimensional data settings with survival outcomes. It also has an embedded feature to identify variables of importance. Therefore, it is an ideal candidate for gene selection in high-dimensional data with survival outcomes. In this paper, we develop a novel method based on the random forests to identify a set of prognostic genes. We compare our method with several machine learning methods and various node split criteria using several real data sets. Our method performed well in both simulations and real data analysis.Additionally, we have shown the advantages of our approach over single-gene-based approaches. Our method incorporates multivariate correlations in microarray data for survival outcomes. The described method allows us to better utilize the information available from microarray data with survival outcomes.

Convergent differential regulation of parvalbumin in the brains of vocal learners.

Spoken language and learned song are complex communication behaviors found in only a few species, including humans and three groups of distantly related birds--songbirds, parrots, and hummingbirds. Despite their large phylogenetic distances, these vocal learners show convergent behaviors and associated brain pathways for vocal communication. However, it is not clear whether this behavioral and anatomical convergence is associated with molecular convergence. Here we used oligo microarrays to screen for genes differentially regulated in brain nuclei necessary for producing learned vocalizations relative to adjacent brain areas that control other behaviors in avian vocal learners versus vocal non-learners. A top candidate gene in our screen was a calcium-binding protein, parvalbumin (PV). In situ hybridization verification revealed that PV was expressed significantly higher throughout the song motor pathway, including brainstem vocal motor neurons relative to the surrounding brain regions of all distantly related avian vocal learners. This differential expression was specific to PV and vocal learners, as it was not found in avian vocal non-learners nor for control genes in learners and non-learners. Similar to the vocal learning birds, higher PV up-regulation was found in the brainstem tongue motor neurons used for speech production in humans relative to a non-human primate, macaques. These results suggest repeated convergent evolution of differential PV up-regulation in the brains of vocal learners separated by more than 65-300 million years from a common ancestor and that the specialized behaviors of learned song and speech may require extra calcium buffering and signaling.

Coordinated activation of candidate proto-oncogenes and cancer testes antigens via promoter demethylation in head and neck cancer and lung cancer.

BACKGROUND: Epigenetic alterations have been implicated in the pathogenesis of solid tumors, however, proto-oncogenes activated by promoter demethylation have been sporadically reported. We used an integrative method to analyze expression in primary head and neck squamous cell carcinoma (HNSCC) and pharmacologically demethylated cell lines to identify aberrantly demethylated and expressed candidate proto-oncogenes and cancer testes antigens in HNSCC. METHODOLOGY/PRINCIPAL FINDINGS: We noted coordinated promoter demethylation and simultaneous transcriptional upregulation of proto-oncogene candidates with promoter homology, and phylogenetic footprinting of these promoters demonstrated potential recognition sites for the transcription factor BORIS. Aberrant BORIS expression correlated with upregulation of candidate proto-oncogenes in multiple human malignancies including primary non-small cell lung cancers and HNSCC, induced coordinated proto-oncogene specific promoter demethylation and expression in non-tumorigenic cells, and transformed NIH3T3 cells. CONCLUSIONS/SIGNIFICANCE: Coordinated, epigenetic unmasking of multiple genes with growth promoting activity occurs in aerodigestive cancers, and BORIS is implicated in the coordinated promoter demethylation and reactivation of epigenetically silenced genes in human cancers.

Sensitive and precise quantification of insulin-like mRNA expression in Caenorhabditis elegans.

Insulin-like signaling regulates developmental arrest, stress resistance and lifespan in the nematode Caenorhabditis elegans. However, the genome encodes 40 insulin-like peptides, and the regulation and function of individual peptides is largely uncharacterized. We used the nCounter platform to measure mRNA expression of all 40 insulin-like peptides as well as the insulin-like receptor daf-2, its transcriptional effector daf-16, and the daf-16 target gene sod-3. We validated the platform using 53 RNA samples previously characterized by high density oligonucleotide microarray analysis. For this set of genes and the standard nCounter protocol, sensitivity and precision were comparable between the two platforms. We optimized conditions of the nCounter assay by varying the mass of total RNA used for hybridization, thereby increasing sensitivity up to 50-fold and reducing the median coefficient of variation as much as 4-fold. We used deletion mutants to demonstrate specificity of the assay, and we used optimized conditions to assay insulin-like gene expression throughout the C. elegans life cycle. We detected expression for nearly all insulin-like genes and find that they are expressed in a variety of distinct patterns suggesting complexity of regulation and specificity of function. We identified insulin-like genes that are specifically expressed during developmental arrest, larval development, adulthood and embryogenesis. These results demonstrate that the nCounter platform provides a powerful approach to analyzing insulin-like gene expression dynamics, and they suggest hypotheses about the function of individual insulin-like genes.

Evaluation of an epithelial plasticity biomarker panel in men with localized prostate cancer.

BACKGROUND: Given the potential importance of epithelial plasticity (EP) to cancer metastasis, we sought to investigate biomarkers related to EP in men with localized prostate cancer (PC) for the association with time to PSA recurrence and other clinical outcomes after surgery. METHODS: Men with localized PC treated with radical prostatectomy at the Durham VA Medical Center and whose prostatectomy tissues were included in a tissue microarray (TMA) linked to long-term outcomes. We performed immunohistochemical studies using validated antibodies against E-cadherin and Ki-67 and mesenchymal biomarkers including N-cadherin, vimentin, SNAIL, ZEB1 and TWIST. Association studies were conducted for each biomarker with baseline clinical/pathologic characteristics an risk of PSA recurrence over time. RESULTS: Two hundred and five men contributed TMA tissue and had long-term follow-up (median 11 years). Forty-three percent had PSA recurrence; three died of PC. The majority had high E-cadherin expression (86%); 14% had low/absent E-cadherin expression. N-cadherin was rarely expressed (<4%) and we were unable to identify an E-to-N-cadherin switch as independently prognostic. No associations with clinical risk group, PSA recurrence or Gleason sum were noted for SNAIL, ZEB1, vimentin or TWIST, despite heterogeneous expression between patients. We observed an association of higher Ki-67 expression with Gleason sum (P=0.043), National Comprehensive Cancer Network risk (P=0.013) and PSA recurrence (hazard ratio 1.07, P=0.016). CONCLUSIONS: The expression of EP biomarkers in this cohort of men with a low risk of PC-specific mortality was not associated with aggressive features or PSA relapse after surgery.

Bordetella pertussis from functional genomics to intranasal vaccination.

Whooping cough still represents a major health problem, despite the use of effective vaccines for several decades. Being classically a typical childhood disease, whooping cough in young adults is now more common than it used to be, suggesting that protection after vaccination wanes during adolescence. As an alternative to the current vaccines, we wish to develop live attenuated vaccines to be delivered by the nasal route, such as to mimic the natural route of infection and to induce long lasting immunity. Bordetella pertussis, the etiological agent of whooping cough, produces a number of virulence factors, including toxins. Its recently determined genome sequence makes it now possible to apply functional genomics, such as transcriptomics and systematic knock-out mutagenesis. The expression of most known B. pertussis virulence genes is controlled by the two-component system BvgA/S. DNA microarray analyses have led to the identification of novel genes in the BvgA/S regulon, some of which are activated by BvgA/S and others are repressed by BvgA/S. In addition, some genes appear to be differentially modulated by nicotinic acid and MgSO4, both known to modulate the expression of BvgA/S-regulated genes. Among others, the functional genomics approach has uncovered two strongly BvgA/S-activated genes, named hotA and hotB (for 'homolog of toxin'), the products of which show high sequence similarities to pertussis toxin subunits. The identification of the full array of virulence factors, as well as an integrated understanding of the bacterial physiology should allow us to design attenuated B. pertussis strains useful for intranasal vaccination. A first generation of attenuated strains has already shown full protection in mice after a single intranasal administration. Such strains may also serve as vaccine carriers for heterologous antigens, in order to vaccinate against several different pathogens simultaneously.

Selection of yeast strains for bioethanol production from UK seaweeds

Macroalgae (seaweeds) are a promising feedstock for the production of third generation bioethanol, since they have high carbohydrate contents, contain little or no lignin and are available in abundance. However, seaweeds typically contain a more diverse array of monomeric sugars than are commonly present in feedstocks derived from lignocellulosic material which are currently used for bioethanol production. Hence, identification of a suitable fermentative microorganism that can utilise the principal sugars released from the hydrolysis of macroalgae remains a major objective. The present study used a phenotypic microarray technique to screen 24 different yeast strains for their ability to metabolise individual monosaccharides commonly found in seaweeds, as well as hydrolysates following an acid pre-treatment of five native UK seaweed species (Laminaria digitata, Fucus serratus, Chondrus crispus, Palmaria palmata and Ulva lactuca). Five strains of yeast (three Saccharomyces spp, one Pichia sp and one Candida sp) were selected and subsequently evaluated for bioethanol production during fermentation of the hydrolysates. Four out of the five selected strains converted these monomeric sugars into bioethanol, with the highest ethanol yield (13 g L−1) resulting from a fermentation using C. crispus hydrolysate with Saccharomyces cerevisiae YPS128. This study demonstrated the novel application of a phenotypic microarray technique to screen for yeast capable of metabolising sugars present in seaweed hydrolysates; however, metabolic activity did not always imply fermentative production of ethanol.

Detection of Emerging and Re-Emerging Pathogens in Surface Waters Close to an Urban Area

Current knowledge about the spread of pathogens in aquatic environments is scarce probably because bacteria, viruses, algae and their toxins tend to occur at low concentrations in water, making them very difficult to measure directly. The purpose of this study was the development and validation of tools to detect pathogens in freshwater systems close to an urban area. In order to evaluate anthropogenic impacts on water microbiological quality, a phylogenetic microarray was developed in the context of the EU project µAQUA to detect simultaneously numerous pathogens and applied to samples from two different locations close to an urban area located upstream and downstream of Rome in the Tiber River. Furthermore, human enteric viruses were also detected. Fifty liters of water were collected and concentrated using a hollow-fiber ultrafiltration approach. The resultant concentrate was further size-fractionated through a series of decreasing pore size filters. RNA was extracted from pooled filters and hybridized to the newly designed microarray to detect pathogenic bacteria, protozoa and toxic cyanobacteria. Diatoms as indicators of the water quality status, were also included in the microarray to evaluate water quality. The microarray results gave positive signals for bacteria, diatoms, cyanobacteria and protozoa. Cross validation of the microarray was performed using standard microbiological methods for the bacteria. The presence of oral-fecal transmitted human enteric-viruses were detected using q-PCR. Significant concentrations of Salmonella, Clostridium, Campylobacter and Staphylococcus as well as Hepatitis E Virus (HEV), noroviruses GI (NoGGI) and GII (NoGII) and human adenovirus 41 (ADV 41) were found in the Mezzocammino site, whereas lower concentrations of other bacteria and only the ADV41 virus was recovered at the Castel Giubileo site. This study revealed that the pollution level in the Tiber River was considerably higher downstream rather than upstream of Rome and the downstream location was contaminated by emerging and re-emerging pathogens.

Elucidation of the interactions between the transcription factor E2A and the transcriptional co-activator CBP/p300

E2A is a transcription factor that plays a particularly critical role in lymphopoiesis. The chromosomal translocation 1;19, disrupts the E2A gene and results in the expression of the fusion oncoprotein E2A-PBX1, which is implicated in acute lymphoblastic leukemia. Both E2A and E2A-PBX1 contain two activation domains, AD1 and AD2, which comprise conserved ΦxxΦΦ motifs where Φ denotes a hydrophobic amino acid. These domains function to recruit transcriptional co-activators and repressors, including the histone acetyl transferase CREB binding protein (CBP) and its paralog p300. The PCET motif within E2A AD1 interacts with the KIX domain of CBP/p300, the disruption of which abrogates the transcriptional activation by E2A and the transformative properties of E2A-PBX1. The generation of a peptide-based inhibitor targeting the PCET:KIX interaction would serve useful in further assessing the role of E2A and E2A-PBX1 in lymphopoiesis and leukemogenesis. An interaction between E2A AD2 and the KIX domain has also been recently identified, and the TAZ domains of CBP/p300 have been shown to interact with several transcription factors that contain ΦxxΦΦ motifs. Thus the design of an inhibitor of the E2A:CBP/p300 interaction requires the full complement of interactions between E2A and the various domains of CBP/p300 to be elucidated. Here, we have used nuclear magnetic resonance (NMR) spectroscopy to determine that AD2 interacts with KIX at the same site as PCET, which indicates that the E2A:KIX interaction can be disrupted by targeting a single binding site. Using an iterative synthetic peptide microarray approach, a peptide with the sequence DKELQDLLDFSLQY was derived from PCET to interact with KIX with higher affinity than the wild type sequence. This peptide now serves as a lead molecule for further development as an inhibitor of the E2A:CBP/p300 interaction. Fluorescence anisotropy, peptide microarray technology, and isothermal titration calorimetry were employed to characterize interactions between both TAZ domains of CBP/p300 and the PCET motif and AD2 of E2A. Alanine substitution of residues within PCET demonstrated that the ΦxxΦΦ motif is a key mediator of these interactions, analogous to the PCET:KIX interaction. These findings now inform future work to establish possible physiological roles for the E2A:TAZ1 and E2A:TAZ2 interactions.

Cybr, a cytokine-inducible protein that binds cytohesin-1 and regulates its activity

Cytokines regulate lymphocyte development and differentiation, but precisely how they control these processes is still poorly understood. By using microarray technology to detect cytokine-induced genes, we identified a cDNA encoding Cybr, which was increased markedly in cells incubated with IL-2 and IL-12. The mRNA was most abundant in hematopoietic cells and tissues. The predicted amino acid sequence is similar to that of GRP-1-associated protein (GRASP), a recently identified retinoic acid-induced cytohesin-binding protein. Physical interaction, dependent on the coiled-coil domains of Cybr and cytohesin-1, was demonstrated by coimmunoprecipitation of the overexpressed proteins from 293T cells. Cytohesin-1, in addition to its role in cell adhesion, is a guanine nucleotide-exchange protein activator of ARF GTPases. Acceleration of guanosine 5'-O-(thiotriphosphate) binding to ARF by cytohesin-1 in vitro was enhanced by Cybr. Because the binding protein modified activation of ADP ribosylation factor by cytohesin-1, we designate this cytokine-inducible protein Cybr (cytohesin binder and regulator).

5-Fluorouracil: mechanisms of action and clinical strategies

5-fluorouracil (5-FU) is widely used in the treatment of cancer. Over the past 20 years, increased understanding of the mechanism of action of 5-FU has led to the development of strategies that increase its anticancer activity. Despite these advances, drug resistance remains a significant limitation to the clinical use of 5-FU. Emerging technologies, such as DNA microarray profiling, have the potential to identify novel genes that are involved in mediating resistance to 5-FU. Such target genes might prove to be therapeutically valuable as new targets for chemotherapy, or as predictive biomarkers of response to 5-FU-based chemotherapy.

BRCA1 Regulates IFN-γ Signaling through a Mechanism Involving the Type I IFNs

BRCA1 encodes a tumor suppressor gene that is mutated in the germ line of women with a genetic predisposition to breast and ovarian cancer. BRCA1 has been implicated in a number of important cellular functions including DNA damage repair, transcriptional regulation, cell cycle control, and ubiquitination. Using an Affymetrix U95A microarray, IRF-7 was identified as a BRCA1 transcriptional target and was also shown to be synergistically up-regulated by BRCA1 specifically in the presence of IFN-gamma, coincident with the synergistic induction of apoptosis. We show that BRCA1, signal transducer and activator of transcription (STAT)-1, and STAT2 are all required for the induction of IRF-7 following stimulation with IFN-gamma. We also show that the induction of IRF-7 by BRCA1 and IFN-gamma is dependent on the type I IFNs, IFN-alpha and IFN-beta. We show that BRCA1 is required for the up-regulation of STAT1, STAT2, and the type I IFNs in response to IFN-gamma. We show that BRCA1 is localized at the promoters of the molecules involved in type I IFN signaling leading to their up-regulation. Blocking this intermediary type I IFN step using specific antisera shows the requirement for IFN-alpha and IFN-beta in the induction of IRF-7 and apoptosis. Finally, we outline a mechanism for the BRCA1/IFN-gamma regulation of target genes involved in the innate immune response, which is dependent on type I IFN signaling.