Increased burden of ultra-rare structural variants localizing to boundaries of topologically associated domains in schizophrenia.

Despite considerable progress in schizophrenia genetics, most findings have been for large rare structural variants and common variants in well-imputed regions with few genes implicated from exome sequencing. Whole genome sequencing (WGS) can potentially provide a more complete enumeration of etiological genetic variation apart from the exome and regions of high linkage disequilibrium. We analyze high-coverage WGS data from 1162 Swedish schizophrenia cases and 936 ancestry-matched population controls. Our main objective is to evaluate the contribution to schizophrenia etiology from a variety of genetic variants accessible to WGS but not by previous technologies. Our results suggest that ultra-rare structural variants that affect the boundaries of topologically associated domains (TADs) increase risk for schizophrenia. Alterations in TAD boundaries may lead to dysregulation of gene expression. Future mechanistic studies will be needed to determine the precise functional effects of these variants on biology.

A geographically matched control population efficiently limits the number of candidate disease-causing variants in an unbiased whole-genome analysis.

Whole-genome sequencing is a promising approach for human autosomal dominant disease studies. However, the vast number of genetic variants observed by this method constitutes a challenge when trying to identify the causal variants. This is often handled by restricting disease studies to the most damaging variants, e.g. those found in coding regions, and overlooking the remaining genetic variation. Such a biased approach explains in part why the genetic causes of many families with dominantly inherited diseases, in spite of being included in whole-genome sequencing studies, are left unsolved today. Here we explore the use of a geographically matched control population to minimize the number of candidate disease-causing variants without excluding variants based on assumptions on genomic position or functional predictions. To exemplify the benefit of the geographically matched control population we apply a typical disease variant filtering strategy in a family with an autosomal dominant form of colorectal cancer. With the use of the geographically matched control population we end up with 26 candidate variants genome wide. This is in contrast to the tens of thousands of candidates left when only making use of available public variant datasets. The effect of the local control population is dual, it (1) reduces the total number of candidate variants shared between affected individuals, and more importantly (2) increases the rate by which the number of candidate variants are reduced as additional affected family members are included in the filtering strategy. We demonstrate that the application of a geographically matched control population effectively limits the number of candidate disease-causing variants and may provide the means by which variants suitable for functional studies are identified genome wide.

Deletion of FtsH11 protease has impact on chloroplast structure and function in Arabidopsis thaliana when grown under continuous light.

The membrane-integrated metalloprotease FtsH11 of Arabidopsis thaliana is proposed to be dual-targeted to mitochondria and chloroplasts. A bleached phenotype was observed in ftsh11 grown at long days or continuous light, pointing to disturbances in the chloroplast. Within the chloroplast, FtsH11 was found to be located exclusively in the envelope. Two chloroplast-located proteins of unknown function (Tic22-like protein and YGGT-A) showed significantly higher abundance in envelope membranes and intact chloroplasts of ftsh11 and therefore qualify as potential substrates for the FtsH11 protease. No proteomic changes were observed in the mitochondria of 6-week-old ftsh11 compared with wild type, and FtsH11 was not immunodetected in these organelles. The abundance of plastidic proteins, especially of photosynthetic proteins, was altered even during standard growth conditions in total leaves of ftsh11. At continuous light, the amount of photosystem I decreased relative to photosystem II, accompanied by a drastic change of the chloroplast morphology and a drop of non-photochemical quenching. FtsH11 is crucial for chloroplast structure and function during growth in prolonged photoperiod.

The Plant Genome Integrative Explorer Resource: PlantGenIE.org.

Accessing and exploring large-scale genomics data sets remains a significant challenge to researchers without specialist bioinformatics training. We present the integrated PlantGenIE.org platform for exploration of Populus, conifer and Arabidopsis genomics data, which includes expression networks and associated visualization tools. Standard features of a model organism database are provided, including genome browsers, gene list annotation, Blast homology searches and gene information pages. Community annotation updating is supported via integration of WebApollo. We have produced an RNA-sequencing (RNA-Seq) expression atlas for Populus tremula and have integrated these data within the expression tools. An updated version of the ComPlEx resource for performing comparative plant expression analyses of gene coexpression network conservation between species has also been integrated. The PlantGenIE.org platform provides intuitive access to large-scale and genome-wide genomics data from model forest tree species, facilitating both community contributions to annotation improvement and tools supporting use of the included data resources to inform biological insight.

Synergy: a web resource for exploring gene regulation in Synechocystis sp. PCC6803.

Despite being a highly studied model organism, most genes of the cyanobacterium Synechocystis sp. PCC 6803 encode proteins with completely unknown function. To facilitate studies of gene regulation in Synechocystis, we have developed Synergy (http://synergy.plantgenie.org), a web application integrating co-expression networks and regulatory motif analysis. Co-expression networks were inferred from publicly available microarray experiments, while regulatory motifs were identified using a phylogenetic footprinting approach. Automatically discovered motifs were shown to be enriched in the network neighborhoods of regulatory proteins much more often than in the neighborhoods of non-regulatory genes, showing that the data provide a sound starting point for studying gene regulation in Synechocystis. Concordantly, we provide several case studies demonstrating that Synergy can be used to find biologically relevant regulatory mechanisms in Synechocystis. Synergy can be used to interactively perform analyses such as gene/motif search, network visualization and motif/function enrichment. Considering the importance of Synechocystis for photosynthesis and biofuel research, we believe that Synergy will become a valuable resource to the research community.

ComPlEx: conservation and divergence of co-expression networks in A. thaliana, Populus and O. sativa.

BACKGROUND: Divergence in gene regulation has emerged as a key mechanism underlying species differentiation. Comparative analysis of co-expression networks across species can reveal conservation and divergence in the regulation of genes. RESULTS: We inferred co-expression networks of A. thaliana, Populus spp. and O. sativa using state-of-the-art methods based on mutual information and context likelihood of relatedness, and conducted a comprehensive comparison of these networks across a range of co-expression thresholds. In addition to quantifying gene-gene link and network neighbourhood conservation, we also applied recent advancements in network analysis to do cross-species comparisons of network properties such as scale free characteristics and gene centrality as well as network motifs. We found that in all species the networks emerged as scale free only above a certain co-expression threshold, and that the high-centrality genes upholding this organization tended to be conserved. Network motifs, in particular the feed-forward loop, were found to be significantly enriched in specific functional subnetworks but where much less conserved across species than gene centrality. Although individual gene-gene co-expression had massively diverged, up to ~80% of the genes still had a significantly conserved network neighbourhood. For genes with multiple predicted orthologs, about half had one ortholog with conserved regulation and another ortholog with diverged or non-conserved regulation. Furthermore, the most sequence similar ortholog was not the one with the most conserved gene regulation in over half of the cases. CONCLUSIONS: We have provided a comprehensive analysis of gene regulation evolution in plants and built a web tool for Comparative analysis of Plant co-Expression networks (ComPlEx, http://complex.plantgenie.org/). The tool can be particularly useful for identifying the ortholog with the most conserved regulation among several sequence-similar alternatives and can thus be of practical importance in e.g. finding candidate genes for perturbation experiments.

Co-expression analysis, proteomic and metabolomic study on the impact of a Deg/HtrA protease triple mutant in Synechocystis sp. PCC 6803 exposed to temperature and high light stress.

Members of the DegP/HtrA protease family are widespread in nature and play an important role in proteolysis of misfolded and damaged proteins. The cyanobacterium Synechocystis sp. PCC 6803 contains three Deg proteases, HhoA (Sll1679), HhoB (Sll1427) and HtrA (Slr1204). Using the proteomic or metabolomic approach we investigated a triple deletion mutant (Δdeg) exposed to light or temperature stress. To cope with the stress conditions the triple mutant reduces its energy metabolism and stress-related proteins are induced to protect the cells. Additionally the co-expression of the genes encoding the three proteases with other genes in Synechocystis sp. PCC 6803 was analyzed. While HhoA seems to be involved in house-keeping processes related to protein (re)folding, protein clearance and signaling, the hhoB expression cluster is dominated by genes encoding periplasmic proteins linked to metabolism or signal transduction pathways. The htrA expression pattern is similar to that of genes encoding proteins of the electron transport chain, iron- and nitrogen metabolism. Our integrative approach indicates significant rearrangements in cells depleted of the Deg/HtrA proteases when exposed to stress, both, in the cytoplasmic and extracytoplasmic space.

Co-swarming and local collapse: quorum sensing conveys resilience to bacterial communities by localizing cheater mutants in Pseudomonas aeruginosa.

BACKGROUND: Members of swarming bacterial consortia compete for nutrients but also use a co-operation mechanism called quorum sensing (QS) that relies on chemical signals as well as other secreted products ("public goods") necessary for swarming. Deleting various genes of this machinery leads to cheater mutants impaired in various aspects of swarming cooperation. METHODOLOGY/PRINCIPAL FINDINGS: Pairwise consortia made of Pseudomonas aeruginosa, its QS mutants as well as B. cepacia cells show that a interspecies consortium can "combine the skills" of its participants so that the strains can cross together barriers that they could not cross alone. In contrast, deleterious mutants are excluded from consortia either by competition or by local population collapse. According to modeling, both scenarios are the consequence of the QS signalling mechanism itself. CONCLUSION/SIGNIFICANCE: The results indirectly explain why it is an advantage for bacteria to maintain QS systems that can cross-talk among different species, and conversely, why certain QS mutants which can be abundant in isolated niches, cannot spread and hence remain localized.

A simple model for the early events of quorum sensing in Pseudomonas aeruginosa: modeling bacterial swarming as the movement of an "activation zone".

BACKGROUND: Quorum sensing (QS) is a form of gene regulation based on cell-density that depends on inter-cellular communication. While there are a variety of models for bacterial colony morphology, there is little work linking QS genes to movement in an open system. RESULTS: The onset of swarming in environmental P. aeruginosa PUPa3 was described with a simplified computational model in which cells in random motion communicate via a diffusible signal (representing N-acyl homoserine lactones, AHL) as well as diffusible, secreted factors (enzymes, biosurfactans, i.e. "public goods") that regulate the intensity of movement and metabolism in a threshold-dependent manner. As a result, an "activation zone" emerges in which nutrients and other public goods are present in sufficient quantities, and swarming is the spontaneous displacement of this high cell-density zone towards nutrients and/or exogenous signals. The model correctly predicts the behaviour of genomic knockout mutants in which the QS genes responsible either for the synthesis (lasI, rhlI) or the sensing (lasR, rhlR) of AHL signals were inactivated. For wild type cells the model predicts sustained colony growth that can however be collapsed by the overconsumption of nutrients. CONCLUSION: While in more complex models include self-orienting abilities that allow cells to follow concentration gradients of nutrients and chemotactic agents, in this model, displacement towards nutrients or environmental signals is an emergent property of the community that results from the action of a few, well-defined QS genes and their products. Still the model qualitatively describes the salient properties of QS bacteria, i.e. the density-dependent onset of swarming as well as the response to exogenous signals or cues.

Evolution of robust and efficient system topologies.

Mutation/selection algorithms were applied to increase the efficiency and the robustness of sparse random networks. Selection for better efficiency leads to the well-known star topology, while selection for robustness only results in a relatively dense core and a small periphery. Concomitant selection for both efficiency and robustness leads to networks with intermittent center/periphery values. Networks evolving under multiple attack regimes develop distinct topologies with larger cores, and are characterized by parameter distributions different from those developing under single-attack regimes.