Center for Computational and Theoretical Biology
Refine
Has Fulltext
- yes (40)
Is part of the Bibliography
- yes (40)
Document Type
- Journal article (40) (remove)
Keywords
- machine learning (5)
- active zone (4)
- Caenorhabditis elegans (2)
- bee decline (2)
- dSTORM (2)
- deep learning (2)
- evolution (2)
- foraging (2)
- juvenile hormone (2)
- metabarcoding (2)
Institute
- Center for Computational and Theoretical Biology (40)
- Theodor-Boveri-Institut für Biowissenschaften (19)
- Pathologisches Institut (4)
- Physiologisches Institut (4)
- Julius-von-Sachs-Institut für Biowissenschaften (3)
- Klinik und Poliklinik für Unfall-, Hand-, Plastische und Wiederherstellungschirurgie (Chirurgische Klinik II) (3)
- Neurochirurgische Klinik und Poliklinik (3)
- Urologische Klinik und Poliklinik (3)
- Comprehensive Cancer Center Mainfranken (2)
- Klinik und Poliklinik für Dermatologie, Venerologie und Allergologie (2)
EU-Project number / Contract (GA) number
- 2020010013 (1)
- 250194-Carnivorom (1)
- 835102) (1)
Tropical peat swamp forests sequester globally significant stores of carbon in deep layers of waterlogged, anoxic, acidic and nutrient-depleted peat. The roles of microbes in supporting these forests through the formation of peat, carbon sequestration and nutrient cycling are virtually unknown. This study investigated physicochemical peat properties and microbial diversity between three dominant tree species: Shorea uliginosa (Dipterocarpaceae), Koompassia malaccensis (legumes associated with nitrogen-fixing bacteria), Eleiodoxa conferta (palm) and depths (surface, 45 and 90 cm) using microbial 16S rRNA gene amplicon sequencing. Water pH, oxygen, nitrogen, phosphorus, total phenolic contents and C/N ratio differed significantly between depths, but not tree species. Depth also strongly influenced microbial diversity and composition, while both depth and tree species exhibited significant impact on the archaeal communities. Microbial diversity was highest at the surface, where fresh leaf litter accumulates, and nutrient supply is guaranteed. Nitrogen was the core parameter correlating to microbial communities, but the interactive effects from various environmental variables displayed significant correlation to relative abundance of major microbial groups. Proteobacteria was the dominant phylum and the most abundant genus, Rhodoplanes, might be involved in nitrogen fixation. The most abundant methanogens and methanotrophs affiliated, respectively, to families Methanomassiliicoccaceae and Methylocystaceae. Our results demonstrated diverse microbial communities and provide valuable insights on microbial ecology in these extreme ecosystems.
Young grapevines (Vitis vinifera) suffer and eventually can die from the crown gall disease caused by the plant pathogen Allorhizobium vitis (Rhizobiaceae). Virulent members of A. vitis harbor a tumor-inducing plasmid and induce formation of crown galls due to the oncogenes encoded on the transfer DNA. The expression of oncogenes in transformed host cells induces unregulated cell proliferation and metabolic and physiological changes. The crown gall produces opines uncommon to plants, which provide an important nutrient source for A. vitis harboring opine catabolism enzymes. Crown galls host a distinct bacterial community, and the mechanisms establishing a crown gall–specific bacterial community are currently unknown. Thus, we were interested in whether genes homologous to those of the tumor-inducing plasmid coexist in the genomes of the microbial species coexisting in crown galls. We isolated 8 bacterial strains from grapevine crown galls, sequenced their genomes, and tested their virulence and opine utilization ability in bioassays. In addition, the 8 genome sequences were compared with 34 published bacterial genomes, including closely related plant-associated bacteria not from crown galls. Homologous genes for virulence and opine anabolism were only present in the virulent Rhizobiaceae. In contrast, homologs of the opine catabolism genes were present in all strains including the nonvirulent members of the Rhizobiaceae and non-Rhizobiaceae. Gene neighborhood and sequence identity of the opine degradation cluster of virulent and nonvirulent strains together with the results of the opine utilization assay support the important role of opine utilization for cocolonization in crown galls, thereby shaping the crown gall community.
Abstract
Cell lineage decisions occur in three-dimensional spatial patterns that are difficult to identify by eye. There is an ongoing effort to replicate such patterns using mathematical modeling. One approach uses long ranging cell-cell communication to replicate common spatial arrangements like checkerboard and engulfing patterns. In this model, the cell-cell communication has been implemented as a signal that disperses throughout the tissue. On the other hand, machine learning models have been developed for pattern recognition and pattern reconstruction tasks. We combined synthetic data generated by the mathematical model with spatial summary statistics and deep learning algorithms to recognize and reconstruct cell fate patterns in organoids of mouse embryonic stem cells. Application of Moran’s index and pair correlation functions for in vitro and synthetic data from the model showed local clustering and radial segregation. To assess the patterns as a whole, a graph neural network was developed and trained on synthetic data from the model. Application to in vitro data predicted a low signal dispersion value. To test this result, we implemented a multilayer perceptron for the prediction of a given cell fate based on the fates of the neighboring cells. The results show a 70% accuracy of cell fate imputation based on the nine nearest neighbors of a cell. Overall, our approach combines deep learning with mathematical modeling to link cell fate patterns with potential underlying mechanisms.
Author summary
Mammalian embryo development relies on organized differentiation of stem cells into different lineages. Particularly at the early stages of embryogenesis, cells of different fates form three-dimensional spatial patterns that are difficult to identify by eye. Pattern quantification and mathematical modeling have produced first insights into potential mechanisms for the cell fate arrangements. However, these approaches have relied on classifications of the patterns such as inside-out or random, or used summary statistics such as pair correlation functions or cluster radii. Deep neural networks allow characterizing patterns directly. Since the tissue context can be readily reproduced by a graph, we implemented a graph neural network to characterize the patterns of embryonic stem cell organoids as a whole. In addition, we implemented a multilayer perceptron model to reconstruct the fate of a given cell based on its neighbors. To train and test the models, we used synthetic data generated by our mathematical model for cell-cell communication. This interplay of deep learning and mathematical modeling in combination with summary statistics allowed us to identify a potential mechanism for cell fate determination in mouse embryonic stem cells. Our results agree with a mechanism with a dispersion of the intercellular signal that links a cell’s fate to those of the local neighborhood.
Summary
Embryos develop in a concerted sequence of spatiotemporal arrangements of cells. In the preimplantation mouse embryo, the distribution of the cells in the inner cell mass evolves from a salt-and-pepper pattern to spatial segregation of two distinct cell types. The exact properties of the salt-and-pepper pattern have not been analyzed so far. We investigate the spatiotemporal distribution of NANOG- and GATA6-expressing cells in the ICM of the mouse blastocysts with quantitative three-dimensional single-cell-based neighborhood analyses. A combination of spatial statistics and agent-based modeling reveals that the cell fate distribution follows a local clustering pattern. Using ordinary differential equations modeling, we show that this pattern can be established by a distance-based signaling mechanism enabling cells to integrate information from the whole inner cell mass into their cell fate decision. Our work highlights the importance of longer-range signaling to ensure coordinated decisions in groups of cells to successfully build embryos.
Highlights
• The local cell neighborhood and global ICM population composition correlate
• ICM cells show characteristics of local clustering in early and mid mouse blastocysts
• ICM patterning requires integration of signals from cells beyond the first neighbors
Understanding the causal relationship between genotype and phenotype is a major objective in biology. Genome-wide association studies (GWAS) correlate genetic polymorphisms with trait variation and have already identified causative variants for various traits in many different organisms, from humans to plants. Importantly, many adaptive traits, like the regulation of flowering time in plants, are not regulated by distinct genetic effects, but by more sophisticated gene regulatory networks.
(1) Background: C-X-C Motif Chemokine Receptor 4 (CXCR4) and Fibroblast Activation Protein Alpha (FAP) are promising theranostic targets. However, it is unclear whether CXCR4 and FAP positivity mark distinct microenvironments, especially in solid tumors. (2) Methods: Using Random Forest (RF) analysis, we searched for entity-independent mRNA and microRNA signatures related to CXCR4 and FAP overexpression in our pan-cancer cohort from The Cancer Genome Atlas (TCGA) database — representing n = 9242 specimens from 29 tumor entities. CXCR4- and FAP-positive samples were assessed via StringDB cluster analysis, EnrichR, Metascape, and Gene Set Enrichment Analysis (GSEA). Findings were validated via correlation analyses in n = 1541 tumor samples. TIMER2.0 analyzed the association of CXCR4 / FAP expression and infiltration levels of immune-related cells. (3) Results: We identified entity-independent CXCR4 and FAP gene signatures representative for the majority of solid cancers. While CXCR4 positivity marked an immune-related microenvironment, FAP overexpression highlighted an angiogenesis-associated niche. TIMER2.0 analysis confirmed characteristic infiltration levels of CD8+ cells for CXCR4-positive tumors and endothelial cells for FAP-positive tumors. (4) Conclusions: CXCR4- and FAP-directed PET imaging could provide a non-invasive decision aid for entity-agnostic treatment of microenvironment in solid malignancies. Moreover, this machine learning workflow can easily be transferred towards other theranostic targets.
Osmotic stress can be detrimental to plants, whose survival relies heavily on proteomic plasticity. Protein ubiquitination is a central post-translational modification in osmotic-mediated stress. In this study, we used the K-Ɛ-GG antibody enrichment method integrated with high-resolution mass spectrometry to compile a list of 719 ubiquitinated lysine (K-Ub) residues from 450 Arabidopsis root membrane proteins (58% of which are transmembrane proteins), thereby adding to the database of ubiquitinated substrates in plants. Although no ubiquitin (Ub) motifs could be identified, the presence of acidic residues close to K-Ub was revealed. Our ubiquitinome analysis pointed to a broad role of ubiquitination in the internalization and sorting of cargo proteins. Moreover, the simultaneous proteome and ubiquitinome quantification showed that ubiquitination is mostly not involved in membrane protein degradation in response to short osmotic treatment but that it is putatively involved in protein internalization, as described for the aquaporin PIP2;1. Our in silico analysis of ubiquitinated proteins shows that two E2 Ub-conjugating enzymes, UBC32 and UBC34, putatively target membrane proteins under osmotic stress. Finally, we revealed a positive role for UBC32 and UBC34 in primary root growth under osmotic stress.
Automatic image reconstruction is critical to cope with steadily increasing data from advanced microscopy. We describe here the Fiji macro 3D ART VeSElecT which we developed to study synaptic vesicles in electron tomograms. We apply this tool to quantify vesicle properties (i) in embryonic Danio rerio 4 and 8 days past fertilization (dpf) and (ii) to compare Caenorhabditis elegans N2 neuromuscular junctions (NMJ) wild-type and its septin mutant (unc-59(e261)). We demonstrate development-specific and mutant-specific changes in synaptic vesicle pools in both models. We confirm the functionality of our macro by applying our 3D ART VeSElecT on zebrafish NMJ showing smaller vesicles in 8 dpf embryos then 4 dpf, which was validated by manual reconstruction of the vesicle pool. Furthermore, we analyze the impact of C. elegans septin mutant unc-59(e261) on vesicle pool formation and vesicle size. Automated vesicle registration and characterization was implemented in Fiji as two macros (registration and measurement). This flexible arrangement allows in particular reducing false positives by an optional manual revision step. Preprocessing and contrast enhancement work on image-stacks of 1nm/pixel in x and y direction. Semi-automated cell selection was integrated. 3D ART VeSElecT removes interfering components, detects vesicles by 3D segmentation and calculates vesicle volume and diameter (spherical approximation, inner/outer diameter). Results are collected in color using the RoiManager plugin including the possibility of manual removal of non-matching confounder vesicles. Detailed evaluation considered performance (detected vesicles) and specificity (true vesicles) as well as precision and recall. We furthermore show gain in segmentation and morphological filtering compared to learning based methods and a large time gain compared to manual segmentation. 3D ART VeSElecT shows small error rates and its speed gain can be up to 68 times faster in comparison to manual annotation. Both automatic and semi-automatic modes are explained including a tutorial.
Ultrastructural analysis of wild-type and RIM1α knockout active zones in a large cortical synapse
(2022)
Rab3A-interacting molecule (RIM) is crucial for fast Ca\(^{2+}\)-triggered synaptic vesicle (SV) release in presynaptic active zones (AZs). We investigated hippocampal giant mossy fiber bouton (MFB) AZ architecture in 3D using electron tomography of rapid cryo-immobilized acute brain slices in RIM1α\(^{−/−}\) and wild-type mice. In RIM1α\(^{−/−}\), AZs are larger with increased synaptic cleft widths and a 3-fold reduced number of tightly docked SVs (0–2 nm). The distance of tightly docked SVs to the AZ center is increased from 110 to 195 nm, and the width of their electron-dense material between outer SV membrane and AZ membrane is reduced. Furthermore, the SV pool in RIM1α\(^{−/−}\) is more heterogeneous. Thus, RIM1α, besides its role in tight SV docking, is crucial for synaptic architecture and vesicle pool organization in MFBs.
Biodiversity loss, as often found in intensively managed agricultural landscapes, correlates with reduced ecosystem functioning, for example, pollination by insects, and with altered plant composition, diversity, and abundance. But how does this change in floral resource diversity and composition relate to occurrence and resource use patterns of trap-nesting solitary bees? To better understand the impact of land-use intensification on communities of trap-nesting solitary bees in managed grasslands, we investigated their pollen foraging, reproductive fitness, and the nutritional quality of larval food along a land-use intensity gradient in Germany. We found bee species diversity to decrease with increasing land-use intensity irrespective of region-specific community compositions and interaction networks. Land use also strongly affected the diversity and composition of pollen collected by bees. Lack of suitable pollen sources likely explains the absence of several bee species at sites of high land-use intensity. The only species present throughout, Osmia bicornis (red mason bee), foraged on largely different pollen sources across sites. In doing so, it maintained a relatively stable, albeit variable nutritional quality of larval diets (i.e., protein to lipid (P:L) ratio). The observed changes in bee–plant pollen interaction patterns indicate that only the flexible generalists, such as O. bicornis, may be able to compensate the strong alterations in floral resource landscapes and to obtain food of sufficient quality through readily shifting to alternative plant sources. In contrast, other, less flexible, bee species disappear.
Introduction
Neurotransmitter release at presynaptic active zones (AZs) requires concerted protein interactions within a dense 3D nano-hemisphere. Among the complex protein meshwork the (M)unc-13 family member Unc-13 of Drosophila melanogaster is essential for docking of synaptic vesicles and transmitter release.
Methods
We employ minos-mediated integration cassette (MiMIC)-based gene editing using GFSTF (EGFP-FlAsH-StrepII-TEV-3xFlag) to endogenously tag all annotated Drosophila Unc-13 isoforms enabling visualization of endogenous Unc-13 expression within the central and peripheral nervous system.
Results and discussion
Electrophysiological characterization using two-electrode voltage clamp (TEVC) reveals that evoked and spontaneous synaptic transmission remain unaffected in unc-13\(^{GFSTF}\) 3rd instar larvae and acute presynaptic homeostatic potentiation (PHP) can be induced at control levels. Furthermore, multi-color structured-illumination shows precise co-localization of Unc-13\(^{GFSTF}\), Bruchpilot, and GluRIIA-receptor subunits within the synaptic mesoscale. Localization microscopy in combination with HDBSCAN algorithms detect Unc-13\(^{GFSTF}\) subclusters that move toward the AZ center during PHP with unaltered Unc-13\(^{GFSTF}\) protein levels.
Candida auris is a globally emerging fungal pathogen responsible for causing nosocomial outbreaks in healthcare associated settings. It is known to cause infection in all age groups and exhibits multi-drug resistance with high potential for horizontal transmission. Because of this reason combined with limited therapeutic choices available, C. auris infection has been acknowledged as a potential risk for causing a future pandemic, and thus seeking a promising strategy for its treatment is imperative. Here, we combined evolutionary information with reverse vaccinology approach to identify novel epitopes for vaccine design that could elicit CD4+ T-cell responses against C. auris. To this end, we extensively scanned the family of proteins encoded by C. auris genome. In addition, a pathogen may acquire substitutions in epitopes over a period of time which could cause its escape from the immune response thus rendering the vaccine ineffective. To lower this possibility in our design, we eliminated all rapidly evolving genes of C. auris with positive selection. We further employed highly conserved regions of multiple C. auris strains and identified two immunogenic and antigenic T-cell epitopes that could generate the most effective immune response against C. auris. The antigenicity scores of our predicted vaccine candidates were calculated as 0.85 and 1.88 where 0.5 is the threshold for prediction of fungal antigenic sequences. Based on our results, we conclude that our vaccine candidates have the potential to be successfully employed for the treatment of C. auris infection. However, in vivo experiments are imperative to further demonstrate the efficacy of our design.
Personalized oncology is a rapidly evolving area and offers cancer patients therapy options that are more specific than ever. However, there is still a lack of understanding regarding transcriptomic similarities or differences of metastases and corresponding primary sites. Applying two unsupervised dimension reduction methods (t-Distributed Stochastic Neighbor Embedding (t-SNE) and Uniform Manifold Approximation and Projection (UMAP)) on three datasets of metastases (n = 682 samples) with three different data transformations (unprocessed, log10 as well as log10 + 1 transformed values), we visualized potential underlying clusters. Additionally, we analyzed two datasets (n = 616 samples) containing metastases and primary tumors of one entity, to point out potential familiarities. Using these methods, no tight link between the site of resection and cluster formation outcome could be demonstrated, or for datasets consisting of solely metastasis or mixed datasets. Instead, dimension reduction methods and data transformation significantly impacted visual clustering results. Our findings strongly suggest data transformation to be considered as another key element in the interpretation of visual clustering approaches along with initialization and different parameters. Furthermore, the results highlight the need for a more thorough examination of parameters used in the analysis of clusters.
After the recent emergence of SARS-CoV-2 infection, unanswered questions remain related to its evolutionary history, path of transmission or divergence and role of recombination. There is emerging evidence on amino acid substitutions occurring in key residues of the receptor-binding domain of the spike glycoprotein in coronavirus isolates from bat and pangolins. In this article, we summarize our current knowledge on the origin of SARS-CoV-2. We also analyze the host ACE2-interacting residues of the receptor-binding domain of spike glycoprotein in SARS-CoV-2 isolates from bats, and compare it to pangolin SARS-CoV-2 isolates collected from Guangdong province (GD Pangolin-CoV) and Guangxi autonomous regions (GX Pangolin-CoV) of South China. Based on our comparative analysis, we support the view that the Guangdong Pangolins are the intermediate hosts that adapted the SARS-CoV-2 and represented a significant evolutionary link in the path of transmission of SARS-CoV-2 virus. We also discuss the role of intermediate hosts in the origin of Omicron.
Abstract
Introgressive hybridization is a process that enables gene flow across species barriers through the backcrossing of hybrids into a parent population. This may make genetic material, potentially including relevant environmental adaptations, rapidly available in a gene pool. Consequently, it has been postulated to be an important mechanism for enabling evolutionary rescue, that is the recovery of threatened populations through rapid evolutionary adaptation to novel environments. However, predicting the likelihood of such evolutionary rescue for individual species remains challenging. Here, we use the example of Zosterops silvanus, an endangered East African highland bird species suffering from severe habitat loss and fragmentation, to investigate whether hybridization with its congener Zosterops flavilateralis might enable evolutionary rescue of its Taita Hills population. To do so, we employ an empirically parameterized individual‐based model to simulate the species' behaviour, physiology and genetics. We test the population's response to different assumptions of mating behaviour and multiple scenarios of habitat change. We show that as long as hybridization does take place, evolutionary rescue of Z. silvanus is likely. Intermediate hybridization rates enable the greatest long‐term population growth, due to trade‐offs between adaptive and maladaptive introgressed alleles. Habitat change did not have a strong effect on population growth rates, as Z. silvanus is a strong disperser and landscape configuration is therefore not the limiting factor for hybridization. Our results show that targeted gene flow may be a promising avenue to help accelerate the adaptation of endangered species to novel environments, and demonstrate how to combine empirical research and mechanistic modelling to deliver species‐specific predictions for conservation planning.
The negative impact of juvenile undernourishment on adult behavior has been well reported for vertebrates, but relatively little is known about invertebrates. In honeybees, nutrition has long been known to affect task performance and timing of behavioral transitions. Whether and how a dietary restriction during larval development affects the task performance of adult honeybees is largely unknown. We raised honeybees in-vitro, varying the amount of a standardized diet (150 µl, 160 µl, 180 µl in total). Emerging adults were marked and inserted into established colonies. Behavioral performance of nurse bees and foragers was investigated and physiological factors known to be involved in the regulation of social organization were quantified. Surprisingly, adult honeybees raised under different feeding regimes did not differ in any of the behaviors observed. No differences were observed in physiological parameters apart from weight. Honeybees were lighter when undernourished (150 µl), while they were heavier under the overfed treatment (180 µl) compared to the control group raised under a normal diet (160 µl). These data suggest that dietary restrictions during larval development do not affect task performance or physiology in this social insect despite producing clear effects on adult weight. We speculate that possible effects of larval undernourishment might be compensated during the early period of adult life.
Background
Localization-based super-resolution microscopy resolves macromolecular structures down to a few nanometers by computationally reconstructing fluorescent emitter coordinates from diffraction-limited spots. The most commonly used algorithms are based on fitting parametric models of the point spread function (PSF) to a measured photon distribution. These algorithms make assumptions about the symmetry of the PSF and thus, do not work well with irregular, non-linear PSFs that occur for example in confocal lifetime imaging, where a laser is scanned across the sample. An alternative method for reconstructing sparse emitter sets from noisy, diffraction-limited images is compressed sensing, but due to its high computational cost it has not yet been widely adopted. Deep neural network fitters have recently emerged as a new competitive method for localization microscopy. They can learn to fit arbitrary PSFs, but require extensive simulated training data and do not generalize well. A method to efficiently fit the irregular PSFs from confocal lifetime localization microscopy combining the advantages of deep learning and compressed sensing would greatly improve the acquisition speed and throughput of this method.
Results
Here we introduce ReCSAI, a compressed sensing neural network to reconstruct localizations for confocal dSTORM, together with a simulation tool to generate training data. We implemented and compared different artificial network architectures, aiming to combine the advantages of compressed sensing and deep learning. We found that a U-Net with a recursive structure inspired by iterative compressed sensing showed the best results on realistic simulated datasets with noise, as well as on real experimentally measured confocal lifetime scanning data. Adding a trainable wavelet denoising layer as prior step further improved the reconstruction quality.
Conclusions
Our deep learning approach can reach a similar reconstruction accuracy for confocal dSTORM as frame binning with traditional fitting without requiring the acquisition of multiple frames. In addition, our work offers generic insights on the reconstruction of sparse measurements from noisy experimental data by combining compressed sensing and deep learning. We provide the trained networks, the code for network training and inference as well as the simulation tool as python code and Jupyter notebooks for easy reproducibility.
In vitro rearing of honeybee larvae is an established method that enables exact control and monitoring of developmental factors and allows controlled application of pesticides or pathogens. However, only a few studies have investigated how the rearing method itself affects the behavior of the resulting adult honeybees. We raised honeybees in vitro according to a standardized protocol: marking the emerging honeybees individually and inserting them into established colonies. Subsequently, we investigated the behavioral performance of nurse bees and foragers and quantified the physiological factors underlying the social organization. Adult honeybees raised in vitro differed from naturally reared honeybees in their probability of performing social tasks. Further, in vitro-reared bees foraged for a shorter duration in their life and performed fewer foraging trips. Nursing behavior appeared to be unaffected by rearing condition. Weight was also unaffected by rearing condition. Interestingly, juvenile hormone titers, which normally increase strongly around the time when a honeybee becomes a forager, were significantly lower in three- and four-week-old in vitro bees. The effects of the rearing environment on individual sucrose responsiveness and lipid levels were rather minor. These data suggest that larval rearing conditions can affect the task performance and physiology of adult bees despite equal weight, pointing to an important role of the colony environment for these factors. Our observations of behavior and metabolic pathways offer important novel insight into how the rearing environment affects adult honeybees.
Osmotic adaptation and accumulation of compatible solutes is a key process for life at high osmotic pressure and elevated salt concentrations. Most important solutes that can protect cell structures and metabolic processes at high salt concentrations are glycine betaine and ectoine. The genome analysis of more than 130 phototrophic bacteria shows that biosynthesis of glycine betaine is common among marine and halophilic phototrophic Proteobacteria and their chemotrophic relatives, as well as in representatives of Pirellulaceae and Actinobacteria, but are also found in halophilic Cyanobacteria and Chloroherpeton thalassium. This ability correlates well with the successful toleration of extreme salt concentrations. Freshwater bacteria in general lack the possibilities to synthesize and often also to take up these compounds. The biosynthesis of ectoine is found in the phylogenetic lines of phototrophic Alpha- and Gammaproteobacteria, most prominent in the Halorhodospira species and a number of Rhodobacteraceae. It is also common among Streptomycetes and Bacilli. The phylogeny of glycine-sarcosine methyltransferase (GMT) and diaminobutyrate-pyruvate aminotransferase (EctB) sequences correlate well with otherwise established phylogenetic groups. Most significantly, GMT sequences of cyanobacteria form two major phylogenetic branches and the branch of Halorhodospira species is distinct from all other Ectothiorhodospiraceae. A variety of transport systems for osmolytes are present in the studied bacteria.
Understanding the genetic architecture of complex traits is a major objective in biology. The standard approach for doing so is genome-wide association studies (GWAS), which aim to identify genetic polymorphisms responsible for variation in traits of interest. In human genetics, consistency across studies is commonly used as an indicator of reliability. However, if traits are involved in adaptation to the local environment, we do not necessarily expect reproducibility. On the contrary, results may depend on where you sample, and sampling across a wide range of environments may decrease the power of GWAS because of increased genetic heterogeneity. In this study, we examine how sampling affects GWAS in the model plant species Arabidopsis thaliana. We show that traits like flowering time are indeed influenced by distinct genetic effects in local populations. Furthermore, using gene expression as a molecular phenotype, we show that some genes are globally affected by shared variants, whereas others are affected by variants specific to subpopulations. Remarkably, the former are essentially all cis-regulated, whereas the latter are predominately affected by trans-acting variants. Our result illustrate that conclusions about genetic architecture can be extremely sensitive to sampling and population structure.