Chemodiversity in flowers of Tanacetum vulgare has consequences on a florivorous beetle

The chemical composition of plant individuals can vary, leading to high intraspecific chemodiversity. Diversity of floral chemistry may impact the responses of flower-feeding insects. Tanacetum vulgare plants vary significantly in their leaf terpenoid composition, forming distinct chemotypes. We investigated the composition of terpenoids and nutrients of flower heads and pollen in plants belonging to three chemotypes - dominated either by β-thujone (BThu), artemisia ketone (Keto) or a mixture of (Z)-myroxide, santolina triene, and artemisyl acetate (Myrox) - using different analytical platforms. We tested the effects of these differences on preferences, weight gain and performance of adults of the shining flower beetle, Olibrus aeneus. The terpenoid composition and diversity of flower heads and pollen significantly differed among individuals belonging to the above chemotypes, while total concentrations of pollen terpenoids, sugars, amino acids, and lipids did not differ. Beetles preferred BThu over the Myrox chemotype in both olfactory and contact choice assays, while the Keto chemotype was marginally repellent according to olfactory assays. The beetles gained the least weight within 48 h and their initial mortality was highest when feeding exclusively on floral tissues of the Myrox chemotype. Short-term weight gain and long-term performance were highest when feeding on the BThu chemotype. In conclusion, the beetles showed chemotype-specific responses towards different T. vulgare chemotypes, which may be attributed to the terpenoid composition in flower heads and pollen rather than to differences in nutrient profiles. Both richness and overall diversity are important factors when determining chemodiversity of individual plants and their consequences on interacting insects.

Cohesin Subunit RAD21 Regulates the Differentiation and Self-Renewal of Hematopoietic Stem and Progenitor Cells

Recent studies suggest that chromosomal cohesin complex proteins are important in regulating hematopoiesis and may contribute to myeloid malignancies. To investigate the effects of perturbing the cohesin subunit protein RAD21 on normal hematopoiesis, we used conditional knockout (cKO) mouse models. While cohesin is vital for hematopoietic stem cell (HSC) function, Rad21 haploinsufficiency (Rad21Δ/+) led to distinct hematopoietic phenotypes. Our findings revealed that Rad21Δ/+ cells exhibited decreased hematopoietic reconstitution in competitive bone marrow transplantation assays. This reduction in peripheral blood chimerism was specifically observed in the lymphoid compartment, while the chimerism in the myeloid compartment remained unaffected. Rad21 haploinsufficiency also resulted in changes in the hematopoietic stem and progenitor cells (HSPC) and myeloid progenitor compartments, with a significant accumulation of granulocyte-macrophage progenitors in the bone marrow. We observed differential gene expression in Rad21Δ/+ LSK (Lin- Sca1-Kit+) cells, including genes required for HSPC function and differentiation, such as Setdb1, Hmga2, Ncor1, and Myb. In addition, we observed a notable decrease in the expression of genes related to the interferon response and a significant reduction in the expression of genes involved in the IL2-STAT5 signaling pathways. Our studies suggest that RAD21 protein and level of its post-translational modifications in the bone marrow cells may play a potential role in hematopoiesis. Overall, Rad21 haploinsufficiency impairs hematopoietic differentiation and increases HSC self-renewal.

Combinatorial molecule screening identified a novel diterpene and the BET inhibitor CPI-203 as differentiation inducers of primary acute myeloid leukemia cells

Combination treatment has proven effective for patients with acute promyelocytic leukemia, exemplifying the importance of therapy targeting multiple components of oncogenic regulation for a successful outcome. However, recent studies have shown that the mutational complexity of acute myeloid leukemia (AML) precludes the translation of molecular targeting into clinical success. Here, as a complement to genetic profiling, we used unbiased, combinatorial in vitro drug screening to identify pathways that drive AML and to develop personalized combinatorial treatments. First, we screened 513 natural compounds on primary AML cells and identified a novel diterpene (H4) that preferentially induced differentiation of FLT3 wild-type AML, while FLT3-ITD/mutations conferred resistance. The samples responding to H4, displayed increased expression of myeloid markers, a clear decrease in the nuclear-cytoplasmic ratio and the potential of re-activation of the monocytic transcriptional program reducing leukemia propagation in vivo. By combinatorial screening using H4 and molecules with defined targets, we demonstrated that H4 induces differentiation by the activation of the protein kinase C (PKC) signaling pathway, and in line with this, activates PKC phosphorylation and translocation of PKC to the cell membrane. Furthermore, the combinatorial screening identified a bromo- and extra-terminal domain (BET) inhibitor that could further improve H4-dependent leukemic differentiation in FLT3 wild-type monocytic AML. These findings illustrate the value of an unbiased, multiplex screening platform for developing combinatorial therapeutic approaches for AML.

Haploinsufficiency of cohesin protease, Separase, promotes regeneration of hematopoietic stem cells in mice

Cohesin recently emerged as a new regulator of hematopoiesis and leukemia. In addition to cohesin, whether proteins that regulate cohesin's function have any direct role in hematopoiesis and hematologic diseases have not been fully examined. Separase, encoded by the ESPL1 gene, is an important regulator of cohesin's function. Canonically, protease activity of Separase resolves sister chromatid cohesion by cleaving cohesin subunit-Rad21 at the onset of anaphase. Using a Separase haploinsufficient mouse model, we have uncovered a novel role of Separase in hematopoiesis. We report that partial disruption of Separase distinctly alters the functional characteristics of hematopoietic stem/progenitor cells (HSPCs). Although analyses of peripheral blood and bone marrow of Espl1+/Hyp mice broadly displayed unperturbed hematopoietic parameters during normal hematopoiesis, further probing of the composition of early hematopoietic cells in Espl1+/Hyp bone marrow revealed a mild reduction in the frequencies of the Lin- Sca1+ Kit- (LSK) or LSK CD48+ CD150- multipotent hematopoietic progenitors population without a significant change in either long-term or short-term hematopoietic stem cells (HSCs) subsets at steady state. Surprisingly, however, we found that Separase haploinsufficiency promotes regeneration activity of HSCs in serial in vivo repopulation assays. In vitro colony formation assays also revealed an enhanced serial replating capacity of hematopoietic progenitors isolated from Espl1+/Hyp mice. Microarray analysis of differentially expressed genes showed that Separase haploinsufficiency in HSCs (SP-KSL) leads to enrichment of gene signatures that are upregulated in HSCs compared to committed progenitors and mature cells. Taken together, our findings demonstrate a key role of Separase in promoting hematopoietic regeneration of HSCs.

A high throughput method for quantifying number and size distribution of Arabidopsis seeds using large particle flow cytometry

BACKGROUND

Seed size and number are important plant traits from an ecological and horticultural/agronomic perspective. However, in small-seeded species such as Arabidopsis thaliana, research on seed size and number is limited by the absence of suitable high throughput phenotyping methods.

RESULTS

We report on the development of a high throughput method for counting seeds and measuring individual seed sizes. The method uses a large-particle flow cytometer to count individual seeds and sort them according to size, allowing an average of 12,000 seeds/hour to be processed. To achieve this high throughput, post harvested seeds are first separated from remaining plant material (dust and chaff) using a rapid sedimentation-based method. Then, classification algorithms are used to refine the separation process in silico. Accurate identification of all seeds in the samples was achieved, with relative errors below 2%.

CONCLUSION

The tests performed reveal that there is no single classification algorithm that performs best for all samples, so the recommended strategy is to train and use multiple algorithms and use the median predictions of seed size and number across all algorithms. To facilitate the use of this method, an R package (SeedSorter) that implements the methodology has been developed and made freely available. The method was validated with seed samples from several natural accessions of Arabidopsis thaliana, but our analysis pipeline is applicable to any species with seed sizes smaller than 1.5 mm.

Transcriptomic changes in Echinochloa colona in response to treatment with the herbicide imazamox

Presented here is the first Echinochloa colona leaf transcriptome. Analysis of gene expression before and after herbicide treatment reveals that E. colona mounts a stress response upon exposure to herbicide. Herbicides are the most frequently used means of controlling weeds. For many herbicides, the target site is known; however, it is considerably less clear how plant gene expression changes in response to herbicide exposure. In this study, changes in gene expression in response to herbicide exposure in imazamox-sensitive (S) and- resistant (R) junglerice (Echinochloa colona L.) biotypes was examined. As no reference genome is available for this weed, a reference leaf transcriptome was generated. Messenger RNA was isolated from imazamox-treated- and untreated R and S plants and the resulting cDNA libraries were sequenced on an Illumina HiSeq2000. The transcriptome was assembled, annotated, and differential gene expression analysis was performed to identify transcripts that were upregulated or downregulated in response to herbicide exposure for both biotypes. Differentially expressed transcripts included transcription factors, protein-modifying enzymes, and enzymes involved in metabolism and signaling. A literature search revealed that members of the families represented in this analysis were known to be involved in abiotic stress response in other plants, suggesting that imazamox exposure induced a stress response. A time course study examining a subset of transcripts showed that expression peaked within 4-12 h and then returned to untreated levels within 48 h of exposure. Testing of plants from two additional biotypes showed a similar change in gene expression 4 h after herbicide exposure compared to the resistant and sensitive biotypes. This study shows that within 48 h junglerice mounts a stress response to imazamox exposure.

An expanded evaluation of protein function prediction methods shows an improvement in accuracy

BACKGROUND

A major bottleneck in our understanding of the molecular underpinnings of life is the assignment of function to proteins. While molecular experiments provide the most reliable annotation of proteins, their relatively low throughput and restricted purview have led to an increasing role for computational function prediction. However, assessing methods for protein function prediction and tracking progress in the field remain challenging.

RESULTS

We conducted the second critical assessment of functional annotation (CAFA), a timed challenge to assess computational methods that automatically assign protein function. We evaluated 126 methods from 56 research groups for their ability to predict biological functions using Gene Ontology and gene-disease associations using Human Phenotype Ontology on a set of 3681 proteins from 18 species. CAFA2 featured expanded analysis compared with CAFA1, with regards to data set size, variety, and assessment metrics. To review progress in the field, the analysis compared the best methods from CAFA1 to those of CAFA2.

CONCLUSIONS

The top-performing methods in CAFA2 outperformed those from CAFA1. This increased accuracy can be attributed to a combination of the growing number of experimental annotations and improved methods for function prediction. The assessment also revealed that the definition of top-performing algorithms is ontology specific, that different performance metrics can be used to probe the nature of accurate predictions, and the relative diversity of predictions in the biological process and human phenotype ontologies. While there was methodological improvement between CAFA1 and CAFA2, the interpretation of results and usefulness of individual methods remain context-dependent.

Medial HOXA genes demarcate haematopoietic stem cell fate during human development

Pluripotent stem cells (PSC) may provide a potential source of haematopoietic stem/progenitor cells (HSPCs) for transplantation; however, unknown molecular barriers prevent the self-renewal of PSC-HSPCs. Using two-step differentiation, human embryonic stem cells (hESCs) differentiated in vitro into multipotent haematopoietic cells that had CD34⁺CD38^−/loCD90⁺CD45⁺GPI-80⁺ foetal liver (FL) HSC immunophenotype, but displayed poor expansion potential and engraftment ability. Transcriptome analysis of immunophenotypic hESC-HSPCs revealed that, despite their molecular resemblance to FL-HSPCs, medial HOXA genes remained suppressed. Knockdown of HOXA7 disrupted FL-HSPC function and caused transcriptome dysregulation that resembled hESC-derived progenitors. Overexpression of medial HOXA genes prolonged FL-HSPC maintenance but was insufficient to confer self-renewal to hESC-HSPCs. Stimulation of retinoic acid signalling during endothelial-to-haematopoietic transition induced the HOXA cluster and other HSC/definitive haemogenic endothelium genes, and prolonged HSPC maintenance in culture. Thus, retinoic acid signalling-induced medial HOXA gene expression marks the establishment of the definitive HSC fate and controls HSC identity and function.

Scl binds to primed enhancers in mesoderm to regulate hematopoietic and cardiac fate divergence

Scl/Tal1 confers hemogenic competence and prevents ectopic cardiomyogenesis in embryonic endothelium by unknown mechanisms. We discovered that Scl binds to hematopoietic and cardiac enhancers that become epigenetically primed in multipotent cardiovascular mesoderm, to regulate the divergence of hematopoietic and cardiac lineages. Scl does not act as a pioneer factor but rather exploits a pre-established epigenetic landscape. As the blood lineage emerges, Scl binding and active epigenetic modifications are sustained in hematopoietic enhancers, whereas cardiac enhancers are decommissioned by removal of active epigenetic marks. Our data suggest that, rather than recruiting corepressors to enhancers, Scl prevents ectopic cardiogenesis by occupying enhancers that cardiac factors, such as Gata4 and Hand1, use for gene activation. Although hematopoietic Gata factors bind with Scl to both activated and repressed genes, they are dispensable for cardiac repression, but necessary for activating genes that enable hematopoietic stem/progenitor cell development. These results suggest that a unique subset of enhancers in lineage-specific genes that are accessible for regulators of opposing fates during the time of the fate decision provide a platform where the divergence of mutually exclusive fates is orchestrated.

Learning from nature: the use of non-model species to identify novel acclimations to flooding stress

Excess water in the form of waterlogged soil or deeper submergence (generically termed 'flooding') influences plant growth, survival and species distribution in many natural ecosystems. It also has a negative impact on crop growth and yield since many agricultural species are flooding intolerant. The often devastating effect of flooding on plant performance is related to its interference with gas exchange between the plant and its environment. This results in energy deficiency and carbohydrate starvation. In the near future, flooding frequency is expected to increase due to global climate change and the human population is expected to increase to ∼9 billion people by 2050. The need for increased agricultural productivity is self-evident and this will require a better mechanistic understanding of the interaction between plants and abiotic stresses such as flooding. We argue that, in seeking this understanding, we should not restrict the research to model species such as rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana). This is because some stress-tolerance mechanisms are not found in these species. Examples are given of how flooding tolerance is achieved by non-model species of Rumex and Rorippa that have evolved to cope with flooding in natural environments. These findings will add usefully to the spread of resources available to plant breeding programmes aimed at improving flooding tolerance in crop plants.

Abstract 356: Scl Represses Cardiogenesis via Distant Enhancers during Hemogenic Endothelium Specification

Understanding the mechanisms directing mesoderm specification holds a great potential to advance the development of cell-based therapies for cardiovascular and blood disorders. The bHLH transcription factor Scl is known as the master regulator of the hematopoietic fate. We recently discovered that, in addition to its critical function in promoting the establishment of hemogenic endothelium during hematopoietic stem/progenitor cell (HS/PC) development, Scl is also required to repress cardiomyogenesis in endothelium in hematopoietic tissues and endocardium in the heart. However, the mechanisms for the cardiac repression have remained unknown.

Using ChIP-sequencing and microarray analysis of Flk+ mesoderm differentiated from mouse ES cells, we show that Scl both directly activates a broad gene regulatory network required for hemogenic endothelium and HS/PC development (e.g. Runx1, cMyb, Lyl1, Mef2C, Sox7 etc.), and directly represses transcriptional regulators required for cardiogenesis (e.g. Gata4, Gata6, Myocd, etc.) and mesoderm development (Eomes, Mixl1, Etv2, etc.). Repression of cardiac and mesodermal programs occurs during a short developmental window through Scl binding to distant enhancers, while binding to hematopoietic regulators extends throughout HS/PC and red blood cell development and encompasses both distant and proximal binding sites. We also discovered that, surprisingly, Scl complex partners Gata 1 and 2 are dispensable for hematopoietic vs. cardiac specification and Scl binding to majority of its target genes. Nevertheless, Gata factors co-operate with Scl to activate selected transcription factors that facilitate HS/PC emergence from hemogenic endothelium. These results denote Scl as a true master regulator of hematopoietic vs. cardiac fate choice and suggest a mechanism by which lineage-specific bHLH factors direct the divergence of competing fates.

Abstract 290: Haemogenic Endocardium Contribute To Definitive Hematopoiesis During Cardiogenesis

The circulatory system is the first functional organ system that develops during mammalian life. Accumulating evidences suggest that cardiac and endocardial cells can arise from a single common progenitor cell during mammalian cardiogenesis. Notably, these early cardiac progenitors express multiple hematopoietic transcription factors, consistent with previous reports. Indeed, a close relationship among cardiac, endocardial and hematopoietic lineages has been suggested in fly, zebrafish, and embryonic stem cell in vitro differentiation models. However, it is unclear when, where and how this hematopoietic gene program is in operation during in vivo mammalian cardiogenesis. Hematopoietic colony assay suggests that mouse heart explants generate myeloids and erythroids in the absence of circulation, suggesting that the heart tube is a de novo site for the definitive hematopoiesis. Lineage tracing revealed that putative cardiac-derived Nkx2-5+/Isl1+ endocardial cells give rise to CD41+ hematopoietic progenitors that contribute to definitive hematopoiesis in vivo and ex vivo during embryogenesis earlier than in the AGM region. Furthermore, Nkx2-5 and Isl1 are both required for the hemogenic activity of the endocardium. Together, identification of Nkx2-5/Isl1-dependent hemogenic endocardial cells (1) adds hematopoietic component in the cardiogenesis lineage tree, (2) changes the long-held dogma that AGM is the only major source of definitive hematopoiesis in the embryo proper, and (3) represents phylogenetically conserved fundamental mechanism of cardio-vasculo-hematopoietic differentiation pathway during the development of circulatory system.

Ethylene--and oxygen signalling--drive plant survival during flooding

Flooding is a widely occurring environmental stress both for natural and cultivated plant species. The primary problems associated with flooding arise due to restricted gas diffusion underwater. This hampers gas exchange needed for the critical processes of photosynthesis and respiration. Plant acclimation to flooding includes the adaptation of a suite of traits that helps alleviate or avoid these stressful conditions and improves or restores exchange of O2 and CO2 . The manifestation of these traits is, however, reliant on the timely perception of signals that convey the underwater status. Flooding-associated reduced gas diffusion imposes a drastic change in the internal gas composition within submerged plant organs. One of the earliest changes is an increase in the levels of the gaseous plant hormone ethylene. Depending on the species, organ, flooding conditions and time of the day, plants will also subsequently experience a reduction in oxygen levels. This review provides a comprehensive overview on the roles of ethylene and oxygen as critical signals of flooding stress. It includes a discussion of the dynamics of these gases in plants when underwater, their interaction, current knowledge of their perception mechanisms and the resulting downstream changes that mediate important acclimative processes that allow endurance and survival under flooded conditions.

Expansion on stromal cells preserves the undifferentiated state of human hematopoietic stem cells despite compromised reconstitution ability

Lack of HLA-matched hematopoietic stem cells (HSC) limits the number of patients with life-threatening blood disorders that can be treated by HSC transplantation. So far, insufficient understanding of the regulatory mechanisms governing human HSC has precluded the development of effective protocols for culturing HSC for therapeutic use and molecular studies. We defined a culture system using OP9M2 mesenchymal stem cell (MSC) stroma that protects human hematopoietic stem/progenitor cells (HSPC) from differentiation and apoptosis. In addition, it facilitates a dramatic expansion of multipotent progenitors that retain the immunophenotype (CD34+CD38-CD90+) characteristic of human HSPC and proliferative potential over several weeks in culture. In contrast, transplantable HSC could be maintained, but not significantly expanded, during 2-week culture. Temporal analysis of the transcriptome of the ex vivo expanded CD34+CD38-CD90+ cells documented remarkable stability of most transcriptional regulators known to govern the undifferentiated HSC state. Nevertheless, it revealed dynamic fluctuations in transcriptional programs that associate with HSC behavior and may compromise HSC function, such as dysregulation of PBX1 regulated genetic networks. This culture system serves now as a platform for modeling human multilineage hematopoietic stem/progenitor cell hierarchy and studying the complex regulation of HSC identity and function required for successful ex vivo expansion of transplantable HSC.

Haemogenic endocardium contributes to transient definitive haematopoiesis

Haematopoietic cells arise from spatiotemporally restricted domains in the developing embryo. Although studies of non-mammalian animal and in vitro embryonic stem cell models suggest a close relationship among cardiac, endocardial and haematopoietic lineages, it remains unknown whether the mammalian heart tube serves as a haemogenic organ akin to the dorsal aorta. Here we examine the haemogenic activity of the developing endocardium. Mouse heart explants generate myeloid and erythroid colonies in the absence of circulation. Haemogenic activity arises from a subset of endocardial cells in the outflow cushion and atria earlier than in the aorta-gonad-mesonephros region, and is transient and definitive in nature. Interestingly, key cardiac transcription factors, Nkx2-5 and Isl1, are expressed in and required for the haemogenic population of the endocardium. Together, these data suggest that a subset of endocardial/endothelial cells serve as a de novo source for transient definitive haematopoietic progenitors.

Scl represses cardiomyogenesis in prospective hemogenic endothelium and endocardium

Endothelium in embryonic hematopoietic tissues generates hematopoietic stem/progenitor cells; however, it is unknown how its unique potential is specified. We show that transcription factor Scl/Tal1 is essential for both establishing the hematopoietic transcriptional program in hemogenic endothelium and preventing its misspecification to a cardiomyogenic fate. Scl(-/-) embryos activated a cardiac transcriptional program in yolk sac endothelium, leading to the emergence of CD31+Pdgfrα+ cardiogenic precursors that generated spontaneously beating cardiomyocytes. Ectopic cardiogenesis was also observed in Scl(-/-) hearts, where the disorganized endocardium precociously differentiated into cardiomyocytes. Induction of mosaic deletion of Scl in Scl(fl/fl)Rosa26Cre-ER(T2) embryos revealed a cell-intrinsic, temporal requirement for Scl to prevent cardiomyogenesis from endothelium. Scl(-/-) endothelium also upregulated the expression of Wnt antagonists, which promoted rapid cardiomyocyte differentiation of ectopic cardiogenic cells. These results reveal unexpected plasticity in embryonic endothelium such that loss of a single master regulator can induce ectopic cardiomyogenesis from endothelial cells.

GFam: a platform for automatic annotation of gene families

We have developed GFam, a platform for automatic annotation of gene/protein families. GFam provides a framework for genome initiatives and model organism resources to build domain-based families, derive meaningful functional labels and offers a seamless approach to propagate functional annotation across periodic genome updates. GFam is a hybrid approach that uses a greedy algorithm to chain component domains from InterPro annotation provided by its 12 member resources followed by a sequence-based connected component analysis of un-annotated sequence regions to derive consensus domain architecture for each sequence and subsequently generate families based on common architectures. Our integrated approach increases sequence coverage by 7.2 percentage points and residue coverage by 14.6 percentage points higher than the coverage relative to the best single-constituent database within InterPro for the proteome of Arabidopsis. The true power of GFam lies in maximizing annotation provided by the different InterPro data sources that offer resource-specific coverage for different regions of a sequence. GFam’s capability to capture higher sequence and residue coverage can be useful for genome annotation, comparative genomics and functional studies. GFam is a general-purpose software and can be used for any collection of protein sequences. The software is open source and can be obtained from http://www.paccanarolab.org/software/gfam/.

Reorganization of the host epigenome by a viral oncogene

Adenovirus small e1a oncoprotein causes ∼70% reduction in cellular levels of histone H3 lysine 18 acetylation (H3K18ac). It is unclear, however, where this dramatic reduction occurs genome-wide. ChIP-sequencing revealed that by 24 h after expression, e1a erases 95% of H3K18ac peaks in normal, contact-inhibited fibroblasts and replaces them with one-third as many at new genomic locations. The H3K18ac peaks at promoters and intergenic regions of genes with fibroblast-related functions are eliminated after infection, and new H3K18ac peaks are established at promoters of highly induced genes that regulate cell cycling and at new putative enhancers. Strikingly, the regions bound by the retinoblastoma family of proteins in contact-inhibited fibroblasts gain new peaks of H3K18ac in the e1a-expressing cells, including 55% of RB1-bound loci. In contrast, over half of H3K9ac peaks are similarly distributed before and after infection, independently of RB1. The strategic redistribution of H3K18ac by e1a highlights the importance of this modification for transcriptional activation and cellular transformation as well as functional differences between the RB-family member proteins.

The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools

The Arabidopsis Information Resource (TAIR, http://arabidopsis.org) is a genome database for Arabidopsis thaliana, an important reference organism for many fundamental aspects of biology as well as basic and applied plant biology research. TAIR serves as a central access point for Arabidopsis data, annotates gene function and expression patterns using controlled vocabulary terms, and maintains and updates the A. thaliana genome assembly and annotation. TAIR also provides researchers with an extensive set of visualization and analysis tools. Recent developments include several new genome releases (TAIR8, TAIR9 and TAIR10) in which the A. thaliana assembly was updated, pseudogenes and transposon genes were re-annotated, and new data from proteomics and next generation transcriptome sequencing were incorporated into gene models and splice variants. Other highlights include progress on functional annotation of the genome and the release of several new tools including Textpresso for Arabidopsis which provides the capability to carry out full text searches on a large body of research literature.

Natural variation of submergence tolerance among Arabidopsis thaliana accessions

• The exploitation of natural variation in Arabidopsis thaliana (Arabidopsis) provides a huge potential for the identification of the molecular mechanisms underlying this variation as a result of the availability of a vast array of genetic and genomic resources for this species. Eighty-six Arabidopsis accessions were screened for natural variation in flooding tolerance. This forms the first step towards the identification and characterization of the role of candidate genes contributing to flooding tolerance. • Arabidopsis accessions at the 10-leaf stage were subjected to complete submergence in the dark. Survival curves were plotted to estimate median lethal times as a measure of tolerance. Flooding-associated survival parameters, such as root and shoot oxygen content, initial carbohydrate content and petiole elongation under water, were also measured. • There was a significant variation in submergence tolerance among Arabidopsis accessions. However, the order of tolerance did not correlate with root and shoot oxygen content or initial amounts of shoot starch and total soluble sugars. A negative correlation was observed between submergence tolerance and underwater petiole elongation. • Arabidopsis accessions show considerable variation in the ability to tolerate complete submergence, making it a good species in which to identify and characterize genes and to study mechanisms that contribute to survival under water.

Integrative analysis of the Caenorhabditis elegans genome by the modENCODE project

We systematically generated large-scale data sets to improve genome annotation for the nematode Caenorhabditis elegans, a key model organism. These data sets include transcriptome profiling across a developmental time course, genome-wide identification of transcription factor-binding sites, and maps of chromatin organization. From this, we created more complete and accurate gene models, including alternative splice forms and candidate noncoding RNAs. We constructed hierarchical networks of transcription factor-binding and microRNA interactions and discovered chromosomal locations bound by an unusually large number of transcription factors. Different patterns of chromatin composition and histone modification were revealed between chromosome arms and centers, with similarly prominent differences between autosomes and the X chromosome. Integrating data types, we built statistical models relating chromatin, transcription factor binding, and gene expression. Overall, our analyses ascribed putative functions to most of the conserved genome.

Comparison and calibration of transcriptome data from RNA-Seq and tiling arrays

Background

Tiling arrays have been the tool of choice for probing an organism's transcriptome without prior assumptions about the transcribed regions, but RNA-Seq is becoming a viable alternative as the costs of sequencing continue to decrease. Understanding the relative merits of these technologies will help researchers select the appropriate technology for their needs.

Results

Here, we compare these two platforms using a matched sample of poly(A)-enriched RNA isolated from the second larval stage of C. elegans. We find that the raw signals from these two technologies are reasonably well correlated but that RNA-Seq outperforms tiling arrays in several respects, notably in exon boundary detection and dynamic range of expression. By exploring the accuracy of sequencing as a function of depth of coverage, we found that about 4 million reads are required to match the sensitivity of two tiling array replicates. The effects of cross-hybridization were analyzed using a "nearest neighbor" classifier applied to array probes; we describe a method for determining potential "black list" regions whose signals are unreliable. Finally, we propose a strategy for using RNA-Seq data as a gold standard set to calibrate tiling array data. All tiling array and RNA-Seq data sets have been submitted to the modENCODE Data Coordinating Center.

Conclusions

Tiling arrays effectively detect transcript expression levels at a low cost for many species while RNA-Seq provides greater accuracy in several regards. Researchers will need to carefully select the technology appropriate to the biological investigations they are undertaking. It will also be important to reconsider a comparison such as ours as sequencing technologies continue to evolve.

Pregnancy and breast cancer

Breast cancer is one of the most commonly diagnosed malignancies during pregnancy. Pregnancy-associated breast cancer (PABC) presents a challenging clinical situation. This article reviews the current evidence around the management of PABC and the safety of pregnancy after breast cancer. The trend towards later age at first childbirth has resulted in an increase in the number of breast cancer cases coexistent with pregnancy. The management of breast cancer during pregnancy requires a multidisciplinary team approach. Breast surgery can be safely performed during any trimester of pregnancy. Radiation therapy, if required, must be delayed until after delivery. The majority of patients with PABC require chemotherapy. The timing of delivery in relation to chemotherapy administration should be carefully considered. There is no evidence to date that pregnancy termination influences overall survival for the mother. To date, there is no clear evidence that subsequent pregnancy after breast cancer is associated with worse maternal survival. There is a suggestion that subsequent pregnancy may in fact be associated with an improved survival. However, the available studies are limited by potential biases.

SCPS: a fast implementation of a spectral method for detecting protein families on a genome-wide scale

Background

An important problem in genomics is the automatic inference of groups of homologous proteins from pairwise sequence similarities. Several approaches have been proposed for this task which are "local" in the sense that they assign a protein to a cluster based only on the distances between that protein and the other proteins in the set. It was shown recently that global methods such as spectral clustering have better performance on a wide variety of datasets. However, currently available implementations of spectral clustering methods mostly consist of a few loosely coupled Matlab scripts that assume a fair amount of familiarity with Matlab programming and hence they are inaccessible for large parts of the research community.

Results

SCPS (Spectral Clustering of Protein Sequences) is an efficient and user-friendly implementation of a spectral method for inferring protein families. The method uses only pairwise sequence similarities, and is therefore practical when only sequence information is available. SCPS was tested on difficult sets of proteins whose relationships were extracted from the SCOP database, and its results were extensively compared with those obtained using other popular protein clustering algorithms such as TribeMCL, hierarchical clustering and connected component analysis. We show that SCPS is able to identify many of the family/superfamily relationships correctly and that the quality of the obtained clusters as indicated by their F-scores is consistently better than all the other methods we compared it with. We also demonstrate the scalability of SCPS by clustering the entire SCOP database (14,183 sequences) and the complete genome of the yeast Saccharomyces cerevisiae (6,690 sequences).

Conclusions

Besides the spectral method, SCPS also implements connected component analysis and hierarchical clustering, it integrates TribeMCL, it provides different cluster quality tools, it can extract human-readable protein descriptions using GI numbers from NCBI, it interfaces with external tools such as BLAST and Cytoscape, and it can produce publication-quality graphical representations of the clusters obtained, thus constituting a comprehensive and effective tool for practical research in computational biology. Source code and precompiled executables for Windows, Linux and Mac OS X are freely available at http://www.paccanarolab.org/software/scps.

An approach to comparing tiling array and high throughput sequencing technologies for genomic transcript mapping

Background

There are two main technologies for transcriptome profiling, namely, tiling microarrays and high-throughput sequencing. Recently there has been a tremendous amount of excitement about the latter because of the advent of next-generation sequencing technologies and its promises. Consequently, the question of the moment is how these two technologies compare. Here we attempt to develop an approach to do a fair comparison of transcripts identified from tiling microarray and MPSS sequencing data.

Findings

This comparison is a challenging task because the sequencing data is discrete while the tiling array data is continuous. We use the published rice and Arabidopsis datasets which provide currently best matched sets of arrays and sequencing experiments using a slightly earlier generation of sequencing, the MPSS tag sequencing technology. After scoring the arrays consistently in both the organisms, a first pass comparison reveals a surprisingly small overlap in transcripts of 22% and 66% respectively, in rice and Arabidopsis. However, when we do the analysis in detail, we find that this is an underestimate. In particular, when we map the probe intensities onto the sequencing tags and then look at their intensity distribution, we see that they are very similar to exons. Furthermore, restricting our comparison to only protein-coding gene loci revealed a very good overlap between the two technologies.

Conclusion

Our approach to compare genome tiling microarray and MPSS sequencing data suggests that there is actually a reasonable overlap in transcripts identified by the two technologies. This overlap is distorted by the scoring and thresholding in the tiling array scoring procedure.

An approach to compare genome tiling microarray and MPSS sequencing data for transcript mapping

We are correcting the abstract of our published article ([1]). The sentence that starts "We observe that 4.5% of MPSS tags...." was not scientifically complete in the original abstract, having only two of the four numbers required to describe a comparison of two technologies in two different organisms. The abstract below more accurately describes our findings, as documented in Figure 1 of the manuscript.

Transmembrane protein oxygen content and compartmentalization of cells

Recently, there was a report that explored the oxygen content of transmembrane proteins over macroevolutionary time scales where the authors observed a correlation between the geological time of appearance of compartmentalized cells with atmospheric oxygen concentration. The authors predicted, characterized and correlated the differences in the structure and composition of transmembrane proteins from the three kingdoms of life with atmospheric oxygen concentrations in geological timescale. They hypothesized that transmembrane proteins in ancient taxa were selectively excluding oxygen and as this constraint relaxed over time with increase in the levels of atmospheric oxygen the size and number of communication-related transmembrane proteins increased. In summary, they concluded that compartmentalized and non-compartmentalized cells can be distinguished by how oxygen is partitioned at the proteome level. They derived this conclusion from an analysis of 19 taxa. We extended their analysis on a larger sample of taxa comprising 309 eubacterial, 34 archaeal, and 30 eukaryotic complete proteomes and observed that one can not absolutely separate the two groups of cells based on partition of oxygen in their membrane proteins. In addition, the origin of compartmentalized cells is likely to have been driven by an innovation than happened 2700 million years ago in the membrane composition of cells that led to the evolution of endocytosis and exocytosis rather than due to the rise in concentration of atmospheric oxygen.

The selection of acceptable protein mutations

We have determined the general constraints that govern sequence divergence in proteins that retain entirely, or very largely, the same structure and function. To do this we collected data from three different groups of orthologous sequences: those found in humans and mice, in humans and chickens, and in Escherichia coli and Salmonella enterica. In total, these organisms have 21,738 suitable pairs of orthologs, and these contain nearly 2 million mutations. The three groups differ greatly in the taxa from which they come and/or in the time that separates them from their last common ancestor. Nevertheless, the results we obtain from the three different groups are strikingly similar. For each group, the orthologous sequence pairs were assigned to six different divergence categories on the basis of their sequence identities. For categories with the same divergence, common accepted mutations have similar frequencies and rank orders in the three groups. With divergence, the width of the range of common mutations grows in the same manner in each group. We examined the distribution of mutations in protein structures. With increasing divergence, mutations increase at different rates in the buried, intermediate, and exposed regions of protein structures in a manner that explains the exponential relationship between the divergence of structure and sequence. This work implies that commonly allowed mutations are selected by a set of general constraints that are well defined and whose nature varies with divergence.

Global identification and characterization of transcriptionally active regions in the rice genome

Genome tiling microarray studies have consistently documented rich transcriptional activity beyond the annotated genes. However, systematic characterization and transcriptional profiling of the putative novel transcripts on the genome scale are still lacking. We report here the identification of 25,352 and 27,744 transcriptionally active regions (TARs) not encoded by annotated exons in the rice (Oryza. sativa) subspecies japonica and indica, respectively. The non-exonic TARs account for approximately two thirds of the total TARs detected by tiling arrays and represent transcripts likely conserved between japonica and indica. Transcription of 21,018 (83%) japonica non-exonic TARs was verified through expression profiling in 10 tissue types using a re-array in which annotated genes and TARs were each represented by five independent probes. Subsequent analyses indicate that about 80% of the japonica TARs that were not assigned to annotated exons can be assigned to various putatively functional or structural elements of the rice genome, including splice variants, uncharacterized portions of incompletely annotated genes, antisense transcripts, duplicated gene fragments, and potential non-coding RNAs. These results provide a systematic characterization of non-exonic transcripts in rice and thus expand the current view of the complexity and dynamics of the rice transcriptome.

Domain insertions in protein structures

Domains are the structural, functional or evolutionary units of proteins. Proteins can comprise a single domain or a combination of domains. In multi-domain proteins, the domains almost always occur end-to-end, i.e., one domain follows the C-terminal end of another domain. However, there are exceptions to this common pattern, where multi-domain proteins are formed by insertion of one domain (insert) into another domain (parent). Here, we provide a quantitative description of known insertions in the Protein Data Bank (PDB). We found that 9% of domain combinations observed in non-redundant PDB are insertions. Although 90% of all insertions involve only one insert, proteins can clearly have multiple (nested, two-domain and three-domain) inserts. We also observed correlations between the structure and function of a domain and its tendency to be found as a parent or an insert. There is a bias in insert position towards the C terminus of parents. We observed that the atomic distance between the N and C terminus of an insert is significantly smaller when compared to the N-to-C distance in a parent context or a single domain context. Insertions are found always to occur in loop regions of parent domains. Our observations regarding the relationship between domain insertions and the structure, function and evolution of proteins have implications for protein engineering.

DomIns: a web resource for domain insertions in known protein structures

Proteins can be formed by single or multiple domains. The process of recombination at the molecular level has generated a wide variety of multi-domain proteins with specific domain organization to cater to the functional requirements of an organism. The functional and structural costs of inserting a domain into another means that multi-domain proteins are usually formed by covalently linking the N-terminus of one domain to the C-terminus of the preceding domain. While this is true in a large proportion of multi-domain proteins, we find a significant fraction of proteins that are the result of domain insertion. The inserted domain breaks the sequence contiguity of the domain into which it is inserted leading to a novel domain organization. This web resource aims to document domain insertions in known protein structures that are classified in the SCOP database. The web server can be accessed from http://stash.mrc-lmb.cam. ac.uk/DomIns/.