Ensembl Genomes 2016: more genomes, more complexity

A rationally designed oral vaccine induces immunoglobulin A in the murine gut that directs the evolution of attenuated Salmonella variants

The ability of gut bacterial pathogens to escape immunity by antigenic variation-particularly via changes to surface-exposed antigens-is a major barrier to immune clearance¹. However, not all variants are equally fit in all environments^2,3. It should therefore be possible to exploit such immune escape mechanisms to direct an evolutionary trade-off. Here, we demonstrate this phenomenon using Salmonella enterica subspecies enterica serovar Typhimurium (S.Tm). A dominant surface antigen of S.Tm is its O-antigen: a long, repetitive glycan that can be rapidly varied by mutations in biosynthetic pathways or by phase variation^4,5. We quantified the selective advantage of O-antigen variants in the presence and absence of O-antigen-specific immunoglobulin A and identified a set of evolutionary trajectories allowing immune escape without an associated fitness cost in naive mice. Through the use of rationally designed oral vaccines, we induced immunoglobulin A responses blocking all of these trajectories. This selected for Salmonella mutants carrying deletions of the O-antigen polymerase gene wzyB. Due to their short O-antigen, these evolved mutants were more susceptible to environmental stressors (detergents or complement) and predation (bacteriophages) and were impaired in gut colonization and virulence in mice. Therefore, a rationally induced cocktail of intestinal antibodies can direct an evolutionary trade-off in S.Tm. This lays the foundations for the exploration of mucosal vaccines capable of setting evolutionary traps as a prophylactic strategy.

Inferring Adaptive Codon Preference to Understand Sources of Selection Shaping Codon Usage Bias

Alternative synonymous codons are often used at unequal frequencies. Classically, studies of such codon usage bias (CUB) attempted to separate the impact of neutral from selective forces by assuming that deviations from a predicted neutral equilibrium capture selection. However, GC-biased gene conversion (gBGC) can also cause deviation from a neutral null. Alternatively, selection has been inferred from CUB in highly expressed genes, but the accuracy of this approach has not been extensively tested, and gBGC can interfere with such extrapolations (e.g., if expression and gene conversion rates covary). It is therefore critical to examine deviations from a mutational null in a species with no gBGC. To achieve this goal, we implement such an analysis in the highly AT rich genome of Dictyostelium discoideum, where we find no evidence of gBGC. We infer neutral CUB under mutational equilibrium to quantify "adaptive codon preference," a nontautologous genome wide quantitative measure of the relative selection strength driving CUB. We observe signatures of purifying selection consistent with selection favoring adaptive codon preference. Preferred codons are not GC rich, underscoring the independence from gBGC. Expression-associated "preference" largely matches adaptive codon preference but does not wholly capture the influence of selection shaping patterns across all genes, suggesting selective constraints associated specifically with high expression. We observe patterns consistent with effects on mRNA translation and stability shaping adaptive codon preference. Thus, our approach to quantifying adaptive codon preference provides a framework for inferring the sources of selection that shape CUB across different contexts within the genome.

Longitudinal study of the scalp microbiome suggests coconut oil to enrich healthy scalp commensals

Dandruff is a recurrent chronic scalp disorder, affecting majority of the population worldwide. Recently a metagenomic study of the Indian scalp microbiome described an imperative role of bacterial commensals in providing essential vitamins and amino acids to the scalp. Coconut oil and its formulations are commonly applied on the scalp in several parts of the world to maintain scalp health. Thus, in this study we examined the effect of topical application of coconut oil on the scalp microbiome (bacterial and fungal) at the taxonomic and functional levels and their correlation with scalp physiological parameters. A 16-weeks-long time-course study was performed including 12-weeks of treatment and 4-weeks of relapse phase on a cohort of 140 (70 healthy and 70 dandruff) Indian women, resulting in ~ 900 metagenomic samples. After the treatment phase, an increase in the abundance of Cutibacterium acnes and Malassezia globosa in dandruff scalp was observed, which were negatively correlated to dandruff parameters. At the functional level, an enrichment of healthy scalp-related bacterial pathways, such as biotin metabolism and decrease in the fungal pathogenesis pathways was observed. The study provides novel insights on the effect of coconut oil in maintaining a healthy scalp and in modulating the scalp microbiome.

Translational adaptation to heat stress is mediated by RNA 5-methylcytosine in Caenorhabditis elegans

Methylation of carbon-5 of cytosines (m5 C) is a post-transcriptional nucleotide modification of RNA found in all kingdoms of life. While individual m5 C-methyltransferases have been studied, the impact of the global cytosine-5 methylome on development, homeostasis and stress remains unknown. Here, using Caenorhabditis elegans, we generated the first organism devoid of m5 C in RNA, demonstrating that this modification is non-essential. Using this genetic tool, we determine the localisation and enzymatic specificity of m5 C sites in the RNome in vivo. We find that NSUN-4 acts as a dual rRNA and tRNA methyltransferase in C. elegans mitochondria. In agreement with leucine and proline being the most frequently methylated tRNA isoacceptors, loss of m5 C impacts the decoding of some triplets of these two amino acids, leading to reduced translation efficiency. Upon heat stress, m5 C loss leads to ribosome stalling at UUG triplets, the only codon translated by an m5 C34-modified tRNA. This leads to reduced translation efficiency of UUG-rich transcripts and impaired fertility, suggesting a role of m5 C tRNA wobble methylation in the adaptation to higher temperatures.

Combination of Ginsenosides Rb2 and Rg3 Promotes Angiogenic Phenotype of Human Endothelial Cells via PI3K/Akt and MAPK/ERK Pathways

Shexiang Baoxin Pill (SBP) is an oral formulation of Chinese materia medica for the treatment of angina pectoris. It displays pleiotropic roles in protecting the cardiovascular system. However, the mode of action of SBP in promoting angiogenesis, and in particular the synergy between its constituents is currently not fully understood. The combination of ginsenosides Rb2 and Rg3 were studied in human umbilical vein endothelial cells (HUVECs) for their proangiogenic effects. To understand the mode of action of the combination in more mechanistic detail, RNA-Seq analysis was conducted, and differentially expressed genes (DEGs), pathway analysis and Weighted Gene Correlation Network Analysis (WGCNA) were applied to further identify important genes that a play pivotal role in the combination treatment. The effects of pathway-specific inhibitors were observed to provide further support for the hypothesized mode of action of the combination. Ginsenosides Rb2 and Rg3 synergistically promoted HUVEC proliferation and tube formation under defined culture conditions. Also, the combination of Rb2/Rg3 rescued cells from homocysteine-induced damage. mRNA expression of CXCL8, CYR61, FGF16 and FGFRL1 was significantly elevated by the Rb2/Rg3 treatment, and representative signaling pathways induced by these genes were found. The increase of protein levels of phosphorylated-Akt and ERK42/44 by the Rb2/Rg3 combination supports the notion that it promotes endothelial cell proliferation via the PI3K/Akt and MAPK/ERK signaling pathways. The present study provides the hypothesis that SBP, via ginsenosides Rb2 and Rg3, involves the CXCR1/2 CXCL8 (IL8)-mediated PI3K/Akt and MAPK/ERK signaling pathways in achieving its proangiogenic effects.

Gene Coexpression Network Reveals Insights into the Origin and Evolution of a Theanine-Associated Regulatory Module in Non-Camellia and Camellia Species

Theanine (thea) is one of the most important plant-derived characteristic secondary metabolites and a major healthcare product because of its beneficial biological activities, such as being an antianxiety agent, promoting memory, and lowering blood pressure. Thea mostly accumulates in Camellia plants and is especially rich in Camellia sinensis (tea plant). Although some functional genes (e.g., TS, GOGAT, and GS) attributed to thea accumulation have been separately well explored in tea plants, the evolution of a regulatory module (highly interacting gene group) related to thea metabolism remains to be elaborated. Herein, a thea-associated regulatory module (TARM) was mined by using a comprehensive analysis of a weighted gene coexpression network in Camellia and non-Camellia species. Comparative genomic analysis of 84 green plant species revealed that TARM originated from the ancestor of green plants (algae) and that TARM genes were recruited from different evolutionary nodes with the most gene duplication events at the early stage. Among the TARM genes, two core transcription factors of NAC080 and LBD38 were deduced, which may play a crucial role in regulating the biosynthesis of thea. Our findings provide the first insights into the origin and evolution of TARM and indicate a promising paradigm for identifying vital regulatory genes involved in thea metabolism.

A network-based comparative framework to study conservation and divergence of proteomes in plant phylogenies

Comparative functional genomics offers a powerful approach to study species evolution. To date, the majority of these studies have focused on the transcriptome in mammalian and yeast phylogenies. Here, we present a novel multi-species proteomic dataset and a computational pipeline to systematically compare the protein levels across multiple plant species. Globally we find that protein levels diverge according to phylogenetic distance but is more constrained than the mRNA level. Module-level comparative analysis of groups of proteins shows that proteins that are more highly expressed tend to be more conserved. To interpret the evolutionary patterns of conservation and divergence, we develop a novel network-based integrative analysis pipeline that combines publicly available transcriptomic datasets to define co-expression modules. Our analysis pipeline can be used to relate the changes in protein levels to different species-specific phenotypic traits. We present a case study with the rhizobia-legume symbiosis process that supports the role of autophagy in this symbiotic association.

Sequence and Structure-Based Analysis of Specificity Determinants in Eukaryotic Protein Kinases

Protein kinases lie at the heart of cell-signaling processes and are often mutated in disease. Kinase target recognition at the active site is in part determined by a few amino acids around the phosphoacceptor residue. However, relatively little is known about how most preferences are encoded in the kinase sequence or how these preferences evolved. Here, we used alignment-based approaches to predict 30 specificity-determining residues (SDRs) for 16 preferences. These were studied with structural models and were validated by activity assays of mutant kinases. Cancer mutation data revealed that kinase SDRs are mutated more frequently than catalytic residues. We have observed that, throughout evolution, kinase specificity has been strongly conserved across orthologs but can diverge after gene duplication, as illustrated by the G protein-coupled receptor kinase family. The identified SDRs can be used to predict kinase specificity from sequence and aid in the interpretation of evolutionary or disease-related genomic variants.

Whole genome sequencing of silver carp (Hypophthalmichthys molitrix) and bighead carp (Hypophthalmichthys nobilis) provide novel insights into their evolution and speciation

The edible silver carp (Hypophthalmichthys molitrix) and bighead carp (H. nobilis), which are two of the "Four Domesticated Fish" of China, are cultivated intensively worldwide. Here, we constructed 837- and 845-Mb draft genome assemblies for the silver carp and the bighead carp, respectively, including 24,571 and 24,229 annotated protein-coding genes. Genetic maps, anchoring 71.7% and 83.8% of all scaffolds, were obtained for the silver and bighead carp, respectively. Phylogenetic analysis showed that the bighead carp formed a clade with the silver carp, with an estimated divergence time of 3.6 million years ago; the time of divergence between the silver carp and zebrafish was 50.7 million years ago. An East Asian cyprinid genome-specific chromosome fusion took place ~9.2 million years after this clade diverged from the clade containing the common carp and Sinocyclocheilus. KEGG and GO analyses indicated that the expanded gene families in the silver and bighead carp were associated with diseases, the immune system and environmental adaptations. Genomic regions differentiating the silver and bighead carp populations were detected based on the whole-genome sequences of 42 individuals. Genes associated with the divergent regions were associated with reproductive system development and the development of primary female sexual characteristics. Thus, our results provided a novel systematic genomic analysis of the East Asian cyprinids, as well as the evolution and speciation of the silver carp and bighead carp.

Global Analysis of Gene Expression Profiles Provides Novel Insights into the Development and Evolution of the Large Crustacean Eriocheir sinensis

Chinese mitten crab (Eriocheir sinensis) is an important aquaculture species in Crustacea. Functional analysis, although essential, has been hindered due to the lack of sufficient genomic or transcriptomic resources. In this study, transcriptome sequencing was conducted on 59 samples representing diverse developmental stages (fertilized eggs, zoea, megalopa, three sub-stages of larvae, juvenile crabs, and adult crabs) and different tissues (eyestalk, hepatopancreas, and muscle from juvenile crabs, and eyestalk, hepatopancreas, muscle, heart, stomach, gill, thoracic ganglia, intestine, ovary, and testis from adult crabs) of E. sinensis. A comprehensive reference transcriptome was assembled, including 19,023 protein-coding genes. Hierarchical clustering based on 128 differentially expressed cuticle-related genes revealed two distinct expression patterns during the early larval developmental stages, demonstrating the distinct roles of these genes in “crab-like” cuticle formation during metamorphosis and cuticle calcification after molting. Phylogenetic analysis of 1406 one-to-one orthologous gene families identified from seven arthropod species and Caenorhabditis elegans strongly supported the hypothesis that Malacostraca and Branchiopoda do not form a monophyletic group. Furthermore, Branchiopoda is more phylogenetically closely related to Hexapoda, and the clade of Hexapoda and Branchiopoda and the clade of Malacostraca belong to the Pancrustacea. This study offers a high-quality transcriptome resource for E. sinensis and demonstrates the evolutionary relationships of major arthropod groups. The differentially expressed genes identified in this study facilitate further investigation of the cuticle-related gene expression networks which are likely associated with “crab-like” cuticle formation during metamorphosis and cuticle calcification after molting.

Using AnABlast for intergenic sORF prediction in the Caenorhabditis elegans genome

Motivation

Short bioactive peptides encoded by small open reading frames (sORFs) play important roles in eukaryotes. Bioinformatics prediction of ORFs is an early step in a genome sequence analysis, but sORFs encoding short peptides, often using non-AUG initiation codons, are not easily discriminated from false ORFs occurring by chance.

Results

AnABlast is a computational tool designed to highlight putative protein-coding regions in genomic DNA sequences. This protein-coding finder is independent of ORF length and reading frame shifts, thus making of AnABlast a potentially useful tool to predict sORFs. Using this algorithm, here, we report the identification of 82 putative new intergenic sORFs in the Caenorhabditis elegans genome. Sequence similarity, motif presence, expression data and RNA interference experiments support that the underlined sORFs likely encode functional peptides, encouraging the use of AnABlast as a new approach for the accurate prediction of intergenic sORFs in annotated eukaryotic genomes.

Availability and implementation

AnABlast is freely available at http://www.bioinfocabd.upo.es/ab/. The C.elegans genome browser with AnABlast results, annotated genes and all data used in this study is available at http://www.bioinfocabd.upo.es/celegans.

Supplementary information

Supplementary data are available at Bioinformatics online.

Algorithms for ribosome traffic engineering and their potential in improving host cells' titer and growth rate

mRNA translation is a fundamental cellular process consuming most of the intracellular energy; thus, it is under extensive evolutionary selection for optimization, and its efficiency can affect the host's growth rate. We describe a generic approach for improving the growth rate (fitness) of any organism by introducing synonymous mutations based on comprehensive computational models. The algorithms introduce silent mutations that may improve the allocation of ribosomes in the cells via the decreasing of their traffic jams during translation respectively. As a result, resources availability in the cell changes leading to improved growth-rate. We demonstrate experimentally the implementation of the method on Saccharomyces cerevisiae: we show that by introducing a few mutations in two computationally selected genes the mutant's titer increased. Our approach can be employed for improving the growth rate of any organism providing the existence of data for inferring models, and with the relevant genomic engineering tools; thus, it is expected to be extremely useful in biotechnology, medicine, and agriculture.

Calculating the most likely intron splicing orders in S. pombe, fruit fly, Arabidopsis thaliana, and humans

Background

Introns have been shown to be spliced in a defined order, and this order influences both alternative splicing regulation and splicing fidelity, but previous studies have only considered neighbouring introns. The detailed intron splicing order remains unknown.

Results

In this work, a method was developed that can calculate the intron splicing orders of all introns in each transcript. A simulation study showed that this method can accurately calculate intron splicing orders. I further applied this method to real S. pombe, fruit fly, Arabidopsis thaliana, and human sequencing datasets and found that intron splicing orders change from gene to gene and that humans contain more not in-order spliced transcripts than S. pombe, fruit fly and Arabidopsis thaliana. In addition, I reconfirmed that the first introns in humans are spliced slower than those in S. pombe, fruit fly, and Arabidopsis thaliana genome-widely. Both the calculated most likely orders and the method developed here are available on the web.

Conclusions

A novel computational method was developed to calculate the intron splicing orders and applied the method to real sequencing datasets. I obtained intron splicing orders for hundreds or thousands of genes in four organisms. I found humans contain more number of not in-order spliced transcripts.

Transcriptomic and Proteomic Insights into Amborella trichopoda Male Gametophyte Functions

Flowering plants (angiosperms) are characterized by pollen tubes (PTs; male gametophytes) carrying two immobile sperm cells that grow over long distances through the carpel toward the ovules, where double fertilization is executed. It is not understood how these reproductive structures evolved, which genes occur de novo in male gametophytes of angiosperms, and to which extent PT functions are conserved among angiosperms. To contribute to a deeper understanding of the evolution of gametophyte functions, we generated RNA sequencing data from seven reproductive and two vegetative control tissues of the basal angiosperm Amborella trichopoda and complemented these with proteomic data of pollen grains (PGs) and PTs. The eudicot model plant Arabidopsis (Arabidopsis thaliana) served as a reference organism for data analysis, as more than 200 genes have been associated with male gametophyte functions in this species. We describe methods to collect bicellular A. trichopoda PGs, to induce their germination in vitro, and to monitor PT growth and germ cell division. Transcriptomic and proteomic analyses indicate that A. trichopoda PGs are prepared for germination requiring lipids, energy, but likely also reactive oxygen species, while PTs are especially characterized by catabolic/biosynthetic and transport processes including cell wall biosynthesis and gene regulation. Notably, a number of pollen-specific genes were lacking in Arabidopsis, and the number of genes involved in pollen signaling is significantly reduced in A. trichopoda In conclusion, we provide insight into male gametophyte functions of the most basal angiosperm and establish a valuable resource for future studies on the evolution of flowering plants.

Comparison of DNA Methylation Profiles of Hemostatic Genes between Liver Tissue and Peripheral Blood within Individuals

DNA methylation has become increasingly recognized in the etiology of complex diseases, including thrombotic disorders. Blood is often collected in epidemiological studies for genotyping and has recently also been used to examine DNA methylation in epigenome-wide association studies. DNA methylation patterns are often tissue-specific, thus, peripheral blood may not accurately reflect the methylation pattern in the tissue of relevance. Here, we collected paired liver and blood samples concurrently from 27 individuals undergoing liver surgery. We performed targeted bisulfite sequencing for a set of 35 hemostatic genes primarily expressed in liver to analyze DNA methylation levels of >10,000 cytosine-phosphate-guanine (CpG) dinucleotides. We evaluated whether DNA methylation in blood could serve as a proxy for DNA methylation in liver at individual CpGs. Approximately 30% of CpGs were nonvariable and were predominantly hypo- (<25%) or hypermethylated (>70%) in both tissues. While blood can serve as a proxy for liver at these CpGs, the low variability renders these unlikely to explain phenotypic differences. We therefore focused on CpG sites with variable methylation levels in liver. The level of blood-liver tissue correlation varied widely across these variable CpGs; moderate correlations (0.5 ≤ r < 0.75) were detected for 6% and strong correlations (r ≥ 0.75) for a further 4%. Our findings indicate that it is essential to study the concordance of DNA methylation between blood and liver at individual CpGs. This paired blood-liver dataset is intended as a resource to aid interpretation of blood-based DNA methylation results.

No Evidence That Homologs of Key Circadian Clock Genes Direct Circadian Programs of Development or mRNA Abundance in Verticillium dahliae

Many organisms harbor circadian clocks that promote their adaptation to the rhythmic environment. While a broad knowledge of the molecular mechanism of circadian clocks has been gained through the fungal model Neurospora crassa, little is known about circadian clocks in other fungi. N. crassa belongs to the same class as many important plant pathogens including the vascular wilt fungus Verticillium dahliae. We identified homologs of N. crassa clock proteins in V. dahliae, which showed high conservation in key protein domains. However, no evidence for an endogenous, free-running and entrainable rhythm was observed in the daily formation of conidia and microsclerotia. In N. crassa the frequency (frq) gene encodes a central clock protein expressed rhythmically and in response to light. In contrast, expression of Vdfrq is not light-regulated. Temporal gene expression profiling over 48 h in constant darkness and temperature revealed no circadian expression of key clock genes. Furthermore, RNA-seq over a 24 h time-course revealed no robust oscillations of clock-associated transcripts in constant darkness. Comparison of gene expression between wild-type V. dahliae and a ΔVdfrq mutant showed that genes involved in metabolism, transport and redox processes are mis-regulated in the absence of Vdfrq. In addition, VdΔfrq mutants display growth defects and reduced pathogenicity in a strain dependent manner. Our data indicate that if a circadian clock exists in Verticillium, it is based on alternative mechanisms such as post-transcriptional interactions of VdFRQ and the WC proteins or the components of a FRQ-less oscillator. Alternatively, it could be that whilst the original functions of the clock proteins have been maintained, in this species the interactions that generate robust rhythmicity have been lost or are only triggered when specific environmental conditions are met. The presence of conserved clock genes in genomes should not be taken as definitive evidence of circadian function.

Genome-wide identification and expression analysis of the calmodulin-binding transcription activator (CAMTA) gene family in wheat (Triticum aestivum L.)

Background

Plant calmodulin-binding transcription activator (CAMTA) proteins play important roles in hormone signal transduction, developmental regulation, and environmental stress tolerance. However, in wheat, the CAMTA gene family has not been systematically characterized.

Results

In this work, 15 wheat CAMTA genes were identified using a genome-wide search method. Their chromosome location, physicochemical properties, subcellular localization, gene structure, protein domain, and promoter cis-elements were systematically analyzed. Phylogenetic analysis classified the TaCAMTA genes into three groups (groups A, B, and C), numbered 7, 6, and 2, respectively. The results showed that most TaCAMTA genes contained stress-related cis-elements. Finally, to obtain tissue-specific and stress-responsive candidates, the expression profiles of the TaCAMTAs in various tissues and under biotic and abiotic stresses were investigated. Tissue-specific expression analysis showed that all of the 15 TaCAMTA genes were expressed in multiple tissues with different expression levels, as well as under abiotic stress, the expressions of each TaCAMTA gene could respond to at least one abiotic stress. It also found that 584 genes in wheat genome were predicted to be potential target genes by CAMTA, demonstrating that CAMTA can be widely involved in plant development and growth, as well as coping with stresses.

Conclusions

This work systematically identified the CAMTA gene family in wheat at the whole-genome-wide level, providing important candidates for further functional analysis in developmental regulation and the stress response in wheat.

miR155 regulation of behavior, neuropathology, and cortical transcriptomics in Alzheimer's disease

MicroRNAs are recognized as important regulators of many facets of physiological brain function while also being implicated in the pathogenesis of several neurological disorders. Dysregulation of miR155 is widely reported across a variety of neurodegenerative conditions, including Alzheimer's disease (AD), Parkinson's disease, amyotrophic lateral sclerosis, and traumatic brain injury. In previous work, we observed that experimentally validated miR155 gene targets were consistently enriched among genes identified as differentially expressed across multiple brain tissue and disease contexts. In particular, we found that human herpesvirus-6A (HHV-6A) suppressed miR155, recapitulating reports of miR155 inhibition by HHV-6A in infected T-cells, thyrocytes, and natural killer cells. In earlier studies, we also reported the effects of constitutive deletion of miR155 on accelerating the accumulation of Aβ deposits in 4-month-old APP/PSEN1 mice. Herein, we complete the cumulative characterization of transcriptomic, electrophysiological, neuropathological, and learning behavior profiles from 4-, 8- and 10-month-old WT and APP/PSEN1 mice in the absence or presence of miR155. We also integrated human post-mortem brain RNA-sequences from four independent AD consortium studies, together comprising 928 samples collected from six brain regions. We report that gene expression perturbations associated with miR155 deletion in mouse cortex are in aggregate observed to be concordant with AD-associated changes across these independent human late-onset AD (LOAD) data sets, supporting the relevance of our findings to human disease. LOAD has recently been formulated as the clinicopathological manifestation of a multiplex of genetic underpinnings and pathophysiological mechanisms. Our accumulated data are consistent with such a formulation, indicating that miR155 may be uniquely positioned at the intersection of at least four components of this LOAD "multiplex": (1) innate immune response pathways; (2) viral response gene networks; (3) synaptic pathology; and (4) proamyloidogenic pathways involving the amyloid β peptide (Aβ).

Diversity in the intrinsic apoptosis pathway of nematodes

Early studies of the free-living nematode C. elegans informed us how BCL-2-regulated apoptosis in humans is regulated. However, subsequent studies showed C. elegans apoptosis has several unique features compared with human apoptosis. To date, there has been no detailed analysis of apoptosis regulators in nematodes other than C. elegans. Here, we discovered BCL-2 orthologues in 89 free-living and parasitic nematode taxa representing four evolutionary clades (I, III, IV and V). Unlike in C. elegans, 15 species possess multiple (two to five) BCL-2-like proteins, and some do not have any recognisable BCL-2 sequences. Functional studies provided no evidence that BAX/BAK proteins have evolved in nematodes, and structural studies of a BCL-2 protein from the basal clade I revealed it lacks a functionally important feature of the C. elegans orthologue. Clade I CED-4/APAF-1 proteins also possess WD40-repeat sequences associated with apoptosome assembly, not present in C. elegans, or other nematode taxa studied.

Wheat Disease Resistance Genes and Their Diversification Through Integrated Domain Fusions

Plants are in a constant evolutionary arms race with their pathogens. At the molecular level, the plant nucleotide-binding leucine-rich repeat receptors (NLRs) family has coevolved with rapidly evolving pathogen effectors. While many NLRs utilize variable leucine-rich repeats (LRRs) to detect effectors, some have gained integrated domains (IDs) that may be involved in receptor activation or downstream signaling. The major objectives of this project were to identify NLR genes in wheat (Triticum aestivum L.) and assess IDs associated with immune signaling (e.g., kinase and transcription factor domains). We identified 2,151 NLR-like genes in wheat, of which 1,298 formed 547 gene clusters. Among the non-toll/interleukin-1 receptor NLR (non-TNL)-like genes, 1,552 encode LRRs, 802 are coiled-coil (CC) domain-encoding (CC-NBS-LRR or CNL) genes, and three encode resistance to powdery mildew 8 (RPW8) domains (RPW8-NBS-LRR or RNL). The expansion of the NLR gene family in wheat is attributable to its origin by recent polyploidy events. Gene clusters were likely formed by tandem duplications, and wheat NLR phylogenetic relationships were similar to those in barley and Aegilops. We also identified wheat NLR-ID fusion proteins as candidates for NLR functional diversification, often as kinase and transcription factor domains. Comparative analyses of the IDs revealed evolutionary conservation of more than 80% amino acid sequence similarity. Homology assessment indicates that these domains originated as functional non-NLR-encoding genes that were incorporated into NLR-encoding genes through duplication events. We also found that many of the NLR-ID genes encode alternative transcripts that include or exclude IDs, a phenomenon that seems to be conserved among species. To verify this, we have analyzed the alternative transcripts that include or exclude an ID of an NLR-ID from another monocotyledon species, rice (Oryza sativa). This indicates that plants employ alternative splicing to regulate IDs, possibly using them as baits, decoys, and functional signaling components. Genomic and expression data support the hypothesis that wheat uses alternative splicing to include and exclude IDs from NLR proteins.

The regulatory landscape of early maize inflorescence development

BACKGROUND

The functional genome of agronomically important plant species remains largely unexplored, yet presents a virtually untapped resource for targeted crop improvement. Functional elements of regulatory DNA revealed through profiles of chromatin accessibility can be harnessed for fine-tuning gene expression to optimal phenotypes in specific environments.

RESULT

Here, we investigate the non-coding regulatory space in the maize (Zea mays) genome during early reproductive development of pollen- and grain-bearing inflorescences. Using an assay for differential sensitivity of chromatin to micrococcal nuclease (MNase) digestion, we profile accessible chromatin and nucleosome occupancy in these largely undifferentiated tissues and classify at least 1.6% of the genome as accessible, with the majority of MNase hypersensitive sites marking proximal promoters, but also 3' ends of maize genes. This approach maps regulatory elements to footprint-level resolution. Integration of complementary transcriptome profiles and transcription factor occupancy data are used to annotate regulatory factors, such as combinatorial transcription factor binding motifs and long non-coding RNAs, that potentially contribute to organogenesis, including tissue-specific regulation between male and female inflorescence structures. Finally, genome-wide association studies for inflorescence architecture traits based solely on functional regions delineated by MNase hypersensitivity reveals new SNP-trait associations in known regulators of inflorescence development as well as new candidates.

CONCLUSIONS

These analyses provide a comprehensive look into the cis-regulatory landscape during inflorescence differentiation in a major cereal crop, which ultimately shapes architecture and influences yield potential.

Short- and Long-Term Transcriptomic Responses of Escherichia coli to Biocides: a Systems Analysis

The mechanisms of the bacterial response to biocides are poorly understood, despite their broad application. To identify the genetic basis and pathways implicated in the biocide stress response, we exposed Escherichia coli populations to 10 ubiquitous biocides. By comparing the transcriptional responses between a short-term exposure (30 min) and a long-term exposure (8 to 12 h) to biocide stress, we established the common gene and pathway clusters that are implicated in general and biocide-specific stress responses. Our analysis revealed a temporal choreography, starting from the upregulation of chaperones to the subsequent repression of motility and chemotaxis pathways and the induction of an anaerobic pool of enzymes and biofilm regulators. A systematic analysis of the transcriptional data identified a zur-regulated gene cluster to be highly active in the stress response against sodium hypochlorite and peracetic acid, presenting a link between the biocide stress response and zinc homeostasis. Susceptibility assays with knockout mutants further validated our findings and provide clear targets for downstream investigation of the implicated mechanisms of action.IMPORTANCE Antiseptics and disinfectant products are of great importance to control and eliminate pathogens, especially in settings such as hospitals and the food industry. Such products are widely distributed and frequently poorly regulated. Occasional outbreaks have been associated with microbes resistant to such compounds, and researchers have indicated potential cross-resistance with antibiotics. Despite that, there are many gaps in knowledge about the bacterial stress response and the mechanisms of microbial resistance to antiseptics and disinfectants. We investigated the stress response of the bacterium Escherichia coli to 10 common disinfectant and antiseptic chemicals to shed light on the potential mechanisms of tolerance to such compounds.

Predicting the optimal growth temperatures of prokaryotes using only genome derived features

MOTIVATION

Optimal growth temperature is a fundamental characteristic of all living organisms. Knowledge of this temperature is central to the study of a prokaryote, the thermal stability and temperature dependent activity of its genes, and the bioprospecting of its genome for thermally adapted proteins. While high throughput sequencing methods have dramatically increased the availability of genomic information, the growth temperatures of the source organisms are often unknown. This limits the study and technological application of these species and their genomes. Here, we present a novel method for the prediction of growth temperatures of prokaryotes using only genomic sequences.

RESULTS

By applying the reverse ecology principle that an organism's genome includes identifiable adaptations to its native environment, we can predict a species' optimal growth temperature with an accuracy of 5.17°C root-mean-square error and a coefficient of determination of 0.835. The accuracy can be further improved for specific taxonomic clades or by excluding psychrophiles. This method provides a valuable tool for the rapid calculation of organism growth temperature when only the genome sequence is known.

AVAILABILITY AND IMPLEMENTATION

Source code, genomes analyzed and features calculated are available at: https://github.com/DavidBSauer/OGT_prediction.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Dicer-like 5 deficiency confers temperature-sensitive male sterility in maize

Small RNAs play important roles during plant development by regulating transcript levels of target mRNAs, maintaining genome integrity, and reinforcing DNA methylation. Dicer-like 5 (Dcl5) is proposed to be responsible for precise slicing in many monocots to generate diverse 24-nt phased, secondary small interfering RNAs (phasiRNAs), which are exceptionally abundant in meiotic anthers of diverse flowering plants. The importance and functions of these phasiRNAs remain unclear. Here, we characterized several mutants of dcl5, including alleles generated by the clustered regularly interspaced short palindromic repeats (CRISPR)–Cas9 system and a transposon-disrupted allele. We report that dcl5 mutants have few or no 24-nt phasiRNAs, develop short anthers with defective tapetal cells, and exhibit temperature-sensitive male fertility. We propose that DCL5 and 24-nt phasiRNAs are critical for fertility under growth regimes for optimal yield.

Small RNAs act to regulate gene or transposon activity during plant development. Here, the authors show that maize Dicer-like 5 is required for 24-nt phased, secondary small interfering RNA production in anthers and that dicer-like 5 mutants show abnormal tapetal development and temperature-sensitive sterility.

PRMdb: A Repository of Predicted RNA Modifications in Plants

Evidence is mounting that RNA modifications play essential roles in posttranscriptional regulation of gene expression. So far, over 150 RNA modifications catalyzed by distinct enzymes have been documented. In plants, genome-wide identification of RNA modifications is largely limited to the model species Arabidopsis thaliana, while lacking in diverse non-model plants. Here, we present PRMdb, a plant RNA modification database, based on the analysis of thousands of RNA-seq, degradome-seq and small RNA-seq data from a wide range of plant species using the well-documented tool HAMR (high-throughput analysis of modified ribonucleotide). PRMdb provides a user-friendly interface that enables easy browsing and searching of the tRNA and mRNA modification data. We show that PRMdb collects high-confidence RNA modifications including novel RNA modification sites that can be validated by genomic PCR and reverse transcription PCR. In summary, PRMdb provides a valuable web resource for deciphering the epitranscriptomes in diverse plant species and will facilitate functional studies of RNA modifications in plants. RPMdb is available via http://www.biosequencing.cn/PRMdb/.

RNAcentral: a hub of information for non-coding RNA sequences

RNAcentral is a comprehensive database of non-coding RNA (ncRNA) sequences, collating information on ncRNA sequences of all types from a broad range of organisms. We have recently added a new genome mapping pipeline that identifies genomic locations for ncRNA sequences in 296 species. We have also added several new types of functional annotations, such as tRNA secondary structures, Gene Ontology annotations, and miRNA-target interactions. A new quality control mechanism based on Rfam family assignments identifies potential contamination, incomplete sequences, and more. The RNAcentral database has become a vital component of many workflows in the RNA community, serving as both the primary source of sequence data for academic and commercial groups, as well as a source of stable accessions for the annotation of genomic and functional features. These examples are facilitated by an improved RNAcentral web interface, which features an updated genome browser, a new sequence feature viewer, and improved text search functionality. RNAcentral is freely available at https://rnacentral.org.

Contribution of bacterial biodiversity on the operational performance of a styrene biotrickling filter

Long-term operational stability of biotrickling filters (BTFs) degrading volatile organic compounds (VOCs) is dependent on both physicochemical as well as biological properties. Effects of increasingly stressful levels of air pollutants on the microbial structure of biofilms within BTFs are not well understood, especially for VOCs such as styrene. To investigate the relationship between biofilm biodiversity and operational stability, the temporal dynamics of a biofilm from a biotrickling filter subjected to stepwise increasing levels of air polluted with styrene was investigated using 16S rDNA pyrosequencing and PCR-denaturing gradient gel electrophoresis (PCR-DGGE). As styrene contaminant loads were increased, microbial community composition was distinctly altered and diversity was initially reduced in early stages but gradually stabilized and increased diversity in later stages, suggesting a recovery and acclimatization period within the microbial community during incremental exposure of the pollutant. Although temporary reductions in known styrene-degrading bacterial genera (Pseudomonas and Rhodococcus) occurred under increased styrene loads, stable BTF performance was maintained due to functional redundancy. New candidate genera for styrene degradation (Azoarcus, Dokdonella) were identified in conditions of high styrene loads, and may have supported the observed stable BTF performance throughout the experiment. Styrene inlet load was found to be important modulator of community composition and may have been partly responsible for the observed temporary reductions of Pseudomonas. Notable differences between dominant genera detected via pyrosequencing compared to species detected by PCR-DGGE suggests that simultaneous implementation of both techniques is valuable for fully characterizing dynamic microbial communities.

MicroScope-an integrated resource for community expertise of gene functions and comparative analysis of microbial genomic and metabolic data

The overwhelming list of new bacterial genomes becoming available on a daily basis makes accurate genome annotation an essential step that ultimately determines the relevance of thousands of genomes stored in public databanks. The MicroScope platform (http://www.genoscope.cns.fr/agc/microscope) is an integrative resource that supports systematic and efficient revision of microbial genome annotation, data management and comparative analysis. Starting from the results of our syntactic, functional and relational annotation pipelines, MicroScope provides an integrated environment for the expert annotation and comparative analysis of prokaryotic genomes. It combines tools and graphical interfaces to analyze genomes and to perform the manual curation of gene function in a comparative genomics and metabolic context. In this article, we describe the free-of-charge MicroScope services for the annotation and analysis of microbial (meta)genomes, transcriptomic and re-sequencing data. Then, the functionalities of the platform are presented in a way providing practical guidance and help to the nonspecialists in bioinformatics. Newly integrated analysis tools (i.e. prediction of virulence and resistance genes in bacterial genomes) and original method recently developed (the pan-genome graph representation) are also described. Integrated environments such as MicroScope clearly contribute, through the user community, to help maintaining accurate resources.

Visualizing and comparing circular genomes using the CGView family of tools

Graphical genome maps are widely used to assess genome features and sequence characteristics. The CGView (Circular Genome Viewer) software family is a popular collection of tools for generating genome maps for bacteria, organelles and viruses. In this review, we describe the capabilities of the original CGView program along with those of subsequent companion applications, including the CGView Server and the CGView Comparison Tool. We also discuss GView, a graphical user interface-enabled rewrite of CGView, and the GView Server, which offers several integrated analyses for identifying shared or unique genome regions relative to a collection of comparison genomes. We conclude with some remarks about our current development efforts related to CGView aimed at adding new functionality while increasing ease of use.

PEATmoss (Physcomitrella Expression Atlas Tool): a unified gene expression atlas for the model plant Physcomitrella patens

Physcomitrella patens is a bryophyte model plant that is often used to study plant evolution and development. Its resources are of great importance for comparative genomics and evo-devo approaches. However, expression data from Physcomitrella patens were so far generated using different gene annotation versions and three different platforms: CombiMatrix and NimbleGen expression microarrays and RNA sequencing. The currently available P. patens expression data are distributed across three tools with different visualization methods to access the data. Here, we introduce an interactive expression atlas, Physcomitrella Expression Atlas Tool (PEATmoss), that unifies publicly available expression data for P. patens and provides multiple visualization methods to query the data in a single web-based tool. Moreover, PEATmoss includes 35 expression experiments not previously available in any other expression atlas. To facilitate gene expression queries across different gene annotation versions, and to access P. patens annotations and related resources, a lookup database and web tool linked to PEATmoss was implemented. PEATmoss can be accessed at https://peatmoss.online.uni-marburg.de.

Metabolic Labeling of RNAs Uncovers Hidden Features and Dynamics of the Arabidopsis Transcriptome

Transcriptome analysis by RNA sequencing (RNA-seq) has become an indispensable research tool in modern plant biology. Virtually all RNA-seq studies provide a snapshot of the steady state transcriptome, which contains valuable information about RNA populations at a given time but lacks information about the dynamics of RNA synthesis and degradation. Only a few specialized sequencing techniques, such as global run-on sequencing, have been used to provide information about RNA synthesis rates in plants. Here, we demonstrate that RNA labeling with the modified, nontoxic uridine analog 5-ethynyl uridine (5-EU) in Arabidopsis (Arabidopsis thaliana) seedlings provides insight into plant transcriptome dynamics. Pulse labeling with 5-EU revealed nascent and unstable RNAs, RNA processing intermediates generated by splicing, and chloroplast RNAs. Pulse-chase experiments with 5-EU allowed us to determine RNA stabilities without the need for chemical transcription inhibitors such as actinomycin and cordycepin. Inhibitor-free, genome-wide analysis of polyadenylated RNA stability via 5-EU pulse-chase experiments revealed RNAs with shorter half-lives than those reported after chemical inhibition of transcription. In summary, our results indicate that the Arabidopsis nascent transcriptome contains unstable RNAs and RNA processing intermediates and suggest that polyadenylated RNAs have low stability in plants. Our technique lays the foundation for easy, affordable, nascent transcriptome analysis and inhibitor-free analysis of RNA stability in plants.

Architecture of Parallel Adaptation in Ten Lacustrine Threespine Stickleback Populations from the White Sea Area

Adaptation of threespine stickleback to freshwater involves parallel recruitment of freshwater alleles in clusters of closely linked sites, or divergence islands (DIs). However, it remains unclear to what extent the DIs and the alleles that constitute them coincide between populations that underwent adaptation to freshwater independently. We examine threespine sticklebacks from ten freshwater lakes that emerged 500–1500 years ago in the White Sea basin, with the emphasis on repeatability of genomic patterns of adaptation among the lake populations and the role of local recombination rate in the distribution and structure of DIs. The 65 detected DIs are clustered in the genome, forming 12 aggregations, and this clustering cannot be explained by the variation of the recombination rate. Only 21 of the DIs are present in all the freshwater populations, likely being indispensable for successful colonization of freshwater environment by the ancestral marine population. Within most DIs, the same set of single nucleotide polymorphisms (SNPs) distinguish marine and freshwater haplotypes in all the lake populations; however, in some DIs, freshwater alleles differ between populations, suggesting that they could have been established by recruitment of different haplotypes in different populations.

Multilayered Tuning of Dosage Compensation and Z-Chromosome Masculinization in the Wood White (Leptidea sinapis) Butterfly

In species with genetic sex determination, dosage compensation can evolve to equal expression levels of sex-linked and autosomal genes. Current knowledge about dosage compensation has mainly been derived from male-heterogametic (XX/XY) model organisms, whereas less is understood about the process in female-heterogametic systems (ZZ/ZW). In moths and butterflies, downregulation of Z-linked expression in males (ZZ) to match the expression level in females (ZW) is often observed. However, little is known about the underlying regulatory mechanisms, or if dosage compensation patterns vary across ontogenetic stages. In this study, we assessed dynamics of Z-linked and autosomal expression levels across developmental stages in the wood white (Leptidea sinapis). We found that although expression of Z-linked genes in general was reduced compared with autosomal genes, dosage compensation was actually complete for some categories of genes, in particular sex-biased genes, but equalization in females was constrained to a narrower gene set. We also observed a noticeable convergence in Z-linked expression between males and females after correcting for sex-biased genes. Sex-biased expression increased successively across developmental stages, and male-biased genes were enriched on the Z-chromosome. Finally, all five core genes associated with the ribonucleoprotein dosage compensation complex male-specific lethal were detected in adult females, in correspondence with a reduction in the expression difference between autosomes and the single Z-chromosome. We show that tuning of gene dosage is multilayered in Lepidoptera and argue that expression balance across chromosomal classes may predominantly be driven by enrichment of male-biased genes on the Z-chromosome and cooption of available dosage regulators.

Transcriptome differences in adipose stromal cells derived from pre- and postmenopausal women

BACKGROUND

As the population ages, an increasing number of postmenopausal women are donors of adipose stromal cells (ASCs) and may benefit from autologous ASC-related treatments. However, the effect of menopausal status on ASCs has not been investigated.

METHODS

RNA sequencing data were downloaded, and differentially expressed genes (DEGs) were identified. Hierarchical clustering, Gene Ontology, and pathway analyses were applied to the DEGs. Two gene coexpression network analysis approaches were applied to the DEGs to provide a holistic view and preserve gene interactions. Hub genes of the gene coexpression network were identified, and their expression profiles were examined with clinical samples. ASCs from pre- and postmenopausal women were co-cultured with monocytes and T cells to determine their immunoregulatory role.

RESULTS

In total, 2299 DEGs were identified and presented distinct expression profiles between pre- and postmenopausal women. Gene Ontology and pathway analyses revealed some fertility-, sex hormone-, immune-, aging-, and angiogenesis-related terms and pathways. Gene coexpression networks were constructed, and the top hub genes, including TIE1, ANGPT2, RNASE1, PLVAP, CA2, and MPZL2, were consistent between the two approaches. Expression profiles of hub genes from the RNA sequencing data and clinical samples were consistent. ASCs from postmenopausal women elicit M1 polarization, while their counterparts facilitate CD3/4+ T cell proliferation.

CONCLUSIONS

The present study reveals the transcriptome differences in ASCs derived from pre- and postmenopausal women and provides holistic views by preserving gene interactions via gene coexpression network analysis. The top hub genes identified by this study could serve as potential targets to enhance the therapeutic potential of ASCs.

MetaOmGraph: a workbench for interactive exploratory data analysis of large expression datasets

The diverse and growing omics data in public domains provide researchers with tremendous opportunity to extract hidden, yet undiscovered, knowledge. However, the vast majority of archived data remain unused. Here, we present MetaOmGraph (MOG), a free, open-source, standalone software for exploratory analysis of massive datasets. Researchers, without coding, can interactively visualize and evaluate data in the context of its metadata, honing-in on groups of samples or genes based on attributes such as expression values, statistical associations, metadata terms and ontology annotations. Interaction with data is easy via interactive visualizations such as line charts, box plots, scatter plots, histograms and volcano plots. Statistical analyses include co-expression analysis, differential expression analysis and differential correlation analysis, with significance tests. Researchers can send data subsets to R for additional analyses. Multithreading and indexing enable efficient big data analysis. A researcher can create new MOG projects from any numerical data; or explore an existing MOG project. MOG projects, with history of explorations, can be saved and shared. We illustrate MOG by case studies of large curated datasets from human cancer RNA-Seq, where we identify novel putative biomarker genes in different tumors, and microarray and metabolomics data from Arabidopsis thaliana. MOG executable and code: http://metnetweb.gdcb.iastate.edu/ and https://github.com/urmi-21/MetaOmGraph/.

RNAdetect: efficient computational detection of novel non-coding RNAs

MOTIVATION

Non-coding RNAs (ncRNAs) are known to play crucial roles in various biological processes, and there is a pressing need for accurate computational detection methods that could be used to efficiently scan genomes to detect novel ncRNAs. However, unlike coding genes, ncRNAs often lack distinctive sequence features that could be used for recognizing them. Although many ncRNAs are known to have a well conserved secondary structure, which provides useful cues for computational prediction, it has been also shown that a structure-based approach alone may not be sufficient for detecting ncRNAs in a single sequence. Currently, the most effective ncRNA detection methods combine structure-based techniques with a comparative genome analysis approach to improve the prediction performance.

RESULTS

In this paper, we propose RNAdetect, a computational method incorporating novel features for accurate detection of ncRNAs in combination with comparative genome analysis. Given a sequence alignment, RNAdetect can accurately detect the presence of functional ncRNAs by incorporating novel predictive features based on the concept of generalized ensemble defect (GED), which assesses the degree of structure conservation across multiple related sequences and the conformation of the individual folding structures to a common consensus structure. Furthermore, n-gram models (NGMs) are used to extract features that can effectively capture sequence homology to known ncRNA families. Utilization of NGMs can enhance the detection of ncRNAs that have sparse folding structures with many unpaired bases. Extensive performance evaluation based on the Rfam database and bacterial genomes demonstrate that RNAdetect can accurately and reliably detect novel ncRNAs, outperforming the current state-of-the-art methods.

AVAILABILITY AND IMPLEMENTATION

The source code for RNAdetect and the benchmark data used in this paper can be downloaded at https://github.com/bjyoontamu/RNAdetect.

Organizing genome engineering for the gigabase scale

Genome-scale engineering holds great potential to impact science, industry, medicine, and society, and recent improvements in DNA synthesis have enabled the manipulation of megabase genomes. However, coordinating and integrating the workflows and large teams necessary for gigabase genome engineering remains a considerable challenge. We examine this issue and recommend a path forward by: 1) adopting and extending existing representations for designs, assembly plans, samples, data, and workflows; 2) developing new technologies for data curation and quality control; 3) conducting fundamental research on genome-scale modeling and design; and 4) developing new legal and contractual infrastructure to facilitate collaboration.

Diversity, structure, and synteny of the cutinase gene of Colletotrichum species

Colletotrichum species complexes are among the top 10 economically important fungal plant pathogens worldwide because they can infect climacteric and nonclimacteric fruit at the pre and/or postharvest stages. C. truncatum is the major pathogen responsible for anthracnose of green and red bell pepper fruit worldwide. C. brevisporum was recently reported to be a minor pathogen of red bell pepper fruit in Trinidad, but has recently been reported as pathogenic to other host species in other countries. The ability of these phytopathogens to produce and secrete cutinase is required for dismantling the cuticle of the host plant and, therefore, crucial to the necrotrophic phase of their infection strategy. In vitro bioassays using different lipid substrates confirmed the ability of C. truncatum and C. brevisporum isolates from green and red bell peppers to secrete cutinase. The diversity, structure and organization and synteny of the cutinase gene were determined among different Colletotrichum species. Cluster analysis indicated a low level of nucleotide variation among C. truncatum sequences. Nucleotide sequences of C. brevisporum were more related to C. truncatum cutinase nucleotide sequences than to C. gloeosporioides. Cluster patterns coincided with haplotype and there was evidence of significant positive selection with no recombination signatures. The structure of the cutinase gene included two exons with one intervening intron and, therefore, one splice variant. Although amino acid sequences were highly conserved among C. truncatum isolates, diversity "hot spots" were revealed when the 66-amino acid coding region of 200 fungal species was compared. Twenty cutinase orthologues were detected among different fungal species, whose common ancestor is Pezizomycotina and it is purported that these orthologues arose through a single gene duplication event prior to speciation. The cutinase domain was retained both in structure and arrangement among 34 different Colletotrichum species. The order of aligned genomic blocks between species and the arrangement of flanking protein domains were also conserved and shared for those domains immediately located at the N- and C-terminus of the cutinase domain. Among these were an RNA recognition motif, translation elongation factor, signal peptide, pentatricopeptide repeat, and Hsp70 family of chaperone proteins, all of which support the expression of the cutinase gene. The findings of this study are important to understanding the evolution of the cutinase gene in C. truncatum as a key component of the biotrophic-necrotrophic switch which may be useful in developing gene-targeting strategies to decrease the pathogenic potential of Colletotrichum species.

proGenomes2: an improved database for accurate and consistent habitat, taxonomic and functional annotations of prokaryotic genomes

Microbiology depends on the availability of annotated microbial genomes for many applications. Comparative genomics approaches have been a major advance, but consistent and accurate annotations of genomes can be hard to obtain. In addition, newer concepts such as the pan-genome concept are still being implemented to help answer biological questions. Hence, we present proGenomes2, which provides 87 920 high-quality genomes in a user-friendly and interactive manner. Genome sequences and annotations can be retrieved individually or by taxonomic clade. Every genome in the database has been assigned to a species cluster and most genomes could be accurately assigned to one or multiple habitats. In addition, general functional annotations and specific annotations of antibiotic resistance genes and single nucleotide variants are provided. In short, proGenomes2 provides threefold more genomes, enhanced habitat annotations, updated taxonomic and functional annotation and improved linkage to the NCBI BioSample database. The database is available at http://progenomes.embl.de/.

The Use of PTI-Marker Genes to Identify Novel Compounds that Establish Induced Resistance in Rice

Compounds that establish induced resistance (IR) in plants are promising alternatives for the pesticides that are progressively being banned worldwide. Screening platforms to identify IR-establishing compounds have been developed, but none were specifically designed for monocot plants. Here, we propose the use of an RT-qPCR screening platform, based on conserved immunity marker genes of rice as proxy for IR induction. Central regulators of biotic stress responses of rice were identified with a weighted gene co-expression network analysis (WGCNA), using more than 350 microarray datasets of rice under various sorts of biotic stress. Candidate genes were narrowed down to six immunity marker genes, based on consistent association with pattern-triggered immunity (PTI), both in rice plants as in rice cell suspension cultures (RCSCs). By monitoring the expression of these genes in RCSCs upon treatment with candidate IR-inducing compounds, we showed that our marker genes can predict IR induction in rice. Diproline, a novel IR-establishing compound for monocots that was detected with these marker genes, was shown to induce rice resistance against root-knot nematodes, without fitness costs. Gene expression profiling of the here-described PTI-marker genes can be executed on fully-grown plants or in RCSCs, providing a novel and versatile tool to predict IR induction.

Identification of New Leaf Rust Resistance Loci in Wheat and Wild Relatives by Array-Based SNP Genotyping and Association Genetics

Leaf rust caused by Puccinia triticina is the most widespread rust disease of wheat. As pathogen populations are constantly evolving, identification of novel sources of resistance is necessary to maintain disease resistance and stay ahead of this plant-pathogen evolutionary arms race. The wild genepool of wheat is a rich source of genetic diversity, accounting for 44% of the Lr genes identified. Here we performed a genome-wide association study (GWAS) on a diverse germplasm of 385 accessions, including 27 different Triticum and Aegilops species. Genetic characterization using the wheat 90 K array and subsequent filtering identified a set of 20,501 single nucleotide polymorphic (SNP) markers. Of those, 9,570 were validated using exome capture and mapped onto the Chinese Spring reference sequence v1.0. Phylogenetic analyses illustrated four major clades, clearly separating the wild species from the T. aestivum and T. turgidum species. GWAS was conducted using eight statistical models for infection types against six leaf rust isolates and leaf rust severity rated in field trials for 3-4 years at 2-3 locations in Canada. Functional annotation of genes containing significant quantitative trait nucleotides (QTNs) identified 96 disease-related loci associated with leaf rust resistance. A total of 21 QTNs were in haplotype blocks or within flanking markers of at least 16 known Lr genes. The remaining significant QTNs were considered loci that putatively harbor new Lr resistance genes. Isolation of these candidate genes will contribute to the elucidation of their role in leaf rust resistance and promote their usefulness in marker-assisted selection and introgression.

Stochastic models for description of structural-statistical properties in DNA sequences

New stochastic models based on a notion of stochastic codon are proposed. These models, presented by special random strings, describe practical structural-statistical properties which are peculiar to coding DNA both from prokaryotic and eukaryotic genomes. In such the case coding regions are considered as the realizations of random strings. The models introduced explain existence of latent profile periodicity with a period which is not only equal to but also multiplied of three in the coding regions. For the sequences with latent profile period multiplied of three, but not equal to three, the proposed models ensure existence of special property of 3-regularity in these sequences which is practically recognized in all coding sequences of the genomes analyzed. Feasibility of the stochastic models proposed was tested in numerical experiments with binary reencoded paragraphs of literary texts (in English and Italian languages), used as analog of DNA coding regions.

PHI-Nets: A Network Resource for Ascomycete Fungal Pathogens to Annotate and Identify Putative Virulence Interacting Proteins and siRNA Targets

Interactions between proteins underlie all aspects of complex biological mechanisms. Therefore, methodologies based on complex network analyses can facilitate identification of promising candidate genes involved in phenotypes of interest and put this information into appropriate contexts. To facilitate discovery and gain additional insights into globally important pathogenic fungi, we have reconstructed computationally inferred interactomes using an interolog and domain-based approach for 15 diverse Ascomycete fungal species, across nine orders, specifically Aspergillus fumigatus, Bipolaris sorokiniana, Blumeria graminis f. sp. hordei, Botrytis cinerea, Colletotrichum gloeosporioides, Colletotrichum graminicola, Fusarium graminearum, Fusarium oxysporum f. sp. lycopersici, Fusarium verticillioides, Leptosphaeria maculans, Magnaporthe oryzae, Saccharomyces cerevisiae, Sclerotinia sclerotiorum, Verticillium dahliae, and Zymoseptoria tritici. Network cartography analysis was associated with functional patterns of annotated genes linked to the disease-causing ability of each pathogen. In addition, for the best annotated organism, namely F. graminearum, the distribution of annotated genes with respect to network structure was profiled using a random walk with restart algorithm, which suggested possible co-location of virulence-related genes in the protein–protein interaction network. In a second ‘use case’ study involving two networks, namely B. cinerea and F. graminearum, previously identified small silencing plant RNAs were mapped to their targets. The F. graminearum phenotypic network analysis implicates eight B. cinerea targets and 35 F. graminearum predicted interacting proteins as prime candidate virulence genes for further testing. All 15 networks have been made accessible for download at www.phi-base.org providing a rich resource for major crop plant pathogens.

Stress-induced changes in the expression of antioxidant system genes for rice (Oryza sativa L.) and bread wheat (Triticum aestivum L.)

BACKGROUND

Plant cell metabolism inevitably forms reactive oxygen species (ROS), which can damage cells or lead to their death. The antioxidant system (AOS) evolved to eliminate a high concentration of ROS. For plants, this system consists of the seven classes of antioxidant enzymes and antioxidant compounds. Each enzymatic class contains a various number of genes which may vary from species to species. In such a multi-copy genetic system, the integration of evolutionary characteristics and expression data makes it possible to effectively predict promising breeding targets for the design of highly-yielding cultivars. In the plant cells, ROS production can increase as a result of abiotic stresses. Accordingly, AOS responds to stress by altering the expression of the genes of its components. Expression profiles of AOS enzymes, including their changes under stress, remains incomplete. A comprehensive study of the system behavior in response to stress for different species gives the key to identify the general mechanisms of AOS regulation. In this article, we studied stress-induced changes in the expression of AOS genes in photosynthetic tissues for rice and bread wheat.

METHODS

A meta-analysis of genome-wide transcriptome data on stress-induced changes in expression profiles of antioxidant genes using microarray and next generation sequencing (NGS) experiments from the GEO NCBI database for rice and bread wheat was carried out. Experimental study of expression changes in short (6 h) and prolonged (24 h) cold stress responses for selected AOS genes of bread wheat cultivars Saratovskaya29 and Yanetzkis Probat was conducted using qPCR.

RESULTS

The large-scale meta-transcriptome and complementary experimental analysis revealed a summary of fold changes in the AOS gene expression in response to cold and water deficiency for rice and bread wheat.

Genome-wide identification and expression analysis of the StSWEET family genes in potato (Solanum tuberosum L.)

BACKGROUND

The sugar will eventually be exported transporter (SWEET) family is a novel type of membrane-embedded sugar transporter that contains seven transmembrane helices with two MtN3/saliva domains. The SWEET family plays crucial roles in multiple processes, including carbohydrate transportation, development, environmental adaptability and host-pathogen interactions. Although SWEET genes, especially those involved in response to biotic stresses, have been extensively characterized in many plants, they have not yet been studied in potato.

OBJECTIVE

The identification of StSWEET genes provides important candidates for further functional analysis and lays the foundation for the production of good quality and high yield potatoes through molecular breeding.

METHODS

In this study, StSWEET genes were identified using a genome-wide search method. A comprehensive analysis of StSWEET family through bioinformatics methods, such as phylogenetic tree, gene structure and promoter prediction analysis. The expression profiles of StSWEET genes in different potato tissues and under P. infestans attack and sugar stress were studied using quantitative real-time polymerase chain reaction (qRT-PCR).

RESULTS

Phylogenetic analysis classified 33 StSWEET genes into four groups containing 12, 5, 12 and 4 genes. Furthermore, the gene structures and conserved motifs found that the StSWEET genes are very conservative during evolution. The chromosomal localization pattern showed that the distribution and density of the StSWEETs on 10 potato chromosomes were uneven and basically clustered. Predictive promoter analysis indicated that StSWEET proteins are associated with cell growth, development, secondary metabolism, and response to biotic and abiotic stresses. Finally, the expression patterns of the StSWEET genes in different tissues and the induction of P. infestans and the process of the sugar stress were investigated to obtain the tissue-specific and stress-responsive candidates.

CONCLUSION

This study systematically identifies the SWEET gene family in potato at the genome-wide level, providing important candidates for further functional analysis and contributing to a better understanding of the molecular basis of development and tolerance in potato.

Duplication and diversification of trehalase confers evolutionary advantages on lepidopteran insects

Gene duplication provides a major source of new genes for evolutionary novelty and ecological adaptation. However, the maintenance of duplicated genes and their relevance to adaptive evolution has long been debated. Insect trehalase (Treh) plays key roles in energy metabolism, growth, and stress recovery. Here, we show that the duplication of Treh in Lepidoptera (butterflies and moths) is linked with their adaptation to various environmental stresses. Generally, two Treh genes are present in insects: Treh1 and Treh2. We report three distinct forms of Treh in lepidopteran insects, where Treh1 was duplicated into two gene clusters (Treh1a and Treh1b). These gene clusters differ in gene expression patterns, enzymatic properties, and subcellular localizations, suggesting that the enzymes probably underwent sub- and/or neofunctionalization in the lepidopteran insects. Interestingly, selective pressure analysis provided significant evidence of positive selection on duplicate Treh1b gene in lepidopteran insect lineages. Most positively selected sites were located in the alpha-helical region, and several sites were close to the trehalose binding and catalytic sites. Subcellular adaptation of duplicate Treh1b driven by positive selection appears to have occurred as a result of selected changes in specific sequences, allowing for rapid reprogramming of duplicated Treh during evolution. Our results suggest that gene duplication of Treh and subsequent functional diversification could increase the survival rate of lepidopteran insects through various regulations of intracellular trehalose levels, facilitating their adaptation to diverse habitats. This study provides evidence regarding the mechanism by which gene family expansion can contribute to species adaptation through gene duplication and subsequent functional diversification.

Identification of a metallothionein gene in honey bee Apis mellifera and its expression profile in response to Cd, Cu and Pb exposure

Metallothioneins are ubiquitous proteins important in metal homeostasis and detoxification. However, they have not previously been identified in honey bees or other Hymenoptera, where metallothioneins could be of ecophysiological and ecotoxicological significance. Better understanding of the molecular responses to stress induced by toxic metals could contribute to honey bee conservation. In addition, honey bee metallothionein could represent a biomarker for monitoring environmental quality. Here we identify and characterize a metallothionein gene in Apis mellifera (AmMT). AmMT is 1,680 bp long and encodes a 48 amino acids protein with 15 cysteines and no aromatic residues. A metal response element upstream of the start codon, coupled with numerous cis-regulatory elements indicate the functional context of AmMT. Molecular modelling predicts several transition metal binding sites, and comparative phylogenetic analysis revealed five putative metallothionein proteins in three other hymenoptera species. AmMT was characterized by cloning the full-length coding sequence of the putative metallothionein. Recombinant AmMT was found to increase metal tolerance upon overexpression in Escherichia coli supplemented with Cd, Cu or Pb. Finally, in laboratory tests on honey bees, gene expression profiles showed a dose-dependant relationship between Cd, Cu and Pb concentrations present in food and AmMT expression, while field experiments showed induction of AmMT in bees from an industrial site compared to those from an urban area. These studies suggest that AmMT has metal binding properties in agreement with a possible role in metal homeostasis. Further functional and structural characterization of metallothionein in honey bees and other Hymenoptera are necessary.

'Why genes in pieces?'-revisited

The alignment between the boundaries of protein domains and the boundaries of exons could provide evidence for the evolution of proteins via domain shuffling, but literature in the field has so far struggled to conclusively show this. Here, on larger data sets than previously possible, we do finally show that this phenomenon is indisputably found widely across the eukaryotic tree. In contrast, the alignment between exons and the boundaries of intrinsically disordered regions of proteins is not a general property of eukaryotes. Most interesting of all is the discovery that domain-exon alignment is much more common in recently evolved protein sequences than older ones.