Whole-genome sequencing and assembly with high-throughput, short-read technologies

While recently developed short-read sequencing technologies may dramatically reduce the sequencing cost and eventually achieve the $1000 goal for re-sequencing, their limitations prevent the de novo sequencing of eukaryotic genomes with the standard shotgun sequencing protocol. We present SHRAP (SHort Read Assembly Protocol), a sequencing protocol and assembly methodology that utilizes high-throughput short-read technologies. We describe a variation on hierarchical sequencing with two crucial differences: (1) we select a clone library from the genome randomly rather than as a tiling path and (2) we sample clones from the genome at high coverage and reads from the clones at low coverage. We assume that 200 bp read lengths with a 1% error rate and inexpensive random fragment cloning on whole mammalian genomes is feasible. Our assembly methodology is based on first ordering the clones and subsequently performing read assembly in three stages: (1) local assemblies of regions significantly smaller than a clone size, (2) clone-sized assemblies of the results of stage 1, and (3) chromosome-sized assemblies. By aggressively localizing the assembly problem during the first stage, our method succeeds in assembling short, unpaired reads sampled from repetitive genomes. We tested our assembler using simulated reads from D. melanogaster and human chromosomes 1, 11, and 21, and produced assemblies with large sets of contiguous sequence and a misassembly rate comparable to other draft assemblies. Tested on D. melanogaster and the entire human genome, our clone-ordering method produces accurate maps, thereby localizing fragment assembly and enabling the parallelization of the subsequent steps of our pipeline. Thus, we have demonstrated that truly inexpensive de novo sequencing of mammalian genomes will soon be possible with high-throughput, short-read technologies using our methodology.     

Sequence similarity-driven proteomics in organisms with unknown genomes by LC-MS/MS and automated de novo sequencing

LC-MS/MS analysis on a linear ion trap LTQ mass spectrometer, combined with data processing, stringent, and sequence-similarity database searching tools, was employed in a layered manner to identify proteins in organisms with unsequenced genomes. Highly specific stringent searches (MASCOT) were applied as a first layer screen to identify either known (i.e. present in a database) proteins, or unknown proteins sharing identical peptides with related database sequences. Once the confidently matched spectra were removed, the remainder was filtered against a nonannotated library of background spectra that cleaned up the dataset from spectra of common protein and chemical contaminants. The rectified spectral dataset was further subjected to rapid batch de novo interpretation by PepNovo software, followed by the MS BLAST sequence-similarity search that used multiple redundant and partially accurate candidate peptide sequences. Importantly, a single dataset was acquired at the uncompromised sensitivity with no need of manual selection of MS/MS spectra for subsequent de novo interpretation. This approach enabled a completely automated identification of novel proteins that were, otherwise, missed by conventional database searches.     

Transglutaminase 2 undergoes a large conformational change upon activation

Human transglutaminase 2 (TG2), a member of a large family of enzymes that catalyze protein crosslinking, plays an important role in the extracellular matrix biology of many tissues and is implicated in the gluten-induced pathogenesis of celiac sprue. Although vertebrate transglutaminases have been studied extensively, thus far all structurally characterized members of this family have been crystallized in conformations with inaccessible active sites. We have trapped human TG2 in complex with an inhibitor that mimics inflammatory gluten peptide substrates and have solved, at 2-Å resolution, its x-ray crystal structure. The inhibitor stabilizes TG2 in an extended conformation that is dramatically different from earlier transglutaminase structures. The active site is exposed, revealing that catalysis takes place in a tunnel, bridged by two tryptophan residues that separate acyl-donor from acyl-acceptor and stabilize the tetrahedral reaction intermediates. Site-directed mutagenesis was used to investigate the acyl-acceptor side of the tunnel, yielding mutants with a marked increase in preference for hydrolysis over transamidation. By providing the ability to visualize this activated conformer, our results create a foundation for understanding the catalytic as well as the non-catalytic roles of TG2 in biology, and for dissecting the process by which the autoantibody response to TG2 is induced in celiac sprue patients.     

DNA cloning demonstrates high genetic heterogeneity in populations of the subaerial green alga Trentepohlia (Trentepohliales,Chlorophyta)

Mats of the green alga Trentepohlia, a genus widely distributed in the tropics as well as temperate regions, have always been perceived as homogeneous (i.e., formed by only one species). As such, their general nature and specific feature play a supportive role in the species delimitation. However, the presence of morphologically dissimilar thalli was observed under the light microscope and when cultivating a piece of a single mat. To address this, we performed DNA cloning of the rbcL gene on mat fragments of T. abietina, T. annulata, T. jolithus and T. umbrina sampled in Europe to reveal if they are composed of one or more species. We revealed that more Trentepohlia haplotypes may coexist in a single mat. In consideration of this, we conclude that the use of material isolated in unialgal culture will be almost mandatory for a taxonomic reassessment of this complicated genus. Another direct implication of this problem is that herbarium specimens consisting of field‐collected material should not be used for direct sequencing. We further hypothesize the reasons why multiple haplotypes of Trentepohlia occur more frequently in the tuft‐like mats. Possibly, fragments and/or cells of other microalgae, including other species of Trentepohlia, might be retained in such mats more easily than in the crustose trentepohlialean mats.     

The effects of 5‐azacytidine and cadmium on global 5‐methylcytosine content and secondary metabolites in the freshwater microalgae Chlamydomonas reinhardtii and Scenedesmus quadricauda

Epigenetic changes are important mechanisms in the regulation of chromatin structure and gene expression. Cytosine methylation is one of the major epigenetic modifications, mediated by DNA methyltransferases, which transfer methyl groups from S‐adenosyl‐L‐methionine (SAM) to the fifth carbon of cytosine. Various external environmental conditions can change the global hypo/hypermethylation pattern of DNA. These alterations may affect the organism's response to stress conditions. In this study, for the first time, we investigated the effects of 5‐azacytidine, a DNA methyltransferase inhibitor, and cadmium, a toxic metal and environmental pollutant, on the growth, biosynthesis of secondary metabolites (phenols, flavonoids, carotenoids), SAM, S‐adenosylhomocysteine, 5′‐methylthioadenosine and global 5‐methylcytosine (5‐mC) in the green microalgae Chlamydomonas reinhardtii and Scenedesmus quadricauda. The studied species showed major differences in 5‐mC content, secondary metabolite content, and antioxidant activity. Cadmium increased GSH (glutathione) content in C. reinhardtii by 60% whereas 5‐azacytidine did not affect GSH. The biosynthesis of GSH in S. quadricauda in response to the stressors was the opposite. Global 5‐mC content of C. reinhardtii was 1%–1.5%, and the content in S. quadricauda was 3.5%. Amount of some investigated methionine cycle metabolites (SAM, S‐adenosyl homocysteine [SAH], methionine) in S. quadricauda distinctly exceeded C. reinhardtii as well. However, chlorophylls a and b, carotenoids, total phenolic content, total flavonoid content and, antioxidant activity were significantly higher in C. reinhardtii than S. quadricauda. Therefore, in further studies it would be advisable to verify whether methylation of cytosine affects the expression of genes encoding certain secondary metabolites.     

The molecular mechanism of induction of unfolded protein response by Chlamydia

The unfolded protein response (UPR) contributes to chlamydial pathogenesis, as a source of lipids and ATP during replication, and for establishing the initial anti-apoptotic state of host cell that ensures successful inclusion development. The molecular mechanism(s) of UPR induction by Chlamydia is unknown. Chlamydia use type III secretion system (T3SS) effector proteins (e.g, the Translocated Actin-Recruiting Phosphoprotein (Tarp) to stimulate host cell's cytoskeletal reorganization that facilitates invasion and inclusion development. We investigated the hypothesis that T3SS effector-mediated assembly of myosin-II complex produces activated non-muscle myosin heavy chain II (NMMHC-II), which then binds the UPR master regulator (BiP) and/or transducers to induce UPR. Our results revealed the interaction of the chlamydial effector proteins (CT228 and Tarp) with components of the myosin II complex and UPR regulator and transducer during infection. These interactions caused the activation and binding of NMMHC-II to BiP and IRE1α leading to UPR induction. In addition, specific inhibitors of myosin light chain kinase, Tarp oligomerization and myosin ATPase significantly reduced UPR activation and Chlamydia replication. Thus, Chlamydia induce UPR through T3SS effector-mediated activation of NMMHC-II components of the myosin complex to facilitate infectivity. The finding provides greater insights into chlamydial pathogenesis with the potential to identify therapeutic targets and formulations.     

Temperature‐dependence of minimum resource requirements alters competitive hierarchies in phytoplankton

Resource competition theory is a conceptual framework that provides mechanistic insights into competition and community assembly of species with different resource requirements. However, there has been little exploration of how resource requirements depend on other environmental factors, including temperature. Changes in resource requirements as influenced by environmental temperature would imply that climate warming can alter the outcomes of competition and community assembly. We experimentally demonstrate that environmental temperature alters the minimum light and nitrogen requirements – as well as other growth parameters – of six widespread phytoplankton species from distinct taxonomic groups. We found that species require the most nitrogen at the highest temperatures while light requirements tend to be lowest at intermediate temperatures, although there are substantial interspecific differences in the exact shape of this relationship. We also experimentally parameterize two competition models, which we use to illustrate how temperature, through its effects on species’ traits, alters competitive hierarchies in multispecies assemblages, determining community dynamics. Developing a mechanistic understanding of how temperature influences the ability to compete for limiting resources is a critical step towards improving forecasts of community dynamics under climate warming.     

When did the ancestor of true bugs become stinky? Disentangling the phylogenomics of Hemiptera–Heteroptera

The phylogeny of true bugs (Hemiptera: Heteroptera), one of the most diverse insect groups in terms of morphology and ecology, has been the focus of attention for decades with respect to several deep nodes between the suborders of Hemiptera and the infraorders of Heteroptera. Here, we assembled a phylogenomic data set of 53 taxa and 3102 orthologous genes to investigate the phylogeny of Hemiptera–Heteroptera, and both concatenation and coalescent methods were used. A binode-control approach for data filtering was introduced to reduce the incongruence between different genes, which can improve the performance of phylogenetic reconstruction. Both hypotheses (Coleorrhyncha + Heteroptera) and (Coleorrhyncha + Auchenorrhyncha) received support from various analyses, in which the former is more consistent with the morphological evidence. Based on a divergence time estimation performed on genes with a strong phylogenetic signal, the origin of true bugs was dated to 290–268 Ma in the Permian, the time in Earth's history with the highest concentration of atmospheric oxygen. During this time interval, at least 1007 apomorphic amino acids were retained in the common ancestor of the extant true bugs. These molecular apomorphies are located in 553 orthologous genes, which suggests the common ancestor of the extant true bugs may have experienced large-scale evolution at the genome level.