Sensitivity of asymmetric rate-dependent critical systems to initial conditions: Insights into cellular decision making

The work reported here aims to address the effects of time-dependent parameters and stochasticity on decision making in biological systems. We achieve this by extending previous studies that resorted to simple bifurcation normal forms, although in the present case we focus primarily on the issue of the system's sensitivity to initial conditions in the presence of two different noise distributions, Gaussian and Lévy. In addition, we also assess the impact of two-way sweeping at different rates through the critical region of a canonical Pitchfork bifurcation with a constant external asymmetry. The parallel with decision making in biocircuits is performed on this simple system since it is equivalent in its available states and dynamics to more complex genetic circuits published previously. Overall we verify that rate-dependent effects, previously reported as being important features of bifurcating systems, are specific to particular initial conditions. Processing of each starting state, which for the normal form underlying this study is akin to a classification task, is affected by the balance between sweeping speed through critical regions and the type of fluctuations added. For the heavy-tailed noise, two-way dynamic bifurcations are more efficient in processing the external signals, here understood to be jointly represented by the critical parameter profile and the external asymmetry amplitude, when compared to the system relying on escape dynamics. This is particular to the case when the system starts at an attractor not favored by the asymmetry and, in conjunction, when the sweeping amplitude is large

Speed-Dependent Cellular Decision Making in Nonequilibrium Genetic Circuits

Despite being governed by the principles of nonequilibrium transitions, gene expression dynamics underlying cell fate decision is poorly understood. In particular, the effect of signaling speed on cellular decision making is still unclear. Here we show that the decision between alternative cell fates, in a structurally symmetric circuit, can be biased depending on the speed at which the system is forced to go through the decision point. The circuit consists of two mutually inhibiting and self-activating genes, forced by two external signals with identical stationary values but different transient times. Under these conditions, slow passage through the decision point leads to a consistently biased decision due to the transient signaling asymmetry, whereas fast passage reduces and eventually eliminates the switch imbalance. The effect is robust to noise and shows that dynamic bifurcations, well known in nonequilibrium physics, are important for the control of genetic circuits

Transposable-Element Associated Small RNAs in Bombyx mori Genome

Small RNAs are a group of regulatory RNA molecules that control gene expression at transcriptional or post-transcriptional levels among eukaryotes. The silkworm, Bombyx mori L., genome harbors abundant repetitive sequences derived from families of retrotransposons and transposons, which together constitute almost half of the genome space and provide ample resource for biogenesis of the three major small RNA families. We systematically discovered transposable-element (TE)-associated small RNAs in B. mori genome based on a deep RNA-sequencing strategy and the effort yielded 182, 788 and 4,990 TE-associated small RNAs in the miRNA, siRNA and piRNA species, respectively. Our analysis suggested that the three small RNA species preferentially associate with different TEs to create sequence and functional diversity, and we also show evidence that a Bombyx non-LTR retrotransposon, bm1645, alone contributes to the generation of TE-associated small RNAs in a very significant way. The fact that bm1645-associated small RNAs partially overlap with each other implies a possibility that this element may be modulated by different mechanisms to generate different products with diverse functions. Taken together, these discoveries expand the small RNA pool in B. mori genome and lead to new knowledge on the diversity and functional significance of TE-associated small RNAs

Association between the cervicovaginal microbiome, BRCA1 mutation status, and risk of ovarian cancer: a case-control study

Summary Background Various factors—including age, family history, inflammation, reproductive factors, and tubal ligation—modulate the risk of ovarian cancer. In this study, our aim was to establish whether women with, or at risk of developing, ovarian cancer have an imbalanced cervicovaginal microbiome. Methods We did a case-control study in two sets of women aged 18–87 years in the Czech Republic, Germany, Italy, Norway, and the UK. The ovarian cancer set comprised women with epithelial ovarian cancer and controls (both healthy controls and those diagnosed with benign gynaecological conditions). The BRCA set comprised women with a BRCA1 mutation but without ovarian cancer and controls who were wild type for BRCA1 and BRCA2 (both healthy controls and those with benign gynaecological conditions). Cervicovaginal samples were gathered from all participants with the ThinPrep system and then underwent 16S rRNA gene sequencing. For each sample, we calculated the proportion of lactobacilli species (ie, Lactobacillus crispatus, Lactobacillus iners, Lactobacillus gasseri, and Lactobacillus jensenii), which are essential for the generation of a protective low vaginal pH, in the cervicovaginal microbiota. We grouped samples into those in which lactobacilli accounted for at least 50% of the species present (community type L) and those in which lactobacilli accounted for less than 50% of the species present (community type O). We assessed the adjusted association between BRCA1 status and ovarian cancer status and cervicovaginal microbiota community type, using a logistic regression model with a bias reduction method. Findings Participants were recruited between Jan 2, 2016, and July 21, 2018. The ovarian cancer set (n=360) comprised 176 women with epithelial ovarian cancer, 115 healthy controls and 69 controls with benign gynaecological conditions. The BRCA set (n=220) included 109 women with BRCA1 mutations, 97 healthy controls wild type for BRCA1 and BRCA2 and 14 controls with a benign gynaecological condition wild type for BRCA1 and BRCA2. On the basis of two-dimensional density plots, receiver–operating characteristic curve analysis, and age thresholds used previously, we divided the cohort into those younger than 50 years and those aged 50 years or older. In the ovarian cancer set, women aged 50 years or older had a higher prevalence of community type O microbiota (81 [61%] of 133 ovarian cancer cases and 84 [59%] of 142 healthy controls) than those younger than 50 years (23 [53%] of 43 cases and 12 [29%] of 42 controls). In the ovarian cancer set, women younger than 50 years with ovarian cancer had a significantly higher prevalence of community type O microbiota than did age-matched controls under a logistic regression model with bias correction (odds ratio [OR] 2·80 [95% CI 1·17–6·94]; p=0·020). In the BRCA set, women with BRCA1 mutations younger than 50 years were also more likely to have community type O microbiota than age-matched controls (OR 2·79 [95% CI 1·25–6·68]; p=0·012), after adjustment for pregnancy (ever). This risk was increased further if more than one first-degree family member was affected by any cancer (OR 5·26 [95% CI 1·83–15·30]; p=0·0022). In both sets, we noted that the younger the participants, the stronger the association between community type O microbiota and ovarian cancer or BRCA1 mutation status (eg, OR for community type O for cases aged <40 years in the ovarian cancer set 7·00 [95% CI 1·27–51·44], p=0·025; OR for community type O for BRCA1 mutation carriers aged <35 years in the BRCA set 4·40 [1·14–24·36], p=0·031). Interpretation The presence of ovarian cancer, or factors known to affect risk for the disease (ie, age and BRCA1 germline mutations), were significantly associated with having a community type O cervicovaginal microbiota. Whether re-instatement of a community type L microbiome by using, for example, vaginal suppositories containing live lactobacilli, would alter the microbiomial composition higher up in the female genital tract and in the fallopian tubes (the site of origin of high-grade serous ovarian cancer), and whether such changes could translate into a reduced incidence of ovarian cancer, needs to be investigated. Funding EU Horizon 2020 Research and Innovation Programme, EU Horizon 2020 European Research Council Programme, and The Eve Appeal

Characterization of binding specificities of bovine leucocyte class I molecules: impacts for rational epitope discovery.

The binding of peptides to classical major histocompatibility complex (MHC) class I proteins is the single most selective step in antigen presentation. However, the peptide-binding specificity of cattle MHC (bovine leucocyte antigen, BoLA) class I (BoLA-I) molecules remains poorly characterized. Here, we demonstrate how a combination of high-throughput assays using positional scanning combinatorial peptide libraries, peptide dissociation, and peptide-binding affinity binding measurements can be combined with bioinformatics to effectively characterize the functionality of BoLA-I molecules. Using this strategy, we characterized eight BoLA-I molecules, and found the peptide specificity to resemble that of human MHC-I molecules with primary anchors most often at P2 and P9, and occasional auxiliary P1/P3/P5/P6 anchors. We analyzed nine reported CTL epitopes from Theileria parva, and in eight cases, stable and high affinity binding was confirmed. A set of peptides were tested for binding affinity to the eight BoLA proteins and used to refine the predictors of peptide–MHC binding NetMHC and NetMHCpan. The inclusion of BoLA-specific peptide-binding data led to a significant improvement in prediction accuracy for reported T. parva CTL epitopes. For reported CTL epitopes with weak or no predicted binding, these refined prediction methods suggested presence of nested minimal epitopes with high-predicted binding affinity. The enhanced affinity of the alternative peptides was in all cases confirmed experimentally. This study demonstrates how biochemical high-throughput assays combined with immunoinformatics can be used to characterize the peptide-binding motifs of BoLA-I molecules, boosting performance of MHC peptide-binding prediction methods, and empowering rational epitope discovery in cattle.Fil: Hansen, Andreas M.. Universidad de Copenhagen; DinamarcaFil: Rasmussen, Michael. Universidad de Copenhagen; DinamarcaFil: Svitek,Nicholas. International Livestock Research Institute; KeniaFil: Harndahl, Mikkel. Universidad de Copenhagen; DinamarcaFil: Golde, William T.. United States Department of Agriculture. Agriculture Research Service.; Estados UnidosFil: Barlow, John. University of Vermont; Estados UnidosFil: Vishvanath, Nene. International Livestock Research Institute; KeniaFil: Buus, Søren. Universidad de Copenhagen; DinamarcaFil: Nielsen, Morten. Technical University Of Denmark; Dinamarca. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Investigaciones Biotecnológicas. Universidad Nacional de San Martín. Instituto de Investigaciones Biotecnológicas; Argentin