Visual Clustering of Transcriptomic Data from Primary and Metastatic Tumors - Dependencies and Novel Pitfalls Personalized Oncology is a rapidly evolving area and offers cancer patients therapy options more specific than ever. Yet, 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 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 site of resection and cluster formation outcome could be demonstrated, neither for datasets consisting of solely metastasis nor 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 interpretation of visual clustering approaches along with initialization and different parameters. Furthermore, the results indicate only minor transcriptional differences for metastases and corresponding primary tumors.
Comprehensive analysis of pyroptosis-associated in molecular classification, immunity and prognostic of glioma Integrative analysis was performed in the Chinese Glioma Genome Atlas and The Cancer Genome Atlas to describe the pyroptosis-associated molecular classification and prognostic signature in glioma. Pyroptosis-related genes were used for consensus clustering and to develop a prognostic signature. The immune statuses, molecular alterations and clinical features of differentially expressed genes were analyzed among different subclasses and risk groups. A lncRNA-miRNA-mRNA network was built, and drug sensitivity analysis was used to identify small molecular drugs for the identified genes. Glioma can be divided into two subclasses using 30 pyroptosis-related genes. Cluster 1 displayed high immune signatures and poor prognosis as well as high immune-related function scores. A prognostic signature based on 15 pyroptosis-related genes of the CGGA cohort can predict the overall survival of glioma and was well validated in the TCGA cohort. Cluster 1 had higher risk scores. The high-risk group had high immune cell and function scores and low DNA methylation of pyroptosis-related genes. The differences in pyroptosis-related gene mutations and somatic copy numbers were significant between the high-risk and low-risk groups. The ceRNA regulatory network uncovered the regulatory patterns of different risk groups in glioma. Nine pairs of target genes and drugs were identified. In vitro, CASP8 promotes the progression of glioma cells. Pyroptosis-related genes can reflect the molecular biological and clinical features of glioma subclasses. The established prognostic signature can predict prognosis and distinguish molecular alterations in glioma patients. Our comprehensive analyses provide valuable guidelines for improving glioma patient management and individualized therapy.
BlenderPhotonics - an integrated open-source software environment for 3-D meshing and photon simulations in complex tissues SignificanceRapid advances in biophotonics techniques require quantitative, model-based computational approaches to obtain functional and structural information from increasingly complex and multi-scaled anatomies. The lack of efficient tools to accurately model tissue structures and subsequently perform quantitative multi-physics modeling greatly impedes the clinical translation of these modalities.

AimWhile the mesh-based Monte Carlo (MMC) method expands our capabilities in simulating complex tissues by using tetrahedral meshes, the generation of such domains often requires specialized meshing tools such as Iso2Mesh. Creating a simplified and intuitive interface for tissue anatomical modeling and optical simulations is essential towards making these advanced modeling techniques broadly accessible to the user community.

ApproachWe responded to the above challenge by combining the powerful, open-source 3-D modeling software, Blender, with state-of-the-art 3-D mesh generation and MC simulation tools, utilizing the interactive graphical user interface (GUI) in Blender as the front-end to allow users to create complex tissue mesh models, and subsequently launch MMC light simulations.

ResultsHere we present a tutorial to our newly developed Python-based Blender add-on - BlenderPhotonics - to interface with Iso2Mesh and MMC, allowing users to create, configure and refine complex simulation domains and run hardware-accelerated 3-D light simulations with only a few clicks. We provide a comprehensive introduction to this new tool and walk readers through 5 examples, ranging from simple shapes to sophisticated realistic tissue models.

ConclusionBlenderPhotonics is user-friendly and open-source, leveraging the vastly rich ecosystem of Blender. It wraps advanced modeling capabilities within an easy-to-use and interactive interface. The latest software can be downloaded at http://mcx.space/bp. Microbial colonization of tannin-rich tropical plants: interplay between degradability, methane production and tannin disappearance in the rumen AO_SCPLOWBSTRACTC_SCPLOWCondensed tannins in plants are found free and attached to protein and fibre but it is not known whether these fractions influence rumen degradation and microbial colonization. This study explored the rumen degradation of tropical tannins-rich plants and the relationship between condensed tannins fractions disappearance and microbial communities colonising plant particles. Leaves from Calliandra calothyrsus, Gliricidia sepium, and Leucaena leucocephala, pods from Acacia nilotica and the leaves of two agricultural by-products: Manihot esculenta and Musa spp. were incubated in situ in the rumen of dairy cows. An in vitro approach was also used to assess the effects of these plants on rumen fermentation parameters. All plants contained more than 100 g/kg of condensed tannins with a large proportion (32 to 61%) bound to proteins. Calliandra calothyrsus had the highest concentration of condensed tannins at 361 g/kg, whereas Acacia nilotica was particularly rich in hydrolysable tannins (350 g/kg). Hydrolysable and free condensed tannins from all plants completely disappeared after 24 h incubation in the rumen. Disappearance of protein-bound condensed tannins was variable with values ranging from 93% for Gliricidia sepium to 21% for Acacia nilotica. In contrast, fibre-bound condensed tannins disappearance averaged ~82% and did not vary between plants. Disappearance of bound fractions of condensed tannins was not associated with degradability of plant fractions. The presence of tannins interfered with the microbial colonisation of plants. Each plant had distinct bacterial and archaeal communities after 3 and 12 h of incubation in the rumen and distinct protozoal communities at 3 h. Adherent communities in tannin-rich plants had a lower relative abundance of fibrolytic microbes, notably Fibrobacter spp. Whereas, archaea diversity was reduced in high tannin-containing Calliandra calothyrsus and Acacia nilotica at 12 h of incubation. Concurrently, in vitro methane production was lower for Calliandra calothyrsus, Acacia nilotica and Leucaena leucocephala although for the latter total volatile fatty acids production was not affected and was similar to control. Here we show that the total amount of hydrolysable and condensed tannins contained in a plant govern the interaction with rumen microbes affecting degradability and fermentation. The effect of protein- and fibre-bound condensed tannins on degradability is less important. Curation, inference, and assessment of a globally reconstructed gene regulatory network for Streptomyces coelicolor Streptomyces coelicolor A3(2) is a model microorganism for the study of Streptomycetes, antibiotic production, and secondary metabolism in general. Even though S. coelicolor has an outstanding variety of regulators among bacteria, little effort to globally study its transcription has been made. We manually curated 29 years of literature and databases to assemble a meta-curated experimentally-validated gene regulatory network (GRN) with 5386 genes and 9707 regulatory interactions ([~]41% of the total expected interactions). This provides the most extensive and up-to-date reconstruction available for the regulatory circuitry of this organism. Only [~]6% (534/9707) are supported by experiments confirming the binding of the transcription factor to the upstream region of the target gene, the so-called "strong" evidence. While for the remaining interactions there is no confirmation of direct binding. To tackle network incompleteness, we performed network inference using several methods (including two proposed here) for motif identification in DNA sequences and GRN inference from transcriptomics. Further, we contrasted the structural properties and functional architecture of the networks to assess the reliability of the predictions, finding the inference from DNA sequence data to be the most trustworthy approach. Finally, we show two applications of the inferred and the curated networks. The inference allowed us to propose novel transcription factors for the key Streptomyces antibiotic regulatory proteins (SARPs). The curated network allowed us to study the conservation of the system-level components between S. coelicolor and Corynebacterium glutamicum. There we identified the basal machinery as the common signature between the two organisms. The curated networks were deposited in Abasy Atlas (https://abasy.ccg.unam.mx/) while the inferences are available as Supplementary Material. The obstacles and potential clues of prime editing applications in tomato, a dicot plant Precision genome editing is highly desired for crop improvement. The recently emerged CRISPR/Cas technology offers great potential applications in precision plant genome engineering. A prime editing (PE) approach combining a reverse transcriptase (RT) with a Cas9 nickase and a "priming" extended guide RNA has shown a high frequency for precise genome modification in mammalian cells and several plant species. Nevertheless, the applications of the PE approach in dicot plants are still limited and inefficient. We designed and tested prime editors for precision editing of a synthetic sequence in a transient assay and for desirable alleles of 10 loci in tomato by stable transformation. Our data obtained by targeted deep sequencing also revealed only low PE efficiencies in both the tobacco and tomato systems. Further assessment of the activities of the PE components uncovered potential reasons for the inefficiency of the PE complexes. Modifying the pegRNA sequences by shortening or introducing mismatches to the primer binding sequences (PBS) in order to reduce their melting temperatures (Tm) did not enhance the PE efficiency at the SlBMP21, SlALC and SlALS1 loci. Our data show challenges of PE approach in tomato, indicating a further improvement of the PE system for the successful applications such as use of improved expression systems. Our work provides an important clue for the successful application of the PE approach in crop improvement. Ecology of active viruses and their bacterial hosts in frozen Arctic peat soil revealed with H218O stable isotope probing metagenomics Winter carbon loss in northern ecosystems is estimated to be greater than the average growing season carbon uptake. However, most ecosystem carbon measurements neglect winter months since carbon losses (primarily driven by microbial decomposers) are assumed to be negligible at low temperatures. We used stable isotope probing (SIP) targeted metagenomics to reveal the genomic potential of active soil microbial populations under winter conditions, with an emphasis on viruses and virus-host dynamics. Peat soils from the Bonanza Creek LTER site in Alaska were incubated under subzero anoxic conditions with H218O for 184 and 370 days. We identified 46 bacterial populations (MAGs; spanning 9 bacterial phyla) and 243 viral populations (vOTUs) that actively took up 18O and produced significant CO2 throughout the incubation. Active hosts, predicted for 33% of the active vOTUs, were some of the most abundant MAGs and capable of fermentation and organic matter degradation. Approximately three-quarters of the active vOTUs carried auxiliary metabolic genes that spanned five functional categories, including carbon utilization, highlighting the potential impact of viruses in this peat soils microbial biogeochemistry. These results illustrate significant bacterial and viral activity and interactions occur in frozen soils, revealing viruses are a major community-structuring agent throughout winter months. Psychedelics flatten the brain's energy landscape: evidence from receptor-informed network control theory Psychedelics like lysergic acid diethylamide (LSD) and psilocybin offer a powerful window into the function of the human brain and mind, by temporarily altering subjective experience through their neurochemical effects. A recent model postulates that serotonin 2a (5-HT2a) receptor agonism allows the brain to explore its dynamic landscape more readily, as reflected by more diverse (entropic) brain activity. We postulate that this increase in entropy may arise in part from a flattening of the brains control energy landscape, which can be observed using network control theory to quantify the energy required to transition between recurrent brain states measured using functional magnetic resonance imaging (fMRI) in individuals under LSD, psilocybin, and placebo conditions. We show that LSD and psilocybin reduce the amount of control energy required for brain state transitions, and, furthermore, that, across individuals, LSDs reduction in control energy correlates with more frequent state transitions and increased entropy of brain state dynamics. Through network control analysis that incorporates the spatial distribution of 5-HT2a receptors from publicly available (non-drug) positron emission tomography (PET) maps, we demonstrate the specific role of this receptor in reducing control energy. Our findings provide evidence that 5-HT2a receptor agonist compounds allow for more facile state transitions and more temporally diverse brain activity. More broadly, by combining receptor-informed network control theory with pharmacological modulation, our work highlights the potential of this approach in studying the impacts of targeted neuropharmacological manipulation on brain activity dynamics.

Significance StatementWe present a multi-modal framework for quantifying the effects of two psychedelic drugs (LSD and psilocybin) on brain dynamics by combining functional magnetic resonance imaging (fMRI), diffusion MRI (dMRI), positron emission tomography (PET) and network control theory. Our findings provide evidence that psychedelics flatten the brains control energy landscape, allowing for more facile state transitions and more temporally diverse brain activity. We also demonstrate that the spatial distribution of serotonin 2a receptors - the main target of LSD and psilocybin - is optimized for generating these effects. This approach could be used to understand how drugs act on different receptors in the brain to influence brain activity dynamics. LPS-TLR4 pathway exaggerates alcohol-induced acute-on-chronic liver failure via provoking NETs formation Backgrounds: Intrahepatic infiltration of neutrophils is a character of alcoholic acute-on-chronic liver failure (AACLF) and neutrophil extracellular traps (NETs) are an important strategy for neutrophils to fix and kill invading microorganisms. Intestinal bacteria and the gut-liver axis have been thought to play a key role in many liver diseases also including AACLF. However, whether NETs appear in AACLF and play a role in AACLF is still unsure. Methods: WT, NE KO, and TLR4 KO mice were used to build the AACLF model, and the intestinal bacteria were eliminated at the same time and LPS was given. Then the formation of NETs and AACLF related markers were detected. Results: The serum MPO-DNA and LPS concentration was increased in AACLF patients and a correlation was revealed between these two indexes. More intrahepatic NETs formed in AACLF mice by testing MPO-DNA, Cit H3, and NE. These markers decreased with gut detergent and restored markers with gut detergent plus LPS supplement. While NETs formation failed to change with gut microbiome or combine LPS supplement in TLR4 KO mice. As we tested AACLF related characters, liver injury, intrahepatic fat deposition, inflammation, and fibrosis alleviated with depletion of NE. These related marks were also attenuated with gut sterilization by antibiotics and recovered with combined treatment with antibiotics plus LPS. But the liver injury, intrahepatic fat, fibro deposition, and liver inflammation-related markers did show a significant difference in TLR4 KO mice when they received the same treatment. Conclusion: Intestinal-derived LPS promotes NETs formation in AACLF through the TLR4 pathway and further accelerates the AACLF process by NETs. Trimeric Photosystem I facilitates energy transfer from phycobilisomes in Synechocystis sp. PCC 6803 In cyanobacteria, phycobilisomes serve as peripheral light-harvesting complexes of the two photosystems, extending their antenna size and the wavelength range of photons available for photosynthesis. The abundance of phycobilisomes, the number of phycobiliproteins they contain, and their light-harvesting function are dynamically adjusted in response to the physiological conditions. Phycobilisomes are also thought to be involved in state transitions that maintain the excitation balance between the two photosystems. Unlike its eukaryotic counterpart, PSI is trimeric in many cyanobacterial species and the physiological significance of this is not well understood. Here we compared the composition and light-harvesting function of phycobilisomes in cells of Synechocystis sp. PCC 6803, which has primarily trimeric PSI, and the{Delta} psaL mutant unable to form trimers. We also investigated a mutant additionally lacking the PsaJ and PsaF subunits of PSI. Both strains with monomeric PSI accumulated significantly more allophycocyanin per chlorophyll, indicating higher abundance of phycobilisomes. On the other hand, a higher phycocyanin:allophycocyanin ratio in WT suggests larger phycobilisomes or the presence of APC-less phycobilisomes (CpcL-type), that are not assembled in cells with monomeric PSI. Steady-state and time-resolved fluorescence spectroscopy at room temperature and 77 K revealed that PSII receives more energy from the phycobilisomes at the expense of PSI in cells with monomeric PSI, regardless of the presence of PsaF. Taken together, these results show that the oligomeric state of PSI has an impact on the excitation energy flow in Synechocystis.

One-sentence summaryCyanobacterial mutants with monomeric PSI show changes in the composition and abundance of phycobilisomes and in the excitation energy transfer to PSII and PSI. The structurome of a Clostridium difficile phage and the remarkable accurate prediction of its novel phage receptor-binding protein As natural bacterial predators, bacteriophages have the potential to be developed to tackle antimicrobial resistance, but our exploitation of them is limited by understanding their vast uncharacterised genetic diversity1,2. Fascinatingly, this genetic diversity reflects many ways that phages can make proteins, performing similar functions that together form the familiar phage particle. Critical to infection are phage receptor-binding proteins (RBPs) that bind bacterial  receptors and initiate bacterial entry3. Here we identified and characterised Gp22, a novel RBP for phage CDHS-1 that infects pathogenic C. difficile, but that had no recognisable RBPs. We showed that Gp22 antibodies neutralised CDHS-1 infection and used immunogold-labelling and transmission electron microscopy to identify their location on the capsid. The Gp22 three-dimensional structure was resolved by X-ray crystallography revealing a new RBP class with an N-terminal L-shaped -helical superhelix domain and a C-terminal Mg2+-binding domain. The findings provide novel insights into C. difficile phage biology and phage-host interactions. This will facilitate optimal phage development and future engineering strategies4,5. Furthermore, the AlphaFold2-predicted Gp22 structure, which was strikingly accurate, paves the way for a structurome based transformation and guidance of future phage studies where many proteins lack sequence homology but have recognisable protein structures. ZIKV disrupts placental ultrastructure and drug transporter expression in mice Congenital Zika virus (ZIKV) infection can induce fetal brain abnormalities. Here, we investigated whether maternal ZIKV infection affects placental physiology and metabolic transport potential and impacts the fetal outcome, regardless of viral presence in the fetus at term. Low (103 PFU-ZIKVPE243; low ZIKV) and high (5x107 PFU-ZIKVPE243; high ZIKV) virus titers were injected into immunocompetent (ICompetent C57BL/6) and immunocompromised (ICompromised A129) mice at gestational day (GD) 12.5 for tissue collection at GD18.5 (term). High ZIKV elicited fetal death rates of 66% and 100%, whereas low ZIKV induced fetal death rates of 0% and 60% in C57BL/6 and A129 dams, respectively. All surviving fetuses exhibited intrauterine growth restriction (IUGR) and decreased placental efficiency. High-ZIKV infection in C57BL/6 and A129 mice resulted in virus detection in maternal spleens and placenta, but only A129 fetuses presented virus RNA in the brain. Nevertheless, pregnancies in both strains produced fetuses with decreased head sizes (p<0.05). Low-ZIKV-A129 dams had higher IL-6 and CXCL1 levels (p<0.05), and their placentas showed increased CCL-2 and CXCL-1 contents (p<0.05). In contrast, low-ZIKV-C57BL/6 dams had an elevated CCL2 serum level and increased type I and II IFN expression in the placenta. Notably, less abundant microvilli and mitochondrial degeneration were evidenced in the placental labyrinth zone (Lz) of ICompromised and high-ZIKV-ICompetent mice but not in low-ZIKV-C57BL/6 mice. In addition, decreased placental expression of the drug transporters P-glycoprotein (P-gp) and breast cancer resistance protein (Bcrp) and the lipid transporter Abca1 was detected in all ZIKV-infected groups, but Bcrp and Abca1 were only reduced in ICompromised and high-ZIKV ICompetent mice. Our data indicate that gestational ZIKV infection triggers specific proinflammatory responses and affects placental turnover and transporter expression in a manner dependent on virus concentration and maternal immune status. Placental damage may impair proper fetal-maternal exchange function and fetal growth/survival, likely contributing to congenital Zika syndrome. Constraints on the evolution of toxin-resistant Na,K-ATPases have limited dependence on sequence divergence A growing body of theoretical and experimental evidence suggests that intramolecular epistasis is a major determinant of rates and patterns of protein evolution and imposes a substantial constraint on the evolution of novel protein functions. Here, we examine the role of intramolecular epistasis in the recurrent evolution of resistance to cardiotonic steroids (CTS) across tetrapods, which occurs via specific amino acid substitutions to the -subunit family of Na,K-ATPases (ATP1A). After identifying a series of recurrent substitutions at two key sites of ATP1A that are predicted to confer CTS resistance in diverse tetrapods, we then performed protein engineering experiments to test the functional consequences of introducing these substitutions onto divergent species backgrounds. In line with previous results, we find that substitutions at these sites can have substantial background-dependent effects on CTS resistance. Globally, however, these substitutions also have pleiotropic effects that are consistent with additive rather than background-dependent effects. Moreover, the magnitude of a substitutions effect on activity does not depend on the overall extent of ATP1A sequence divergence between species. Our results suggest that epistatic constraints on the evolution of CTS-resistant forms of Na,K-ATPase likely depend on a small number of sites, with little dependence on overall levels of protein divergence. We propose that dependence on a limited number sites may account for the observation of convergent CTS resistance substitutions observed among taxa with highly divergent Na,K-ATPases.

Significance StatementIndividual amino acid residues within a protein work in concert to produce a functionally coherent structure that must be maintained even as orthologous proteins in different species diverge over time. Given this dependence, we expect identical mutations to have more similar effects on protein function in more closely related species. We tested this hypothesis by performing protein-engineering experiments on ATP1A, an enzyme mediating target-site insensitivity to cardiotonic steroids (CTS) in diverse animals. These experiments reveal that the phenotypic effects of substitutions can sometimes be background-dependent, but also that the magnitude of these phenotypic effects does not correlate with overall levels of ATP1A sequence divergence. Our results suggest that epistatic constraints are determined by states at a small number of sites, potentially explaining the frequent convergent CTS resistance substitutions among Na,K-ATPases of highly divergent taxa. Using Obviously-Related Instrumental Variables to Increase the Predictive Power of Polygenic Scores The conventional way of boosting the predictive power of polygenic scores is to increase the GWAS sample size by meta-analyzing GWAS results of multiple cohorts. In this paper, we challenge this convention. Through simulations, we show that Instrumental Variable (IV) regression using two polygenic scores constructed from independent GWAS summary statistics outperforms the typical Ordinary Least Squares (OLS) model employing a single meta-analysis based polygenic score in terms of bias, root mean squared error, and statistical power. We verify the empirical validity of the simulations by predicting educational attainment (EA) and height in a sample of siblings from the UK Biobank. We show that IV regression between-families approaches the SNP-based heritability, and improves the predictive power of polygenic scores by 12% (height) to 22% (EA). Furthermore, IV regression within-families provides the tightest lower bound for the direct genetic effect, increasing the lower bound for EA from 2.0% to 3.4%, and for height from 28.9% to 37.7%. A terrain treadmill to study animal locomotion through large obstacles A major challenge to understanding locomotion in complex 3-D terrain with large obstacles is to create tools for controlled, systematic lab experiments. Existing terrain arenas only allow observations at small spatiotemporal scales ([~]10 body length, [~]10 stride cycles). Here, we create a terrain treadmill to enable high-resolution observations of animal locomotion through large obstacles over large spatiotemporal scales. An animal moves through modular obstacles on an inner sphere, while a rigidly-attached, concentric, transparent outer sphere rotated with the opposite velocity via closed-loop feedback to keep the animal on top. During sustained locomotion, a discoid cockroach moved through pillar obstacles for 25 minutes ({approx} 2500 strides) over 67 m ({approx} 1500 body lengths), and was contained within a radius of 4 cm (0.9 body length) for 83% of the duration, even at speeds of up to 10 body length/s. The treadmill enabled observation of diverse locomotor behaviors and quantification of animal-obstacle interaction.

Summary StatementA device keeps a small animal on top of a sphere while it traverses obstacles, creating a "terrain treadmill" to study locomotion over a long time and distance at high-resolution. Distinct SARS-CoV-2 antibody reactivity patterns elicited by natural infection and mRNA Vaccination We analyzed data from two ongoing COVID-19 longitudinal serological surveys in Orange County, CA., between April 2020 and March 2021. A total of 8,476 finger stick blood specimens were collected before and after an aggressive mRNA vaccination campaign. IgG levels were determined using a multiplex antigen microarray containing 10 SARS-CoV-2 antigens, 4 SARS, 3 MERS, 12 Common CoV, and 8 Influenza antigens. Twenty-six percent of 3,347 specimens from unvaccinated Orange County residents in December 2020 were SARS-CoV-2 seropositive. The Ab response was predominantly against nucleocapsid (NP), full length spike and the spike S2 domain. Anti-receptor binding domain (RBD) reactivity was low and there was no cross-reactivity against SARS S1 or SARS RBD. An aggressive mRNA vaccination campaign at the UCI Medical Center started on December 16, 2020 and 6,724 healthcare workers were vaccinated within 3 weeks. Seroprevalence increased from 13% in December to 79% in January, 93% in February and 99% in March. mRNA vaccination induced much higher Ab levels especially against the RBD domain and significant cross-reactivity against SARS RBD and S1 was also observed. Nucleocapsid protein Abs can be used to distinguish individuals in a population of vaccinees to classify those who have been previously infected and those who have not, because nucleocapsid is not in the vaccine. Previously infected individuals developed higher Ab titers to the vaccine than those who have not been previously exposed. These results indicate that mRNA vaccination rapidly induces a much stronger and broader Ab response than SARS-CoV-2 infection. Modeling microRNA-driven post-transcriptional regulation by using exon-intron split analysis (EISA) in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses. Mechanisms of splicing substrate definition by modulation of pre-mRNA 3D structural scaffolds Coordination of different serine-arginine-rich (SR) proteins - a class of critical splicing activators - facilitates recognition of the highly degenerate mammalian splice signal sequences against the background sequences. Yet, the mechanistic details of their actions remain unclear. Here we show that cooperative binding of SR proteins to exonic and intronic motifs remodels the pre-mRNA 3D structural scaffold. The scaffold generated by the pre-mRNA-specific combinations of different SR proteins in an appropriate stoichiometry is recognized by U1 snRNP - the 5' splice site recognizing factor. A large excess of U1 snRNP particles displaces majority of the bound SR protein molecules from the remodeled pre-mRNA. A higher than optimal stoichiometry of SR proteins occludes the binding sites on the pre-mRNA, raising the U1 snRNP levels required for SR protein displacement and potentially impeding spliceosome assembly. This novel splicing-related function of U1 snRNP is important for distinguishing the substrate and the non-substrate by U2AF65 - the primary 3' splice site-recognizing factor. Overall, this work elucidates early regulatory steps of mammalian splicing substrate definition by SR proteins. 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. The chirality of the mitotic spindle provides a mechanical response to forces and depends on microtubule motors and augmin Forces produced by motor proteins and microtubule dynamics within the mitotic spindle are crucial for proper chromosome segregation. In addition to linear forces, rotational forces or torques are present in the spindle, reflected in the left-handed twisted shapes of microtubule bundles that make the spindle chiral. However, the biological role and molecular origins of spindle chirality are unknown. By developing methods for measuring spindle twist, we show that spindles are most chiral near the metaphase-to-anaphase transition. To assess the role of chirality in maintaining spindle robustness under force, we compressed the spindles along their axis. This resulted in stronger left-handed twist, suggesting that the twisted shape allows for a mechanical response to forces. Inhibition or depletion of motor proteins that perform chiral stepping, Eg5/kinesin-5, Kif18A/kinesin-8, MKLP1/kinesin-6, and dynein, decreased the left-handed twist or led to right-handed twist, implying that these motors regulate the twist by rotating microtubules within their antiparallel overlaps or at the spindle pole. Right-handed twist was also observed after the depletion of the microtubule nucleator augmin, indicating its contribution to the twist through the nucleation of antiparallel bridging microtubules. The uncovered switch from left-handed to right-handed twist reveals the existence of competing mechanisms that promote twisting in opposite directions. As round spindles were more twisted than elongated ones, we infer that bending and twisting moments are generated by similar molecular mechanisms and propose a physiological role for spindle chirality in allowing the spindle to absorb mechanical load. AutoGGN: A Gene Graph Network AutoML Tool for Multi-Omics Research Omics data identifies biological characteristics from genetic to phenotypic levels during the life span. Molecular interaction networks have a fundamental impact on life activities. Integrating omics data and molecular interaction networks will help researchers delve into comprehensive information underlying the data. Here, we proposed a new multimodal method called AutoGGN to aggregate multi-omics data and molecular interaction networks based on graph convolutional neural networks. We evaluated AutoGGN using two different tasks: cancer type classification and single-cell stage classification. On both tasks, AutoGGN showed better performance compared to other methods, the trend is relevant to the ability of utilizing much more information from biological data. The phenomenon indicated AutoGGN has the potential to incorporate valuable information from molecular interaction networks and multi-omics data effectively. Furthermore, in order to provide a better understanding of the mechanism of prediction results, we assessed the explanation using SHAP module and identified the key genes contributing to the prediction of classification, which will provide insights for the downstream design of biological experiments. Sugar Signaling Induces Dynamic Changes during Meristem Development to Regulate Flowering in Arabidopsis Aerial parts of plants originate from pluripotent stem cells in the shoot apical meristem. Their population is maintained via the maintenance regulators WUSCHEL and CLAVATA3 in a negative feed-back loop. Meristem size is dynamic and undergoes a more than 2-fold expansion upon the transition to reproductive growth. The mechanism controlling this doming is largely unknown, but coinciding increased trehalose 6-phosphate and changed meristem size in overexpressing or knockdown lines of TREHALOSE PHOSPHATE SYNTHASE1 suggest a participation of sugar signaling. Here we show that TREHALOSE PHOSPHATE PHOSPHATASEJ is directly regulated by WUSCHEL. Plants with reduced levels of TREHALOSE PHOSPHATE PHOSPHATASEJ in the outer meristem layer are flowering early and its reduction in the late flowering clavata3 mutant, restores wild-type flowering. This is caused by a reduction of mature miR156 and increased expression of SQUAMOSA PROMOTER-BINDING PROTEIN-LIKE genes. We demonstrate that these are important for age pathway-induced flowering, in a negative feed-back loop with WUSCHEL downstream of the trehalose 6-phosphate pathway. In summary, our findings demonstrate a dynamic feed-back regulation between central maintenance and flowering time regulators with sugar signaling.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=169 SRC="FIGDIR/small/439483v2_ufig1.gif" ALT="Figure 1">
View larger version (34K):
org.highwire.dtl.DTLVardef@988d6eorg.highwire.dtl.DTLVardef@16d469borg.highwire.dtl.DTLVardef@1366ce5org.highwire.dtl.DTLVardef@2753f0_HPS_FORMAT_FIGEXP  M_FIG O_FLOATNOSynopsis:C_FLOATNO Schematic illustration of dynamic feedback-regulations between sugar signaling, meristem maintenance and the age pathway of the flowering network in the transition SAM as compared to the vegetative SAM.During vegetative growth SAM maintenance is controlled by spatially separated WUS and CLV3 expression in a negative feedback-loop. At floral transition increased activity of the T6P pathway in the SAM uncouples this regulation. This results in an unreported special relocation of WUS expression and involves a transient negative feedback-loop between WUS and SPL4.

C_FIG 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. Cellular feedback to organotelluranes displays modulation of antioxidant proteins gene expression Organotelluranes RT3 and RT4 are thiol reagents that induce mitochondrial transition pore (MTP) opening in a sensitive and insensitive manner to cyclosporin A. Although RT3 and RT4 promote glutathione depletion, paradoxically, they are also an efficient antioxidant for membrane lipids. These compounds antagonistic effects elicited the challenging question of how the gene expression of antioxidant enzymes would respond to treatment with these compounds. The influence of RT3 and RT4 on antioxidant enzyme expression was investigated in cultured aortic smooth muscle cells (ASMC). RT3 and RT4 promoted disruption of ionic calcium homeostasis, mitochondrial transmembrane potential ({Delta}{Psi}), and cell death in a dose-dependent manner. The cell death mechanisms responded qualitatively to the increase of the organotellurane concentration and changed from apoptosis to necrosis. RT3 and RT4 increased the expression of thioredoxin significantly. RT3 also increased the expression of glutaredoxin and glutathione peroxidase, slightly the catalase expression without significant effects on SOD expression. The results are consistent with GSH and protein thiol depletion and discussed based on the cell toxicity mechanism exhibited by these compounds.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=150 SRC="FIGDIR/small/425411v4_ufig1.gif" ALT="Figure 1">
View larger version (73K):
org.highwire.dtl.DTLVardef@52bb82org.highwire.dtl.DTLVardef@19476d9org.highwire.dtl.DTLVardef@6aa99org.highwire.dtl.DTLVardef@1aa605e_HPS_FORMAT_FIGEXP  M_FIG C_FIG Life history adaptations to fluctuating environments: Combined effects of demographic buffering and lability of demographic parameters Demographic buffering and lability have been identified as adaptive strategies to optimise fitness in a fluctuating environment. These are not mutually exclusive, however we lack efficient methods to measure their relative importance for a given life history. Here, we decompose the stochastic growth rate (fitness) into components arising from nonlinear responses and variance-covariance of demographic parameters to an environmental driver, which allows studying joint effects of buffering and lability. We apply this decomposition for 154 animal matrix population models under different scenarios, to explore how these main fitness components vary across life histories. Faster-living species appear more responsive to environmental fluctuations, either positively or negatively. They have the highest potential for strong adaptive demographic lability, while demographic buffering is a main strategy in slow-living species. Our decomposition provides a comprehensive framework to study how organisms adapt to variability through buffering and lability, and to predict species responses to climate change. Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations.