Repressive Gene Regulation Synchronizes Development with Cellular Metabolism.

Metabolic conditions affect the developmental tempo of animals. Developmental gene regulatory networks (GRNs) must therefore synchronize their dynamics with a variable timescale. We find that layered repression of genes couples GRN output with variable metabolism. When repressors of transcription or mRNA and protein stability are lost, fewer errors in Drosophila development occur when metabolism is lowered. We demonstrate the universality of this phenomenon by eliminating the entire microRNA family of repressors and find that development to maturity can be largely rescued when metabolism is reduced. Using a mathematical model that replicates GRN dynamics, we find that lowering metabolism suppresses the emergence of developmental errors by curtailing the influence of auxiliary repressors on GRN output. We experimentally show that gene expression dynamics are less affected by loss of repressors when metabolism is reduced. Thus, layered repression provides robustness through error suppression and may provide an evolutionary route to a shorter reproductive cycle.

Energy metabolism modulates the regulatory impact of activators on gene expression.

Gene expression is a regulated process fueled by ATP consumption. Therefore, regulation must be coupled to constraints imposed by the level of energy metabolism. Here, we explore this relationship both theoretically and experimentally. A stylized mathematical model predicts that activators of gene expression have variable impact depending on metabolic rate. Activators become less essential when metabolic rate is reduced and more essential when metabolic rate is enhanced. We find that, in the Drosophila eye, expression dynamics of the yan gene are less affected by loss of EGFR-mediated activation when metabolism is reduced, and the opposite effect is seen when metabolism is enhanced. The effects are also seen at the level of pattern regularity in the adult eye, where loss of EGFR-mediated activation is mitigated by lower metabolism. We propose that gene activation is tuned by energy metabolism to allow for faithful expression dynamics in the face of variable metabolic conditions.

Ratiometric sensing of Pnt and Yan transcription factor levels confers ultrasensitivity to photoreceptor fate transitions in Drosophila.

Cell state transitions are often triggered by large changes in the concentrations of transcription factors and therefore large differences in their stoichiometric ratios. Whether cells can elicit transitions using modest changes in the ratios of co-expressed factors is unclear. Here, we investigate how cells in the Drosophila eye resolve state transitions by quantifying the expression dynamics of the ETS transcription factors Pnt and Yan. Eye progenitor cells maintain a relatively constant ratio of Pnt/Yan protein, despite expressing both proteins with pulsatile dynamics. A rapid and sustained twofold increase in the Pnt/Yan ratio accompanies transitions to photoreceptor fates. Genetic perturbations that modestly disrupt the Pnt/Yan ratio produce fate transition defects consistent with the hypothesis that transitions are normally driven by a twofold shift in the ratio. A biophysical model based on cooperative Yan-DNA binding coupled with non-cooperative Pnt-DNA binding illustrates how twofold ratio changes could generate ultrasensitive changes in target gene transcription to drive fate transitions. Thus, coupling cell state transitions to the Pnt/Yan ratio sensitizes the system to modest fold-changes, conferring robustness and ultrasensitivity to the developmental program.

The characteristics of early-stage research into human genes are substantially different from subsequent research.

Throughout the last 2 decades, several scholars observed that present day research into human genes rarely turns toward genes that had not already been extensively investigated in the past. Guided by hypotheses derived from studies of science and innovation, we present here a literature-wide data-driven meta-analysis to identify the specific scientific and organizational contexts that coincided with early-stage research into human genes throughout the past half century. We demonstrate that early-stage research into human genes differs in team size, citation impact, funding mechanisms, and publication outlet, but that generalized insights derived from studies of science and innovation only partially apply to early-stage research into human genes. Further, we demonstrate that, presently, genome biology accounts for most of the initial early-stage research, while subsequent early-stage research can engage other life sciences fields. We therefore anticipate that the specificity of our findings will enable scientists and policymakers to better promote early-stage research into human genes and increase overall innovation within the life sciences.

A new approach for extracting information from protein dynamics.

Increased ability to predict protein structures is moving research focus towards understanding protein dynamics. A promising approach is to represent protein dynamics through networks and take advantage of well-developed methods from network science. Most studies build protein dynamics networks from correlation measures, an approach that only works under very specific conditions, instead of the more robust inverse approach. Thus, we apply the inverse approach to the dynamics of protein dihedral angles, a system of internal coordinates, to avoid structural alignment. Using the well-characterized adhesion protein, FimH, we show that our method identifies networks that are physically interpretable, robust, and relevant to the allosteric pathway sites. We further use our approach to detect dynamical differences, despite structural similarity, for Siglec-8 in the immune system, and the SARS-CoV-2 spike protein. Our study demonstrates that using the inverse approach to extract a network from protein dynamics yields important biophysical insights.

The first step is recognizing there is a problem: a methodology for adjusting for variability in disease severity when estimating clinician performance.

BACKGROUND: Adoption of innovations in the field of medicine is frequently hindered by a failure to recognize the condition targeted by the innovation. This is particularly true in cases where recognition requires integration of patient information from different sources, or where disease presentation can be heterogeneous and the recognition step may be easier for some patients than for others. METHODS: We propose a general data-driven metric for clinician recognition that accounts for the variability in patient disease severity and for institutional standards. As a case study, we evaluate the ventilatory management of 362 patients with acute respiratory distress syndrome (ARDS) at a large academic hospital, because clinician recognition of ARDS has been identified as a major barrier to adoption to evidence-based ventilatory management. We calculate our metric for the 48 critical care physicians caring for these patients and examine the relationships between differences in ARDS recognition performance from overall institutional levels and provider characteristics such as demographics, social network position, and self-reported barriers and opinions. RESULTS: Our metric was found to be robust to patient characteristics previously demonstrated to affect ARDS recognition, such as disease severity and patient height. Training background was the only factor in this study that showed an association with physician recognition. Pulmonary and critical care medicine (PCCM) training was associated with higher recognition (ß = 0.63, 95% confidence interval 0.46-0.80, p < 7 × 10- 5). Non-PCCM physicians recognized ARDS cases less frequently and expressed greater satisfaction with the ability to get the information needed for making an ARDS diagnosis (p < 5 × 10- 4), suggesting that lower performing clinicians may be less aware of institutional barriers. CONCLUSIONS: We present a data-driven metric of clinician disease recognition that accounts for variability in patient disease severity and for institutional standards. Using this metric, we identify two unique physician populations with different intervention needs. One population consistently recognizes ARDS and reports barriers vs one does not and reports fewer barriers.

Conformational stability of the bacterial adhesin, FimH, with an inactivating mutation.

Allostery governing two conformational states is one of the proposed mechanisms for catch-bond behavior in adhesive proteins. In FimH, a catch-bond protein expressed by pathogenic bacteria, separation of two domains disrupts inhibition by the pilin domain. Thus, tensile force can induce a conformational change in the lectin domain, from an inactive state to an active state with high affinity. To better understand allosteric inhibition in two-domain FimH (H2 inactive), we use molecular dynamics simulations to study the lectin domain alone, which has high affinity (HL active), and also the lectin domain stabilized in the low-affinity conformation by an Arg-60-Pro mutation (HL mutant). Because ligand-binding induces an allostery-like conformational change in HL mutant, this more experimentally tractable version has been proposed as a "minimal model" for FimH. We find that HL mutant has larger backbone fluctuations than both H2 inactive and HL active, at the binding pocket and allosteric interdomain region. We use an internal coordinate system of dihedral angles to identify protein regions with differences in backbone and side chain dynamics beyond the putative allosteric pathway sites. By characterizing HL mutant dynamics for the first time, we provide additional insight into the transmission of allosteric information across the lectin domain and build upon structural and thermodynamic data in the literature to further support the use of HL mutant as a "minimal model." Understanding how to alter protein dynamics to prevent the allosteric conformational change may guide drug development to prevent infection by blocking FimH adhesion.

Large-scale investigation of the reasons why potentially important genes are ignored.

Biomedical research has been previously reported to primarily focus on a minority of all known genes. Here, we demonstrate that these differences in attention can be explained, to a large extent, exclusively from a small set of identifiable chemical, physical, and biological properties of genes. Together with knowledge about homologous genes from model organisms, these features allow us to accurately predict the number of publications on individual human genes, the year of their first report, the levels of funding awarded by the National Institutes of Health (NIH), and the development of drugs against disease-associated genes. By explicitly identifying the reasons for gene-specific bias and performing a meta-analysis of existing computational and experimental knowledge bases, we describe gene-specific strategies for the identification of important but hitherto ignored genes that can open novel directions for future investigation.

Reply to "Far away from the lamppost".

In this Formal Comment, the authors of the recent publication "Large-scale investigation of the reasons why potentially important genes are ignored" maintain that it can be read as an opportunity to explore the unknown.

Fly-QMA: Automated analysis of mosaic imaginal discs in Drosophila.

Mosaic analysis provides a means to probe developmental processes in situ by generating loss-of-function mutants within otherwise wildtype tissues. Combining these techniques with quantitative microscopy enables researchers to rigorously compare RNA or protein expression across the resultant clones. However, visual inspection of mosaic tissues remains common in the literature because quantification demands considerable labor and computational expertise. Practitioners must segment cell membranes or cell nuclei from a tissue and annotate the clones before their data are suitable for analysis. Here, we introduce Fly-QMA, a computational framework that automates each of these tasks for confocal microscopy images of Drosophila imaginal discs. The framework includes an unsupervised annotation algorithm that incorporates spatial context to inform the genetic identity of each cell. We use a combination of real and synthetic validation data to survey the performance of the annotation algorithm across a broad range of conditions. By contributing our framework to the open-source software ecosystem, we aim to contribute to the current move toward automated quantitative analysis among developmental biologists.

Multidimensional Assessment of the Host Response in Mechanically Ventilated Patients with Suspected Pneumonia.

Rationale: The identification of informative elements of the host response to infection may improve the diagnosis and management of bacterial pneumonia. Objectives: To determine whether the absence of alveolar neutrophilia can exclude bacterial pneumonia in critically ill patients with suspected infection and to test whether signatures of bacterial pneumonia can be identified in the alveolar macrophage transcriptome. Methods: We determined the test characteristics of alveolar neutrophilia for the diagnosis of bacterial pneumonia in three cohorts of mechanically ventilated patients. In one cohort, we also isolated macrophages from alveolar lavage fluid and used the transcriptome to identify signatures of bacterial pneumonia. Finally, we developed a humanized mouse model of Pseudomonas aeruginosa pneumonia to determine if pathogen-specific signatures can be identified in human alveolar macrophages. Measurements and Main Results: An alveolar neutrophil percentage less than 50% had a negative predictive value of greater than 90% for bacterial pneumonia in both the retrospective (n = 851) and validation cohorts (n = 76 and n = 79). A transcriptional signature of bacterial pneumonia was present in both resident and recruited macrophages. Gene signatures from both cell types identified patients with bacterial pneumonia with test characteristics similar to alveolar neutrophilia. Conclusions: The absence of alveolar neutrophilia has a high negative predictive value for bacterial pneumonia in critically ill patients with suspected infection. Macrophages can be isolated from alveolar lavage fluid obtained during routine care and used for RNA-Seq analysis. This novel approach may facilitate a longitudinal and multidimensional assessment of the host response to bacterial pneumonia.

Single-Cell Transcriptomic Analysis of Human Lung Provides Insights into the Pathobiology of Pulmonary Fibrosis.

Rationale: The contributions of diverse cell populations in the human lung to pulmonary fibrosis pathogenesis are poorly understood. Single-cell RNA sequencing can reveal changes within individual cell populations during pulmonary fibrosis that are important for disease pathogenesis. Objectives: To determine whether single-cell RNA sequencing can reveal disease-related heterogeneity within alveolar macrophages, epithelial cells, or other cell types in lung tissue from subjects with pulmonary fibrosis compared with control subjects. Methods: We performed single-cell RNA sequencing on lung tissue obtained from eight transplant donors and eight recipients with pulmonary fibrosis and on one bronchoscopic cryobiospy sample from a patient with idiopathic pulmonary fibrosis. We validated these data using in situ RNA hybridization, immunohistochemistry, and bulk RNA-sequencing on flow-sorted cells from 22 additional subjects. Measurements and Main Results: We identified a distinct, novel population of profibrotic alveolar macrophages exclusively in patients with fibrosis. Within epithelial cells, the expression of genes involved in Wnt secretion and response was restricted to nonoverlapping cells. We identified rare cell populations including airway stem cells and senescent cells emerging during pulmonary fibrosis. We developed a web-based tool to explore these data. Conclusions: We generated a single-cell atlas of pulmonary fibrosis. Using this atlas, we demonstrated heterogeneity within alveolar macrophages and epithelial cells from subjects with pulmonary fibrosis. These results support the feasibility of discovery-based approaches using next-generation sequencing technologies to identify signaling pathways for targeting in the development of personalized therapies for patients with pulmonary fibrosis.

Within-Gene Shine-Dalgarno Sequences Are Not Selected for Function.

The Shine-Dalgarno (SD) sequence motif facilitates translation initiation and is frequently found upstream of bacterial start codons. However, thousands of instances of this motif occur throughout the middle of protein coding genes in a typical bacterial genome. Here, we use comparative evolutionary analysis to test whether SD sequences located within genes are functionally constrained. We measure the conservation of SD sequences across Enterobacteriales, and find that they are significantly less conserved than expected. Further, the strongest SD sequences are the least conserved whereas we find evidence of conservation for the weakest possible SD sequences given amino acid constraints. Our findings indicate that most SD sequences within genes are likely to be deleterious and removed via selection. To illustrate the origin of these deleterious costs, we show that ATG start codons are significantly depleted downstream of SD sequences within genes, highlighting the constraint that these sequences impose on the surrounding nucleotides to minimize the potential for erroneous translation initiation.

Diversity of Translation Initiation Mechanisms across Bacterial Species Is Driven by Environmental Conditions and Growth Demands.

The Shine-Dalgarno (SD) sequence motif is frequently found upstream of protein coding genes and is thought to be the dominant mechanism of translation initiation used by bacteria. Experimental studies have shown that the SD sequence facilitates start codon recognition and enhances translation initiation by directly interacting with the highly conserved anti-SD sequence on the 30S ribosomal subunit. However, the proportion of SD-led genes within a genome varies across species and the factors governing this variation in translation initiation mechanisms remain largely unknown. Here, we conduct a phylogenetically informed analysis and find that species capable of rapid growth contain a higher proportion of SD-led genes throughout their genomes. We show that SD sequence utilization covaries with a suite of genomic features that are important for efficient translation initiation and elongation. In addition to these endogenous genomic factors, we further show that exogenous environmental factors may influence the evolution of translation initiation mechanisms by finding that thermophilic species contain significantly more SD-led genes than mesophiles. Our results demonstrate that variation in translation initiation mechanisms across bacterial species is predictable and is a consequence of differential life-history strategies related to maximum growth rate and environmental-specific constraints.

Differences in Collaboration Patterns across Discipline, Career Stage, and Gender.

Collaboration plays an increasingly important role in promoting research productivity and impact. What remains unclear is whether female and male researchers in science, technology, engineering, and mathematical (STEM) disciplines differ in their collaboration propensity. Here, we report on an empirical analysis of the complete publication records of 3,980 faculty members in six STEM disciplines at select U.S. research universities. We find that female faculty have significantly fewer distinct co-authors over their careers than males, but that this difference can be fully accounted for by females' lower publication rate and shorter career lengths. Next, we find that female scientists have a lower probability of repeating previous co-authors than males, an intriguing result because prior research shows that teams involving new collaborations produce work with higher impact. Finally, we find evidence for gender segregation in some sub-disciplines in molecular biology, in particular in genomics where we find female faculty to be clearly under-represented.

Cross-evaluation of metrics to estimate the significance of creative works.

In a world overflowing with creative works, it is useful to be able to filter out the unimportant works so that the significant ones can be identified and thereby absorbed. An automated method could provide an objective approach for evaluating the significance of works on a universal scale. However, there have been few attempts at creating such a measure, and there are few "ground truths" for validating the effectiveness of potential metrics for significance. For movies, the US Library of Congress's National Film Registry (NFR) contains American films that are "culturally, historically, or aesthetically significant" as chosen through a careful evaluation and deliberation process. By analyzing a network of citations between 15,425 United States-produced films procured from the Internet Movie Database (IMDb), we obtain several automated metrics for significance. The best of these metrics is able to indicate a film's presence in the NFR at least as well or better than metrics based on aggregated expert opinions or large population surveys. Importantly, automated metrics can easily be applied to older films for which no other rating may be available. Our results may have implications for the evaluation of other creative works such as scientific research.

NullSeq: A Tool for Generating Random Coding Sequences with Desired Amino Acid and GC Contents.

The existence of over- and under-represented sequence motifs in genomes provides evidence of selective evolutionary pressures on biological mechanisms such as transcription, translation, ligand-substrate binding, and host immunity. In order to accurately identify motifs and other genome-scale patterns of interest, it is essential to be able to generate accurate null models that are appropriate for the sequences under study. While many tools have been developed to create random nucleotide sequences, protein coding sequences are subject to a unique set of constraints that complicates the process of generating appropriate null models. There are currently no tools available that allow users to create random coding sequences with specified amino acid composition and GC content for the purpose of hypothesis testing. Using the principle of maximum entropy, we developed a method that generates unbiased random sequences with pre-specified amino acid and GC content, which we have developed into a python package. Our method is the simplest way to obtain maximally unbiased random sequences that are subject to GC usage and primary amino acid sequence constraints. Furthermore, this approach can easily be expanded to create unbiased random sequences that incorporate more complicated constraints such as individual nucleotide usage or even di-nucleotide frequencies. The ability to generate correctly specified null models will allow researchers to accurately identify sequence motifs which will lead to a better understanding of biological processes as well as more effective engineering of biological systems.

Impact of heterogeneity and socioeconomic factors on individual behavior in decentralized sharing ecosystems.

Tens of millions of individuals around the world use decentralized content distribution systems, a fact of growing social, economic, and technological importance. These sharing systems are poorly understood because, unlike in other technosocial systems, it is difficult to gather large-scale data about user behavior. Here, we investigate user activity patterns and the socioeconomic factors that could explain the behavior. Our analysis reveals that (i) the ecosystem is heterogeneous at several levels: content types are heterogeneous, users specialize in a few content types, and countries are heterogeneous in user profiles; and (ii) there is a strong correlation between socioeconomic indicators of a country and users behavior. Our findings open a research area on the dynamics of decentralized sharing ecosystems and the socioeconomic factors affecting them, and may have implications for the design of algorithms and for policymaking.

The role of mentorship in protégé performance.

The role of mentorship in protégé performance is a matter of importance to academic, business and governmental organizations. Although the benefits of mentorship for protégés, mentors and their organizations are apparent, the extent to which protégés mimic their mentors' career choices and acquire their mentorship skills is unclear. The importance of a science, technology, engineering and mathematics workforce to economic growth and the role of effective mentorship in maintaining a 'healthy' such workforce demand the study of the role of mentorship in academia. Here we investigate one aspect of mentor emulation by studying mentorship fecundity-the number of protégés a mentor trains-using data from the Mathematics Genealogy Project, which tracks the mentorship record of thousands of mathematicians over several centuries. We demonstrate that fecundity among academic mathematicians is correlated with other measures of academic success. We also find that the average fecundity of mentors remains stable over 60 years of recorded mentorship. We further discover three significant correlations in mentorship fecundity. First, mentors with low mentorship fecundities train protégés that go on to have mentorship fecundities 37% higher than expected. Second, in the first third of their careers, mentors with high fecundities train protégés that go on to have fecundities 29% higher than expected. Finally, in the last third of their careers, mentors with high fecundities train protégés that go on to have fecundities 31% lower than expected.

Quantifying position-dependent codon usage bias.

Although the mapping of codon to amino acid is conserved across nearly all species, the frequency at which synonymous codons are used varies both between organisms and between genes from the same organism. This variation affects diverse cellular processes including protein expression, regulation, and folding. Here, we mathematically model an additional layer of complexity and show that individual codon usage biases follow a position-dependent exponential decay model with unique parameter fits for each codon. We use this methodology to perform an in-depth analysis on codon usage bias in the model organism Escherichia coli. Our methodology shows that lowly and highly expressed genes are more similar in their codon usage patterns in the 5'-gene regions, but that these preferences diverge at distal sites resulting in greater positional dependency (pD, which we mathematically define later) for highly expressed genes. We show that position-dependent codon usage bias is partially explained by the structural requirements of mRNAs that results in increased usage of A/T rich codons shortly after the gene start. However, we also show that the pD of 4- and 6-fold degenerate codons is partially related to the gene copy number of cognate-tRNAs supporting existing hypotheses that posit benefits to a region of slow translation in the beginning of coding sequences. Lastly, we demonstrate that viewing codon usage bias through a position-dependent framework has practical utility by improving accuracy of gene expression prediction when incorporating positional dependencies into the Codon Adaptation Index model.