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.
Exploratory analysis of multiple traits co-adaptations in the population history During the history of range expansion, the populations encounter with variety of environments. They respond to the local environments by modifying the mutually interacting traits. Therefore, to understand the whole life history of the populations, it is ideal to capture the history of their range expansion with reference to the series of surrounding environments and to infer the coadaptation of the multiple traits. Toward this end, we provide an exploratory analysis based on the features of populations: site frequency spectra of populations, population-specific FST, association between genes and environments, positive selections on traits mapped on the admixture graph, and GWAS results. Correspondence analysis of genes, environments, and traits provides a birds-eye view of the history of population differentiation and range expansion and various types of environmental selections at the times. Principal component analysis of the estimated trait-specific polygenic adaptations mapped on the admixture graph enables to understand the coadaptation of multiple traits. The potential usefulness was confirmed by analyzing a public dataset of wild poplar in northwestern America. In response to the northern cold temperature and longer daylength, the populations increased the photosynthetic activity and nutrient use efficiency at the expense of the risk of pathogen invasion, and in response to warm temperature, they increased the growth. At higher altitude, they shifted the maximum activity to earlier period in spring to reduce the activity in dry summer. The R codes for our representation method and simulations of population colonization used in this study are available as supplementary script.
Liver environment-imposed constraints diversify movement strategies of liver-localized CD8 T cells Pathogen-specific CD8 T cells face the problem of finding rare cells that present their cognate antigen either in the lymph node or infected tissue. While quantitative details of T cell movement strategies in some tissues such lymph nodes or skin have been relatively well characterized, we still lack quantitative understanding of T cell movement in many other important tissues such as the spleen, lung, liver, and gut. Furthermore, how environmental constraints influence movement patterns of T cells in tissues remains incompletely characterized. Liver is one of the major organs in mice and humans and is the target site for many pathogens such as hepatitis B and C viruses in humans and Plasmodium parasites in multiple mammalian species. We developed a protocol to generate stable numbers of liver-located CD8 T cells, used intravital microscopy to record movement patterns of CD8 T cells in livers of live mice, and analyzed these and previously published data using well-established statistical and computational methods. We show that in most of our experiments Plasmodium-specific liver-localized CD8 T cells perform correlated random walks characterized by transiently superdiffusive displacement with persistence times of 10-15 min that exceed those observed for T cells in lymph nodes. Liver-localized CD8 T cells typically crawl on the lumenal side of liver sinusoids (i.e., are in the blood); simulating T cell movement in digital structures derived from the liver sinusoids illustrates that liver structure alone is sufficient to explain the relatively long superdiffusive displacement of T cells. In experiments when CD8 T cells in the liver poorly attach to the sinusoids (e.g., one week post immunization with radiation-attenuated Plasmodium sporozoites) T cells also undergo Levy flights - large displacements occurring due to cell deattaching from the endothelium, floating with the blood flow, and re-attaching at another location. Our analysis, thus, provides quantitative details of movement patterns of liver-localized CD8 T cells and illustrates how structural and physiological details of the tissue may impact T cell movement patterns. "Adopt a Microorganism" methodology enhances the discourse richness of high school students in both hybrid learning and remote education The Internet has changed the way teachers and students access information and build knowledge. The recent COVID-19 pandemic has created challenges for both teachers and students and a demand for new methodologies of remote learning. In the life sciences, mixing online content with practical activities represents an even greater challenge. In microbiology, the implementation of an active teaching methodology, the #Adopt project, based on the social network Facebook(R), represents an excellent option for connecting remote education with classroom activities. In 2020, the version applied in high school, "Adopt a Microorganism", was adapted to meet the demands of emergency remote education owing to the suppression of face-to-face activities caused by the pandemic. In the present study, we assessed how the change in methodology impacted the discourse richness of students from high school integrated with technical education in the Business Administration program of the Federal Institute of Sao Paulo, Sorocaba Campus. Three questionnaires related to the groups of microorganisms (Archaea, Bacteria, Virus, Fungi, and Protozoan) were applied. The students responses in the 2019 and 2020 classes were compared concerning content richness and multiplicity of concepts through the application of the Shannon diversity index, an approach that is generally used to assess biodiversity in different environments. The observed results suggest that remote learning provided students with a conceptual basis and richness of content equivalent to that achieved by students subjected to the hybrid teaching model. In conclusion, this study suggests that the #Adopt project methodology increases students discourse richness in microbiology even without face-to-face traditional classes. Distal-less and spalt are distal organisers of pierid wing patterns Two genes, Distal-less (Dll) and spalt (sal), are known to be involved in establishing nymphalid butterfly wing patterns. They function in several ways: in the differentiation of the eyespots central signaling cells, or foci; in the differentiation of the surrounding black disc; in overall scale melanisation (Dll); and in elaborating marginal patterns, such as parafocal elements. However, little is known about the functions of these genes in the development of wing patterns in other butterfly families. Here, we study the expression and function of Dll and sal in the development of spots and other melanic wing patterns of the Indian cabbage white, Pieris canidia, a pierid butterfly. In P. canidia, both Dll and Sal proteins are expressed in the scale-building cells at the wing tips, in chevron patterns along the pupal wing margins, and in areas of future scale melanisation. Additionally, Sal alone is expressed in the future black spots. CRISPR knockouts of Dll and sal showed that each gene is required for the development of melanic wing pattern elements, and repressing pteridine granule formation, in the areas where they are expressed. We conclude that both genes likely play ancestral roles in organising distal butterfly wing patterns, across pierid and nymphalid butterflies, but are unlikely to be differentiating signalling centers in pierids black spots. The genetic and developmental mechanisms that set up the location of spots and eyespots are likely distinct in each lineage. Probing the rules of cell coordination in live tissues by interpretable machine learning based on graph neural networks Robustness in developing and homeostatic tissues is supported by various types of spatiotemporal cell-to-cell interactions. Although live imaging and cell tracking are powerful in providing direct evidence of cell coordination rules, extracting and comparing these rules across many tissues with potentially different length and timescales of coordination requires a versatile framework of analysis. Here we demonstrate that graph neural network (GNN) models are suited for this purpose, by showing how they can be applied to predict cell fate in tissues and utilized to infer the cell interactions governing the multicellular dynamics. Analyzing the live mammalian epidermis data, where spatiotemporal graphs constructed from cell tracks and cell contacts are given as inputs, GNN discovers distinct neighbor cell fate coordination rules that depend on the region of the body. This approach demonstrates how the GNN framework is powerful in inferring general cell interaction rules from live data without prior knowledge of the signaling involved. Universal protein misfolding intermediates can bypass the proteostasis network and remain soluble and less-functional Misfolded protein conformations with decreased functionality can bypass the proteostasis machinery and remain soluble in vivo. This is an unexpected phenomenon as several cellular quality control mechanisms have evolved to rid cells of misfolded proteins. Three questions, then, are: how is it structurally possible for long-lived, soluble, misfolded proteins to bypass the proteostasis machinery and processes? How widespread are these soluble, misfolded states across the proteome? And how long do they persist for? Here, we address these questions using coarse-grain molecular dynamics simulations of the synthesis, termination, and post-translational dynamics of a representative set of cytosolic E. coli proteins. We predict that half of all proteins exhibit subpopulations of misfolded conformations that are likely to bypass molecular chaperones, avoid aggregation, and not be rapidly degraded. These misfolded states may persist for months or longer for some proteins. Structurally characterizing these misfolded states, we observe they have a large amount of native structure, but also contain localized misfolded regions from non-native changes in entanglement, in which a protein segment threads through a loop formed by another portion of the protein that is not found in the native state. The surface properties of these misfolded states are native like, suggesting they may bypass the proteostasis machinery and its regulatory processes to remain soluble, while their entanglements make these states long-lived kinetic traps, as disentanglement requires unfolding of already folded portions of the protein. In terms of function, we predict that one-third of proteins have subpopulations that misfold into less-functional states that have structurally perturbed functional sites yet remain soluble. Data from limited-proteolysis mass spectrometry experiments, which interrogate the misfolded conformations populated by proteins upon unfolding and refolding, are consistent with the structural changes seen in the entangled states of glycerol-3-phosphate dehydrogenase upon misfolding. These results provide an explanation for how proteins can misfold into soluble conformations with reduced functionality that can bypass cellular quality controls, and indicate, unexpectedly, this may be a wide-spread phenomenon in proteomes. Such entanglements are observed in many native structures, suggesting the non-native entanglements we observe are plausible. More broadly, these near-native entangled structures suggest a hypothesis for how synonymous mutations can modulate downstream protein structure and function, with these mutations partitioning nascent proteins between these kinetically trapped states. Modulations in motor unit discharge are related to changes in fascicle length during isometric contractions The integration of electromyography (EMG) and ultrasound imaging has provided important information about the mechanisms of muscle activation and contraction. Unfortunately, EMG does not allow an accurate assessment of the interplay between the neural drive received by muscles, changes in fascicle length (FL) and the force/torque produced. We aimed to assess the relationship between modulations in tibialis anterior (TA) motor unit (MU) firing rate, FL and dorsiflexion torque (DT) using ultrasound-transparent high-density EMG electrodes. EMG and ultrasound images were recorded simultaneously from TA, using a 32-electrode silicon matrix, while performing isometric dorsiflexion, at diverse ankle joint positions (0{degrees} and 30{degrees} plantar flexion) and torques (20% and 40% of maximum). EMG signals were decomposed into individual MUs and changes in FL were assessed with a fascicle-tracking algorithm. MU firings were converted into a cumulative spike train (CST) that was cross-correlated with DT (CST-DT) and FL (CST-FL). High cross-correlations were found for CST-FL, 0.60 (range: 0.31-0.85) and CST-DT 0.71 (range: 0.31-0.88). Cross-correlation lags revealed that the delay between CST-FL (~75ms) was significantly smaller than CST-DT (~150ms, p<0.001). These delays affected the interpretation of MU recruitment/de-recruitment thresholds, with FL showing similar lengths for both recruitment and de-recruitment. This study is the first to demonstrate that changes in TA FL are closely related to both modulations in MU firing frequency and DT. These relationships allow assessment of the interplay between neural drive, muscle contraction and resultant torque, thereby providing a better understanding of the mechanisms responsible for the generation of muscle force.

NEW AND NOTEWORTHYBy employing ultrasound-transparent high-density surface EMG electrodes, we show that modulations in tibialis anterior motor unit discharge rate were closely related to both changes in its fascicle length and resultant torque. These relationships allowed quantifying delays between neural drive and muscle shortening as well as muscle shortening and torque during submaximal isometric contractions, providing an accurate estimation of the time required to generate muscle force and subsequent production of torque via the tendon. Optimization of avian perching manoeuvres Perching at speed is amongst the most challenging flight behaviours that birds perform1,2, and beyond the capability of current autonomous vehicles. Smaller birds may touchdown by hovering3-8, but larger birds typically swoop upward to perch1,2 - presumably because the adverse scaling of their power margin prohibits slow flapping flight9, and because swooping transfers excess kinetic to potential energy1,2,10. Perching is risky in larger birds6,11, demanding precise control of velocity and pose12-15, but it is unknown how they optimize this challenging manoeuvre. More generally, whereas cruising flight behaviours such as migration and commuting are adapted to minimize cost-of-transport or time-of-flight16, the optimization of unsteady flight manoeuvres remains largely unexplored7,17. Here we show that swooping minimizes neither the time nor energy required to perch safely in Harris hawks Parabuteo unicinctus, but instead minimizes the distance flown under hazardous post-stall conditions. By combining motion capture data from 1,563 flights with flight dynamics modelling, we found that the birds choice of where to transition from powered dive to unpowered climb minimizes the distance from the landing perch over which very high lift coefficients are required. Time and energy are therefore invested to maintain the control authority needed to execute a safe landing, rather than being minimized continuously as in technical applications of autonomous perching under nonlinear feedback control13 and deep reinforcement learning18,19. Naive birds learn this behaviour on-the-fly, so our findings suggest an alternative reward function for reinforcement learning of autonomous perching in air vehicles. Large protein complex interfaces have evolved to promote cotranslational assembly Assembly pathways of protein complexes should be precise and efficient to minimise misfolding and unwanted interactions with other proteins in the cell. One way to achieve this efficiency is by seeding assembly pathways during translation via the cotranslational assembly of subunits. While recent evidence suggests that such cotranslational assembly is widespread, little is known about the properties of protein complexes associated with the phenomenon. Here, using a combination of proteome-specific protein complex structures and publicly available ribosome profiling data, we show that cotranslational assembly is particularly common between subunits that form large intermolecular interfaces. To test whether large interfaces have evolved to promote cotranslational assembly, as opposed to cotranslational assembly being a non-adaptive consequence of large interfaces, we compared the sizes of first and last translated interfaces of heteromers in bacterial, yeast, and human complexes. We detect a clear preference in N-terminal interfaces, which will be translated first and thus more likely to form cotranslationally, to be larger than C-terminal interfaces, suggesting that large interfaces have evolved as a means to maximise the chance of successful cotranslational subunit binding. Polycystins are upregulated in fibrotic kidneys and promote renal fibrosis Autosomal polycystic kidney disease (ADPKD) is caused by mutations in PKD1or PKD2 gene, which encodes polycystin-1 or polycystin-2 respectively. Paradoxically, polycystins are up-regulated in human ADPKD kidneys and overexpression of polycystin-1 even induced renal cyst formation in mice, suggesting that up-regulation of polycystin is pathogenic. In this study, we determined the role of polycystin in renal fibrosis and explored its underlying mechanisms. We found that the expression of polycystin-1 or polycystin-2 was positively correlated with the expression of collagen-I in the kidney of three different fibrotic mouse models. Pharmaceutical inhibition or genetic knockout of polycystin-2 reduced the expression of pro-fibrotic markers in fibrotic kidneys of two animal models. Moreover, conditional knockout of Pkd1 gene attenuated renal fibrosis in two mouse models. Enhancer of Zeste homolog 2 (EZH2) is a methyltransferase which plays important role in renal fibrosis. We found that EZH2 was positively correlated with the expression of polycystins in fibrotic kidneys. Conditional knockout of the Ezh2 attenuated the anti-fibrotic responses induced by Pkd1 deletion in fibrotic kidneys.

In conclusion, polycystins are up-regulated in injured kidneys and polycystins promote renal fibrosis through up-regulation of EZH2 in fibrotic kidneys, suggesting that renal fibrosis in polycystic kidneys is not the direct effect of PKD gene mutations.

SignificanceMutations of polycystins cause autosomal polycystic kidney disease (ADPKD) which is associated with increased renal fibrosis in cystic kidneys. However, the role of polycystin in renal fibrosis remains unclear. Here, we showed that the expression of polycystins is positively correlated with the development of renal fibrosis in injured kidneys. Pharmaceutical inhibition or genetic deletion of polycystin-2 attenuated renal fibrosis in fibrotic kidneys. Moreover, genetic deletion of polycystin-1 inhibited renal fibrosis in two different animal models. Finally, we showed that ablation of the methyltransferase EZH2 mitigated the anti-fibrotic effect of polycystin-1 deletion in fibrotic kidneys. We conclude that polycystin promotes through EZH2 and renal fibrosis in cystic kidney is not the direct effect of PKD gene mutations. Genomic regions associate with major axes of variation driven by gas exchange and leaf construction traits in cultivated sunflower (Helianthus annuus L.) Stomata and leaf veins play an essential role in transpiration and the movement of water throughout leaves. These traits are thus thought to play a key role in the adaptation of plants to drought and a better understanding of the genetic basis of their variation and coordination could inform efforts to improve drought tolerance. Here, we explore patterns of variation and covariation in leaf anatomical traits and analyze their genetic architecture via genome-wide association (GWA) analyses in cultivated sunflower (Helianthus annuus L.). Traits related to stomatal density and morphology as well as lower order veins were manually measured from digital images while the density of minor veins was estimated using a novel deep learning approach. Leaf, stomatal, and vein traits exhibited numerous significant correlations that generally followed expectations based on functional relationships. Correlated suites of traits could further be separated along three major principal component (PC) axes that were heavily influenced by variation in traits related to gas exchange, leaf hydraulics, and leaf construction. While there was limited evidence of colocalization when individual traits were subjected to GWA analyses, major multivariate PC axes that were most strongly influenced by several traits related to gas exchange or leaf construction did exhibit significant genomic associations. These results provide insight into the genetic basis of leaf trait covariation and showcase potential targets for future efforts aimed at modifying leaf anatomical traits in sunflower.

Significance StatementUsing traditional and automated/high-throughput (using a novel deep learning approach) phenotyping methods we studied leaf anatomical variation in sunflower. Genome-wide association (GWA) analyses identified numerous genomic regions underlying individual trait variation and regions underlying major multivariate axes of phenotypic variation. These results illustrate the value of employing a multivariate approach to GWA analyses and shed light on the extent to which leaf trait (co-)variation can be genetically decoupled to explore novel phenotypic space. First Report of Lactoria cornuta from Karnataka, West coast of India, Role of tropical cyclones in Range Extension India is known for its rich biodiversity and is fortunate to have several endemic species from different classes of vertebrates. India is home to 7.5% of the global fish diversity, with 91 endemic species of ray finned fishes, the actinopterygians. Some fish species have never been reported until recently, one such example is the long-horned Cowfish (Ostracidae), best known for inhabiting only certain areas of the Indo-Pacific oceans. However, this fish has been reported from Bay of Bengal, Costs of Goa and Lakshadweep. This is the first time this fish has been found on the Coasts of Karnataka, west coast of India. Given its morphology, the migration seems highly improbable since it can only endure lethargic swimming. In this article, we have explored the role of other external forces that could have contributed to its range extension in the east and west coast of India. In this paper we are exploring the role of Super Cyclones as a vital force in determining the expansion of its range. The fact of collection longhorn cowfish for the first time from the coast of Karwar coincides with the passage of Cyclone Ockhi on the west coast of India this suggests the possible role of super cyclone in its range extension. A year later, another specimen was recovered from a location very close to the first, which indicates the successful establishment of Lacturia cornuta in its new environment. Related morphometric and meristic studies of our specimens are consistent with previous studies reported in the Bay of Bengal. A humanized nanobody phage display library yields potent binders of SARS CoV-2 spike Neutralizing antibodies targeting the SARS-CoV-2 spike protein have shown a great preventative/therapeutic potential. Here, we report a rapid and efficient strategy for the development and design of SARS-CoV-2 neutralizing humanized nanobody constructs with sub-nanomolar affinities and nanomolar potencies. CryoEM-based structural analysis of the nanobodies in complex with spike revealed two distinct binding modes. The most potent nanobody, RBD-1-2G(NCATS-BL8125), tolerates the N501Y RBD mutation and remains capable of neutralizing the B.1.1.7 (Alpha) variant. Molecular dynamics simulations provide a structural basis for understanding the neutralization process of nanobodies exclusively focused on the spike-ACE2 interface with and without the N501Y mutation on RBD. A primary human airway air-lung interface (ALI) ex vivo model showed that RBD-1-2G-Fc antibody treatment was effective at reducing viral burden following WA1 and B.1.1.7 SARS-CoV-2 infections. Therefore, this presented strategy will serve as a tool to mitigate the threat of emerging SARS-CoV-2 variants. The 55K SNP-Based Exploration of QTL for Spikelet Number Per Spike in a Tetraploid Wheat (Triticum turgidum L.) Population: Chinese Landrace Ailanmai X Wild Emmer The spikelet number per spike (SNS) is a primary factor determining wheat yield. Common wheat breeding reduces the genetic diversity among elite germplasm resources, leading to a detrimental effect on future wheat production. It is, therefore, necessary to explore new genetic resources for SNS to increase wheat yield. A tetraploid landrace  Ailanmai x wild emmer wheat recombinant inbred line (RIL) population was used to construct a genetic map using the wheat 55K single nucleotide polymorphism (SNP) array. The linkage map containing 1150 bin markers with a total genetic distance of 2411.8 cM was obtained. Based on phenotypic data from eight environments and best linear unbiased prediction (BLUP) values, five quantitative trait loci (QTL) for SNS were identified, explaining 6.71-29.40% of the phenotypic variation. Two of them, QSns.sau-AM-2B.3 and QSns.sau-AM-3B.2, were detected as the major and novel QTL. Their effects were further validated in two additional F2 populations using tightly linked Kompetitive Allele-Specific PCR (KASP) markers. Potential candidate genes within the physical intervals of corresponding QTL were predicted to participate in inflorescence development and spikelet formation. Genetic associations between SNS and other agronomic traits were also detected and analyzed. This study demonstrates the feasibility of wheat 55K SNP array developed for common wheat in genetic mapping of tetraploid population and shows the potential application of wheat related species in wheat improvement programs. Entorhinal cortex epigenome-wide association study highlights four novel loci showing differential methylation in Alzheimer's disease BackgroundStudies on DNA methylation (DNAm) in Alzheimers disease (AD) have recently highlighted several genomic loci showing association with disease onset and progression.

MethodsHere, we conducted an epigenome-wide association study (EWAS) using DNAm profiles in entorhinal cortex (EC) from 149 AD patients and control brains and combined these with two previously published EC datasets by meta-analysis (total n=337).

ResultsWe identified 12 cytosine-phosphate-guanine (CpG) sites showing epigenome-wide significant association with either case-control status or Braaks tau-staging. Four of these CpGs, located in proximity to CNFN/LIPE, TENT5A, PALD1/PRF1, and DIRAS1, represent novel findings. Integrating DNAm levels with RNA sequencing-based mRNA expression data generated in the same individuals showed significant DNAm-mRNA correlations for 6 of the 12 significant CpGs. Lastly, by calculating rates of epigenetic age acceleration using two recently proposed "epigenetic clock" estimators we found a significant association with accelerated epigenetic aging in AD patients vs. controls.

ConclusionIn summary, our study represents the hitherto most comprehensive EWAS in AD using EC and highlights several novel differentially methylated loci with potential effects on gene expression. 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. Analysis of alternative mRNA splicing in vemurafenib-resistant melanoma cells Alternative mRNA splicing is common in cancers. In BRAF V600E mutated malignant melanoma a frequent mechanism of acquired resistance to BRAF inhibitors involves alternative splicing (AS) of BRAF. The resulting shortened BRAF protein constitutively dimerizes and conveys drug resistance. Here, we have analysed AS in SKMEL-239 melanoma cells and a BRAF inhibitor (vemurafenib) resistant derivative that expresses an AS, shortened BRAF V600E transcript. Transcriptome analysis showed differential expression of spliceosome components between the two cell lines. As there is no consensus approach to analysing AS events, we used and compared four common AS softwares based on different principles, DEXSeq, rMATS, ASpli, and LeafCutter. Two of them correctly identified the BRAF V600E AS in the vemurafenib resistant cells. Only 12 AS events were identified by all four softwares. Testing the AS predictions experimentally showed that these overlapping predictions are highly accurate. Interestingly, they identified AS caused alterations in the expression of melanin synthesis and cell migration genes in the vemurafenib resistant cells. This analysis shows that combining different AS analysis approaches produce reliable results and meaningful, biologically testable hypotheses.

Graphic Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=70 SRC="FIGDIR/small/484656v2_ufig1.gif" ALT="Figure 1">
View larger version (24K):
org.highwire.dtl.DTLVardef@9f4920org.highwire.dtl.DTLVardef@179e9bcorg.highwire.dtl.DTLVardef@e05632org.highwire.dtl.DTLVardef@3b1c41_HPS_FORMAT_FIGEXP  M_FIG C_FIG How much does it take to be old? Modelling the time since the last harvesting to infer the distribution of overmature forests in France AimThe distribution of overmature forests in metropolitan France is poorly known, with only a few well-studied prominent sites, and has never been evaluated countrywide. Here, we modelled French forest reserves time since the last harvesting operation - a proxy for forest maturity - then inferred the current statistical distribution of overmature forests (i.e. forests over 50 years without harvesting) in France.

LocationMetropolitan France

MethodsWe used inventories from forest reserves and managed forests to calibrate a generalised linear mixed model explaining the time since the last harvesting with selected structural attributes and environmental variables. We then projected this model on the independent National Forest Inventory dataset. We thus obtained an updated estimation of the proportion and a rough distribution of overmature forest stands in metropolitan France.

ResultsWe found that high basal area of very large trees, high volumes of standing and downed deadwood, high diversity of tree-related microhabitats and more marginally diversity of decay stages best characterized the time since the last harvesting. Volumes of stumps and high density of coppices translating legacy of past forest management also distinguished more overmature plots. Our projection yielded an estimated 3% of French forests over 50 years without harvesting mostly located in more inaccessible areas (i.e. mountainous areas) and a promising proportion of future overmature forests if left unharvested.

Main conclusionsOur study showed that the time since the last harvesting is a good proxy for a combination of stand structure attributes key in characterising overmature temperate forests. It gives the first robust statistical estimate of the proportion of overmature forests and may serve to report on their status in metropolitan France. Our method could be implemented at a larger spatial scale, notably in countries with accessible National Forest Inventory and calibration data, to produce indicators at international level. ZAF -- The First Open Source Fully Automated Feeder for Aquatic Facilities In the past few decades, aquatic animals have become popular model organisms in biology, spurring a growing need for establishing aquatic facilities. Zebrafish are widely studied and relatively easy to culture using commercial systems. However, a challenging aspect of maintaining aquatic facilities is animal feeding, which is both time- and resource-consuming. We have developed an open-source fully automatic daily feeding system, Zebrafish Automatic Feeder (ZAF). ZAF is reliable, provides a standardized amount of food to every tank, is cost-efficient and easy to build. The advanced version, ZAF+, allows for the precise control of food distribution as a function of fish density per tank, and has a user-friendly interface. Both ZAF and ZAF+ are adaptable to any laboratory environment and facilitate the implementation of aquatic colonies. Here we provide all blueprints and instructions for building the mechanics, electronics, fluidics, as well as to setup the control software and its user-friendly graphical interface. Importantly, the design is modular and can be scaled to meet different user needs. Furthermore, our results show that ZAF and ZAF+ do not adversely affect zebrafish culture, enabling fully automatic feeding for any aquatic facility. Production and characterization of virus-free, CRISPR-CAR T cells capable of inducing solid tumor regression BackgroundChimeric antigen receptor (CAR) T cells traditionally harbor viral vector-based sequences that encode the CAR transgene in the genome. These T cell products have yet to show consistent anti-tumor activity in patients with solid tumors. Further, viral vector manufacturing is resource intensive, suffers from batch-to-batch variability, and includes several animal components, adding regulatory and supply chain pressures.

MethodsAnti-GD2 CAR T cells were generated using CRISPR/Cas9 within nine days using recombinant Cas9 protein and nucleic acids, without any viral vectors or animal components. The CAR was specifically targeted to the T Cell Receptor Alpha Constant gene (TRAC). T cell products were characterized at the level of the genome, transcriptome, proteome, and secretome using CHANGE-seq, scRNA-seq, spectral cytometry, and ELISA assays. Functionality was evaluated in vivo in an NSG xenograft neuroblastoma model.

ResultsIn comparison to traditional retroviral CAR T cells, virus-free CRISPR CAR (VFC-CAR) T cells exhibit TRAC-targeted genomic integration of the CAR transgene, elevation of transcriptional and protein characteristics associated with a memory phenotype, and low tonic signaling prior to infusion arising in part from the the knockout of the TCR. Upon exposure to the GD2 target antigen, anti-GD2 VFC-CAR T cells exhibited specific cytotoxicity against GD2+ cells in vitro and induced solid tumor regression in vivo, with robust homing, persistence, and low exhaustion against a human neuroblastoma xenograft model.

ConclusionsThis proof-of-principle study leveraging virus-free genome editing technology could enable flexible manufacturing of clinically relevant, high-quality CAR T cells to treat cancers, including solid tumors. 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. Beneficial effects of novel strains of Aureobasidium pullulans produced 1,3-1,6 β-glucans on non-esterified fatty acid levels in diabetic KKAy mice Obesity, metabolic syndrome, associated lipotoxicity and its cascade of events contribute to the majority of the burden related to non-communicable diseases globally. Preventive lifestyle changes aside, several beneficial effects have been reported in type II diabetes mellitus and dyslipidaemia patients with biological response modifier glucans (BRMG) produced as an exopolysaccharide by Aureobasidium pullulans. In this study, we compared two strains (AFO-202 and N-163) that produce beta glucans in alleviating lipotoxicity. This study was performed in obese diabetic mice model of KK-Ay mice, in four groups with six subjects in each group - Group 1: sacrificed on Day 0 for baseline values; Group 2: control (drinking water); Group 3: AFO-202 beta glucan--200 mg/kg/day; Group 4: N-163 beta glucan--300 mg/kg/day. The animals in groups 2-4 had the test solutions administered by gavage once daily for 28 consecutive days. Biochemical analyses were conducted of blood glucose, triglycerides, total cholesterol, LDL cholesterol, HDL cholesterol and non-esterified fatty acids (NEFA). Group 4 (N-163) had the lowest NEFA levels, as compared to the other groups, and marginally decreased triglyceride levels. The groups had no significant differences in blood glucose, HbA1c, triglycerides, or LDL and HDL cholesterol. N-163 produced by A. pullulans decreased NEFA in a diabetic mice model in 28 days. These results, although modest, warrant further in-depth research into lipotoxicity and associated inflammatory cascades in both healthy and disease affected subjects to develop novel strategies for prevention and management. Fast ATP-dependent Subunit Rotation in reconstituted FoF1-ATP synthase trapped in Solution FoF1-ATP synthases are ubiquitous membrane-bound, rotary motor enzymes that can catalyze ATP synthesis and hydrolysis. Their enzyme kinetics are controlled by internal subunit rotation, by substrate and product concentrations, by mechanical inhibitory mechanisms, but also by the electrochemical potential of protons across the membrane. Single-molecule Forster resonance energy transfer (smFRET) has been used to detect subunit rotation within FoF1-ATP synthases embedded in freely diffusing liposomes. We now report that kinetic monitoring of functional rotation can be prolonged from milliseconds to seconds by utilizing an Anti-Brownian electrokinetic trap (ABEL trap). These extended observation times allowed us to observe fluctuating rates of functional rotation for individual FoF1-liposomes in solution. Broad distributions of ATP-dependent catalytic rates were revealed. The buildup of an electrochemical potential of protons was confirmed to limit the maximum rate of ATP hydrolysis. In the presence of ionophores or uncouplers, the fastest subunit rotation speeds measured in single reconstituted FoF1-ATP synthases were 180 full rounds per second. This was much faster than measured by biochemical ensemble averaging, but not as fast as the maximum rotational speed reported previously for isolated single F1 complexes uncoupled from the membrane-embedded Fo complex. Further application of ABEL trap measurements should help resolve the mechanistic causes of such fluctuating rates of subunit rotation.

TOC GRAPHIC

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=108 SRC="FIGDIR/small/436299v5_ufig1.gif" ALT="Figure 1">
View larger version (35K):
org.highwire.dtl.DTLVardef@1bd800org.highwire.dtl.DTLVardef@228aeorg.highwire.dtl.DTLVardef@1bce39forg.highwire.dtl.DTLVardef@3f01a_HPS_FORMAT_FIGEXP  M_FIG O_FLOATNOFigure 1.C_FLOATNO C_FIG Phylogenomic evidence for the Origin of Obligately Anaerobic Anammox Bacteria around the Great Oxidation Event The anaerobic ammonium oxidation (anammox) bacteria could transform ammonium and nitrite to dinitrogen gas, and this obligate anaerobic process accounts for up to half of the global nitrogen loss in surface environments. Yet its origin and evolution, which may give important insights into the biogeochemistry of early Earth, remains enigmatic. Here, we performed comprehensive phylogenomic analysis and showed a single origin of anammox bacteria within the phylum Planctomycetes. After accommodating the uncertainties and factors influencing time estimates, which includes implementing both a traditional cyanobacteria-based and a recently developed mitochondria-based approach, we estimated that anammox bacteria originated at early Proterozoic and most likely around the so-called Great Oxidation Event (GOE; 2.32 to 2.5 billion years ago [Ga]) which fundamentally changed global biogeochemical cycles. We further showed that during the origin of anammox bacteria, genes involved in oxidative stress, bioenergetics and anammox granules formation were recruited, which might have contributed to their survival on an increasingly oxic Earth. Our findings suggest the rising levels of atmospheric oxygen, which made nitrite increasingly available, was a potential driving force for the emergence of anammox bacteria. This is one of the first studies that link the GOE to the evolution of obligate anaerobic bacteria. An enhanced inverted encoding model for neural reconstructions Here we present a more interpretable and flexible approach for reconstructing the contents of perception, attention, and memory from neuroimaging data. Our enhanced inverted encoding model (eIEM) incorporates methodological improvements including proper accounting of population-level tuning functions and a trial-by-trial prediction error-based metric where reconstruction quality is measured in meaningful units. Improved flexibility is further gained via eIEMs novel goodness-of-fit feature: for trial-by-trial reconstructions, goodness-of-fits are obtained independently (non-circularly) to prediction error and can be applied to any IEM procedure or decoding metric, resulting in improved reconstruction quality and brain-behavior correlations. We validate eIEM from methodological principles, simulated neuroimaging datasets, and three pre-existing fMRI datasets spanning perception, attention, and working memory. Notably, eIEM is easy to apply and broadly accessible - our Python package (https://pypi.org/project/inverted-encoding) implements eIEM in one line of code - and is easily modifiable to compare performance metrics and/or scale up to more complex models.