Systematic functional characterization of antisense eRNA of protocadherin α composite enhancer.

Enhancers generate bidirectional noncoding enhancer RNAs (eRNAs) that may regulate gene expression. At present, the eRNA function remains enigmatic. Here, we report a 5' capped antisense eRNA PEARL (Pcdh eRNA associated with R-loop formation) that is transcribed from the protocadherin (Pcdh) α HS5-1 enhancer region. Through loss- and gain-of-function experiments with CRISPR/Cas9 DNA fragment editing, CRISPRi, and CRISPRa, as well as locked nucleic acid strategies, in conjunction with ChIRP, MeDIP, DRIP, QHR-4C, and HiChIP experiments, we found that PEARL regulates Pcdhα gene expression by forming local RNA-DNA duplexes (R-loops) in situ within the HS5-1 enhancer region to promote long-distance chromatin interactions between distal enhancers and target promoters. In particular, increased levels of eRNA PEARL via perturbing transcription elongation factor SPT6 lead to strengthened local three-dimensional chromatin organization within the Pcdh superTAD. These findings have important implications regarding molecular mechanisms by which the HS5-1 enhancer regulates stochastic Pcdhα promoter choice in single cells in the brain.

Unraveling the mesoscale organization induced by network-driven processes.

Complex systems are characterized by emergent patterns created by the nontrivial interplay between dynamical processes and the networks of interactions on which these processes unfold. Topological or dynamical descriptors alone are not enough to fully embrace this interplay in all its complexity, and many times one has to resort to dynamics-specific approaches that limit a comprehension of general principles. To address this challenge, we employ a metric-that we name Jacobian distance-which captures the spatiotemporal spreading of perturbations, enabling us to uncover the latent geometry inherent in network-driven processes. We compute the Jacobian distance for a broad set of nonlinear dynamical models on synthetic and real-world networks of high interest for applications from biological to ecological and social contexts. We show, analytically and computationally, that the process-driven latent geometry of a complex network is sensitive to both the specific features of the dynamics and the topological properties of the network. This translates into potential mismatches between the functional and the topological mesoscale organization, which we explain by means of the spectrum of the Jacobian matrix. Finally, we demonstrate that the Jacobian distance offers a clear advantage with respect to traditional methods when studying human brain networks. In particular, we show that it outperforms classical network communication models in explaining functional communities from structural data, therefore highlighting its potential in linking structure and function in the brain.

Prevention efficacy of the broadly neutralizing antibody VRC01 depends on HIV-1 envelope sequence features.

In the Antibody Mediated Prevention (AMP) trials (HVTN 704/HPTN 085 and HVTN 703/HPTN 081), prevention efficacy (PE) of the monoclonal broadly neutralizing antibody (bnAb) VRC01 (vs. placebo) against HIV-1 acquisition diagnosis varied according to the HIV-1 Envelope (Env) neutralization sensitivity to VRC01, as measured by 80% inhibitory concentration (IC80). Here, we performed a genotypic sieve analysis, a complementary approach to gaining insight into correlates of protection that assesses how PE varies with HIV-1 sequence features. We analyzed HIV-1 Env amino acid (AA) sequences from the earliest available HIV-1 RNA-positive plasma samples from AMP participants diagnosed with HIV-1 and identified Env sequence features that associated with PE. The strongest Env AA sequence correlate in both trials was VRC01 epitope distance that quantifies the divergence of the VRC01 epitope in an acquired HIV-1 isolate from the VRC01 epitope of reference HIV-1 strains that were most sensitive to VRC01-mediated neutralization. In HVTN 704/HPTN 085, the Env sequence-based predicted probability that VRC01 IC80 against the acquired isolate exceeded 1 µg/mL also significantly associated with PE. In HVTN 703/HPTN 081, a physicochemical-weighted Hamming distance across 50 VRC01 binding-associated Env AA positions of the acquired isolate from the most VRC01-sensitive HIV-1 strain significantly associated with PE. These results suggest that incorporating mutation scoring by BLOSUM62 and weighting by the strength of interactions at AA positions in the epitope:VRC01 interface can optimize performance of an Env sequence-based biomarker of VRC01 prevention efficacy. Future work could determine whether these results extend to other bnAbs and bnAb combinations.

Incorporating network diffusion and peak location information for better single-cell ATAC-seq data analysis.

Single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) data provided new insights into the understanding of epigenetic heterogeneity and transcriptional regulation. With the increasing abundance of dataset resources, there is an urgent need to extract more useful information through high-quality data analysis methods specifically designed for scATAC-seq. However, analyzing scATAC-seq data poses challenges due to its near binarization, high sparsity and ultra-high dimensionality properties. Here, we proposed a novel network diffusion-based computational method to comprehensively analyze scATAC-seq data, named Single-Cell ATAC-seq Analysis via Network Refinement with Peaks Location Information (SCARP). SCARP formulates the Network Refinement diffusion method under the graph theory framework to aggregate information from different network orders, effectively compensating for missing signals in the scATAC-seq data. By incorporating distance information between adjacent peaks on the genome, SCARP also contributes to depicting the co-accessibility of peaks. These two innovations empower SCARP to obtain lower-dimensional representations for both cells and peaks more effectively. We have demonstrated through sufficient experiments that SCARP facilitated superior analyses of scATAC-seq data. Specifically, SCARP exhibited outstanding cell clustering performance, enabling better elucidation of cell heterogeneity and the discovery of new biologically significant cell subpopulations. Additionally, SCARP was also instrumental in portraying co-accessibility relationships of accessible regions and providing new insight into transcriptional regulation. Consequently, SCARP identified genes that were involved in key Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to diseases and predicted reliable cis-regulatory interactions. To sum up, our studies suggested that SCARP is a promising tool to comprehensively analyze the scATAC-seq data.

Spatial proximity and gene function: a new dimension in prokaryotic gene association network analysis with 3D-GeneNet.

Understanding the biological functions and processes of genes, particularly those not yet characterized, is crucial for advancing molecular biology and identifying therapeutic targets. The hypothesis guiding this study is that the 3D proximity of genes correlates with their functional interactions and relevance in prokaryotes. We introduced 3D-GeneNet, an innovative software tool that utilizes high-throughput sequencing data from chromosome conformation capture techniques and integrates topological metrics to construct gene association networks. Through a series of comparative analyses focused on spatial versus linear distances, we explored various dimensions such as topological structure, functional enrichment levels, distribution patterns of linear distances among gene pairs, and the area under the receiver operating characteristic curve by utilizing model organism Escherichia coli K-12. Furthermore, 3D-GeneNet was shown to maintain good accuracy compared to multiple algorithms (neighbourhood, co-occurrence, coexpression, and fusion) across multiple bacteria, including E. coli, Brucella abortus, and Vibrio cholerae. In addition, the accuracy of 3D-GeneNet's prediction of long-distance gene interactions was identified by bacterial two-hybrid assays on E. coli K-12 MG1655, where 3D-GeneNet not only increased the accuracy of linear genomic distance tripled but also achieved 60% accuracy by running alone. Finally, it can be concluded that the applicability of 3D-GeneNet will extend to various bacterial forms, including Gram-negative, Gram-positive, single-, and multi-chromosomal bacteria through Hi-C sequencing and analysis. Such findings highlight the broad applicability and significant promise of this method in the realm of gene association network. 3D-GeneNet is freely accessible at https://github.com/gaoyuanccc/3D-GeneNet.

DDK-Linker: a network-based strategy identifies disease signals by linking high-throughput omics datasets to disease knowledge.

The high-throughput genomic and proteomic scanning approaches allow investigators to measure the quantification of genome-wide genes (or gene products) for certain disease conditions, which plays an essential role in promoting the discovery of disease mechanisms. The high-throughput approaches often generate a large gene list of interest (GOIs), such as differentially expressed genes/proteins. However, researchers have to perform manual triage and validation to explore the most promising, biologically plausible linkages between the known disease genes and GOIs (disease signals) for further study. Here, to address this challenge, we proposed a network-based strategy DDK-Linker to facilitate the exploration of disease signals hidden in omics data by linking GOIs to disease knowns genes. Specifically, it reconstructed gene distances in the protein-protein interaction (PPI) network through six network methods (random walk with restart, Deepwalk, Node2Vec, LINE, HOPE, Laplacian) to discover disease signals in omics data that have shorter distances to disease genes. Furthermore, benefiting from the establishment of knowledge base we established, the abundant bioinformatics annotations were provided for each candidate disease signal. To assist in omics data interpretation and facilitate the usage, we have developed this strategy into an application that users can access through a website or download the R package. We believe DDK-Linker will accelerate the exploring of disease genes and drug targets in a variety of omics data, such as genomics, transcriptomics and proteomics data, and provide clues for complex disease mechanism and pharmacological research. DDK-Linker is freely accessible at http://ddklinker.ncpsb.org.cn/.

MGCNSS: miRNA-disease association prediction with multi-layer graph convolution and distance-based negative sample selection strategy.

Identifying disease-associated microRNAs (miRNAs) could help understand the deep mechanism of diseases, which promotes the development of new medicine. Recently, network-based approaches have been widely proposed for inferring the potential associations between miRNAs and diseases. However, these approaches ignore the importance of different relations in meta-paths when learning the embeddings of miRNAs and diseases. Besides, they pay little attention to screening out reliable negative samples which is crucial for improving the prediction accuracy. In this study, we propose a novel approach named MGCNSS with the multi-layer graph convolution and high-quality negative sample selection strategy. Specifically, MGCNSS first constructs a comprehensive heterogeneous network by integrating miRNA and disease similarity networks coupled with their known association relationships. Then, we employ the multi-layer graph convolution to automatically capture the meta-path relations with different lengths in the heterogeneous network and learn the discriminative representations of miRNAs and diseases. After that, MGCNSS establishes a highly reliable negative sample set from the unlabeled sample set with the negative distance-based sample selection strategy. Finally, we train MGCNSS under an unsupervised learning manner and predict the potential associations between miRNAs and diseases. The experimental results fully demonstrate that MGCNSS outperforms all baseline methods on both balanced and imbalanced datasets. More importantly, we conduct case studies on colon neoplasms and esophageal neoplasms, further confirming the ability of MGCNSS to detect potential candidate miRNAs. The source code is publicly available on GitHub https://github.com/15136943622/MGCNSS/tree/master.

Scoring alignments by embedding vector similarity.

Sequence similarity is of paramount importance in biology, as similar sequences tend to have similar function and share common ancestry. Scoring matrices, such as PAM or BLOSUM, play a crucial role in all bioinformatics algorithms for identifying similarities, but have the drawback that they are fixed, independent of context. We propose a new scoring method for amino acid similarity that remedies this weakness, being contextually dependent. It relies on recent advances in deep learning architectures that employ self-supervised learning in order to leverage the power of enormous amounts of unlabelled data to generate contextual embeddings, which are vector representations for words. These ideas have been applied to protein sequences, producing embedding vectors for protein residues. We propose the E-score between two residues as the cosine similarity between their embedding vector representations. Thorough testing on a wide variety of reference multiple sequence alignments indicate that the alignments produced using the new $E$-score method, especially ProtT5-score, are significantly better than those obtained using BLOSUM matrices. The new method proposes to change the way alignments are computed, with far-reaching implications in all areas of textual data that use sequence similarity. The program to compute alignments based on various $E$-scores is available as a web server at e-score.csd.uwo.ca. The source code is freely available for download from github.com/lucian-ilie/E-score.

scDOT: enhancing single-cell RNA-Seq data annotation and uncovering novel cell types through multi-reference integration.

The proliferation of single-cell RNA-seq data has greatly enhanced our ability to comprehend the intricate nature of diverse tissues. However, accurately annotating cell types in such data, especially when handling multiple reference datasets and identifying novel cell types, remains a significant challenge. To address these issues, we introduce Single Cell annotation based on Distance metric learning and Optimal Transport (scDOT), an innovative cell-type annotation method adept at integrating multiple reference datasets and uncovering previously unseen cell types. scDOT introduces two key innovations. First, by incorporating distance metric learning and optimal transport, it presents a novel optimization framework. This framework effectively learns the predictive power of each reference dataset for new query data and simultaneously establishes a probabilistic mapping between cells in the query data and reference-defined cell types. Secondly, scDOT develops an interpretable scoring system based on the acquired probabilistic mapping, enabling the precise identification of previously unseen cell types within the data. To rigorously assess scDOT's capabilities, we systematically evaluate its performance using two diverse collections of benchmark datasets encompassing various tissues, sequencing technologies and diverse cell types. Our experimental results consistently affirm the superior performance of scDOT in cell-type annotation and the identification of previously unseen cell types. These advancements provide researchers with a potent tool for precise cell-type annotation, ultimately enriching our understanding of complex biological tissues.

Genetic factors predict hybrid formation in the British flora.

Natural hybridization can have a profound evolutionary impact, with consequences ranging from the extinction of rare taxa to the origin of new species. Natural hybridization is particularly common in plants; however, our understanding of the general factors that promote or prevent hybridization is hampered by the highly variable outcomes in different lineages. Here, we quantify the influence of different predictors on hybrid formation across species from an entire flora. We combine estimates of hybridization with ecological attributes and a new species-level phylogeny for over 1,100 UK flowering plant species. Our results show that genetic factors, particularly parental genetic distance, as well as phylogenetic position and ploidy, are key determinants of hybrid formation, whereas many other factors such as range overlap and genus size explain much less variation in hybrid formation. Overall, intrinsic genetic factors shape the evolutionary and ecological consequences of natural hybridization across species in a flora.

Estimating perceptions of the relative COVID risk of different social-distancing behaviors from respondents' pairwise assessments.

How do people compare the effectiveness of different social-distancing behaviors in avoiding the spread of viral infection? During the COVID pandemic, we showed 676 online respondents in the United States, United Kingdom, and Israel 30 pairs of brief videos of acquaintances meeting. We asked respondents to indicate which video from each pair depicted greater risk of COVID infection. Their choices imply that on average, respondents considered talking 14 min longer to be as risky as standing 1 foot closer, being indoors as standing 3 feet closer, being exposed to coughs or sneezes as 3 to 4 ft closer, greeting with a hug as 7 ft closer, and with a handshake as 5 ft closer. Respondents considered properly masking as protecting the wearer and interlocutor equally, removing the mask entirely or only when talking as standing 4 to 5 ft closer but wearing it under the nose as only 1 to 2 ft closer. We provide weaker evidence on beliefs about the interaction effects of different behaviors. In a more limited, ex post analysis, we find little evidence of differences in beliefs across subpopulations.

Rapid Adaptation and Interspecific Introgression in the North American Crop Pest Helicoverpa zea.

Insect crop pests threaten global food security. This threat is amplified through the spread of nonnative species and through adaptation of native pests to control measures. Adaptations such as pesticide resistance can result from selection on variation within a population, or through gene flow from another population. We investigate these processes in an economically important noctuid crop pest, Helicoverpa zea, which has evolved resistance to a wide range of pesticides. Its sister species Helicoverpa armigera, first detected as an invasive species in Brazil in 2013, introduced the pyrethroid-resistance gene CYP337B3 to South American H. zea via adaptive introgression. To understand whether this could contribute to pesticide resistance in North America, we sequenced 237 H. zea genomes across 10 sample sites. We report H. armigera introgression into the North American H. zea population. Two individuals sampled in Texas in 2019 carry H. armigera haplotypes in a 4 Mbp region containing CYP337B3. Next, we identify signatures of selection in the panmictic population of nonadmixed H. zea, identifying a selective sweep at a second cytochrome P450 gene: CYP333B3. We estimate that its derived allele conferred a ∼5% fitness advantage and show that this estimate explains independently observed rare nonsynonymous CYP333B3 mutations approaching fixation over a ∼20-year period. We also detect putative signatures of selection at a kinesin gene associated with Bt resistance. Overall, we document two mechanisms of rapid adaptation: the introduction of fitness-enhancing alleles through interspecific introgression, and selection on intraspecific variation.

Inter-domain distance prediction based on deep learning for domain assembly.

AlphaFold2 achieved a breakthrough in protein structure prediction through the end-to-end deep learning method, which can predict nearly all single-domain proteins at experimental resolution. However, the prediction accuracy of full-chain proteins is generally lower than that of single-domain proteins because of the incorrect interactions between domains. In this work, we develop an inter-domain distance prediction method, named DeepIDDP. In DeepIDDP, we design a neural network with attention mechanisms, where two new inter-domain features are used to enhance the ability to capture the interactions between domains. Furthermore, we propose a data enhancement strategy termed DPMSA, which is employed to deal with the absence of co-evolutionary information on targets. We integrate DeepIDDP into our previously developed domain assembly method SADA, termed SADA-DeepIDDP. Tested on a given multi-domain benchmark dataset, the accuracy of SADA-DeepIDDP inter-domain distance prediction is 11.3% and 21.6% higher than trRosettaX and trRosetta, respectively. The accuracy of the domain assembly model is 2.5% higher than that of SADA. Meanwhile, we reassemble 68 human multi-domain protein models with TM-score ≤ 0.80 from the AlphaFold protein structure database, where the average TM-score is improved by 11.8% after the reassembly by our method. The online server is at http://zhanglab-bioinf.com/DeepIDDP/.

RNet: a network strategy to predict RNA binding preferences.

Determining the RNA binding preferences remains challenging because of the bottleneck of the binding interactions accompanied by subtle RNA flexibility. Typically, designing RNA inhibitors involves screening thousands of potential candidates for binding. Accurate binding site information can increase the number of successful hits even with few candidates. There are two main issues regarding RNA binding preference: binding site prediction and binding dynamical behavior prediction. Here, we propose one interpretable network-based approach, RNet, to acquire precise binding site and binding dynamical behavior information. RNetsite employs a machine learning-based network decomposition algorithm to predict RNA binding sites by analyzing the local and global network properties. Our research focuses on large RNAs with 3D structures without considering smaller regulatory RNAs, which are too small and dynamic. Our study shows that RNetsite outperforms existing methods, achieving precision values as high as 0.701 on TE18 and 0.788 on RB9 tests. In addition, RNetsite demonstrates remarkable robustness regarding perturbations in RNA structures. We also developed RNetdyn, a distance-based dynamical graph algorithm, to characterize the interface dynamical behavior consequences upon inhibitor binding. The simulation testing of competitive inhibitors indicates that RNetdyn outperforms the traditional method by 30%. The benchmark testing results demonstrate that RNet is highly accurate and robust. Our interpretable network algorithms can assist in predicting RNA binding preferences and accelerating RNA inhibitor design, providing valuable insights to the RNA research community.

Detection of PatIent-Level distances from single cell genomics and pathomics data with Optimal Transport (PILOT).

Although clinical applications represent the next challenge in single-cell genomics and digital pathology, we still lack computational methods to analyze single-cell or pathomics data to find sample-level trajectories or clusters associated with diseases. This remains challenging as single-cell/pathomics data are multi-scale, i.e., a sample is represented by clusters of cells/structures, and samples cannot be easily compared with each other. Here we propose PatIent Level analysis with Optimal Transport (PILOT). PILOT uses optimal transport to compute the Wasserstein distance between two individual single-cell samples. This allows us to perform unsupervised analysis at the sample level and uncover trajectories or cellular clusters associated with disease progression. We evaluate PILOT and competing approaches in single-cell genomics or pathomics studies involving various human diseases with up to 600 samples/patients and millions of cells or tissue structures. Our results demonstrate that PILOT detects disease-associated samples from large and complex single-cell or pathomics data. Moreover, PILOT provides a statistical approach to find changes in cell populations, gene expression, and tissue structures related to the trajectories or clusters supporting interpretation of predictions.

Genetic distance and ancestry proportion modify the association between maternal genetic risk score of type 2 diabetes and fetal growth.

BACKGROUND: Maternal genetic risk of type 2 diabetes (T2D) has been associated with fetal growth, but the influence of genetic ancestry is not yet fully understood. We aimed to investigate the influence of genetic distance (GD) and genetic ancestry proportion (GAP) on the association of maternal genetic risk score of T2D (GRST2D) with fetal weight and birthweight. METHODS: Multi-ancestral pregnant women (n = 1,837) from the NICHD Fetal Growth Studies - Singletons cohort were included in the current analyses. Fetal weight (in grams, g) was estimated from ultrasound measurements of fetal biometry, and birthweight (g) was measured at delivery. GRST2D was calculated using T2D-associated variants identified in the latest trans-ancestral genome-wide association study and was categorized into quartiles. GD and GAP were estimated using genotype data of four reference populations. GD was categorized into closest, middle, and farthest tertiles, and GAP was categorized as highest, medium, and lowest. Linear regression analyses were performed to test the association of GRST2D with fetal weight and birthweight, adjusted for covariates, in each GD and GAP category. RESULTS: Among women with the closest GD from African and Amerindigenous ancestries, the fourth and third GRST2D quartile was significantly associated with 5.18 to 7.48 g (weeks 17-20) and 6.83 to 25.44 g (weeks 19-27) larger fetal weight compared to the first quartile, respectively. Among women with middle GD from European ancestry, the fourth GRST2D quartile was significantly associated with 5.73 to 21.21 g (weeks 18-26) larger fetal weight. Furthermore, among women with middle GD from European and African ancestries, the fourth and second GRST2D quartiles were significantly associated with 117.04 g (95% CI = 23.88-210.20, p = 0.014) and 95.05 g (95% CI = 4.73-185.36, p = 0.039) larger birthweight compared to the first quartile, respectively. The absence of significant association among women with the closest GD from East Asian ancestry was complemented by a positive significant association among women with the highest East Asian GAP. CONCLUSIONS: The association between maternal GRST2D and fetal growth began in early-second trimester and was influenced by GD and GAP. The results suggest the use of genetic GD and GAP could improve the generalizability of GRS.

Enhanced group-level dorsolateral prefrontal cortex subregion parcellation through functional connectivity-based distance-constrained spectral clustering with application to autism spectrum disorder.

The dorsolateral prefrontal cortex (DLPFC) assumes a central role in cognitive and behavioral control, emerging as a crucial target region for interventions in autism spectrum disorder neuroregulation. Consequently, we endeavor to unravel the functional subregions within the DLPFC to shed light on the intricate functions of the brain. We introduce a distance-constrained spectral clustering (SC-DW) methodology that leverages functional connection to identify distinctive functional subregions within the DLPFC. Furthermore, we verify the relationship between the functional characteristics of these subregions and their clinical implications. Our methodology begins with principal component analysis to extract the salient features. Subsequently, we construct an adjacency matrix, which is constrained by the spatial properties of the brain, by linearly combining the distance matrix and a similarity matrix. The quality of spectral clustering is further optimized through multiple cluster evaluation coefficient. The results from SC-DW revealed four uniform and contiguous subregions within the bilateral DLPFC. Notably, we observe a substantial positive correlation between the functional characteristics of the third and fourth subregions in the left DLPFC with clinical manifestations. These findings underscore the unique insights offered by our proposed methodology in the realms of brain subregion delineation and therapeutic targeting.

Locust density shapes energy metabolism and oxidative stress resulting in divergence of flight traits.

Flight ability is essential for the enormous diversity and evolutionary success of insects. The migratory locusts exhibit flight capacity plasticity in gregarious and solitary individuals closely linked with different density experiences. However, the differential mechanisms underlying flight traits of locusts are largely unexplored. Here, we investigated the variation of flight capacity by using behavioral, physiological, and multiomics approaches. Behavioral assays showed that solitary locusts possess high initial flight speeds and short-term flight, whereas gregarious locusts can fly for a longer distance at a relatively lower speed. Metabolome-transcriptome analysis revealed that solitary locusts have more active flight muscle energy metabolism than gregarious locusts, whereas gregarious locusts show less evidence of reactive oxygen species production during flight. The repression of metabolic activity by RNA interference markedly reduced the initial flight speed of solitary locusts. Elevating the oxidative stress by paraquat injection remarkably inhibited the long-distance flight of gregarious locusts. In respective crowding and isolation treatments, energy metabolic profiles and flight traits of solitary and gregarious locusts were reversed, indicating that the differentiation of flight capacity depended on density and can be reshaped rapidly. The density-dependent flight traits of locusts were attributed to the plasticity of energy metabolism and degree of oxidative stress production but not energy storage. The findings provided insights into the mechanism underlying the trade-off between velocity and sustainability in animal locomotion and movement.

Linguistic measures of psychological distance track symptom levels and treatment outcomes in a large set of psychotherapy transcripts.

SignificanceUsing language to "distance" ourselves from distressing situations (i.e., by talking less about ourselves and the present moment) can help us manage emotions. Here, we translate this basic research to discover that such "linguistic distancing" is a replicable measure of mental health in a large set of therapy transcripts (N = 6,229). Additionally, clustering techniques showed that language alone could identify participants who differed on both symptom severity and treatment outcomes. These findings lay the foundation for 1) tools that can rapidly identify people in need of psychological services based on language alone and 2) linguistic interventions that can improve mental health.

Spatial turnover of soil viral populations and genotypes overlain by cohesive responses to moisture in grasslands.

Viruses shape microbial communities, food web dynamics, and carbon and nutrient cycling in diverse ecosystems. However, little is known about the patterns and drivers of viral community composition, particularly in soil, precluding a predictive understanding of viral impacts on terrestrial habitats. To investigate soil viral community assembly processes, here we analyzed 43 soil viromes from a rainfall manipulation experiment in a Mediterranean grassland in California. We identified 5,315 viral populations (viral operational taxonomic units [vOTUs] with a representative sequence ≥10 kbp) and found that viral community composition exhibited a highly significant distance-decay relationship within the 200-m2 field site. This pattern was recapitulated by the intrapopulation microheterogeneity trends of prevalent vOTUs (detected in ≥90% of the viromes), which tended to exhibit negative correlations between spatial distance and the genomic similarity of their predominant allelic variants. Although significant spatial structuring was also observed in the bacterial and archaeal communities, the signal was dampened relative to the viromes, suggesting differences in local assembly drivers for viruses and prokaryotes and/or differences in the temporal scales captured by viromes and total DNA. Despite the overwhelming spatial signal, evidence for environmental filtering was revealed in a protein-sharing network analysis, wherein a group of related vOTUs predicted to infect actinobacteria was shown to be significantly enriched in low-moisture samples distributed throughout the field. Overall, our results indicate a highly diverse, dynamic, active, and spatially structured soil virosphere capable of rapid responses to changing environmental conditions.