RESUMO
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.
Assuntos
Adenocarcinoma , Neoplasias Colorretais , Humanos , Adenocarcinoma/patologia , Neoplasias Colorretais/genética , Neoplasias Colorretais/imunologia , Neoplasias Colorretais/patologia , Processamento de Imagem Assistida por Computador , Oncogenes , Microambiente TumoralRESUMO
Microbes drive most ecosystems and are modulated by viruses that impact their lifespan, gene flow, and metabolic outputs. However, ecosystem-level impacts of viral community diversity remain difficult to assess due to classification issues and few reference genomes. Here, we establish an â¼12-fold expanded global ocean DNA virome dataset of 195,728 viral populations, now including the Arctic Ocean, and validate that these populations form discrete genotypic clusters. Meta-community analyses revealed five ecological zones throughout the global ocean, including two distinct Arctic regions. Across the zones, local and global patterns and drivers in viral community diversity were established for both macrodiversity (inter-population diversity) and microdiversity (intra-population genetic variation). These patterns sometimes, but not always, paralleled those from macro-organisms and revealed temperate and tropical surface waters and the Arctic as biodiversity hotspots and mechanistic hypotheses to explain them. Such further understanding of ocean viruses is critical for broader inclusion in ecosystem models.
Assuntos
Organismos Aquáticos/genética , Biodiversidade , Vírus de DNA/genética , DNA Viral/genética , Metagenoma , Microbiologia da ÁguaRESUMO
The neocortex is a complex neurobiological system with many interacting regions. How these regions work together to subserve flexible behavior and cognition has become increasingly amenable to rigorous research. Here, I review recent experimental and theoretical work on the modus operandi of a multiregional cortex. These studies revealed several general principles for the neocortical interareal connectivity, low-dimensional macroscopic gradients of biological properties across cortical areas, and a hierarchy of timescales for information processing. Theoretical work suggests testable predictions regarding differential excitation and inhibition along feedforward and feedback pathways in the cortical hierarchy. Furthermore, modeling of distributed working memory and simple decision-making has given rise to a novel mathematical concept, dubbed bifurcation in space, that potentially explains how different cortical areas, with a canonical circuit organization but gradients of biological heterogeneities, are able to subserve their respective (e.g., sensory coding versus executive control) functions in a modularly organized brain.
Assuntos
Neocórtex , Cognição/fisiologia , Função Executiva , Memória de Curto Prazo/fisiologia , Neocórtex/fisiologia , Rede Nervosa/fisiologiaRESUMO
Most single animal cells have an internal vector that determines where recycling membrane is added to the cell's surface. Because of the specific molecular composition of this added membrane, a dynamic asymmetry is formed on the surface of the cell. The consequences of this dynamic asymmetry are discussed, together with what they imply for how cells move. The polarity of a single-celled embryo, such as that of the nematode Caenorhabditis elegans, is explored in a similar framework.
Assuntos
Membrana Celular/química , Endocitose , Animais , Biologia/métodos , Caenorhabditis elegans , Movimento Celular , Citoesqueleto/metabolismo , Dictyostelium , Fibroblastos/metabolismo , Células HeLa , Hemaglutininas/química , Humanos , Proteínas de Membrana/química , Proteínas/químicaRESUMO
Orientational order, encoded in anisotropic fields, plays an important role during the development of an organism. A striking example of this is the freshwater polyp Hydra, where topological defects in the muscle fiber orientation have been shown to localize to key features of the body plan. This body plan is organized by morphogen concentration gradients, raising the question how muscle fiber orientation, morphogen gradients and body shape interact. Here, we introduce a minimal model that couples nematic orientational order to the gradient of a morphogen field. We show that on a planar surface, alignment to a radial concentration gradient can induce unbinding of topological defects, as observed during budding and tentacle formation in Hydra, and stabilize aster/vortex-like defects, as observed at a Hydra's mouth. On curved surfaces mimicking the morphologies of Hydra in various stages of development-from spheroid to adult-our model reproduces the experimentally observed reorganization of orientational order. Our results suggest how gradient alignment and curvature effects may work together to control orientational order during development and lay the foundations for future modeling efforts that will include the tissue mechanics that drive shape deformations.
Assuntos
Hydra , Animais , Anisotropia , Morfogênese , Hydra/fisiologia , Regeneração/fisiologia , Padronização CorporalRESUMO
The COVID-19 pandemic struck societies directly and indirectly, not just challenging population health but disrupting many aspects of life. Different effects of the spreading virus-and the measures to fight it-are reported and discussed in different scientific fora, with hard-to-compare methods and metrics from different traditions. While the pandemic struck some groups more than others, it is difficult to assess the comprehensive impact on social inequalities. This paper gauges social inequalities using individual-level administrative data for Sweden's entire population. We describe and analyze the relative risks for different social groups in four dimensions-gender, education, income, and world region of birth-to experience three types of COVID-19 incidence, as well as six additional negative life outcomes that reflect general health, access to medical care, and economic strain. During the pandemic, the overall population faced severe morbidity and mortality from COVID-19 and saw higher all-cause mortality, income losses and unemployment risks, as well as reduced access to medical care. These burdens fell more heavily on individuals with low income or education and on immigrants. Although these vulnerable groups experienced larger absolute risks of suffering the direct and indirect consequences of the pandemic, the relative risks in pandemic years (2020 and 2021) were conspicuously similar to those in prepandemic years (2016 to 2019).
Assuntos
COVID-19 , Humanos , COVID-19/epidemiologia , Pandemias , Suécia/epidemiologia , Risco , Classe SocialRESUMO
Recent work has recognized a gradient-like organization in cortical function, spanning from primary sensory to transmodal cortices. It has been suggested that this axis is aligned with regional differences in neurotransmitter expression. Given the abundance of dopamine D1-receptors (D1DR), and its importance for modulation and neural gain, we tested the hypothesis that D1DR organization is aligned with functional architecture, and that inter-regional relationships in D1DR co-expression modulate functional cross talk. Using the world's largest dopamine D1DR-PET and MRI database (N = 180%, 50% female), we demonstrate that D1DR organization follows a unimodal-transmodal hierarchy, expressing a high spatial correspondence to the principal gradient of functional connectivity. We also demonstrate that individual differences in D1DR density between unimodal and transmodal regions are associated with functional differentiation of the apices in the cortical hierarchy. Finally, we show that spatial co-expression of D1DR primarily modulates couplings within, but not between, functional networks. Together, our results show that D1DR co-expression provides a biomolecular layer to the functional organization of the brain.
Assuntos
Encéfalo , Dopamina , Feminino , Humanos , Masculino , Imageamento por Ressonância Magnética/métodosRESUMO
Axonal projections from layer V neurons of distinct neocortical areas are topographically organized into discrete clusters within the pontine nuclei during the establishment of voluntary movements. However, the molecular determinants controlling corticopontine connectivity are insufficiently understood. Here, we show that an intrinsic cortical genetic program driven by Nr2f1 graded expression is directly implicated in the organization of corticopontine topographic mapping. Transgenic mice lacking cortical expression of Nr2f1 and exhibiting areal organization defects were used as model systems to investigate the arrangement of corticopontine projections. By combining three-dimensional digital brain atlas tools, Cre-dependent mouse lines and axonal tracing, we show that Nr2f1 expression in postmitotic neurons spatially and temporally controls somatosensory topographic projections, whereas expression in progenitor cells influences the ratio between corticopontine and corticospinal fibres passing the pontine nuclei. We conclude that cortical gradients of area-patterning genes are directly implicated in the establishment of a topographic somatotopic mapping from the cortex onto pontine nuclei.
Assuntos
Mapeamento Encefálico , Ponte , Animais , Axônios , Córtex Cerebral , Camundongos , Vias Neurais/fisiologia , Neurônios , Ponte/fisiologiaRESUMO
The importance of regulatory features in health and disease is increasing, making it crucial to identify the hallmarks of these features. Self-attention networks (SAN) have given rise to numerous models for the prediction of complex phenomena. But the potential of SANs in biological models was limited because of high memory requirement proportional to input token length and lack of interpretability of self-attention scores. To overcome these constraints, we propose a deep learning model named Interpretable Self-Attention Network for REGulatory interactions (ISANREG) that combines both block self-attention and attention-attribution mechanisms. This model predicts transcription factor-bound motif instances and DNA-mediated TF-TF interactions using self-attention attribution scores derived from the network, overcoming the limitations of previous deep learning models. ISANREG will serve as a framework for other biological models in interpreting the contribution of the input with single-nucleotide resolution.
Assuntos
Modelos Biológicos , Nucleotídeos , Fatores de TranscriçãoRESUMO
Network neuroscience offers a unique framework to understand the organizational principles of the human brain. Despite recent progress, our understanding of how the brain is modulated by focal lesions remains incomplete. Resection of the temporal lobe is the most effective treatment to control seizures in pharmaco-resistant temporal lobe epilepsy (TLE), making this syndrome a powerful model to study lesional effects on network organization in young and middle-aged adults. Here, we assessed the downstream consequences of a focal lesion and its surgical resection on the brain's structural connectome, and explored how this reorganization relates to clinical variables at the individual patient level. We included adults with pharmaco-resistant TLE (n = 37) who underwent anterior temporal lobectomy between two imaging time points, as well as age- and sex-matched healthy controls who underwent comparable imaging (n = 31). Core to our analysis was the projection of high-dimensional structural connectome data-derived from diffusion MRI tractography from each subject-into lower-dimensional gradients. We then compared connectome gradients in patients relative to controls before surgery, tracked surgically-induced connectome reconfiguration from pre- to postoperative time points, and examined associations to patient-specific clinical and imaging phenotypes. Before surgery, individuals with TLE presented with marked connectome changes in bilateral temporo-parietal regions, reflecting an increased segregation of the ipsilateral anterior temporal lobe from the rest of the brain. Surgery-induced connectome reorganization was localized to this temporo-parietal subnetwork, but primarily involved postoperative integration of contralateral regions with the rest of the brain. Using a partial least-squares analysis, we uncovered a latent clinical imaging signature underlying this pre- to postoperative connectome reorganization, showing that patients who displayed postoperative integration in bilateral fronto-occipital cortices also had greater preoperative ipsilateral hippocampal atrophy, lower seizure frequency and secondarily generalized seizures. Our results bridge the effects of focal brain lesions and their surgical resections with large-scale network reorganization and interindividual clinical variability, thus offering new avenues to examine the fundamental malleability of the human brain.
Assuntos
Lobectomia Temporal Anterior , Conectoma , Epilepsia do Lobo Temporal , Lobo Temporal , Humanos , Feminino , Masculino , Adulto , Epilepsia do Lobo Temporal/cirurgia , Epilepsia do Lobo Temporal/fisiopatologia , Epilepsia do Lobo Temporal/diagnóstico por imagem , Epilepsia do Lobo Temporal/patologia , Lobo Temporal/patologia , Lobo Temporal/cirurgia , Lobo Temporal/diagnóstico por imagem , Lobectomia Temporal Anterior/métodos , Pessoa de Meia-Idade , Adulto Jovem , Imagem de Tensor de Difusão , Rede Nervosa/diagnóstico por imagem , Rede Nervosa/patologia , Epilepsia Resistente a Medicamentos/cirurgia , Epilepsia Resistente a Medicamentos/diagnóstico por imagem , Epilepsia Resistente a Medicamentos/fisiopatologia , Epilepsia Resistente a Medicamentos/patologiaRESUMO
Declarative memory encompasses episodic and semantic divisions. Episodic memory captures singular events with specific spatiotemporal relationships, while semantic memory houses context-independent knowledge. Behavioural and functional neuroimaging studies have revealed common and distinct neural substrates of both memory systems, implicating mesiotemporal lobe (MTL) regions such as the hippocampus and distributed neocortices. Here, we explored declarative memory system reorganization in patients with unilateral temporal lobe epilepsy (TLE) as a human disease model to test the impact of variable degrees of MTL pathology on memory function. Our cohort included 31 patients with TLE as well as 60 age and sex-matched healthy controls, and all participants underwent episodic and semantic retrieval tasks during a multimodal MRI session. The functional MRI tasks were closely matched in terms of stimuli and trial design. Capitalizing on non-linear connectome gradient mapping techniques, we derived task-based functional topographies during episodic and semantic memory states, both in the MTL and in neocortical networks. Comparing neocortical and hippocampal functional gradients between TLE patients and healthy controls, we observed a marked topographic reorganization of both neocortical and MTL systems during episodic memory states. Neocortical alterations were characterized by reduced functional differentiation in TLE across lateral temporal and midline parietal cortices in both hemispheres. In the MTL, on the other hand, patients presented with a more marked functional differentiation of posterior and anterior hippocampal segments ipsilateral to the seizure focus and pathological core, indicating perturbed intrahippocampal connectivity. Semantic memory reorganization was also found in bilateral lateral temporal and ipsilateral angular regions, while hippocampal functional topographies were unaffected. Leveraging MRI proxies of MTL pathology, we furthermore observed alterations in hippocampal microstructure and morphology that were associated with TLE-related functional reorganization during episodic memory. Moreover, correlation analysis and statistical mediation models revealed that these functional alterations contributed to behavioural deficits in episodic, but again not semantic memory in patients. Altogether, our findings suggest that semantic processes rely on distributed neocortical networks, while episodic processes are supported by a network involving both the hippocampus and neocortex. Alterations of such networks can provide a compact signature of state-dependent reorganization in conditions associated with MTL damage, such as TLE.
RESUMO
Visual hallucinations in Parkinson's disease can be viewed from a systems-level perspective, whereby dysfunctional communication between brain networks responsible for perception predisposes a person to hallucinate. To this end, abnormal functional interactions between higher-order and primary sensory networks have been implicated in the pathophysiology of visual hallucinations in Parkinson's disease, however the precise signatures remain to be determined. Dimensionality reduction techniques offer a novel means for simplifying the interpretation of multidimensional brain imaging data, identifying hierarchical patterns in the data that are driven by both within- and between-functional network changes. Here, we applied two complementary non-linear dimensionality reduction techniques-diffusion-map embedding and t-distributed stochastic neighbour embedding (t-SNE)-to resting state functional MRI data, in order to characterize the altered functional hierarchy associated with susceptibility to visual hallucinations. Our study involved 77 people with Parkinson's disease (31 with hallucinations; 46 without hallucinations) and 19 age-matched healthy control subjects. In patients with visual hallucinations, we found compression of the unimodal-heteromodal gradient consistent with increased functional integration between sensory and higher order networks. This was mirrored in a traditional functional connectivity analysis, which showed increased connectivity between the visual and default mode networks in the hallucinating group. Together, these results suggest a route by which higher-order regions may have excessive influence over earlier sensory processes, as proposed by theoretical models of hallucinations across disorders. By contrast, the t-SNE analysis identified distinct alterations in prefrontal regions, suggesting an additional layer of complexity in the functional brain network abnormalities implicated in hallucinations, which was not apparent in traditional functional connectivity analyses. Together, the results confirm abnormal brain organization associated with the hallucinating phenotype in Parkinson's disease and highlight the utility of applying convergent dimensionality reduction techniques to investigate complex clinical symptoms. In addition, the patterns we describe in Parkinson's disease converge with those seen in other conditions, suggesting that reduced hierarchical differentiation across sensory-perceptual systems may be a common transdiagnostic vulnerability in neuropsychiatric disorders with perceptual disturbances.
Assuntos
Doença de Parkinson , Humanos , Doença de Parkinson/complicações , Imageamento por Ressonância Magnética/métodos , Alucinações/etiologia , Encéfalo/diagnóstico por imagem , Mapeamento EncefálicoRESUMO
This study aimed to explore the topographic features of thalamic subregions, functional connectomes and hierarchical organizations between thalamus and cortex in poststroke fatigue patients. We consecutively recruited 121 acute ischemic stroke patients (mean age: 59 years) and 46 healthy controls matched for age, sex, and educational level. The mean age was 59 years (range 19-80) and 38% of acute stroke patients were females. Resting-state functional and structural magnetic resonance imaging were conducted on all participants. The fatigue symptoms were measured using the Fatigue Severity Scale. The thalamic functional subdivisions corresponding to the canonical functional network were defined using the winner-take-all parcellation method. Thalamic functional gradients were derived using the diffusion embedding analysis. The results suggested abnormal functional connectivity of thalamic subregions primarily located in the temporal lobe, posterior cingulate gyrus, parietal lobe, and precuneus. The thalamus showed a gradual increase from the medial to the lateral in all groups, but the right thalamus shifted more laterally in poststroke fatigue patients than in non- poststroke fatigue patients. Poststroke fatigue patients also had higher gradient scores in the somatomotor network and the right medial prefrontal and premotor thalamic regions, but lower values in the right lateral prefrontal thalamus. The findings suggested that poststroke fatigue patients had altered functional connectivity and thalamocortical hierarchical organizations, providing new insights into the neural mechanisms of the thalamus.
Assuntos
Conectoma , AVC Isquêmico , Acidente Vascular Cerebral , Feminino , Humanos , Adulto Jovem , Adulto , Pessoa de Meia-Idade , Idoso , Idoso de 80 Anos ou mais , Masculino , Conectoma/métodos , AVC Isquêmico/patologia , Tálamo/patologia , Imageamento por Ressonância Magnética/métodos , Acidente Vascular Cerebral/complicações , Acidente Vascular Cerebral/diagnóstico por imagem , Acidente Vascular Cerebral/patologia , Fadiga/diagnóstico por imagem , Fadiga/etiologiaRESUMO
SignificanceThe temperature difference between low and high latitudes is one measure of the efficiency of the global climate system in redistributing heat and is used to test the ability of models to represent the climate system through time. Here, we show that the latitudinal temperature gradient has exhibited a consistent inverse relationship with global mean sea-surface temperature for at least the past 95 million years. Our results help reduce conflicts between climate models and empirical estimates of temperature and argue for a fundamental consistency in the dynamics of heat transport and radiative transfer across vastly different background states.
RESUMO
Adolescence is a time of profound changes in the physical wiring and function of the brain. Here, we analyzed structural and functional brain network development in an accelerated longitudinal cohort spanning 14 to 25 y (n = 199). Core to our work was an advanced in vivo model of cortical wiring incorporating MRI features of corticocortical proximity, microstructural similarity, and white matter tractography. Longitudinal analyses assessing age-related changes in cortical wiring identified a continued differentiation of multiple corticocortical structural networks in youth. We then assessed structure-function coupling using resting-state functional MRI measures in the same participants both via cross-sectional analysis at baseline and by studying longitudinal change between baseline and follow-up scans. At baseline, regions with more similar structural wiring were more likely to be functionally coupled. Moreover, correlating longitudinal structural wiring changes with longitudinal functional connectivity reconfigurations, we found that increased structural differentiation, particularly between sensory/unimodal and default mode networks, was reflected by reduced functional interactions. These findings provide insights into adolescent development of human brain structure and function, illustrating how structural wiring interacts with the maturation of macroscale functional hierarchies.
Assuntos
Desenvolvimento do Adolescente , Encéfalo , Conectoma , Adolescente , Encéfalo/fisiologia , Encéfalo/ultraestrutura , Estudos Transversais , Humanos , Imageamento por Ressonância Magnética , Rede Nervosa/fisiologia , Rede Nervosa/ultraestruturaRESUMO
To rapidly process temporal information at a low metabolic cost, biological neurons integrate inputs as an analog sum, but communicate with spikes, binary events in time. Analog neuromorphic hardware uses the same principles to emulate spiking neural networks with exceptional energy efficiency. However, instantiating high-performing spiking networks on such hardware remains a significant challenge due to device mismatch and the lack of efficient training algorithms. Surrogate gradient learning has emerged as a promising training strategy for spiking networks, but its applicability for analog neuromorphic systems has not been demonstrated. Here, we demonstrate surrogate gradient learning on the BrainScaleS-2 analog neuromorphic system using an in-the-loop approach. We show that learning self-corrects for device mismatch, resulting in competitive spiking network performance on both vision and speech benchmarks. Our networks display sparse spiking activity with, on average, less than one spike per hidden neuron and input, perform inference at rates of up to 85,000 frames per second, and consume less than 200 mW. In summary, our work sets several benchmarks for low-energy spiking network processing on analog neuromorphic hardware and paves the way for future on-chip learning algorithms.
Assuntos
Redes Neurais de Computação , Potenciais de Ação/fisiologia , Algoritmos , Encéfalo/fisiologia , Computadores , Modelos Biológicos , Modelos Neurológicos , Modelos Teóricos , Neurônios/fisiologiaRESUMO
The gut microbiota features important genetic diversity, and the specific spatial features of the gut may shape evolution within this environment. We investigate the fixation probability of neutral bacterial mutants within a minimal model of the gut that includes hydrodynamic flow and resulting gradients of food and bacterial concentrations. We find that this fixation probability is substantially increased, compared with an equivalent well-mixed system, in the regime where the profiles of food and bacterial concentration are strongly spatially dependent. Fixation probability then becomes independent of total population size. We show that our results can be rationalized by introducing an active population, which consists of those bacteria that are actively consuming food and dividing. The active population size yields an effective population size for neutral mutant fixation probability in the gut.
Assuntos
Bactérias , Biodiversidade , Microbioma Gastrointestinal , Hidrodinâmica , Bactérias/genética , Evolução Biológica , Alimentos , Microbiologia de Alimentos , Humanos , Densidade Demográfica , RNA Ribossômico 16S/genéticaRESUMO
Host specialization plays a critical role in the ecology and evolution of plant-microbe symbiosis. Theory predicts that host specialization is associated with microbial genome streamlining and is influenced by the abundance of host species, both of which can vary across latitudes, leading to a latitudinal gradient in host specificity. Here, we quantified the host specificity and composition of plant-bacteria symbioses on leaves across 329 tree species spanning a latitudinal gradient. Our analysis revealed a predominance of host-specialized leaf bacteria. The degree of host specificity was negatively correlated with bacterial genome size and the local abundance of host plants. Additionally, we found an increased host specificity at lower latitudes, aligning with the high prevalence of small bacterial genomes and rare host species in the tropics. These findings underscore the importance of genome streamlining and host abundance in the evolution of host specificity in plant-associated bacteria along the latitudinal gradient.
Assuntos
Tamanho do Genoma , Especificidade de Hospedeiro , Folhas de Planta , Simbiose , Folhas de Planta/microbiologia , Bactérias/genética , Bactérias/classificação , Genoma Bacteriano , Árvores/microbiologiaRESUMO
AbstractA key question about macroevolutionary speciation rates is whether they are controlled by microevolutionary processes operating at the population level. For example, does spatial variation in population genetic differentiation underlie geographical gradients in speciation rates? Previous work suggests that speciation rates increase with elevation in Neotropical birds, but underlying population-level gradients remain unexplored. Here, we characterize elevational phylogeographic diversity between montane and lowland birds in the megadiverse Andes-Amazonian system and assess its relationship to speciation rates to evaluate the link between population-level differentiation and species-level diversification. We aggregated and georeferenced nearly 7,000 mitochondrial DNA sequences across 103 species or species complexes in the Andes and Amazonia and used these sequences to describe phylogeographic differentiation across both regions. Our results show increased levels of both discrete and continuous metrics of population structure in the Andean mountains compared with the Amazonian lowlands. However, higher levels of population differentiation do not predict higher rates of speciation in our dataset. Multiple potential factors may lead to our observed decoupling of initial population divergence and speciation rates, including the ephemerality of incipient species and the multifaceted nature of the speciation process, as well as methodological challenges associated with estimating rates of population differentiation and speciation.
Assuntos
Aves , DNA Mitocondrial , Animais , Filogenia , Filogeografia , Aves/genética , DNA Mitocondrial/genética , Deriva Genética , Especiação GenéticaRESUMO
AbstractUnderstanding and predicting the evolutionary responses of complex morphological traits to selection remains a major challenge in evolutionary biology. Because traits are genetically correlated, selection on a particular trait produces both direct effects on the distribution of that trait and indirect effects on other traits in the population. The correlations between traits can strongly impact evolutionary responses to selection and may thus impose constraints on adaptation. Here, we used museum specimens and comparative quantitative genetic approaches to investigate whether the covariation among cranial traits facilitated or constrained the response to selection during the major dietary transitions in one of the world's most ecologically diverse mammalian families-the phyllostomid bats. We reconstructed the set of net selection gradients that would have acted on each cranial trait during the major transitions to feeding specializations and decomposed the selection responses into their direct and indirect components. We found that for all transitions, most traits capturing craniofacial length evolved toward adaptive directions owing to direct selection. Additionally, we showed instances of dietary transitions in which the complex interaction between the patterns of covariation among traits and the strength and direction of selection either constrained or facilitated evolution. Our work highlights the importance of considering the within-species covariation estimates to quantify evolvability and to disentangle the relative contribution of variational constraints versus selective causes for observed patterns.