Search

Single-molecule imaging in living cells is an effective tool for elucidating the mechanisms of cellular phenomena at the molecular level. However, the analysis was not designed for throughput and requires high expertise, preventing it from reaching large scale, which is necessary when searching for rare cells that induce singularity phenomena. To overcome this limitation, we have automated the imaging procedures by combining our own focusing device, artificial intelligence, and robotics. The apparatus, called automated in-cell single-molecule imaging system (AiSIS), achieves a throughput that is a hundred-fold higher than conventional manual imaging operations, enabling the analysis of molecular events by individual cells across a large population. Here, using AiSIS, we demonstrate the single-molecule imaging of molecular behaviors and reactions related to tau protein aggregation, which is considered a singularity phenomenon in neurological disorders. Changes in the dynamics and kinetics of molecular events were observed inside and on the basal membrane of cells after the induction of aggregation. Additionally, to detect rare cells based on the molecular behavior, we developed a method to identify the state of individual cells defined by the quantitative distribution of molecular mobility and clustering. Using this method, cellular variations in receptor behavior were shown to decrease following ligand stimulation. This cell state analysis based on large-scale single-molecule imaging by AiSIS will advance the study of molecular mechanisms causing singularity phenomena., Competing Interests: The authors declare no competing interests., (2024 THE BIOPHYSICAL SOCIETY OF JAPAN.)

Over large regions exposed to natural disasters, cascading effects resulting from complex or concatenated natural processes may represent a large portion of total risk. Populated high-mountain environments are a major concern, and methods for large-scale quantitative risk analyses are urgently required to improve risk mitigation. This article presents a comprehensive quantitative rockfall risk assessment over a large archetypal valley of the Andean mountains, in Central Chile, which integrates a wide spectrum of elements at risk. Risk is expressed as an expected damage both in monetary terms and casualties, at different scales relevant for decision making. Notably, total rockfall risk is divided into its main drivers, which allows quantifying seismically induced rockfall risk. For this purpose, the local seismic hazard is quantified and the yield acceleration, that is, acceleration required to initiate rockfall, is determined at the regional scale. The probability of failure is thereafter derived in terms of annual frequency of rockfall initiation and integrated in the quantitative risk assessment (QRA) process. Our results show the significant role of seismic activity as the triggering mechanism of rockfalls, and highlight elements at risk that have a major contribution to the total risk. Eventually a sensitivity analysis is conducted to (i) assess the robustness of obtained risk estimates to the data and modeling choices and (ii) identify the most influential assumptions. Our approach evidences the feasibility of large-scale QRAs in sensitive environments and opens perspectives for refining QRAs in similar territories significantly affected by cascading effects and multihazards., (© 2023 Society for Risk Analysis.)

Zea mays L is a crucial crop for Brazil, ranking second in terms of production and sixth in terms of exports. In Brazil, the second season, or off-season, accounts for 80 % of the overall maize output, which primarily occurs after the soybean main season. A maize yield forecast model for the off-season was developed and implemented throughout Brazilian territory due to its importance to the country's economy and food security. The model was built using multiple linear regressions that connected outputs simulated from a land surface model used in large-scale analysis for agriculture (JULES-crop), to agrometeorological indicators. The application of the developed model occurred every 10 days from the sowing until the maturation. A comparison of the forecasting model was verified with the official off-season maize yields for the years 2003-2016. Agrometeorological indicators during the reproductive phase accounted for 60 % of the interannual variability in maize production. When outputs simulated by JULES-crop were included, the forecasting model achieved Nash-Sutcliffe modeling efficiency (EF) of 0.77 in the maturation and EF = 0.72 in the filling-grain stage, suggesting that this approach can generate useful predictions for final maize yield beginning on the 80th day of the cycle. Outputs of JULES crop enhanced modeling performance during the vegetative stage, reducing the standard deviation error in prediction from 0.59 to 0.49 Mg ha -1 ., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2024 The Authors.)

The auxin response factor (ARF) and auxin/indole-3-acetic acid (Aux/IAA) family of genes are central components of the auxin signaling pathway and play essential roles in plant growth and development. Their large-scale analysis and evolutionary trajectory of origin are currently not known. Here, we identified the corresponding ARF and Aux/IAA family members and performed a large-scale analysis by scanning 406 plant genomes. The results showed that the ARF and Aux/IAA gene families originated from charophytes. The ARF family sequences were more conserved than the Aux/IAA family sequences. Dispersed duplications were the common expansion mode of ARF and Aux/IAA families in bryophytes, ferns, and gymnosperms; however, whole-genome duplication was the common expansion mode of the ARF and Aux/IAA families in basal angiosperms, magnoliids, monocots, and dicots. Expression and regulatory network analyses revealed that the Arabidopsis thaliana ARF and Aux/IAA families responded to multiple hormone, biotic, and abiotic stresses. The APETALA2 and serum response factor-transcription factor gene families were commonly enriched in the upstream and downstream genes of the ARF and Aux/IAA gene families. Our study provides a comprehensive overview of the evolutionary trajectories, structural functions, expansion mechanisms, expression patterns, and regulatory networks of these two gene families., (© 2024. The Author(s).)

State-space models are widely employed across various research disciplines to study unobserved dynamics. Conventional estimation techniques, such as Kalman filtering and expectation maximisation, offer valuable insights but incur high computational costs in large-scale analyses. Sparse inverse covariance estimators can mitigate these costs, but at the expense of a trade-off between enforced sparsity and increased estimation bias, necessitating careful assessment in low signal-to-noise ratio (SNR) situations. To address these challenges, we propose a three-fold solution: (1) Introducing multiple penalised state-space (MPSS) models that leverage data-driven regularisation; (2) Developing novel algorithms derived from backpropagation, gradient descent, and alternating least squares to solve MPSS models; (3) Presenting a K-fold cross-validation extension for evaluating regularisation parameters. We validate this MPSS regularisation framework through lower and more complex simulations under varying SNR conditions, including a large-scale synthetic magneto- and electro-encephalography (MEG/EEG) data analysis. In addition, we apply MPSS models to concurrently solve brain source localisation and functional connectivity problems for real event-related MEG/EEG data, encompassing thousands of sources on the cortical surface. The proposed methodology overcomes the limitations of existing approaches, such as constraints to small-scale and region-of-interest analyses. Thus, it may enable a more accurate and detailed exploration of cognitive brain functions., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or relations that could have appeared to influence the work reported in this paper., (Copyright © 2023. Published by Elsevier Inc.)

Transfer RNA-derived fragments (tRFs) are noncoding RNAs that arise from either mature transfer RNAs (tRNAs) or their precursors. One important category of tRFs comprises the tRNA halves, which are generated through cleavage at the anticodon. A given tRNA typically gives rise to several co-expressed 5'-tRNA halves (5'-tRHs) that differ in the location of their 3' ends. These 5'-tRHs, even though distinct, have traditionally been treated as indistinguishable from one another due to their near-identical sequences and lengths. We focused on co-expressed 5'-tRHs that arise from the same tRNA and systematically examined their exact sequences and abundances across 10 different human tissues. To this end, we manually curated and analyzed several hundred human RNA-seq datasets from NCBI's Sequence Run Archive (SRA). We grouped datasets from the same tissue into their own collection and examined each group separately. We found that a given tRNA produces different groups of co-expressed 5'-tRHs in different tissues, different cell lines, and different diseases. Importantly, the co-expressed 5'-tRHs differ in their sequences, absolute abundances, and relative abundances, even among tRNAs with near-identical sequences from the same isodecoder or isoacceptor group. The findings suggest that co-expressed 5'-tRHs that are produced from the same tRNA or closely related tRNAs have distinct, context-dependent roles. Moreover, our analyses show that cell lines modeling the same tissue type and disease may not be interchangeable when it comes to experimenting with tRFs.

The SARS-CoV-2 pandemic had huge impacts on global urban populations, activity and health, yet little is known about attendant consequences for urban river ecosystems. We detected significant changes in occurrence and risks from contaminants of emerging concern (CECs) in waterways across Greater London (UK) during the pandemic. We were able to rapidly identify and monitor large numbers of CECs in n = 390 samples across 2019-2021 using novel direct-injection liquid chromatography-mass spectrometry methods for scalable targeted analysis, suspect screening and prioritisation of CEC risks. A total of 10,029 measured environmental concentrations (MECs) were obtained for 66 unique CECs. Pharmaceutical MECs decreased during lockdown in 2020 in the R. Thames (p ≤ 0.001), but then increased significantly in 2021 (p ≤ 0.01). For the tributary rivers, the R. Lee, Beverley Brook, R. Wandle and R. Hogsmill were the most impacted, primarily via wastewater treatment plant effluent and combined sewer overflows. In the R. Hogsmill in particular, pharmaceutical MEC trends were generally correlated with NHS prescription statistics, likely reflecting limited wastewater dilution. Suspect screening of ∼ 1,200 compounds tentatively identified 25 additional CECs at the five most impacted sites, including metabolites such as O-desmethylvenlafaxine, an EU Watch List compound. Lastly, risk quotients (RQs) ≥ 0.1 were calculated for 21 compounds across the whole Greater London freshwater catchment, of which seven were of medium risk (RQ ≥ 1.0) and three were in the high-risk category (RQ ≥ 10), including imidacloprid (RQ = 19.6), azithromycin (15.7) and diclofenac (10.5). This is the largest spatiotemporal dataset of its kind for any major capital city globally and the first for Greater London, representing ∼ 16 % of the population of England, and delivering a foundational One-Health case study in the third largest city in Europe across a global pandemic., Competing Interests: Declaration of Competing Interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Leon Barron reports financial support was provided by Imperial College London., (Copyright © 2023 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Loss of solubility usually leads to the detrimental elimination of protein function. In some cases, the protein aggregation is also required for beneficial functions. Given the duality of this phenomenon, it remains a fundamental question how natural selection controls the aggregation. The exponential growth of genomic sequence data and recent progress with in silico predictors of the aggregation allows approaching this problem by a large-scale bioinformatics analysis. Most of the aggregation-prone regions are hidden within the 3D structure, rendering them inaccessible for the intermolecular interactions responsible for aggregation. Thus, the most realistic census of the aggregation-prone regions requires crossing aggregation prediction with information about the location of the natively unfolded regions. This allows us to detect so-called 'exposed aggregation-prone regions' (EARs). Here, we analyzed the occurrence and distribution of the EARs in 76 reference proteomes from the three kingdoms of life. For this purpose, we used a bioinformatics pipeline, which provides a consensual result based on several predictors of aggregation. Our analysis revealed a number of new statistically significant correlations about the presence of EARs in different organisms, their dependence on protein length, cellular localizations, co-occurrence with short linear motifs and the level of protein expression. We also obtained a list of proteins with the conserved aggregation-prone sequences for further experimental tests. Insights gained from this work led to a deeper understanding of the relationship between protein evolution and aggregation., (© The Author(s) 2023. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Declines in European bird populations are reported for decades but the direct effect of major anthropogenic pressures on such declines remains unquantified. Causal relationships between pressures and bird population responses are difficult to identify as pressures interact at different spatial scales and responses vary among species. Here, we uncover direct relationships between population time-series of 170 common bird species, monitored at more than 20,000 sites in 28 European countries, over 37 y, and four widespread anthropogenic pressures: agricultural intensification, change in forest cover, urbanisation and temperature change over the last decades. We quantify the influence of each pressure on population time-series and its importance relative to other pressures, and we identify traits of most affected species. We find that agricultural intensification, in particular pesticides and fertiliser use, is the main pressure for most bird population declines, especially for invertebrate feeders. Responses to changes in forest cover, urbanisation and temperature are more species-specific. Specifically, forest cover is associated with a positive effect and growing urbanisation with a negative effect on population dynamics, while temperature change has an effect on the dynamics of a large number of bird populations, the magnitude and direction of which depend on species' thermal preferences. Our results not only confirm the pervasive and strong effects of anthropogenic pressures on common breeding birds, but quantify the relative strength of these effects stressing the urgent need for transformative changes in the way of inhabiting the world in European countries, if bird populations shall have a chance of recovering.

Activity cliffs (AC) are formed by pairs of structural analogues that are active against the same target but have a large difference in potency. While much of our knowledge about ACs has originated from the analysis and comparison of compounds and activity data, several studies have reported AC predictions over the past decade. Different from typical compound classification tasks, AC predictions must be carried out at the level of compound pairs representing ACs or nonACs. Most AC predictions reported so far have focused on individual methods or comparisons of two or three approaches and only investigated a few compound activity classes (from 2 to 10). Although promising prediction accuracy has been reported in most cases, different system set-ups, AC definitions, methods, and calculation conditions were used, precluding direct comparisons of these studies. Therefore, we have carried out a large-scale AC prediction campaign across 100 activity classes comparing machine learning methods of greatly varying complexity, ranging from pair-based nearest neighbor classifiers and decision tree or kernel methods to deep neural networks. The results of our systematic predictions revealed the level of accuracy that can be expected for AC predictions across many different compound classes. In addition, prediction accuracy did not scale with methodological complexity but was significantly influenced by memorization of compounds shared by different ACs or nonACs. In many instances, limited training data were sufficient for building accurate models using different methods and there was no detectable advantage of deep learning over simpler approaches for AC prediction. On a global scale, support vector machine models performed best, by only small margins compared to others including simple nearest neighbor classifiers., (© 2023. The Author(s).)

Alternative splicing is an important means of generating the protein diversity necessary for cellular functions. Hence, there is a growing interest in assessing the structural and functional impact of alternative protein isoforms. Typically, experimental studies are used to determine the structures of the canonical proteins ignoring the other isoforms. Therefore, there is still a large gap between abundant sequence information and meager structural data on these isoforms. During the last decade, significant progress has been achieved in the development of bioinformatics tools for structural and functional annotations of proteins. Moreover, the appearance of the AlphaFold program opened up the possibility to model a large number of high-confidence structures of the isoforms. In this study, using state-of-the-art tools, we performed in silico analysis of 58 eukaryotic proteomes. The evaluated structural states included structured domains, intrinsically disordered regions, aggregation-prone regions, and tandem repeats. Among other things, we found that the isoforms have fewer signal peptides, transmembrane regions, or tandem repeat regions in comparison with their canonical counterparts. This could change protein function and/or cellular localization. The AlphaFold modeling demonstrated that frequently isoforms, having differences with the canonical sequences, still can fold in similar structures though with significant structural rearrangements which can lead to changes of their functions. Based on the modeling, we suggested classification of the structural differences between canonical proteins and isoforms. Altogether, we can conclude that a majority of isoforms, similarly to the canonical proteins are under selective pressure for the functional roles.

Glycans are involved in many fundamental cellular processes such as growth, differentiation, and morphogenesis. However, their broad structural diversity makes analysis difficult. Glycomics via mass spectrometry has focused on the composition of glycans, but informatics analysis has not kept pace with the development of instrumentation and measurement techniques. We developed Toolbox Accelerating Glycomics (TAG), in which glycans can be added manually to the glycan list that can be freely designed with labels and sialic acid modifications, and fast processing is possible. In the present work, we improved TAG for large-scale analysis such as cohort analysis of serum samples. The sialic acid linkage-specific alkylamidation (SALSA) method converts differences in linkages such as α2,3- and α2,6-linkages of sialic acids into differences in mass. Glycans modified by SALSA and several structures discovered in recent years were added to the glycan list. A routine to generate calibration curves has been implemented to explore quantitation. These improvements are based on redefinitions of residues and glycans in the TAG List to incorporate information on glycans that could not be attributed because it was not assumed in the previous version of TAG. These functions were verified through analysis of purchased sera and 74 spectra with linearity at the level of R2 > 0.8 with 81 estimated glycan structures obtained including some candidate of rare glycans such as those with the N,N’-diacetyllactosediamine structure, suggesting they can be applied to large-scale analyses.

Tigecycline is one of important antimicrobial agents for the treatment of infections caused by multidrug-resistant (MDR) Gram-negative bacteria. However, the emergence and prevalence of plasmid-mediated tigecycline resistance gene tet (X4) are threatening human and animal health. Fitness cost elicited by resistance plasmids is a key factor affecting the maintenance and transmission of antibiotic resistance genes (ARGs) in the host. A comparative analysis of the fitness cost of different types of tet (X4)-positive plasmids is helpful to understand and predict the prevalence of dominant plasmids. In this study, we performed a large-scale analysis of fitness cost of tet (X4)-positive plasmids origin from clinical isolates. These plasmids were successfully electroporated into a reference strain Escherichia coli TOP10, and a series of transformants carrying the tet (X) gene were obtained. The effects of tet (X4)-positive plasmids on the growth rate, plasmid stability, relative fitness, biofilm formation, and virulence in a Galleria mellonella model were evaluated. Consequently, we found that these plasmids resulted in varying degrees of fitness cost on TOP10, including delayed bacterial growth and attenuated virulence. Out of these plasmids, tet (X4)-harboring IncFII plasmids showed the lowest fitness cost on the host. Furthermore, by means of experimental evolution in the presence of commonly used drugs in clinic, the fitness cost of tet (X4)-positive plasmids was substantially alleviated, accompanied by increased plasmid stability. Collectively, our data reveal the differential fitness cost caused by different types of tet (X4)-positive plasmids and suggest that the wide use of tetracycline antibiotics may promote the evolution of plasmids., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Tang, Cai, Jiang, Wang and Liu.)

X chromosome inactivation is a mammalian dosage compensation mechanism, where one of two X chromosomes is randomly inactivated in female cells. Previous studies have suggested that primed human embryonic stem cells (hESCs) maintain an eroded state of the X chromosome and do not express XIST , while in naïve transition, both XIST and the eroded X chromosome are reactivated. However, the pattern of chromosome X reactivation in naïve hESCs remains mainly unknown. In this study, we examine the variations in the status of X chromosome between primed and naïve hESCs by analyzing RNA sequencing samples from different studies. We show that most samples of naïve hESCs indeed reactivate XIST and there is an increase in gene expression levels on chromosome X. However, most of the naïve samples do not fully activate chromosome X in a uniform manner and present a distinct eroded pattern, probably as a result of XIST reactivation and initiation of re-inactivation of chromosome X. This large-scale analysis provides a higher-resolution description of the changes occurring in chromosome X during primed-to-naïve transition and emphasizes the importance of taking these variations into consideration when studying X inactivation in embryonic development.

Intensification of agriculture and increased insecticide use have been implicated in global losses of farmland biodiversity and ecosystem services. We hypothesized that increased insecticide applications (proportion of area treated with insecticides) in Canada's expansive agricultural landscapes are due, in part, to shifts toward more simplified landscapes. To assess this relationship, we analyzed data from the Canadian Census of Agriculture spanning 20 years including five census periods (1996-2016) and across 225 census units within the four major agricultural regions of Pacific, Prairie, Central, and Atlantic Canada. Generalized mixed effects models were used to evaluate if changes in landscape simplification - defined as the proportion of farmland in crops (cereals, oilseeds, pulses and fruit/vegetables) - alongside other farming and climatic variables, influenced insecticide applications over time. Bayesian spatial-temporal models were further used to estimate the strength of the relationship with landscape simplification over time. We found that landscape simplification increased in 89% and insecticide applications increased in 70% of the Census Division spatial units during the 1996-2016 period. Nationally, significant increases in landscape simplification were observed in the two most agriculturally intensive regions of Prairie (from 55% to 63%) and Central (from 51% to 60%) Canada. For both regions, landscape simplification was a strong and significant predictor of higher insecticide applications, even after accounting for other factors such as climate, farm economics, farm size and land use practices (e.g., area in cash crops and tillage). If current trends continue, we estimated that insecticide applications will increase another 10%-20% by 2036 as a result of landscape simplification alone. To avoid increased reliance on toxic insecticides, agri-environmental policies need to consider that losing diverse natural habitat can increase insect pest pressure and resistance with negative environmental consequences extending beyond the field., (© 2022 The Ecological Society of America.)

Australia responded to the emergence of the COVID-19 global pandemic in 2020 by initiating a border and bio-security zone lockdown and policies emphasising social distancing and hand hygiene. To understand the public response to this, Southern Queensland Rural Health commenced a two-phase research project exploring attitudes and practices towards the COVID-19 pandemic in Australia. An initial online survey foreshadowed 90 qualitative interviews with respondents to explore what the pandemic meant for everyday life. This paper details use of a qualitative approach by a national collaborative of investigators from 9 rural university departments in Australia who came together to research the qualitative phase of the project. Our methodological approach aligned with extant literature describing the management of large-scale interviewing and coding in the context of unfolding and dynamic contexts. The 'RITA' model (Rapid Identification of Themes from Audio recordings) entails a five-step process designed to progress from identifying research foci, through deductive and iterative coding to identify key concepts. We used a combination of coding templates, organisation and tagging of field notes and real-time sharing through a secure cloud drive to create a data set for immersive analysis and generation of ideas. Use of this method has added to the collective knowledge about successful rapid research investigations, recognising the inherent tension between speed and rigour. This is not a binary but a dialectic; trustworthiness is integral to qualitative research. However, use of fresh approaches is accommodated by new technologies and can preserve adequate rigour while enabling collaboration, research capacity building and increasing the pace of data collection and analysis. This project has presented methodological challenges and highlights some strengths of such an approach. It is hoped that reporting our approach and experiences is useful for the broader health and research community considering large-scale qualitative research., Competing Interests: Declaration of Conflicting Interests: The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article., (© The Author(s) 2021.)

Diffuse gliomas are a heterogeneous group of tumors with aggressive biological behavior and a lack of effective treatment methods. Despite new molecular findings, the differences between pathohistological types still require better understanding. In this in silico analysis, we investigated AKT1 , AKT2 , AKT3 , CHUK , GSK3β , EGFR , PTEN , and PIK3AP1 as participants of EGFR-PI3K-AKT-mTOR signaling using data from the publicly available cBioPortal platform. Integrative large-scale analyses investigated changes in copy number aberrations (CNA), methylation, mRNA transcription and protein expression within 751 samples of diffuse astrocytomas, anaplastic astrocytomas and glioblastomas. The study showed a significant percentage of CNA in PTEN (76%), PIK3AP1 and CHUK (75% each), EGFR (74%), AKT2 (39%), AKT1 (32%), AKT3 (19%) and GSK3β (18%) in the total sample. Comprehensive statistical analyses show how genomics and epigenomics affect the expression of examined genes differently across various pathohistological types and grades, suggesting that genes AKT3 , CHUK and PTEN behave like tumor suppressors, while AKT1 , AKT2 , EGFR , and PIK3AP1 show oncogenic behavior and are involved in enhanced activity of the EGFR-PI3K-AKT-mTOR signaling pathway. Our findings contribute to the knowledge of the molecular differences between pathohistological types and ultimately offer the possibility of new treatment targets and personalized therapies in patients with diffuse gliomas.

Antibody repertoire sequencing enables researchers to acquire millions of B cell receptors and investigate these molecules at the single-nucleotide level. This power and resolution in studying humoral responses have led to its wide applications. However, most of these studies were conducted with a limited number of samples. Given the extraordinary diversity, assessment of these key features with a large sample set is demanded. Thus, we collect and systematically analyze 2,152 high-quality heavy-chain antibody repertoires. Our study reveals that 52 core variable genes universally contribute to more than 99% of each individual's repertoire; a distal interspersed preferences characterize V gene recombination; the number of public clones between two repertoires follows a linear model, and the positive selection dominates at RGYW motif in somatic hypermutations. Thus, this population-level analysis resolves some critical features of the antibody repertoire and may have significant value to the large cadre of scientists., Competing Interests: Declaration of interests An institutional affiliation (Nanfang Hospital, Southern Medical University, Guangzhou 510515, China) has patents on library preparation methods for HTS sequencing (China patents CN201910468844.1, CN201910468845.6, and CN201910468857.9). The authors declare no other competing interests., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

High-speed analysis of large (prote)omics sample sets at the rate of thousands or millions of samples per day on a single platform has been a challenge since the beginning of proteomics. For many years, ESI-based MS methods have dominated proteomics because of their high sensitivity and great depth in analyzing complex proteomes. However, despite improvements in speed, ESI-based MS methods are fundamentally limited by their sample introduction, which excludes off-line sample preparation/fractionation because of the time required to switch between individual samples/sample fractions, and therefore being dependent on the speed of on-line sample preparation methods such as liquid chromatography. Laser-based ionization methods have the advantage of moving from one sample to the next without these limitations, being mainly restricted by the speed of modern sample stages, i.e. 10 ms or less between samples. This speed matches the data acquisition speed of modern high-performing mass spectrometers whereas the pulse repetition rate of the lasers (>1 kHz) provides a sufficient number of desorption/ionization events for successful ion signal detection from each sample at the above speed of the sample stages. Other advantages of laser-based ionization methods include the generally higher tolerance to sample additives and contamination compared with ESI MS, and the contact-less and pulsed nature of the laser used for desorption, reducing the risk of cross-contamination. Furthermore, new developments in MALDI have expanded its analytical capabilities, now being able to fully exploit high-performing hybrid mass analyzers and their strengths in sensitivity and MS/MS analysis by generating an ESI-like stable yield of multiply charged analyte ions. Thus, these new developments and the intrinsically high speed of laser-based methods now provide a good basis for tackling extreme sample analysis speed in the omics., Competing Interests: Conflict of interest—Authors declare no competing interests., (© 2020 Cramer.)

Single-molecule imaging analysis has been applied to study the dynamics and kinetics of molecular behaviors and interactions in living cells. In spite of its high potential as a technique to investigate the molecular mechanisms of cellular phenomena, single-molecule imaging analysis has not been extended to a large scale of molecules in cells due to the low measurement throughput as well as required expertise. To overcome these problems, we have automated the imaging processes by using computer operations, robotics and artificial intelligence (AI). AI is an ideal substitute for expertise to obtain high-quality images for quantitative analysis. Our automated in-cell single-molecule imaging system, AiSIS, could analyze 1600 cells in 1 day, which corresponds to ∼ 100-fold higher efficiency than manual analysis. The large-scale analysis revealed cell-to-cell heterogeneity in the molecular behavior, which had not been recognized in previous studies. An analysis of the receptor behavior and downstream signaling was accomplished within a significantly reduced time frame and revealed the detailed activation scheme of signal transduction, advancing cell biology research. Furthermore, by combining the high-throughput analysis with our previous finding that a receptor changes its behavioral dynamics depending on the presence of a ligand/agonist or inhibitor/antagonist, we show that AiSIS is applicable to comprehensive pharmacological analysis such as drug screening. This AI-aided automation has wide applications for single-molecule analysis., (© The Author(s) 2020. Published by Oxford University Press on behalf of The Japanese Society of Microscopy. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Background: Endothelial cells (ECs) make up the innermost layer throughout the entire vasculature. Their phenotypes and physiological functions are initially regulated by developmental signals and extracellular stimuli. The underlying molecular mechanisms responsible for the diverse phenotypes of ECs from different organs are not well understood., Results: To characterize the transcriptomic and epigenomic landscape in the vascular system, we cataloged gene expression and active histone marks in nine types of human ECs (generating 148 genome-wide datasets) and carried out a comprehensive analysis with chromatin interaction data. We developed a robust procedure for comparative epigenome analysis that circumvents variations at the level of the individual and technical noise derived from sample preparation under various conditions. Through this approach, we identified 3765 EC-specific enhancers, some of which were associated with disease-associated genetic variations. We also identified various candidate marker genes for each EC type. We found that the nine EC types can be divided into two subgroups, corresponding to those with upper-body origins and lower-body origins, based on their epigenomic landscape. Epigenomic variations were highly correlated with gene expression patterns, but also provided unique information. Most of the deferentially expressed genes and enhancers were cooperatively enriched in more than one EC type, suggesting that the distinct combinations of multiple genes play key roles in the diverse phenotypes across EC types. Notably, many homeobox genes were differentially expressed across EC types, and their expression was correlated with the relative position of each organ in the body. This reflects the developmental origins of ECs and their roles in angiogenesis, vasculogenesis and wound healing., Conclusions: This comprehensive analysis of epigenome characterization of EC types reveals diverse transcriptional regulation across human vascular systems. These datasets provide a valuable resource for understanding the vascular system and associated diseases.

In the living cells, the majority of proteins does not work alone, but interact with other proteins or other biomolecules to maintain the cellular function, constituting a "protein community". Previous efforts on mass spectroscopy-based protein interaction networks, interactomes, have provided a picture on the protein community. However, these were static information after cells were disrupted. For a better understanding of the protein community in cells, it is important to know the properties of intracellular dynamics and interactions. Since hydrodynamic size and mobility of proteins are related into such properties, direct measurement of diffusional motion of proteins in single living cells will be helpful for uncovering the properties. Here we completed measurement of the diffusion and homo-oligomeric properties of 369 cytoplasmic GFP-fusion proteins in living yeast Saccharomyces cerevisiae cells using fluorescence correlation spectroscopy (FCS). The large-scale analysis showed that the motions of majority of proteins obeyed a two-component (i.e. slow and fast components) diffusion model. Remarkably, both of the two components diffused more slowly than expected monomeric states. In addition, further analysis suggested that more proteins existed as homo-oligomeric states in living cells than previously expected. Our study, which characterizes the dynamics of proteins in living cells on a large-scale, provided a global view on intracellular protein dynamics to understand the protein community., (Copyright © 2019 Elsevier Inc. All rights reserved.)

River impoundment constitutes one of the most important anthropogenic impacts on the World's rivers. An increasing number of studies have tried to quantify the effects of river impoundment on riverine ecosystems over the past two decades, often focusing on the effects of individual large reservoirs. This study is one of the first to use a large-scale, multi-year diatom dataset from a routine biomonitoring network to analyse sample sites downstream of a large number of water supply reservoirs ( n  = 77) and to compare them with paired unregulated control sites. We analysed benthic diatom assemblage structure and a set of derived indices, including ecological guilds, in tandem with multiple spatio-temporal variables to disclose patterns of ecological responses to reservoirs beyond the site-specific scale. Diatom assemblage structure at sites downstream of water supply reservoirs was significantly different to control sites, with the effect being most evident at the regional scale. We found that regional influences were important drivers of differences in assemblage structure at the national scale, although this effect was weaker at downstream sites, indicating the homogenising effect of river impoundment on diatom assemblages. Sites downstream of reservoirs typically exhibited a higher taxonomic richness, with the strongest increases found within the motile guild. In addition, Trophic Diatom Index (TDI) values were typically higher at downstream sites. Water quality gradients appeared to be an important driver of diatom assemblages, but the influence of other abiotic factors could not be ruled out and should be investigated further. Our results demonstrate the value of diatom assemblage data from national-scale biomonitoring networks to detect the effects of water supply reservoirs on instream communities at large spatial scales. This information may assist water resource managers with the future implementation of mitigation measures such as setting environmental flow targets., Competing Interests: Michael J. Dunbar and Libby Greenway are employed by the Environment Agency of England., (©2019 Krajenbrink et al.)

Inferring transcriptional direction (orientation) of the CRISPR array is essential for many applications, including systematically investigating non-canonical CRISPR/Cas functions. The standard method, CRISPRDirection (embedded within CRISPRCasFinder), fails to predict the orientation (ND predictions) for ∼37% of the classified CRISPR arrays (>2200 loci); this goes up to >70% for the II-B subtype where non-canonical functions were first experimentally discovered. Alternatively, Potential Orientation (also embedded within CRISPRCasFinder), has a much smaller frequency of ND predictions but might have significantly lower accuracy. We propose a novel simple criterion, where the CRISPR array direction is assigned according to the direction of its associated cas genes (Cas Orientation). We systematically assess the performance of the three methods (Cas Orientation, CRISPRDirection, and Potential Orientation) across all CRISPR/Cas subtypes, by a mutual crosscheck of their predictions, and by comparing them to the experimental dataset. Interestingly, CRISPRDirection agrees much better with Cas Orientation than with Potential Orientation, despite CRISPRDirection and Potential Orientation being mutually related - Potential Orientation corresponding to one of six (heterogeneous) predictors employed by CRISPRDirection - and being unrelated to Cas Orientation. We find that Cas Orientation has much higher accuracy compared to Potential Orientation and comparable accuracy to CRISPRDirection - while accurately assigning an orientation to ∼95% of the CRISPR arrays that are non-determined by CRISPRDirection. Cas Orientation is, at the same time, simple to employ, requiring only (routine for prokaryotes) the prediction of the associated protein coding gene direction.

A multiple reaction monitoring (MRM) assay was developed for precise quantitation of 87 plasma proteins including the three isoforms of apolipoprotein E (APOE) associated with cardiovascular diseases using nanoscale liquid chromatography separation and stable isotope dilution strategy. The analytical performance of the assay was evaluated and we found an average technical variation of 4.7% in 4-5 orders of magnitude dynamic range (≈0.2 mg/L to 4.5 g/L) from whole plasma digest. Here, we report a complete workflow, including sample processing adapted to 96-well plate format and normalization strategy for large-scale studies. To further investigate the MS-based quantitation the amount of six selected proteins was measured by routinely used clinical chemistry assays as well and the two methods showed excellent correlation with high significance (p-value < 10e-5) for the six proteins, in addition for the cardiovascular predictor factor, APOB: APOA1 ratio (r = 0.969, p-value < 10e-5). Moreover, we utilized the developed assay for screening of biobank samples from patients with myocardial infarction and performed the comparative analysis of patient groups with STEMI (ST- segment elevation myocardial infarction), NSTEMI (non ST- segment elevation myocardial infarction) and type-2 AMI (type-2 myocardial infarction) patients.

Background: Promiscuity in molecular interactions between small-molecules, including drugs, and proteins is widespread. Such unintended interactions can be exploited to suggest drug repurposing possibilities as well as to identify potential molecular mechanisms responsible for observed side-effects., Methods: We perform a large-scale analysis to detect binding-site molecular interaction field similarities between the binding-sites of the primary target of 400 drugs against a dataset of 14082 cavities within 7895 different proteins representing a non-redundant dataset of all proteins with known structure. Statistically-significant cases with high levels of similarities represent potential cases where the drugs that bind the original target may in principle bind the suggested off-target. Such cases are further analysed with docking simulations to verify if indeed the drug could, in principle, bind the off-target. Diverse sources of data are integrated to associated potential cross-reactivity targets with side-effects., Results: We observe that promiscuous binding-sites tend to display higher levels of hydrophobic and aromatic similarities. Focusing on the most statistically significant similarities (Z-score ≥ 3.0) and corroborating docking results (RMSD < 2.0 Å), we find 2923 cases involving 140 unique drugs and 1216 unique potential cross-reactivity protein targets. We highlight a few cases with a potential for drug repurposing (acetazolamide as a chorismate pyruvate lyase inhibitor, raloxifene as a bacterial quorum sensing inhibitor) as well as to explain the side-effects of zanamivir and captopril. A web-interface permits to explore the detected similarities for each of the 400 binding-sites of the primary drug targets and visualise them for the most statistically significant cases., Conclusions: The detection of molecular interaction field similarities provide the opportunity to suggest drug repurposing opportunities as well as to identify potential molecular mechanisms responsible for side-effects. All methods utilized are freely available and can be readily applied to new query binding-sites. All data is freely available and represents an invaluable source to identify further candidates for repurposing and suggest potential mechanisms responsible for side-effects.

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) analysis can detect protein/DNA-binding and histone-modification sites across an entire genome. Recent advances in sequencing technologies and analyses enable us to compare hundreds of samples simultaneously; such large-scale analysis has potential to reveal the high-dimensional interrelationship level for regulatory elements and annotate novel functional genomic regions de novo. Because many experimental considerations are relevant to the choice of a method in a ChIP-seq analysis, the overall design and quality management of the experiment are of critical importance. This review offers guiding principles of computation and sample preparation for ChIP-seq analyses, highlighting the validity and limitations of the state-of-the-art procedures at each step. We also discuss the latest challenges of single-cell analysis that will encourage a new era in this field., (© The Author 2016. Published by Oxford University Press.)

A wide variety of large-scale data have been produced in bioinformatics. In response, the need for efficient handling of biomedical big data has been partly met by parallel computing. However, the time demand of many bioinformatics programs still remains high for large-scale practical uses because of factors that hinder acceleration by parallelization. Recently, new generations of storage devices have emerged, such as NAND flash-based solid-state drives (SSDs), and with the renewed interest in near-data processing, they are increasingly becoming acceleration methods that can accompany parallel processing. In certain cases, a simple drop-in replacement of hard disk drives by SSDs results in dramatic speedup. Despite the various advantages and continuous cost reduction of SSDs, there has been little review of SSD-based profiling and performance exploration of important but time-consuming bioinformatics programs. For an informative review, we perform in-depth profiling and analysis of 23 key bioinformatics programs using multiple types of devices. Based on the insight we obtain from this research, we further discuss issues related to design and optimize bioinformatics algorithms and pipelines to fully exploit SSDs. The programs we profile cover traditional and emerging areas of importance, such as alignment, assembly, mapping, expression analysis, variant calling and metagenomics. We explain how acceleration by parallelization can be combined with SSDs for improved performance and also how using SSDs can expedite important bioinformatics pipelines, such as variant calling by the Genome Analysis Toolkit and transcriptome analysis using RNA sequencing. We hope that this review can provide useful directions and tips to accompany future bioinformatics algorithm design procedures that properly consider new generations of powerful storage devices., (© The Author 2015. Published by Oxford University Press. For Permissions, please email: journals.permissions@oup.com.)

A microRNA (miRNA) is a small non-coding RNA molecule, functioning in transcriptional and post-transcriptional regulation of gene expression. The human genome may encode over 1000 miRNAs. Albeit poorly characterized, miRNAs are widely deemed as important regulators of biological processes. Aberrant expression of miRNAs has been observed in many cancers and other disease states, indicating that they are deeply implicated with these diseases, particularly in carcinogenesis. Therefore, it is important for both basic research and miRNA-based therapy to discriminate the real pre-miRNAs from the false ones (such as hairpin sequences with similar stem-loops). Particularly, with the avalanche of RNA sequences generated in the post-genomic age, it is highly desired to develop computational sequence-based methods for effectively identifying the human pre-miRNAs. Here, we propose a predictor called "iMiRNA-PseDPC", in which the RNA sequences are formulated by a novel feature vector called "pseudo distance-pair composition" (PseDPC) with 10 types of structure statuses. Rigorous cross-validations on a much larger and more stringent newly constructed benchmark data-set showed that our approach has remarkably outperformed the existing ones in either prediction accuracy or efficiency, indicating the new predictor is quite promising or at least may become a complementary tool to the existing predictors in this area. For the convenience of most experimental scientists, a user-friendly web server for the new predictor has been established at http://bioinformatics.hitsz.edu.cn/iMiRNA-PseDPC/, by which users can easily get their desired results without the need to go through the mathematical details. It is anticipated that the new predictor may become a useful high throughput tool for genome analysis particularly in dealing with large-scale data.

Proteins must fold into their native structures in the crowded cellular environment, to perform their functions. Although such macromolecular crowding has been considered to affect the folding properties of proteins, large-scale experimental data have so far been lacking. Here, we individually translated 142 Escherichia coli cytoplasmic proteins using a reconstituted cell-free translation system in the presence of macromolecular crowding reagents (MCRs), Ficoll 70 or dextran 70, and evaluated the aggregation propensities of 142 proteins. The results showed that the MCR effects varied depending on the proteins, although the degree of these effects was modest. Statistical analyses suggested that structural parameters were involved in the effects of the MCRs. Our dataset provides a valuable resource to understand protein folding and aggregation inside cells.