The cyclic dinucleotide 2'3'-cGAMP induces a broad anti-bacterial and anti-viral response in the sea anemone Nematostella vectensis In mammals, cyclic dinucleotides (CDNs) bind and activate STING to initiate an anti-viral type I interferon response. CDNs and STING originated in bacteria and are present in most animals. By contrast, interferons are believed to have emerged in vertebrates; thus, the function of CDN signaling in invertebrates is unclear. Here, we use a CDN, 23-cGAMP, to activate immune responses in a model cnidarian invertebrate, the starlet sea anemone Nematostella vectensis. Using RNA-Seq, we found that 23-cGAMP induces robust transcription of both anti-viral and anti-bacterial genes in N. vectensis. Many of the anti-viral genes induced by 23-cGAMP are homologs of vertebrate interferon stimulated genes, implying that the interferon response predates the evolution of interferons. Knockdown experiments identified a role for NF-{kappa}B in specifically inducing anti-bacterial genes downstream of 23-cGAMP. Some of these putative anti-bacterial genes were also found to be induced during Pseudomonas aeruginosa infection. We characterized the protein product of one of the putative anti-bacterial genes, the N. vectensis homolog of Dae4, and found that it has conserved anti-bacterial activity. This work suggests that a broad anti-bacterial and anti-viral transcriptional response is an evolutionarily ancestral output of 23-cGAMP signaling in animals.

Significance statementCyclic dinucleotides are signaling molecules that originated in bacteria and were subsequently acquired and co-opted by animals for immune signaling. The major cyclic dinucleotide signaling pathway in mammals results in the production of anti-viral molecules called interferons. Invertebrates such as sea anemones lack interferons, and thus it was unclear whether cyclic dinucleotide signaling would play a role in immunity in these animals. Here we report that in the anemone Nematostella vectensis, cyclic dinucleotides activate both anti-viral and anti-bacterial immune responses, and do so through a conserved pathway. These results provide insights into the evolutionary origins of innate immunity, and suggest a broader ancestral role for cyclic dinucleotide signaling that evolved toward more specialized anti-viral functions in mammals.
Metabolic regulation of ILC2 differentiation into ILC1-like cells during Mycobacterium tuberculosis infection Tissue-resident innate lymphoid cells (ILCs) regulate tissue homeostasis, protect against pathogens at mucosal surfaces and are key players at the interface of innate and adaptive immunity. How ILCs adapt their phenotype and function to environmental cues within tissues remains to be fully understood. Here, we show that Mycobacterium tuberculosis infection alters the phenotype and function of immature lung ILC2 toward a protective interferon-{gamma}-producing ILC1-like population. This differentiation is controlled by type 1 cytokines and is associated with a glycolytic program involving the transcription factor HIF1. Collectively, our data reveal how tissue-resident ILCs adapt to type 1 inflammation toward a pathogen tailored immune response.
DNA supercoiling-mediated collective behavior of co-transcribing RNA polymerases Multiple RNA polymerases (RNAPs) transcribing a gene have been known to exhibit collective group behavior, causing the transcription elongation rate to increase with the rate of transcription initiation. Such behavior has long been believed to be driven by a physical interaction or "push" between closely spaced RNAPs. However, recent studies have posited that RNAPs separated by longer distances may cooperate via the DNA segment under transcription. Here, we present a theoretical model incorporating the mechanical coupling between RNAP translocation and the torsional response of supercoiled DNA. Using stochastic simulations, we demonstrate long-range cooperation between co-transcribing RNAPs mediated by DNA supercoiling. We find that inhibiting transcription initiation can slow down the already recruited RNAPs, in agreement with recent experimental observations, and predict that the average transcription elongation rate varies non-monotonically with the rate of transcription initiation. We further show that while RNAPs transcribing neighboring genes oriented in tandem can cooperate, those transcribing genes in divergent or convergent orientations can act antagonistically, and that such behavior holds over a large range of intergenic separations. Our model makes testable predictions, revealing how the mechanical interplay between RNAPs and the DNA they transcribe can govern a key cellular process. Bamboos flower after the return of almost the same sun-moon phasing as at seedling emergence All Asian woody bamboo species of economic importance are semelparous. They remain vegetative during time intervals that are specific to each species and range from three to 120 years, with notable concentrations around a series of values (3, 7-8, 14-17, 29-36, 42-48, 61-64, and 120 years). Then, they flower gregariously within a short period. Entire forests temporarily disappear during these periods, and the dates of these dramatic events have been recorded over the last 200 years.

While searching for a correlated environmental cue, I have found that the concentrations of flowering cycles were highly correlated with the series of successive returns of almost the same sun-moon phasing as at seedling emergence. On basis of knowledge on plant photoperiod sensitivity, I hypothesize that bamboo plants i) run a lunar cellular clock that is set at the full moon, ii) retain in their cellular memory the exact sun-moon phasing of the year of their emergence as seedlings, and iii) inhibit flowering until the occurrence of a unique, species-specific sun-moon phasing that is shifted by a precise amount from the sun-moon phasing at their emergence. Recent evidence of plant responses to lunar cycles supports this hypothesis, for which experimental evidence is now anticipated.

HighlightThe long-standing enigma of the triggering of the gregarious flowering of bamboos is possibly solved. Flowering would be inhibited until the occurrence of the species-specific sun-moon phasing. SERCA function is impaired in skeletal and cardiac muscles from young DBA/2J mdx mice The DBA/2J (D2) mdx mouse has emerged as a more severe model of Duchenne muscular dystrophy when compared to the traditional C57BL/10 (C57) mdx mouse. Here, we questioned whether sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) function would differ in muscles from young D2 and C57 mdx mice. In gastrocnemius muscles, both D2- and C57 mdx mice exhibited signs of impaired Ca2+ uptake, however, this was more severe in D2 mdx mice. Maximal SERCA activity was lowered only in D2 mdx gastrocnemius muscles and not C57 mdx muscles. Furthermore, in the left ventricle and diaphragm, Ca2+ uptake was impaired in C57 mdx muscles with lowered rates of Ca2+ uptake compared with C57 WT mice, whereas in muscles from D2 mdx mice, rates of Ca2+ uptake were unattainable due to the severe impairments in their ability to transport Ca2+. Overall, our study demonstrates that SERCA function is drastically impaired in young D2 mdx mice. Immunomodulatory Contribution of Mast Cells to the Regenerative Biomaterial Microenvironment Bioactive immunomodulatory biomaterials have shown promise for influencing the immune response to promote tissue repair and regeneration. Macrophages and T cells have been associated with this response; however, other immune cell types have been traditionally overlooked. In this study, we investigated the role of mast cells in the regulation of the immune response to decellularized biomaterial scaffolds using a subcutaneous implant model. In mast cell-deficient mice, there was dysregulation of the expected M1 to M2 macrophage transition typically induced by the biomaterial scaffold. Polarization progression deviated in a sex specific manner with an early transition to an M2 profile in female mice, while the male response was unable to properly transition past a pro-inflammatory M1 state. Both were reversed with adoptive mast cell transfer. Further investigation of the later stage immune response in male mice determined a greater sustained pro-inflammatory gene expression profile including the IL-1 cytokine family, IL-6, alarmins, and chemokines. These results highlight mast cells as another important cell type that influences the immune response to pro-regenerative biomaterials. Inhibition of BCMA signaling alters protein translation and Multiple Myeloma growth Recent T Cell therapies have been effective in the treatment of hematological cancers. However, immunotoxicity and treatment relapse pose significant clinical challenges. Here, we revealed distinctive requirement for neutralizing TNF receptor ligands APRIL and BAFF in MM and DLBCL. Furthermore, we investigated the use of BCMA decoy receptor (sBCMA-Fc) as a therapeutic inhibitor of ARPIL and BAFF. While wild-type sBCMA-Fc successfully blocked APRIL signaling with picomolar binding affinity, inhibiting tumor growth in MM models, it lacked efficacy in inhibiting DLBCL progression due to its weak binding for BAFF. To expand the therapeutic utility of sBCMA-Fc, using a protein engineering approach, we generated an affinity enhanced mutant sBCMA-Fc fusion molecule (sBCMA-Fc V3) with 4-folds and 500-folds enhancement in binding to APRIL and BAFF respectively. The sBCMA-Fc V3 clone significantly enhanced antitumor activity against both MM and DLBCL. Importantly, sBCMA-Fc V3 was proven to be a viable clinical candidate by showing adequate toxicity profile and on-target mechanism of action in nonhuman primates.

SUMMARYThis study demonstrates the dichotomous function of APRIL and BAFF in MM and DLBCL, that can be safely targeted by an engineered fusion protein designed to trap APRIL and BAFF with ultra-high binding affinity. Bespoke data augmentation and network construction enable developmental morphological classification on limited microscopy datasets Deep learning is a powerful tool for image classification, yet training deep learning models often requires large datasets, which are not always available in biological research. Here we use a small dataset to train an accurate classifier by employing data augmentation regimes and Bayesian optimisation. We establish a new staging system for the Hamburger-Hamilton stage 10 chick embryo brain, based on morphology, identifying three sub-stages uncoupled from somite number, the usual metric used to stage chick embryos. Incorporating biologically informed, data-driven preprocessing, we train a deep convolutional neural network to classify our sub-stages. We find that augmenting our images with a combination of rotation, blur, shear, and cutout transformations is effective in training an accurate classifier on a small microscopy dataset. We then apply our approach to an unrelated small microscopy dataset comprising images of limb buds, finding good agreement with the results from brain classification. Through saliency analysis, we ensure that relevant features are being used by our classifier, gain insight into the efficacy of different data augmentations, and identify sub-stage specific features. In summary, we apply biological domain expertise to design effective data augmentation regimes for training a neural network-based image classifier on small, specialised datasets, classifying the developing brains and limbs with up to 90.9% and 94.4% accuracy respectively. These reliable classifiers can inform future experiments. Allosteric activation of SARS-CoV-2 RdRp by remdesivir triphosphate and other phosphorylated nucleotides The catalytic subunit of SARS-CoV-2 RNA-dependent RNA polymerase (RdRp), Nsp12, has a unique NiRAN domain that transfers nucleoside monophosphates to the Nsp9 protein. The NiRAN and RdRp modules form a dynamic interface distant from their catalytic sites and both activities are essential for viral replication. We report that codon-optimized (for the pause-free translation) Nsp12 exists in inactive state in which NiRAN/RdRp interactions are broken, whereas translation by slow ribosomes and incubation with accessory Nsp7/8 subunits or NTPs partially rescue RdRp activity. Our data show that adenosine and remdesivir triphosphates promote synthesis of A-less RNAs, as does ppGpp, while amino acid substitutions at the NiRAN/RdRp interface augment activation, suggesting that ligand binding to the NiRAN catalytic site modulates RdRp activity. The existence of allosterically-linked nucleotidyl transferase sites that utilize the same substrates has important implications for understanding the mechanism of SARS-CoV-2 replication and design of its inhibitors.

HighlightsO_LICodon-optimization of Nsp12 triggers misfolding and activity loss
C_LIO_LISlow translation, accessory Nsp7 and Nsp8 subunits, and NTPs rescue Nsp12
C_LIO_LINon-substrate nucleotides activate RNA chain synthesis, likely via NiRAN domain
C_LIO_LICrosstalk between two Nsp12 active sites that bind the same ligands
C_LI Adolescent development of multiscale structural wiring and functional interactions in the human connectome AO_SCPLOWBSTRACTC_SCPLOWAdolescence is a time of profound changes in the physical wiring and function of the brain. Here, we analyzed structural and functional brain network development in an accelerated longitudinal cohort spanning 14-25 years (n = 199). Core to our work was an advanced in vivo model of cortical wiring incorporating MRI features of cortico-cortical proximity, microstructural similarity, and white matter tractography. Longitudinal analyses assessing age-related changes in cortical wiring identified a continued differentiation of multiple cortico-cortical structural networks in youth. Studying resting-state functional MRI measures in the same participants at baseline, we found that regions with more similar structural wiring were more likely to be functionally coupled. Moreover, longitudinal structural wiring changes, particularly between sensory/unimodal and default mode networks, were reflected in tendencies for increased differentiation in brain function. These longitudinal findings provide new insights into adolescent development of human brain structure and function, illustrating how structural wiring interacts with the maturation of macroscale functional hierarchies. A novel 5 prime DNA binding site in RFC facilitates PCNA loading for gap DNA repair RFC uses ATP to assemble PCNA onto primed sites for replicative DNA polymerases {delta} and {gamma}. The RFC pentamer forms a central chamber that binds 3' ss/ds DNA junctions to load PCNA onto DNA during replication. We show here five structures revealing a novel 2nd DNA binding site in RFC that binds a 5' duplex. This 5' DNA site is located in the N-terminal BRCT domain and AAA+ module of the large Rfc1 subunit. It appears 5' DNA binds after 3' DNA binding to RFC. Our structures reveal ideal binding to a 7-nt gap, which includes 2 bp unwound by the clamp loader. Biochemical studies show enhanced binding to 5 and 10 nt gaps, consistent with the structural results. We propose that the 5' site facilitates RFCs PCNA loading activity at a DNA damage-induced gap to recruit gap-filling polymerases. These findings are consistent with genetic studies showing that base excision repair of gaps greater than 1 base requires PCNA and involves the 5' DNA binding domain of Rfc1. We further observe that a 5' end facilitates PCNA loading at an RPA coated 50-nt gap, suggesting a potential role of the RFC 5' DNA site in longer gaps including during lagging strand DNA synthesis. Diurnal activity budget and conservation implication of Hippopotamus amphibius at Bui National Park, Ghana Animals apportion time for their daily behavioral activities. We studied the activity budget of Hippopotamus amphibius at the Black Volta River in the Bui National Park, Ghana. We performed instantaneous scan sample using ground count survey between August 2020 to July 2021. We observed that H. amphibius spent most of their day time resting (54.75 %), followed by feeding (22.93 %), walking (19.2 %) and touching (3.12 %). We found significant difference in the time spent between the daily activities (H = 41.67; p < 0.0001). Except activities involving touching, we found no significant seasonal differences in feeding, resting and walking activities by H. amphibius. Understanding the behavioral activities of H. amphibius will enhance management and conservation of the animals. Gene loss during a transition to multicellularity Multicellular evolution is a major transition associated with momentous diversification of multiple lineages and increased developmental complexity. The volvocine algae comprise a valuable system for the study of this transition, as they span from unicellular to undifferentiated and differentiated multicellular morphologies despite their genomes being highly similar, suggesting multicellular evolution requires few genetic changes to undergo dramatic shifts in developmental complexity. Here, the evolutionary dynamics of six volvocine genomes were examined, where a gradual loss of genes was observed in parallel to the co-option of a few key genes. Protein complexes in the six species exhibited a high degree of novel interactions, suggesting that gene loss plays a role in evolutionary novelty. This finding was supported by gene network modeling, where gene loss outpaces gene gain in generating novel stable network states. These results suggest developmental complexity may be driven by gene loss, and not just by gene gain.

Significance StatementIncreased developmental complexity is thought to evolve mainly through genetic innovation and co-option. Comparison of the genomes of the closely related volvocines revealed that despite large changes in developmental complexity their genomes are exhibiting significant gene loss. Further, a burst of gene loss is shown to occur at the transition to undifferentiated multicellularity by gene inactivation and decay. This likely results in changes to protein-protein interactions within the cell, suggesting evolutionary novelty does not always require gene gain. Using empirical and modelling approaches we demonstrate gene loss can more easily produce protein network novelty than gene gain. We propose that gene loss can be a driver of biological innovation, allowing for reconfiguration of gene networks and differential use of existing functional repertoires. Trade-offs in biosensor optimization for dynamic pathway engineering Recent progress in synthetic biology allows the construction of dynamic control circuits for metabolic engineering. This technology promises to overcome many challenges encountered in traditional pathway engineering, thanks to their ability to self-regulate gene expression in response to bioreactor perturbations. The central components in these control circuits are metabolite biosensors that read out pathway signals and actuate enzyme expression. However, the construction of metabolite biosensors is a major bottleneck for strain design, and a key challenge is to understand the relation between biosensor dose-response curves and pathway performance. Here we employ multiobjective optimization to quantify performance trade-offs that arise in the design and calibration of metabolite biosensors. Our approach reveals strategies for tuning dose-response curves along an optimal trade-off between production flux and the cost of an increased expression burden on the host. We explore properties of control architectures built in the literature, and identify their advantages and caveats in terms of performance and robustness to growth conditions and leaky promoters. We demonstrate the optimality of a control circuit for glucaric acid production in Escherichia coli, which has been shown to increase titer by 2.5-fold as compared to static designs. Our results lay the groundwork for the automated design of control circuits for pathway engineering, with applications in the food, energy and pharmaceutical sectors. Intermittent compression increases flow through a gel: a mechanical model perhaps relevant to flow of cerebrospinal fluid Exchange of molecules between cerebrospinal fluid (CSF) and brain cells contributes to brain function and protection from dementia. Despite widespread acceptance of the  glymphatic theory, the route by which CSF is brought close enough to the neural tissue for solutes to be exchanged by extracellular diffusion is not entirely clear. Exogenous molecules injected into CSF are observed to reach the basement lamina that surrounds the dense capillary network. Transport of solutes by diffusion along the basement lamina, a gel of macromolcules about 100 nm thick, would be too slow; bulk flow in a static geometry would require unphysiologically high pressures. However, it is known that the pulsation of blood aids transport of CSF, and we hypothesized that this is because the pulsation intermittently squeezes the pericapillary lamina. In a primitive mimicry, we have tested whether intermittent squeezing increases flow through an agar gel. In all but one of 216 tests, pulsation caused a reversible increase, sometimes by a factor of 100 or more. The enhancement was greatest for frequencies 5-11 Hz and, over the tested range of pressure heads (20 - 50 cmH2O), was greatest for the lowest pressure. The results suggest a reason why exercise slows the aging of the brain. Herbivore life histories are altered by drought stress in their hosts plants Drought associated with climate change can stress plants, altering their interactions with phytophagous arthropods. Dryness not only impacts cultivated plants but also their parasites, which in some cases are favoured by drought. Herbivorous arthropods feeding on drought-stressed plants typically produce bigger offspring and develop faster. However, how much responses to abiotic factors differ among populations of a species remains poorly documented. Here, we evaluate variability among populations of a major agricultural pest, the two spotted spider mite, Tetranychus urticae, in response to drought-stress. We compare key life-history parameters of twelve populations that originate from climates ranging from wet and cool Atlantic locations to medium to dry hot Mediterranean locations. We evaluated how plant drought stress affects four life history traits: development time, fecundity, sex-ratio and emigration rate in an experiment comparing well-watered and drought-stressed bean plants. Mites feeding on drought-stressed plants developed faster and attempted to leave leaves less often, and young females were more fecund. The mites from wet temperate climates exhibited greater plasticity between the two water regimes than mites originating from dryer and hot climates. Thus, climate in the area of origin influences mite response to drought. Affordable Caenorhabditis elegans tracking system for classroom use For decades the nematode C. elegans has served as an outstanding research organism owing to its unsurpassed experimental amenability. This advantage has also made this tiny worm an attractive vehicle for science instruction across higher learning institutions. However, the prohibitive cost associated with the automated behavioral assessment of these animals remains an obstacle preventing their full adoption in undergraduate and high school settings. To improve this situation, we developed an inexpensive worm tracking system for use by high school interns and undergraduate students. Over the past two years this tracker has been successfully used by undergraduate students in our introductory Cell and Molecular lab (BSC220) at Illinois State University. Here we describe and demonstrate the use of our inexpensive worm tracking system. A chromosome-level assembly of the "Cascade" hop (Humulus lupulus) genome uncovers signatures of molecular evolution and improves time of divergence estimates for the Cannabaceae family O_LIWe present a chromosome-level assembly of the Cascade hop (Humulus lupulus L. var. lupulus) genome. The hop genome is large (2.8 Gb) and complex, and early attempts at assembly resulted in fragmented assemblies. Recent advances have made assembly of the hop genome more tractable, transforming the extent of investigation that can occur.
C_LIO_LIThe chromosome-level assembly of Cascade was developed by scaffolding the previously-reported Cascade assembly generated with PacBio long-read sequencing, and polishing with Illumina short-read DNA sequencing. We developed gene models and repeat annotations, and used a controlled bi-parental mapping population to identify significant sex-associated markers. We assess molecular evolution in gene sequences, gene family expansion and contraction, and time divergence using Bayesian inference.
C_LIO_LIWe identified the putative sex chromosome in the female genome based on significant sex-associated markers from the bi-parental mapping population. While the estimate of repeat content (~64%) is similar to the hemp genome, syntenic blocks in hop contain a greater percentage of LTRs. Hop is enriched for disease resistance-associated genes in syntenic gene blocks and expanded gene families.
C_LIO_LIThe Cascade chromosome-level assembly will inform cultivation strategies and serve to deepen our understanding of the hop genomic landscape, benefiting hop researchers and the Cannabaceae genomics community.
C_LI Induction of exaggerated cytokine production in human peripheral blood mononuclear cells by a recombinant SARS-CoV-2 spike glycoprotein S1 and its inhibition by dexamethasone An understanding of the pathological inflammatory mechanisms involved in SARS-CoV-2 virus infection is necessary in order to discover new molecular pharmacological targets for SARS-CoV-2 cytokine storm. In this study, the effects of a recombinant SARS-CoV-2 spike glycoprotein S1 was investigated in human peripheral blood mononuclear cells (PBMCs). Stimulation of PBMCs with spike glycoprotein S1 (100 ng/mL) resulted in significant elevation in the production of TNF, IL-6, IL-1{beta} and IL-8. However, pre-treatment with dexamethasone (100 nM) caused significant reduction in the release of these cytokines. Further experiments revealed that S1 stimulation of PBMCs increased phosphorylation of NF-{kappa}B p65 and I{kappa}B, and I{kappa}B degradation. DNA binding of NF-{kappa}B p65 was also significantly increased following stimulation with spike glycoprotein S1. Treatment of PBMCs with dexamethasone (100 nM) or BAY11-7082 (1 {micro}M) resulted in inhibition of spike glycoprotein S1-induced NF-{kappa}B activation. Activation of p38 MAPK by S1 was blocked in the presence of dexamethasone and SKF 86002. CRID3, but not dexamethasone pre-treatment produced significant inhibition of S1-induced activation of NLRP3/caspase-1. Further experiments revealed that S1-induced increase in the production of TNF, IL-6, IL-1{beta} and IL-8 was reduced in the presence of BAY11-7082 and SKF 86002, while CRID3 pre-treatment resulted in the reduction of IL-1{beta} production. These results suggest that SARS-CoV-2 spike glycoprotein S1 stimulated PBMCs to release pro-inflammatory cytokines through mechanisms involving activation of NF-{kappa}B, p38 MAPK and NLRP3 inflammasome. It is proposed that the clinical benefits of dexamethasone in COVID-19 is possibly due to its anti-inflammatory activity in reducing SARS-CoV-2 cytokine storm. Genetic architecture of natural variation of cardiac performance in flies Deciphering the genetic architecture of human cardiac disorders is of fundamental importance but their underlying complexity is a major hurdle. We investigated the natural variation of cardiac performance in the sequenced inbred lines of the Drosophila Genetic Reference Panel (DGRP)1. Genome Wide Associations Studies (GWAS) identified genetic networks associated with natural variation of cardiac traits which were used to gain insights as to the molecular and cellular processes affected. Non-coding variants that we identified were used to map potential regulatory non-coding regions, which in turn were employed to predict Transcription Factors (TFs) binding sites. Cognate TFs, many of which themselves bear polymorphisms associated with variations of cardiac performance, were also validated by heart specific knockdown. Additionally, we showed that the natural variations associated with variability in cardiac performance affect a set of genes overlapping those associated with average traits but through different variants in the same genes. Furthermore, we showed that phenotypic variability was also associated with natural variation of gene regulatory networks. More importantly, we documented correlations between genes associated with cardiac phenotypes in both flies and humans, which supports a conserved genetic architecture regulating adult cardiac function from arthropods to mammals. Specifically, roles for PAX9 and EGR2 in the regulation of the cardiac rhythm were established in both models, illustrating that the characteristics of natural variations in cardiac function identified in Drosophila can accelerate discovery in humans. Protein Complex Structure Prediction Powered by Multiple Sequence Alignments of Interologs from Multiple Taxonomic Ranks and AlphaFold2 AlphaFold2 is expected to be able to predict protein complex structures as long as a multiple sequence alignment (MSA) of the interologs of the target protein-protein interaction (PPI) can be provided. However, preparing the MSA of protein-protein interologs is a non-trivial task due to the existence of paralogs. In this study, a simplified phylogeny-based approach was applied to generate the MSA of interologs, which was then used as the input to AlphaFold2 for protein complex structure prediction. Extensively benchmarked this protocol on non-redundant PPI dataset including 107 bacterial PPIs and 442 eukaryotic PPIs, we show complex structures of 79.5% of the bacterial PPIs and 49.8% of the eukaryotic PPIs can be successfully predicted. Considering PPIs may not be conserved in species with long evolutionary distances, we further restricted interologs in the MSA to different taxonomic ranks of the species of the target PPI in protein complex structure prediction. We found the success rates can be increased to 87.9% for the bacterial PPIs and 56.3% for the eukaryotic PPIs if interologs in the MSA are restricted to a specific taxonomic rank of the species of each target PPI. Finally, we show the optimal taxonomic ranks for protein complex structure prediction can be selected with the application of the recalculated predicted TM-scores of the output models. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. Snowball Earths, population bottlenecks, and Prochlorococcus evolution Prochlorococcus are the most abundant photosynthetic organisms in the modern ocean. A massive DNA loss event occurred in their early evolutionary history, leading to highly reduced genomes in nearly all lineages, as well as enhanced efficiency in both nutrient uptake and light absorption. The environmental landscape that shaped this ancient genome reduction, however, remained unknown. Through careful molecular clock analyses, we established that this Prochlorococcus genome reduction occurred during the Neoproterozoic Snowball Earth climate catastrophe. The lethally low temperature and exceedingly dim light during the Snowball Earth event would have inhibited Prochlorococcus growth and proliferation and caused severe population bottlenecks. These bottlenecks are recorded as an excess of deleterious mutations that accumulated across genomic regions in the descendant lineages. Prochlorococcus adaptation to extreme environmental conditions during Snowball Earth intervals can be inferred by tracing the evolutionary paths of genes that encode key metabolic potential. This metabolic potential includes modified lipopolysaccharide structure, strengthened peptidoglycan biosynthesis, the replacement of a sophisticated circadian clock with an hourglass-like mechanism that resets daily for dim light adaption, and the adoption of ammonia diffusion as an efficient membrane transporter-independent mode of nitrogen acquisition. In this way, the Neoproterozoic Snowball Earth event altered the physiological characters of Prochlorococcus, shaping their ecologically vital role as the most abundant primary producers in the modern oceans. Global analysis of human-to-mouse contact-dependent intercellular mRNA and lncRNA transfer in cell culture Full-length mRNAs can transfer between adjacent mammalian cells via direct cell-to-cell connections called tunneling nanotubes (TNTs). However, the extent of mRNA transfer at the transcriptome-wide level (the transferome) is unknown. Here, we analyzed whole transcriptome mRNA and lncRNA transfer between heterogeneous human-mouse cell populations in in vitro co-culture using RNA-sequencing. Our data indicate that mRNA transfer is non-selective, prevalent across the human transcriptome, and that the amount of transfer to mouse embryonic fibroblasts (MEFs) strongly correlates with the endogenous level of gene expression in donor human breast cancer cells (MCF7). These results were validated by both quantitative RT-PCR and in situ hybridization, and analysis shows that typically <1% of endogenous mRNAs and lncRNAs undergo transfer. Non-selective expression-dependent RNA transfer was further validated using synthetic RNA reporters. Notably, significant differential changes in the native MEF transcriptome were observed in response to co-culture, including the upregulation of multiple cancer and cancer-associated fibroblast-related genes and pathways. Together, these results lead us to suggest that TNT-mediated RNA transfer could be a phenomenon of physiological importance under both normal and pathogenic conditions.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=183 SRC="FIGDIR/small/470233v2_ufig1.gif" ALT="Figure 1">
View larger version (44K):
org.highwire.dtl.DTLVardef@13a7176org.highwire.dtl.DTLVardef@1ef1c75org.highwire.dtl.DTLVardef@1355eeforg.highwire.dtl.DTLVardef@ebc275_HPS_FORMAT_FIGEXP  M_FIG C_FIG Fungal community turnover along steep environmental gradients from geothermal soils in Yellowstone National Park Geothermal soils offer unique insight into the way extreme environmental factors shape communities of organisms. However, little is known about the fungi growing in these environments and in particular how localized steep abiotic gradients affect fungal diversity. We used metabarcoding to characterize soil fungi surrounding a hot spring-fed thermal creek with water up to 84 {degrees}C and pH 10 in Yellowstone National Park. We found a significant association between fungal communities and soil variable principal components, and we identify the key trends in co-varying soil variables that explain the variation in fungal community. Saprotrophic and ectomycorrhizal fungi community profiles followed, and were significantly associated with, different soil variable principal components, highlighting potential differences in the factors that structure these different fungal trophic guilds. In addition, in vitro growth experiments in four target fungal species revealed a wide range of tolerances to pH levels but not to heat. Overall, our results documenting turnover in fungal species within a few hundred meters suggest many co-varying environmental factors structure the diverse fungal communities found in the soils of Yellowstone National Park. Extended Graphical Lasso for Multiple Interaction Networks for High Dimensional Omics Data There has been a spate of interest in association networks in biological and medical research, for example, genetic interaction networks. In this paper, we propose a novel method, the extended joint hub graphical lasso (EDOHA), to estimate multiple related interaction networks for high dimensional omics data across multiple distinct classes. To be specific, we construct a convex penalized log likelihood optimization problem and solve it with an alternating direction method of multipliers (ADMM) algorithm. The proposed method can also be adapted to estimate interaction networks for high dimensional compositional data such as microbial interaction networks. The performance of the proposed method in the simulated studies shows that EDOHA has remarkable advantages in recognizing class-specific hubs than the existing comparable methods. We also present three applications of real datasets. Biological interpretations of our results confirm those of previous studies and offer a more comprehensive understanding of the underlying mechanism in disease.

Author summaryReconstruction of multiple association networks from high dimensional omics data is an important topic, especially in biology. Previous studies focused on estimating different networks and detecting common hubs among all classes. Integration of information over different classes of data while allowing difference in the hub nodes is also biologically plausible. Therefore, we propose a method, EDOHA, to jointly construct multiple interaction networks with capacity in finding different hub networks for each class of data. Simulation studies show the better performance over conventional methods. The method has been demonstrated in three real world data.