Background

A central question for ecologists is the extent to which anthropogenic disturbances (e.g. tourism) might impact wildlife and affect the systems under study. From a research perspective, identifying the effects of human disturbance caused by research-related activities is crucial in order to understand and account for potential biases and derive appropriate conclusions from the data.

Results

Here, we document a case of biological adjustment to chronic human disturbance in a colonial seabird, the king penguin (Aptenodytes patagonicus), breeding on remote and protected islands of the Southern ocean. Using heart rate (HR) as a measure of the stress response, we show that, in a colony with areas exposed to the continuous presence of humans (including scientists) for over 50 years, penguins have adjusted to human disturbance and habituated to certain, but not all, types of stressors. When compared to birds breeding in relatively undisturbed areas, birds in areas of high chronic human disturbance were found to exhibit attenuated HR responses to acute anthropogenic stressors of low-intensity (i.e. sounds or human approaches) to which they had been subjected intensely over the years. However, such attenuation was not apparent for high-intensity stressors (i.e. captures for scientific research) which only a few individuals experience each year.

Conclusions

Habituation to anthropogenic sounds/approaches could be an adaptation to deal with chronic innocuous stressors, and beneficial from a research perspective. Alternately, whether penguins have actually habituated to anthropogenic disturbances over time or whether human presence has driven the directional selection of human-tolerant phenotypes, remains an open question with profound ecological and conservation implications, and emphasizes the need for more knowledge on the effects of human disturbance on long-term studied populations.

There is a need for improved therapies for severe asthma. Lebrikizumab, a humanized monoclonal antibody that binds to interleukin (IL)-13, is under development for the treatment of poorly controlled asthma. This article reviews the potential role of IL-13 in the pathogenesis of asthma, the efficacy and safety of lebrikizumab in humans, and progress in patient selection for lebrikizumab therapy. IL-13 is a T-helper (Th2) cell-derived cytokine implicated in inflammatory responses in asthma, including serum immunoglobulin-E synthesis, mucus hypersecretion, and subepithelial fibrosis. Blocking the pro-inflammatory effects of IL-13 with lebrikizumab has the potential to improve asthma control. Published data on the efficacy and safety of lebrikizumab in the treatment of asthma are relatively limited. The late asthmatic response after inhaled allergen challenge is reduced by almost 50%, following treatment with lebrikizumab. In a Phase II study performed in 219 adults with poorly controlled asthma despite inhaled corticosteroids (MILLY trial), lebrikizumab produced an improvement in prebronchodilator forced expiratory volume in 1 second of 5.5% compared with placebo at 12 weeks, but had no effects on other efficacy end points. Adverse effects were similar to placebo, except that musculoskeletal side effects occurred slightly more often with lebrikizumab. Stratifying patients into a high Th2 phenotype using serum periostin, which is upregulated in lung epithelial cells by IL-13, may identify individuals responsive to blockade of IL-13. In the MILLY trial, lebrikizumab treatment was associated with greater improvement in lung function in patients with elevated serum periostin levels compared with those with low periostin levels. Two large Phase III randomized controlled trials in patients with uncontrolled asthma are underway to establish the safety and efficacy of lebrikizumab when administered over a 52-week period. These studies will also help to determine whether identifying patients with a Th2 high inflammatory phenotype using serum periostin allows a personalized approach to the treatment of asthma.

SUMMARY

DNA methylation has been implicated as an epigenetic component of mechanisms that stabilize cell-fate decisions. Here, we have characterized the methylomes of human female hematopoietic stem/progenitor cells (HSPCs) and mature cells from the myeloid and lymphoid lineages. Hypomethylated regions (HMRs) associated with lineage-specific genes were often methylated in the opposing lineage. In HSPCs, these sites tended to show intermediate, complex patterns that resolve to uniformity upon differentiation, by increased or decreased methylation. Promoter HMRs shared across diverse cell types typically display a constitutive core that expands and contracts in a lineage-specific manner to fine-tune the expression of associated genes. Many newly identified intergenic HMRs, both constitutive and lineage specific, were enriched for factor binding sites with an implied role in genome organization and regulation of gene expression, respectively. Overall, our studies represent an important reference data set and provide insights into directional changes in DNA methylation as cells adopt terminal fates.

Translation regulation plays important roles in both normal physiological conditions and diseases states. This regulation requires cis-regulatory elements located mostly in 5′ and 3′ UTRs and trans-regulatory factors (e.g., RNA binding proteins (RBPs)) which recognize specific RNA features and interact with the translation machinery to modulate its activity. In this paper, we discuss important aspects of 5′ UTR-mediated regulation by providing an overview of the characteristics and the function of the main elements present in this region, like uORF (upstream open reading frame), secondary structures, and RBPs binding motifs and different mechanisms of translation regulation and the impact they have on gene expression and human health when deregulated.

Summary

During germ cell and preimplantation development, mammalian cells undergo nearly complete reprogramming of DNA methylation patterns. We profiled the methylomes of human and chimp sperm as a basis for comparison to methylation patterns of ES cells. While the majority of promoters escape methylation in both ES cells and sperm, the corresponding hypomethylated regions show substantial structural differences. Repeat elements are heavily methylated in both germ and somatic cells; however, retrotransposons from several subfamilies evade methylation more effectively during male germ cell development, while other subfamilies show the opposite trend. Comparing methylomes of human and chimp sperm revealed a subset of differentially methylated promoters and strikingly divergent methylation in retrotransposon subfamilies, with an evolutionary impact that is apparent in the underlying genomic sequence. Thus, the features that determine DNA methylation patterns differ between male germ cells and somatic cells, and elements of these features have diverged between humans and chimpanzees.

Motivation: Post-translational modifications to histones have several well known associations with regulation of gene expression. While some modifications appear concentrated narrowly, covering promoters or enhancers, others are dispersed as epigenomic domains. These domains mark contiguous regions sharing an epigenomic property, such as actively transcribed or poised genes, or heterochromatically silenced regions. While high-throughput methods like ChIP-Seq have led to a flood of high-quality data about these epigenomic domains, there remain important analysis problems that are not adequately solved by current analysis tools.

Results: We present the RSEG method for identifying epigenomic domains from ChIP-Seq data for histone modifications. In contrast with other methods emphasizing the locations of ‘peaks’ in read density profiles, our method identifies the boundaries of domains. RSEG is also able to incorporate a control sample and find genomic regions with differential histone modifications between two samples.

Availability: RSEG, including source code and documentation, is freely available at http://smithlab.cmb.usc.edu/histone/rseg/.

Contact: anrewds@usc.edu

Supplementary information: Supplementary data are available at Bioinformatics online.

The lifespan of a mammalian mRNA is determined, in part, by the binding of regulatory proteins and small RNA-guided complexes. The conserved endonuclease activity of Argonaute2 requires extensive complementarity between a small RNA and its target and is not used by animal microRNAs, which pair with their targets imperfectly. Here, we investigate the endonucleolytic function of Ago2 and other nucleases by transcriptome-wide profiling of mRNA cleavage products retaining 5′-phosphate groups in mouse ES. We detect a prominent signature of Ago2-dependent cleavage events and validate several such targets. Unexpectedly, a broader class of Ago2-independent cleavage sites is also observed, indicating participation of additional nucleases in site-specific mRNA cleavage. Within this class, we identify a cohort of Drosha-dependent mRNA cleavage events that functionally regulate mRNA levels in mES cells, including one in the Dgcr8 mRNA. Together, these results highlight the underappreciated role of endonucleolytic cleavage in controlling mRNA fates in mammals.

Summary: We report on a major new version of the RMAP software for mapping reads from short-read sequencing technology. General improvements to accuracy and space requirements are included, along with novel functionality. Included in the RMAP software package are tools for mapping paired-end reads, mapping using more sophisticated use of quality scores, collecting ambiguous mapping locations and mapping bisulfite-treated reads.

Availability: The applications described in this note are available for download at http://www.cmb.usc.edu/people/andrewds/rmap and are distributed as Open Source software under the GPLv3.0. The software has been tested on Linux and OS X platforms.

Contact: andrewds@usc.edu; mzhang@cshl.edu

The RMAP algorithm was introduced by (Smith et al., 2008) as one of the earliest available programs for mapping reads from the Illumina second-generation sequencing technology. One important contribution of RMAP was to incorporate the use of quality scores directly into the mapping process: read positions with too low a quality score were not considered while mapping, and that quality score cutoff could be adjusted by the user. Subsequently, numerous mapping algorithm have appeared (Langmead et al., 2009; Li,H. et al., 2008; Li,R. et al., 2008; Lin et al., 2008; Schatz, 2009; Yanovsky et al., 2008), with improvements in both efficiency and breadth of functionality (e.g. ability to map paired-end reads; integrated SNP calling). Investigators requiring solutions to mapping problems now have many options. As new applications of short-read sequencing emerge, many variations on the analysis task of read mapping emerge. Diversity in performance characteristics of existing mapping tools becomes potentially valuable.

We report the first major update to RMAP. The basic algorithmic framework in RMAP is still to preprocess reads and scan the genome, but several modifications have been made and much additional functionality has been included. Importantly, RMAP has a memory footprint that depends on the number of reads being mapped. This feature allows RMAP to be used effectively in cluster environments with commodity nodes, because partitioning the reads allows natural parallelizations with linear reduction in memory requirements per processor core used.

Included in this release of the RMAP software package is functionality for mapping paired-end reads, making more sophisticated use of quality scores, collecting mapping locations for ambiguously mapping reads and mapping bisulfite-treated reads.

The transcription factor Hypoxia-inducible factor 1 (HIF-1) plays a central role in the transcriptional response to oxygen flux. To gain insight into the molecular pathways regulated by HIF-1, it is essential to identify the downstream-target genes. We report here a strategy to identify HIF-1-target genes based on an integrative genomic approach combining computational strategies and experimental validation. To identify HIF-1-target genes microarrays data sets were used to rank genes based on their differential response to hypoxia. The proximal promoters of these genes were then analyzed for the presence of conserved HIF-1-binding sites. Genes were scored and ranked based on their response to hypoxia and their HIF-binding site score. Using this strategy we recovered 41% of the previously confirmed HIF-1-target genes that responded to hypoxia in the microarrays and provide a catalogue of predicted HIF-1 targets. We present experimental validation for ANKRD37 as a novel HIF-1-target gene. Together these analyses demonstrate the potential to recover novel HIF-1-target genes and the discovery of mammalian-regulatory elements operative in the context of microarray data sets.

Background

Although microarray-based studies have revealed global view of gene expression in cancer cells, we still have little knowledge about regulatory mechanisms underlying the transcriptome. Several computational methods applied to yeast data have recently succeeded in identifying expression modules, which is defined as co-expressed gene sets under common regulatory mechanisms. However, such module discovery methods are not applied cancer transcriptome data.

Results

In order to decode oncogenic regulatory programs in cancer cells, we developed a novel module discovery method termed EEM by extending a previously reported module discovery method, and applied it to breast cancer expression data. Starting from seed gene sets prepared based on cis-regulatory elements, ChIP-chip data, and gene locus information, EEM identified 10 principal expression modules in breast cancer based on their expression coherence. Moreover, EEM depicted their activity profiles, which predict regulatory programs in each subtypes of breast tumors. For example, our analysis revealed that the expression module regulated by the Polycomb repressive complex 2 (PRC2) is downregulated in triple negative breast cancers, suggesting similarity of transcriptional programs between stem cells and aggressive breast cancer cells. We also found that the activity of the PRC2 expression module is negatively correlated to the expression of EZH2, a component of PRC2 which belongs to the E2F expression module. E2F-driven EZH2 overexpression may be responsible for the repression of the PRC2 expression modules in triple negative tumors. Furthermore, our network analysis predicts regulatory circuits in breast cancer cells.

Conclusion

These results demonstrate that the gene set-based module discovery approach is a powerful tool to decode regulatory programs in cancer cells.

Insulator elements affect gene expression by preventing the spread of heterochromatin and restricting transcriptional enhancers from activation of unrelated promoters. In vertebrates, insulator’s function requires association with the CCCTC-binding factor (CTCF), a protein that recognizes long and diverse nucleotide sequences. While insulators are critical in gene regulation, only a few have been reported. Here, we describe 13,804 CTCF binding sites in potential insulators of the human genome, discovered experimentally in primary human fibroblasts. Most of these sequences are located far from the transcriptional start sites, with their distribution strongly correlated with genes. The majority of them fit to a consensus motif highly conserved and suitable for predicting possible insulators driven by CTCF in other vertebrate genomes. In addition, CTCF localization is largely invariant across different cell types. Our results provide resource for investigating insulator function and possible other general and evolutionarily conserved activities of CTCF sites.

Background

Microarray technology has unveiled transcriptomic differences among tumors of various phenotypes, and, especially, brought great progress in molecular understanding of phenotypic diversity of breast tumors. However, compared with the massive knowledge about the transcriptome, we have surprisingly little knowledge about regulatory mechanisms underling transcriptomic diversity.

Results

To gain insights into the transcriptional programs that drive tumor progression, we integrated regulatory sequence data and expression profiles of breast cancer into a Bayesian Network, and searched for cis-regulatory motifs statistically associated with given histological grades and prognosis. Our analysis found that motifs bound by ELK1, E2F, NRF1 and NFY are potential regulatory motifs that positively correlate with malignant progression of breast cancer.

Conclusion

The results suggest that these 4 motifs are principal regulatory motifs driving malignant progression of breast cancer. Our method offers a more concise description about transcriptome diversity among breast tumors with different clinical phenotypes.

Background

Second-generation sequencing has the potential to revolutionize genomics and impact all areas of biomedical science. New technologies will make re-sequencing widely available for such applications as identifying genome variations or interrogating the oligonucleotide content of a large sample (e.g. ChIP-sequencing). The increase in speed, sensitivity and availability of sequencing technology brings demand for advances in computational technology to perform associated analysis tasks. The Solexa/Illumina 1G sequencer can produce tens of millions of reads, ranging in length from ~25–50 nt, in a single experiment. Accurately mapping the reads back to a reference genome is a critical task in almost all applications. Two sources of information that are often ignored when mapping reads from the Solexa technology are the 3' ends of longer reads, which contain a much higher frequency of sequencing errors, and the base-call quality scores.

Results

To investigate whether these sources of information can be used to improve accuracy when mapping reads, we developed the RMAP tool, which can map reads having a wide range of lengths and allows base-call quality scores to determine which positions in each read are more important when mapping. We applied RMAP to analyze data re-sequenced from two human BAC regions for varying read lengths, and varying criteria for use of quality scores. RMAP is freely available for downloading at .

Conclusion

Our results indicate that significant gains in Solexa read mapping performance can be achieved by considering the information in 3' ends of longer reads, and appropriately using the base-call quality scores. The RMAP tool we have developed will enable researchers to effectively exploit this information in targeted re-sequencing projects.

Background

It is becoming increasingly important for researchers to be able to scan through large genomic regions for transcription factor binding sites or clusters of binding sites forming cis-regulatory modules. Correspondingly, there has been a push to develop algorithms for the rapid detection and assessment of cis-regulatory modules. While various algorithms for this purpose have been introduced, most are not well suited for rapid, genome scale scanning.

Results

We introduce methods designed for the detection and statistical evaluation of cis-regulatory modules, modeled as either clusters of individual binding sites or as combinations of sites with constrained organization. In order to determine the statistical significance of module sites, we first need a method to determine the statistical significance of single transcription factor binding site matches. We introduce a straightforward method of estimating the statistical significance of single site matches using a database of known promoters to produce data structures that can be used to estimate p-values for binding site matches. We next introduce a technique to calculate the statistical significance of the arrangement of binding sites within a module using a max-gap model. If the module scanned for has defined organizational parameters, the probability of the module is corrected to account for organizational constraints. The statistical significance of single site matches and the architecture of sites within the module can be combined to provide an overall estimation of statistical significance of cis-regulatory module sites.

Conclusion

The methods introduced in this paper allow for the detection and statistical evaluation of single transcription factor binding sites and cis-regulatory modules. The features described are implemented in the Search Tool for Occurrences of Regulatory Motifs (STORM) and MODSTORM software.

Transcription factor-binding sites and the cis-regulatory modules they compose are central determinants of gene expression. We previously showed that binding site motifs and modules in proximal promoters can be used to predict a significant portion of mammalian tissue-specific transcription. Here, we report on a systematic analysis of promoters controlling tissue-specific expression in heart, kidney, liver, pancreas, skeletal muscle, testis and CD4 T cells, for both human and mouse. We integrated multiple sources of expression data to compile sets of transcripts with strong evidence for tissue-specific regulation. The analysis of the promoters corresponding to these sets produced a catalog of predicted tissue-specific motifs and modules, and cis-regulatory elements. Predicted regulatory interactions are supported by statistical evidence, and provide a foundation for targeted experiments that will improve our understanding of tissue-specific regulatory networks. In a broader context, methods used to construct the catalog provide a model for the analysis of genomic regions that regulate differentially expressed genes.

Background

General protein evolution models help determine the baseline expectations for the evolution of sequences, and they have been extensively useful in sequence analysis and for the computer simulation of artificial sequence data sets.

Results

We have developed a new method of simulating protein sequence evolution, including insertion and deletion (indel) events in addition to amino-acid substitutions. The simulation generates both the simulated sequence family and a true sequence alignment that captures the evolutionary relationships between amino acids from different sequences. Our statistical model for indel evolution is based on the empirical indel distribution determined by Qian and Goldstein. We have parameterized this distribution so that it applies to sequences diverged by varying evolutionary times and generalized it to provide flexibility in simulation conditions. Our method uses a Monte-Carlo simulation strategy, and has been implemented in a C++ program named Simprot.

Conclusion

Simprot will be useful for testing methods of analysis of protein sequence families particularly alignment methods, phylogenetic tree building, detection of recombination and horizontal gene transfer, and homology detection, where knowing the true course of sequence evolution is essential.

Recombinant adeno-associated viruses (rAAVs) have attracted considerable interest as gene delivery systems because they show long-term expression in vivo and transduce numerous cell types. Limitations to successful gene transduction from rAAVs have prompted investigations of a variety of treatments to enhance transgene expression from rAAV vectors. Tyrphostin-1, an epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor, dramatically enhances rAAV transgene expression. Elegant studies have demonstrated that a single-strand D-sequence-binding protein (ssDBP) is phosphorylated by EGFR and binds to the D sequence element in the AAV terminal repeat (TR). Binding of the Tyr-phosphorylated ssDBP prevents conversion of single-stranded vector DNA to a double-strand conformation. We observed dramatic increases in transgene expression in lung epithelial cells (IB3) with tyrphostin treatment. Gel shift analysis of ssDBP revealed that its DNA binding characteristics were unchanged after tyrphostin treatment or adenovirus infection. Tyrphostin stimulated rAAV transgene expression to a greater extent than adenovirus coinfection. Southern hybridizations revealed that the vector DNA remained in the single-strand conformation in tyrphostin-treated cells but double-stranded replicative form monomer DNA was most abundant in adenovirus-infected cells. Northern analyses revealed that tyrphostin treatment enhanced mRNA accumulation more than in adenovirus-infected cultures even though replicative form DNA was undetectable. Analysis of the JNK, ERK, and p38K mitogen-activated protein kinase pathways revealed that tyrphostin treatment stimulated the activity of JNK and p38K. Our data suggest that tyrphostin-induced alteration of stress response pathways results in dramatic enhancement of transcription on linear vector DNA templates in the IB3 cell line. These results expand the downstream targets of the EGFR in regulating rAAV transduction.

Background

Exhaled nitric oxide (FENO) measurements are used as a surrogate marker for eosinophilic airway inflammation. However, many constitutional and environmental factors affect FENO, making it difficult to devise reference values. Our aim was to evaluate the relative importance of factors affecting FENO in a well characterised adult population.

Methods

Data were obtained from 895 members of the Dunedin Multidisciplinary Health and Development Study at age 32. The effects of sex, height, weight, lung function indices, smoking, atopy, asthma and rhinitis on FENO were explored by unadjusted and adjusted linear regression analyses.

Results

The effect of sex on FENO was both statistically and clinically significant, with FENO levels approximately 25% less in females. Overall, current smoking reduced FENO up to 50%, but this effect occurred predominantly in those who smoked on the day of the FENO measurement. Atopy increased FENO by 60%. The sex-related differences in FENO remained significant (p < 0.001) after controlling for all other significant factors affecting FENO.

Conclusion

Even after adjustment, FENO values are significantly different in males and females. The derivation of reference values and the interpretation of FENO in the clinical setting should be stratified by sex. Other common factors such as current smoking and atopy also require to be taken into account.