tropiTree: An NGS-Based EST-SSR Resource for 24 Tropical Tree Species

The development of genetic tools for non-model organisms has been hampered by cost, but advances in next-generation sequencing (NGS) have created new opportunities. In ecological research, this raises the prospect for developing molecular markers to simultaneously study important genetic processes such as gene flow in multiple non-model plant species within complex natural and anthropogenic landscapes. Here, we report the use of bar-coded multiplexed paired-end Illumina NGS for the de novo development of expressed sequence tag-derived simple sequence repeat (EST-SSR) markers at low cost for a range of 24 tree species. Each chosen tree species is important in complex tropical agroforestry systems where little is currently known about many genetic processes. An average of more than 5,000 EST-SSRs was identified for each of the 24 sequenced species, whereas prior to analysis 20 of the species had fewer than 100 nucleotide sequence citations. To make results available to potential users in a suitable format, we have developed an open-access, interactive online database, tropiTree (http://bioinf.hutton.ac.uk/tropiTree), which has a range of visualisation and search facilities, and which is a model for the efficient presentation and application of NGS data.     

NCI Workshop Report: Clinical and Computational Requirements for Correlating Imaging Phenotypes with Genomics Signatures

The National Cancer Institute (NCI) Cancer Imaging Program organized two related workshops on June 26–27, 2013, entitled “Correlating Imaging Phenotypes with Genomics Signatures Research” and “Scalable Computational Resources as Required for Imaging-Genomics Decision Support Systems.” The first workshop focused on clinical and scientific requirements, exploring our knowledge of phenotypic characteristics of cancer biological properties to determine whether the field is sufficiently advanced to correlate with imaging phenotypes that underpin genomics and clinical outcomes, and exploring new scientific methods to extract phenotypic features from medical images and relate them to genomics analyses. The second workshop focused on computational methods that explore informatics and computational requirements to extract phenotypic features from medical images and relate them to genomics analyses and improve the accessibility and speed of dissemination of existing NIH resources. These workshops linked clinical and scientific requirements of currently known phenotypic and genotypic cancer biology characteristics with imaging phenotypes that underpin genomics and clinical outcomes. The group generated a set of recommendations to NCI leadership and the research community that encourage and support development of the emerging radiogenomics research field to address short-and longer-term goals in cancer research.     

Homologues of potato chromosome 5 show variable collinearity in the euchromatin,but dramatic absence of sequence similarity in the pericentromeric heterochromatin

Background

In flowering plants it has been shown that de novo genome assemblies of different species and genera show a significant drop in the proportion of alignable sequence. Within a plant species, however, it is assumed that different haplotypes of the same chromosome align well. In this paper we have compared three de novo assemblies of potato chromosome 5 and report on the sequence variation and the proportion of sequence that can be aligned.

Results

For the diploid potato clone RH89-039-16 (RH) we produced two linkage phase controlled and haplotype-specific assemblies of chromosome 5 based on BAC-by-BAC sequencing, which were aligned to each other and compared to the 52 Mb chromosome 5 reference sequence of the doubled monoploid clone DM 1–3 516 R44 (DM). We identified 17.0 Mb of non-redundant sequence scaffolds derived from euchromatic regions of RH and 38.4 Mb from the pericentromeric heterochromatin. For 32.7 Mb of the RH sequences the correct position and order on chromosome 5 was determined, using genetic markers, fluorescence in situ hybridisation and alignment to the DM reference genome. This ordered fraction of the RH sequences is situated in the euchromatic arms and in the heterochromatin borders. In the euchromatic regions, the sequence collinearity between the three chromosomal homologs is good, but interruption of collinearity occurs at nine gene clusters. Towards and into the heterochromatin borders, absence of collinearity due to structural variation was more extensive and was caused by hemizygous and poorly aligning regions of up to 450 kb in length. In the most central heterochromatin, a total of 22.7 Mb sequence from both RH haplotypes remained unordered. These RH sequences have very few syntenic regions and represent a non-alignable region between the RH and DM heterochromatin haplotypes of chromosome 5.

Conclusions

Our results show that among homologous potato chromosomes large regions are present with dramatic loss of sequence collinearity. This stresses the need for more de novo reference assemblies in order to capture genome diversity in this crop. The discovery of three highly diverged pericentric heterochromatin haplotypes within one species is a novelty in plant genome analysis. The possible origin and cytogenetic implication of this heterochromatin haplotype diversity are discussed.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1578-1) contains supplementary material, which is available to authorized users.     

A new H-2.1-PositiveD region allele,D dx,controlling two molecules,H-2Ddx and H-2Ldx

The inbred strains GRS/A and LIS/A carry the haplotypeH-2 ^dx , which had earlier been shown to have theK ^d ,I ^f ,S ^f , andG ^f alleles and a previously unknownD region allele,D ^dx . We show here that theD ^dx allele determines a new private specificity, H-2.63, is H-2.28 negative, and determines at least one public specificity of the H-2.1 family. It is thus a second example (afterD ^k ) of a H-2.1-positive H-2.28-negativeD region allele. Capping experiments show that the D^dx product comprises two molecules: H-2D^dx bearing the private specificity H-2.63, and H-2L^dx, which is H-2.63-negative but reacts with sera against the H-2.1 family of specificities. SDS gel electrophoresis of detergent-solubilized immunoprecipitated D^dx products shows that the H-2L^dx antigen has a molecular weight of approximately 45,000 daltons and is associated with a smaller polypeptide (mol. wt. 12,000).     

Aphid populations showing differential levels of virulence on Capsicum accessions

The green peach aphid,Myzus persicae,is one of the most threatening pests in pepper cultivation and growers would benefit from resistant varietices.Previously,we identified two Capsicum acessions as susceptible and three as resistant to M.persicae using an aphid population originating from the Netherlands(NL).Later on we identified an aphid population originating from a diferent gcographical region(Switserland,SW)that was virulent on all tested Capsicum acessions.The objeetive of the current work is to describe in detail diferent aspects of the interaction between two aphid populations and two sclected Capsicum acessions(one that was susceptible[PB2013046]and one that was resistant[PB2013071]to population NL),including biochemical processes involved.Electrical penetration graph(EPG)recordings showed similar feeding activities for both aphid populations on PB2013046.On acession PB2013071 the aphid population sw was able to devote significantly more time to phloem ingestion than population NL.We also studied plant defense response and found that plants of acession PB2013046 could not induce an accumulation of reactive oxygen species and callose formation after infestation with either aphid population.However,plants of PB2013071 induced a stronger defense response after infestation by population NL than after infestation by population SW.Based on these results,population SW of M.persicae seems to have overcome the resistance of PB2013071 that prevented feeding of aphids from NL population.The potential mechanism by which SW population overcomes the resistance is discussed.     

Mutability of bacteriophage M13 by ultraviolet light: Role of pyrimidine dimers

Summary The role of pyrimidine dimers in mutagenesis by ultraviolet light was examined by measuring the UV-induced reversion of six different bacteriophage M13 amber mutants for which the neighboring DNA sequences are known. The mutational response at amber (TAG) codons preceded by a guanine or adenine (where no pyrimidine dimer can be formed) were compared with those preceded by thymine or cytosine (where dimer formation is possible). Equivalent levels of UV-induced mutagenesis were observed at both kinds of sites. This observation demonstrates that there is no requirement for a pyrimidine dimer directly at the site of UV-induced mutation in this single-stranded DNA phage. UV irradiation of the phage was also performed in the presence of Ag⁺ ions, which specifically sensitize the DNA to dimer formation. The two methods of irradiation, when compared at equal survival levels (and presumably equal dimer frequencies), produced equivalent frequencies of reversion of the amber phage. We believe these results indicate that while the presence of pyrimidine dimers may be a prerequisite for UV mutagenesis, the actual mutagenic event can occur at a site some distance removed from a dimer.     

Mechanical Cell-Matrix Feedback Explains Pairwise and Collective Endothelial Cell Behavior In Vitro

In vitro cultures of endothelial cells are a widely used model system of the collective behavior of endothelial cells during vasculogenesis and angiogenesis. When seeded in an extracellular matrix, endothelial cells can form blood vessel-like structures, including vascular networks and sprouts. Endothelial morphogenesis depends on a large number of chemical and mechanical factors, including the compliancy of the extracellular matrix, the available growth factors, the adhesion of cells to the extracellular matrix, cell-cell signaling, etc. Although various computational models have been proposed to explain the role of each of these biochemical and biomechanical effects, the understanding of the mechanisms underlying in vitro angiogenesis is still incomplete. Most explanations focus on predicting the whole vascular network or sprout from the underlying cell behavior, and do not check if the same model also correctly captures the intermediate scale: the pairwise cell-cell interactions or single cell responses to ECM mechanics. Here we show, using a hybrid cellular Potts and finite element computational model, that a single set of biologically plausible rules describing (a) the contractile forces that endothelial cells exert on the ECM, (b) the resulting strains in the extracellular matrix, and (c) the cellular response to the strains, suffices for reproducing the behavior of individual endothelial cells and the interactions of endothelial cell pairs in compliant matrices. With the same set of rules, the model also reproduces network formation from scattered cells, and sprouting from endothelial spheroids. Combining the present mechanical model with aspects of previously proposed mechanical and chemical models may lead to a more complete understanding of in vitro angiogenesis.     

Genomic determinants of protein folding thermodynamics in prokaryotic organisms

Here we investigate how thermodynamic properties of orthologous proteins are influenced by the genomic environment in which they evolve. We performed a comparative computational study of 21 protein families in 73 prokaryotic species and obtained the following main results. (i) Protein stability with respect to the unfolded state and with respect to misfolding are anticorrelated. There appears to be a trade-off between these two properties, which cannot be optimized simultaneously. (ii) Folding thermodynamic parameters are strongly correlated with two genomic features, genome size and G+C composition. In particular, the normalized energy gap, an indicator of folding efficiency in statistical mechanical models of protein folding, is smaller in proteins of organisms with a small genome size and a compositional bias towards A+T. Such genomic features are characteristic for bacteria with an intracellular lifestyle. We interpret these correlations in light of mutation pressure and natural selection. A mutational bias toward A+T at the DNA level translates into a mutational bias toward more hydrophobic (and in general more interactive) proteins, a consequence of the structure of the genetic code. Increased hydrophobicity renders proteins more stable against unfolding but less stable against misfolding. Proteins with high hydrophobicity and low stability against misfolding occur in organisms with reduced genomes, like obligate intracellular bacteria. We argue that they are fixed because these organisms experience weaker purifying selection due to their small effective population sizes. This interpretation is supported by the observation of a high expression level of chaperones in these bacteria. Our results indicate that the mutational spectrum of a genome and the strength of selection significantly influence protein folding thermodynamics.