Structure-based discovery of CFTR potentiators and inhibitors.

The cystic fibrosis transmembrane conductance regulator (CFTR) is a crucial ion channel whose loss of function leads to cystic fibrosis, whereas its hyperactivation leads to secretory diarrhea. Small molecules that improve CFTR folding (correctors) or function (potentiators) are clinically available. However, the only potentiator, ivacaftor, has suboptimal pharmacokinetics and inhibitors have yet to be clinically developed. Here, we combine molecular docking, electrophysiology, cryo-EM, and medicinal chemistry to identify CFTR modulators. We docked ∼155 million molecules into the potentiator site on CFTR, synthesized 53 test ligands, and used structure-based optimization to identify candidate modulators. This approach uncovered mid-nanomolar potentiators, as well as inhibitors, that bind to the same allosteric site. These molecules represent potential leads for the development of more effective drugs for cystic fibrosis and secretory diarrhea, demonstrating the feasibility of large-scale docking for ion channel drug discovery.

Unexplored microbial diversity from 2,500 food metagenomes and links with the human microbiome.

Complex microbiomes are part of the food we eat and influence our own microbiome, but their diversity remains largely unexplored. Here, we generated the open access curatedFoodMetagenomicData (cFMD) resource by integrating 1,950 newly sequenced and 583 public food metagenomes. We produced 10,899 metagenome-assembled genomes spanning 1,036 prokaryotic and 108 eukaryotic species-level genome bins (SGBs), including 320 previously undescribed taxa. Food SGBs displayed significant microbial diversity within and between food categories. Extension to >20,000 human metagenomes revealed that food SGBs accounted on average for 3% of the adult gut microbiome. Strain-level analysis highlighted potential instances of food-to-gut transmission and intestinal colonization (e.g., Lacticaseibacillus paracasei) as well as SGBs with divergent genomic structures in food and humans (e.g., Streptococcus gallolyticus and Limosilactobabillus mucosae). The cFMD expands our knowledge on food microbiomes, their role in shaping the human microbiome, and supports future uses of metagenomics for food quality, safety, and authentication.

Long-term mesoscale imaging of 3D intercellular dynamics across a mammalian organ.

A comprehensive understanding of physio-pathological processes necessitates non-invasive intravital three-dimensional (3D) imaging over varying spatial and temporal scales. However, huge data throughput, optical heterogeneity, surface irregularity, and phototoxicity pose great challenges, leading to an inevitable trade-off between volume size, resolution, speed, sample health, and system complexity. Here, we introduce a compact real-time, ultra-large-scale, high-resolution 3D mesoscope (RUSH3D), achieving uniform resolutions of 2.6 × 2.6 × 6 µm3 across a volume of 8,000 × 6,000 × 400 µm3 at 20 Hz with low phototoxicity. Through the integration of multiple computational imaging techniques, RUSH3D facilitates a 13-fold improvement in data throughput and an orders-of-magnitude reduction in system size and cost. With these advantages, we observed premovement neural activity and cross-day visual representational drift across the mouse cortex, the formation and progression of multiple germinal centers in mouse inguinal lymph nodes, and heterogeneous immune responses following traumatic brain injury-all at single-cell resolution, opening up a horizon for intravital mesoscale study of large-scale intercellular interactions at the organ level.

Multiplexed 3D atlas of state transitions and immune interaction in colorectal cancer.

Advanced solid cancers are complex assemblies of tumor, immune, and stromal cells characterized by high intratumoral variation. We use highly multiplexed tissue imaging, 3D reconstruction, spatial statistics, and machine learning to identify cell types and states underlying morphological features of known diagnostic and prognostic significance in colorectal cancer. Quantitation of these features in high-plex marker space reveals recurrent transitions from one tumor morphology to the next, some of which are coincident with long-range gradients in the expression of oncogenes and epigenetic regulators. At the tumor invasive margin, where tumor, normal, and immune cells compete, T cell suppression involves multiple cell types and 3D imaging shows that seemingly localized 2D features such as tertiary lymphoid structures are commonly interconnected and have graded molecular properties. Thus, while cancer genetics emphasizes the importance of discrete changes in tumor state, whole-specimen imaging reveals large-scale morphological and molecular gradients analogous to those in developing tissues.

High-precision coding in visual cortex.

Individual neurons in visual cortex provide the brain with unreliable estimates of visual features. It is not known whether the single-neuron variability is correlated across large neural populations, thus impairing the global encoding of stimuli. We recorded simultaneously from up to 50,000 neurons in mouse primary visual cortex (V1) and in higher order visual areas and measured stimulus discrimination thresholds of 0.35° and 0.37°, respectively, in an orientation decoding task. These neural thresholds were almost 100 times smaller than the behavioral discrimination thresholds reported in mice. This discrepancy could not be explained by stimulus properties or arousal states. Furthermore, behavioral variability during a sensory discrimination task could not be explained by neural variability in V1. Instead, behavior-related neural activity arose dynamically across a network of non-sensory brain areas. These results imply that perceptual discrimination in mice is limited by downstream decoders, not by neural noise in sensory representations.

Changes of Cell Biochemical States Are Revealed in Protein Homomeric Complex Dynamics.

We report here a simple and global strategy to map out gene functions and target pathways of drugs, toxins, or other small molecules based on "homomer dynamics" protein-fragment complementation assays (hdPCA). hdPCA measures changes in self-association (homomerization) of over 3,500 yeast proteins in yeast grown under different conditions. hdPCA complements genetic interaction measurements while eliminating the confounding effects of gene ablation. We demonstrate that hdPCA accurately predicts the effects of two longevity and health span-affecting drugs, the immunosuppressant rapamycin and the type 2 diabetes drug metformin, on cellular pathways. We also discovered an unsuspected global cellular response to metformin that resembles iron deficiency and includes a change in protein-bound iron levels. This discovery opens a new avenue to investigate molecular mechanisms for the prevention or treatment of diabetes, cancers, and other chronic diseases of aging.

Target-Specific Precision of CRISPR-Mediated Genome Editing.

The CRISPR-Cas9 system has successfully been adapted to edit the genome of various organisms. However, our ability to predict the editing outcome at specific sites is limited. Here, we examined indel profiles at over 1,000 genomic sites in human cells and uncovered general principles guiding CRISPR-mediated DNA editing. We find that precision of DNA editing (i.e., recurrence of a specific indel) varies considerably among sites, with some targets showing one highly preferred indel and others displaying numerous infrequent indels. Editing precision correlates with editing efficiency and a preference for single-nucleotide homologous insertions. Precise targets and editing outcome can be predicted based on simple rules that mainly depend on the fourth nucleotide upstream of the protospacer adjacent motif (PAM). Indel profiles are robust, but they can be influenced by chromatin features. Our findings have important implications for clinical applications of CRISPR technology and reveal general patterns of broken end joining that can provide insights into DNA repair mechanisms.

Multifractality in spin glasses.

We unveil the multifractal behavior of Ising spin glasses in their low-temperature phase. Using the Janus II custom-built supercomputer, the spin-glass correlation function is studied locally. Dramatic fluctuations are found when pairs of sites at the same distance are compared. The scaling of these fluctuations, as the spin-glass coherence length grows with time, is characterized through the computation of the singularity spectrum and its corresponding Legendre transform. A comparatively small number of site pairs controls the average correlation that governs the response to a magnetic field. We explain how this scenario of dramatic fluctuations (at length scales smaller than the coherence length) can be reconciled with the smooth, self-averaging behavior that has long been considered to describe spin-glass dynamics.

Multiscale Flow for robust and optimal cosmological analysis.

We propose Multiscale Flow, a generative Normalizing Flow that creates samples and models the field-level likelihood of two-dimensional cosmological data such as weak lensing. Multiscale Flow uses hierarchical decomposition of cosmological fields via a wavelet basis and then models different wavelet components separately as Normalizing Flows. The log-likelihood of the original cosmological field can be recovered by summing over the log-likelihood of each wavelet term. This decomposition allows us to separate the information from different scales and identify distribution shifts in the data such as unknown scale-dependent systematics. The resulting likelihood analysis can not only identify these types of systematics, but can also be made optimal, in the sense that the Multiscale Flow can learn the full likelihood at the field without any dimensionality reduction. We apply Multiscale Flow to weak lensing mock datasets for cosmological inference and show that it significantly outperforms traditional summary statistics such as power spectrum and peak counts, as well as machine learning-based summary statistics such as scattering transform and convolutional neural networks. We further show that Multiscale Flow is able to identify distribution shifts not in the training data such as baryonic effects. Finally, we demonstrate that Multiscale Flow can be used to generate realistic samples of weak lensing data.

Space wandering in the rodent default mode network.

The default mode network (DMN) is a large-scale brain network known to be suppressed during a wide range of cognitive tasks. However, our comprehension of its role in naturalistic and unconstrained behaviors has remained elusive because most research on the DMN has been conducted within the restrictive confines of MRI scanners. Here, we use multisite GCaMP (a genetically encoded calcium indicator) fiber photometry with simultaneous videography to probe DMN function in awake, freely exploring rats. We examined neural dynamics in three core DMN nodes-the retrosplenial cortex, cingulate cortex, and prelimbic cortex-as well as the anterior insula node of the salience network, and their association with the rats' spatial exploration behaviors. We found that DMN nodes displayed a hierarchical functional organization during spatial exploration, characterized by stronger coupling with each other than with the anterior insula. Crucially, these DMN nodes encoded the kinematics of spatial exploration, including linear and angular velocity. Additionally, we identified latent brain states that encoded distinct patterns of time-varying exploration behaviors and found that higher linear velocity was associated with enhanced DMN activity, heightened synchronization among DMN nodes, and increased anticorrelation between the DMN and anterior insula. Our findings highlight the involvement of the DMN in collectively and dynamically encoding spatial exploration in a real-world setting. Our findings challenge the notion that the DMN is primarily a "task-negative" network disengaged from the external world. By illuminating the DMN's role in naturalistic behaviors, our study underscores the importance of investigating brain network function in ecologically valid contexts.

A Guide to Emerging Technologies for Large-Scale and Whole-Brain Optical Imaging of Neuronal Activity.

The mammalian brain is a densely interconnected network that consists of millions to billions of neurons. Decoding how information is represented and processed by this neural circuitry requires the ability to capture and manipulate the dynamics of large populations at high speed and high resolution over a large area of the brain. Although the use of optical approaches by the neuroscience community has rapidly increased over the past two decades, most microscopy approaches are unable to record the activity of all neurons comprising a functional network across the mammalian brain at relevant temporal and spatial resolutions. In this review, we survey the recent development in optical technologies for Ca2+ imaging in this regard and provide an overview of the strengths and limitations of each modality and its potential for scalability. We provide guidance from the perspective of a biological user driven by the typical biological applications and sample conditions. We also discuss the potential for future advances and synergies that could be obtained through hybrid approaches or other modalities.

Biolinguistic graph fusion model for circRNA-miRNA association prediction.

Emerging clinical evidence suggests that sophisticated associations with circular ribonucleic acids (RNAs) (circRNAs) and microRNAs (miRNAs) are a critical regulatory factor of various pathological processes and play a critical role in most intricate human diseases. Nonetheless, the above correlations via wet experiments are error-prone and labor-intensive, and the underlying novel circRNA-miRNA association (CMA) has been validated by numerous existing computational methods that rely only on single correlation data. Considering the inadequacy of existing machine learning models, we propose a new model named BGF-CMAP, which combines the gradient boosting decision tree with natural language processing and graph embedding methods to infer associations between circRNAs and miRNAs. Specifically, BGF-CMAP extracts sequence attribute features and interaction behavior features by Word2vec and two homogeneous graph embedding algorithms, large-scale information network embedding and graph factorization, respectively. Multitudinous comprehensive experimental analysis revealed that BGF-CMAP successfully predicted the complex relationship between circRNAs and miRNAs with an accuracy of 82.90% and an area under receiver operating characteristic of 0.9075. Furthermore, 23 of the top 30 miRNA-associated circRNAs of the studies on data were confirmed in relevant experiences, showing that the BGF-CMAP model is superior to others. BGF-CMAP can serve as a helpful model to provide a scientific theoretical basis for the study of CMA prediction.

Farmland practices are driving bird population decline across Europe.

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.

Large-scale land acquisitions exacerbate local farmland inequalities in Tanzania.

Land inequality stalls economic development, entrenches poverty, and is associated with environmental degradation. Yet, rigorous assessments of land-use interventions attend to inequality only rarely. A land inequality lens is especially important to understand how recent large-scale land acquisitions (LSLAs) affect smallholder and indigenous communities across as much as 100 million hectares around the world. This paper studies inequalities in land assets, specifically landholdings and farm size, to derive insights into the distributional outcomes of LSLAs. Using a household survey covering four pairs of land acquisition and control sites in Tanzania, we use a quasi-experimental design to characterize changes in land inequality and subsequent impacts on well-being. We find convincing evidence that LSLAs in Tanzania lead to both reduced landholdings and greater farmland inequality among smallholders. Households in proximity to LSLAs are associated with 21.1% (P = 0.02) smaller landholdings while evidence, although insignificant, is suggestive that farm sizes are also declining. Aggregate estimates, however, hide that households in the bottom quartiles of farm size suffer the brunt of landlessness and land loss induced by LSLAs that combine to generate greater farmland inequality. Additional analyses find that land inequality is not offset by improvements in other livelihood dimensions, rather farm size decreases among households near LSLAs are associated with no income improvements, lower wealth, increased poverty, and higher food insecurity. The results demonstrate that without explicit consideration of distributional outcomes, land-use policies can systematically reinforce existing inequalities.

Data-driven large-scale genomic analysis reveals an intricate phylogenetic and functional landscape in J-domain proteins.

The 70-kD heat shock protein (Hsp70) chaperone system is a central hub of the proteostasis network that helps maintain protein homeostasis in all organisms. The recruitment of Hsp70 to perform different and specific cellular functions is regulated by the J-domain protein (JDP) co-chaperone family carrying the small namesake J-domain, required to interact and drive the ATPase cycle of Hsp70s. Besides the J-domain, prokaryotic and eukaryotic JDPs display a staggering diversity in domain architecture, function, and cellular localization. Very little is known about the overall JDP family, despite their essential role in cellular proteostasis, development, and its link to a broad range of human diseases. In this work, we leverage the exponentially increasing number of JDP gene sequences identified across all kingdoms owing to the advancements in sequencing technology and provide a broad overview of the JDP repertoire. Using an automated classification scheme based on artificial neural networks (ANNs), we demonstrate that the sequences of J-domains carry sufficient discriminatory information to reliably recover the phylogeny, localization, and domain composition of the corresponding full-length JDP. By harnessing the interpretability of the ANNs, we find that many of the discriminatory sequence positions match residues that form the interaction interface between the J-domain and Hsp70. This reveals that key residues within the J-domains have coevolved with their obligatory Hsp70 partners to build chaperone circuits for specific functions in cells.

Census of exposed aggregation-prone regions in proteomes.

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.

CIForm as a Transformer-based model for cell-type annotation of large-scale single-cell RNA-seq data.

Single-cell omics technologies have made it possible to analyze the individual cells within a biological sample, providing a more detailed understanding of biological systems. Accurately determining the cell type of each cell is a crucial goal in single-cell RNA-seq (scRNA-seq) analysis. Apart from overcoming the batch effects arising from various factors, single-cell annotation methods also face the challenge of effectively processing large-scale datasets. With the availability of an increase in the scRNA-seq datasets, integrating multiple datasets and addressing batch effects originating from diverse sources are also challenges in cell-type annotation. In this work, to overcome the challenges, we developed a supervised method called CIForm based on the Transformer for cell-type annotation of large-scale scRNA-seq data. To assess the effectiveness and robustness of CIForm, we have compared it with some leading tools on benchmark datasets. Through the systematic comparisons under various cell-type annotation scenarios, we exhibit that the effectiveness of CIForm is particularly pronounced in cell-type annotation. The source code and data are available at https://github.com/zhanglab-wbgcas/CIForm.

Empirically adjusted fixed-effects meta-analysis methods in genomic studies.

In recent years, meta-analyzing summary results from multiple studies has become a common practice in genomic research, leading to a significant improvement in the power of statistical detection compared to an individual genomic study. Meta analysis methods that combine statistical estimates across studies are known to be statistically more powerful than those combining statistical significance measures. An approach combining effect size estimates based on a fixed-effects model, called METAL, has gained extreme popularity to perform the former type of meta-analysis. In this article, we discuss the limitations of METAL due to its dependence on the theoretical null distribution, leading to incorrect significance testing results. Through various simulation studies and real genomic data application, we show how modifying the z-scores in METAL, using an empirical null distribution, can significantly improve the results, especially in presence of hidden confounders. For the estimation of the null distribution, we consider two different approaches, and we highlight the scenarios when one null estimation approach outperforms the other. This article will allow researchers to gain an insight into the importance of using an empirical null distribution in the fixed-effects meta-analysis as well as in choosing the appropriate empirical null distribution estimation approach.

Large-scale parameters framework with large convolutional kernel for encoding visual fMRI activity information.

Visual encoding models often use deep neural networks to describe the brain's visual cortex response to external stimuli. Inspired by biological findings, researchers found that large receptive fields built with large convolutional kernels improve convolutional encoding model performance. Inspired by scaling laws in recent years, this article investigates the performance of large convolutional kernel encoding models on larger parameter scales. This paper proposes a large-scale parameters framework with a sizeable convolutional kernel for encoding visual functional magnetic resonance imaging activity information. The proposed framework consists of three parts: First, the stimulus image feature extraction module is constructed using a large-kernel convolutional network while increasing channel numbers to expand the parameter size of the framework. Second, enlarging the input data during the training stage through the multi-subject fusion module to accommodate the increase in parameters. Third, the voxel mapping module maps from stimulus image features to functional magnetic resonance imaging signals. Compared to sizeable convolutional kernel visual encoding networks with base parameter scale, our visual encoding framework improves by approximately 7% on the Natural Scenes Dataset, the dedicated dataset for the Algonauts 2023 Challenge. We further analyze that our encoding framework made a trade-off between encoding performance and trainability. This paper confirms that expanding parameters in visual coding can bring performance improvements.

Facilely tuning the intrinsic catalytic sites of the spinel oxide for peroxymonosulfate activation: From fundamental investigation to pilot-scale demonstration.

Heterogeneous peroxymonosulfate (PMS)-based advanced oxidation processes (AOPs) have shown a great potential for pollutant degradation, but their feasibility for large-scale water treatment application has not been demonstrated. Herein, we develop a facile coprecipitation method for the scalable production (∼10 kg) of the Cu-Fe-Mn spinel oxide (CuFeMnO). Such a catalyst has rich oxygen vacancies and symmetry-breaking sites, which endorse it with a superior PMS-catalytic capacity. We find that the working reactive species and their contributions are highly dependent on the properties of target organic pollutants. For the organics with electron-donating group (e.g., -OH), high-valent metal species are mainly responsible for the pollutant degradation, whereas for the organics with electron-withdrawing group (e.g., -COOH and -NO2), hydroxyl radical (•OH) as the secondary oxidant also plays an important role. We demonstrate that the CuFeMnO-PMS system is able to achieve efficient and stable removal of the pollutants in the secondary effluent from a municipal wastewater plant at both bench and pilot scales. Moreover, we explore the application prospect of this PMS-based AOP process for large-scale wastewater treatment. This work describes an opportunity to scalably prepare robust spinel oxide catalysts for water purification and is beneficial to the practical applications of the heterogeneous PMS-AOPs.