Global topology of human connectome is insensitive to early life environments - A prospective longitudinal study of the general population.

The widely acknowledged detrimental impact of early adversity on child development has driven efforts to understand the underlying mechanisms that may mediate these effects within the developing brain. Recent efforts have begun to move beyond associating adversity with the morphology of individual brain regions towards determining if and how adversity might shape their interconnectivity. However, whether adversity effects a global shift in the organisation of whole-brain networks remains unclear. In this study, we assessed this possibility using parental questionnaire and diffusion imaging data from The Avon Longitudinal Study of Parents and Children (ALSPAC, N = 913), a prospective longitudinal study spanning more than 20 years. We tested whether a wide range of adversities-including experiences of abuse, domestic violence, physical and emotional cruelty, poverty, neglect, and parental separation-measured by questionnaire within the first seven years of life were significantly associated with the tractography-derived connectome in young adulthood. We tested this across multiple measures of organisation and using a computational model that simulated the wiring economy of the brain. We found no significant relationships between early exposure to any form of adversity and the global organisation of the structural connectome in young adulthood. We did detect local differences in the medial prefrontal cortex, as well as an association between weaker brain wiring constraints and greater externalising behaviour in adolescence. Our results indicate that further efforts are necessary to delimit the magnitude and functional implications of adversity-related differences in connectomic organization. RESEARCH HIGHLIGHTS: Diverse prospective measures of the early-life environment do not predict the organisation of the DTI tractography-derived connectome in young adulthood Wiring economy of the connectome is weakly associated with externalising in adolescence, but not internalising or cognitive ability Further work is needed to establish the scope and significance of global adversity-related differences in the structural connectome.

Mechanisms underlying divergent relationships between Ca2+ and YAP/TAZ signalling.

Yes-associated protein (YAP) and its homologue TAZ are transducers of several biochemical and biomechanical signals, integrating multiplexed inputs from the microenvironment into higher level cellular functions such as proliferation, differentiation and migration. Emerging evidence suggests that Ca2+ is a key second messenger that connects microenvironmental input signals and YAP/TAZ regulation. However, studies that directly modulate Ca2+ have reported contradictory YAP/TAZ responses: in some studies, a reduction in Ca2+ influx increases the activity of YAP/TAZ, while in others, an increase in Ca2+ influx activates YAP/TAZ. Importantly, Ca2+ and YAP/TAZ exhibit distinct spatiotemporal dynamics, making it difficult to unravel their connections from a purely experimental approach. In this study, we developed a network model of Ca2+ -mediated YAP/TAZ signalling to investigate how temporal dynamics and crosstalk of signalling pathways interacting with Ca2+ can alter the YAP/TAZ response, as observed in experiments. By including six signalling modules (e.g. GPCR, IP3-Ca2+ , kinases, RhoA, F-actin and Hippo-YAP/TAZ) that interact with Ca2+ , we investigated both transient and steady-state cell response to angiotensin II and thapsigargin stimuli. The model predicts that stimuli, Ca2+ transients and frequency-dependent relationships between Ca2+ and YAP/TAZ are primarily mediated by cPKC, DAG, CaMKII and F-actin. Simulation results illustrate the role of Ca2+ dynamics and CaMKII bistable response in switching the direction of changes in Ca2+ -induced YAP/TAZ activity. A frequency-dependent YAP/TAZ response revealed the competition between upstream regulators of LATS1/2, leading to the YAP/TAZ non-monotonic response to periodic GPCR stimulation. This study provides new insights into underlying mechanisms responsible for the controversial Ca2+ -YAP/TAZ relationship observed in experiments. KEY POINTS: YAP/TAZ integrates biochemical and biomechanical inputs to regulate cellular functions, and Ca2+ acts as a key second messenger linking cellular inputs to YAP/TAZ. Studies have reported contradictory Ca2+ -YAP/TAZ relationships for different cell types and stimuli. A network model of Ca2+ -mediated YAP/TAZ signalling was developed to investigate the underlying mechanisms of divergent Ca2+ -YAP/TAZ relationships. The model predicts context-dependent Ca2+ transient, CaMKII bistable response and frequency-dependent activation of LATS1/2 upstream regulators as mechanisms governing the Ca2+ -YAP/TAZ relationship. This study provides new insights into the underlying mechanisms of the controversial Ca2+ -YAP/TAZ relationship to better understand the dynamics of cellular functions controlled by YAP/TAZ activity.

A common molecular logic determines embryonic stem cell self-renewal and reprogramming.

During differentiation and reprogramming, new cell identities are generated by reconfiguration of gene regulatory networks. Here, we combined automated formal reasoning with experimentation to expose the logic of network activation during induction of naïve pluripotency. We find that a Boolean network architecture defined for maintenance of naïve state embryonic stem cells (ESC) also explains transcription factor behaviour and potency during resetting from primed pluripotency. Computationally identified gene activation trajectories were experimentally substantiated at single-cell resolution by RT-qPCR Contingency of factor availability explains the counterintuitive observation that Klf2, which is dispensable for ESC maintenance, is required during resetting. We tested 124 predictions formulated by the dynamic network, finding a predictive accuracy of 77.4%. Finally, we show that this network explains and predicts experimental observations of somatic cell reprogramming. We conclude that a common deterministic program of gene regulation is sufficient to govern maintenance and induction of naïve pluripotency. The tools exemplified here could be broadly applied to delineate dynamic networks underlying cell fate transitions.

A symptom level perspective on reactive and proactive aggressive behaviours and ADHD symptoms in childhood.

OBJECTIVE: Attention-Deficit/Hyperactivity Disorder (ADHD) is one of the most prevalent childhood disorders, affecting around 3.4% of children worldwide. A common and impairing correlate of ADHD is aggressive behaviour. ADHD symptoms and aggression are both heterogeneous and it has been speculated that certain symptoms of ADHD might be more important in aggressive behaviours of different types than others. This study uses a symptom-level analysis to investigate the concurrent and temporal links between ADHD symptoms and aggressive behaviours. METHODS: Using Gaussian Graphical Models and Graphical Vector Autoregression Models, longitudinal and cross-sectional networks of ADHD symptoms and aggressive behaviours, measured using parent-reported Social Behaviour Questionnaires, were estimated. Participants included 1,246 children taking part in the longitudinal Swiss z-proso cohort study at ages 7, 9 and 11. RESULTS: The longitudinal network highlighted that ADHD symptoms and aggressive behaviours share a multitude of reciprocal temporal relations, with inattentive ADHD symptoms preceding both reactive and proactive aggression. Cross-sectional networks suggested that hyperactive/impulsive symptoms were predominantly connected to reactive aggressive behaviours but also to a form of proactive aggression, namely dominating other children. CONCLUSION: Findings provide preliminary evidence which specific symptoms are the most promising targets for reducing aggressive behaviours in children with ADHD. They also highlight the potential importance of targeting feedback loops resulting from aggressive behaviours. Future research is needed to better understand the mechanisms through which ADHD and aggressive behaviours become linked.

Hierarchical modelling of functional brain networks in population and individuals from big fMRI data.

A major goal of large-scale brain imaging datasets is to provide resources for investigating heterogeneous populations. Characterisation of functional brain networks for individual subjects from these datasets will have an enormous potential for prediction of cognitive or clinical traits. We propose for the first time a technique, Stochastic Probabilistic Functional Modes (sPROFUMO), that is scalable to UK Biobank (UKB) with expected 100,000 participants, and hierarchically estimates functional brain networks in individuals and the population, while allowing for bidirectional flow of information between the two. Using simulations, we show the model's utility, especially in scenarios that involve significant cross-subject variability, or require delineation of fine-grained differences between the networks. Subsequently, by applying the model to resting-state fMRI from 4999 UKB subjects, we mapped resting state networks (RSNs) in single subjects with greater detail than has been possible previously in UKB (>100 RSNs), and demonstrate that these RSNs can predict a range of sensorimotor and higher-level cognitive functions. Furthermore, we demonstrate several advantages of the model over independent component analysis combined with dual-regression (ICA-DR), particularly with respect to the estimation of the spatial configuration of the RSNs and the predictive power for cognitive traits. The proposed model and results can open a new door for future investigations into individualised profiles of brain function from big data.

The importance of transdiagnostic symptom level assessment to understanding prognosis for depressed adults: analysis of data from six randomised control trials.

BACKGROUND: Depression is commonly perceived as a single underlying disease with a number of potential treatment options. However, patients with major depression differ dramatically in their symptom presentation and comorbidities, e.g. with anxiety disorders. There are also large variations in treatment outcomes and associations of some anxiety comorbidities with poorer prognoses, but limited understanding as to why, and little information to inform the clinical management of depression. There is a need to improve our understanding of depression, incorporating anxiety comorbidity, and consider the association of a wide range of symptoms with treatment outcomes. METHOD: Individual patient data from six RCTs of depressed patients (total n = 2858) were used to estimate the differential impact symptoms have on outcomes at three post intervention time points using individual items and sum scores. Symptom networks (graphical Gaussian model) were estimated to explore the functional relations among symptoms of depression and anxiety and compare networks for treatment remitters and those with persistent symptoms to identify potential prognostic indicators. RESULTS: Item-level prediction performed similarly to sum scores when predicting outcomes at 3 to 4 months and 6 to 8 months, but outperformed sum scores for 9 to 12 months. Pessimism emerged as the most important predictive symptom (relative to all other symptoms), across these time points. In the network structure at study entry, symptoms clustered into physical symptoms, cognitive symptoms, and anxiety symptoms. Sadness, pessimism, and indecision acted as bridges between communities, with sadness and failure/worthlessness being the most central (i.e. interconnected) symptoms. Connectivity of networks at study entry did not differ for future remitters vs. those with persistent symptoms. CONCLUSION: The relative importance of specific symptoms in association with outcomes and the interactions within the network highlight the value of transdiagnostic assessment and formulation of symptoms to both treatment and prognosis. We discuss the potential for complementary statistical approaches to improve our understanding of psychopathology.

Incorporating RNA-Seq transcriptomics into glycosylation-integrating metabolic network modelling kinetics: Multiomic Chinese hamster ovary (CHO) cell bioreactors.

In this work, the kinetic model based on the previously developed metabolic and glycan reaction networks of the ovarian cells of the Chinese hamster ovary (CHO) cell line was improved by the inclusion of transcriptomic data that took into account the values of the RPKM gene (Reads per Kilobase of Exon per Million Reads Mapped). The transcriptomic (RNASeq) data were obtained together with metabolic and glycan data from the literature, and the concentrations with RPKM values were collected at several points in time from two fed-batch processes. First, the fluxes were determined by regression analysis of the metabolic data, then these fluxes were corrected by using the fold change in gene expression as a measure of enzyme concentrations. Next, the corrected fluxes in the kinetic model were used to calculate the concentration profiles of the metabolites, and literature data were used to evaluate the predicted results of the model. Compared to other studies where the concentration profiles of CHO cell metabolites were described using a kinetic model without consideration of RNA-Seq data to correct the fluxes, this model is unique. The additional integration of transcriptomic data led to better predictions of metabolic concentrations in the fed-batch process, which is a significant improvement of the modelling technique used.

Effects of social structure and management on risk of disease establishment in wild pigs.

Contact heterogeneity among hosts determines invasion and spreading dynamics of infectious disease, thus its characterization is essential for identifying effective disease control strategies. Yet, little is known about the factors shaping contact networks in many wildlife species and how wildlife management actions might affect contact networks. Wild pigs in North America are an invasive, socially structured species that pose a health concern for domestic swine given their ability to transmit numerous devastating diseases such as African swine fever (ASF). Using proximity loggers and GPS data from 48 wild pigs in Florida and South Carolina, USA, we employed a probabilistic framework to estimate weighted contact networks. We determined the effects of sex, social group and spatial distribution (monthly home-range overlap and distance) on wild pig contact. We also estimated the impacts of management-induced perturbations on contact and inferred their effects on ASF establishment in wild pigs with simulation. Social group membership was the primary factor influencing contacts. Between-group contacts depended primarily on space use characteristics, with fewer contacts among groups separated by >2 km and no contacts among groups >4 km apart within a month. Modelling ASF dynamics on the contact network demonstrated that indirect contacts resulting from baiting (a typical method of attracting wild pigs or game species to a site to enhance recreational hunting) increased the risk of disease establishment by ~33% relative to direct contact. Low-intensity population reduction (<5.9% of the population) had no detectable impact on contact structure but reduced predicted ASF establishment risk relative to no population reduction. We demonstrate an approach for understanding the relative role of spatial, social and individual-level characteristics in shaping contact networks and predicting their effects on disease establishment risk, thus providing insight for optimizing disease control in spatially and socially structured wildlife species.

Online network monitoring.

An important problem in network analysis is the online detection of anomalous behaviour. In this paper, we introduce a network surveillance method bringing together network modelling and statistical process control. Our approach is to apply multivariate control charts based on exponential smoothing and cumulative sums in order to monitor networks generated by temporal exponential random graph models (TERGM). The latter allows us to account for temporal dependence while simultaneously reducing the number of parameters to be monitored. The performance of the considered charts is evaluated by calculating the average run length and the conditional expected delay for both simulated and real data. To justify the decision of using the TERGM to describe network data, some measures of goodness of fit are inspected. We demonstrate the effectiveness of the proposed approach by an empirical application, monitoring daily flights in the United States to detect anomalous patterns.

Network analysis reveals underlying syntactic features in a vocally learnt mammalian display, humpback whale song.

Vocal communication systems have a set of rules that govern the arrangement of acoustic signals, broadly defined as 'syntax'. However, there is a limited understanding of potentially shared or analogous rules across vocal displays in different taxa. Recent work on songbirds has investigated syntax using network-based modelling. This technique quantifies features such as connectivity (adjacent signals in a sequence) and recurring patterns. Here, we apply network-based modelling to the complex, hierarchically structured songs of humpback whales (Megaptera novaeangliae) from east Australia. Given the song's annual evolving pattern and the cultural conformity of males within a population, network modelling captured the patterns of multiple song types over 13 consecutive years. Song arrangements in each year displayed clear 'small-world' network structure, characterized by clusters of highly connected sounds. Transitions between these connected sounds further suggested a combination of both structural stability and variability. Small-world network structure within humpback songs may facilitate the characteristic and persistent vocal learning observed. Similar small-world structures and transition patterns are found in several birdsong displays, indicating common syntactic patterns among vocal learning in multiple taxa. Understanding the syntactic rules governing vocal displays in multiple, independently evolving lineages may indicate what rules or structural features are important to the evolution of complex communication, including human language.

Distinct resting-state functional connections associated with episodic and visuospatial memory in older adults.

Episodic and spatial memory are commonly impaired in ageing and Alzheimer's disease. Volumetric and task-based functional magnetic resonance imaging (fMRI) studies suggest a preferential involvement of the medial temporal lobe (MTL), particularly the hippocampus, in episodic and spatial memory processing. The present study examined how these two memory types were related in terms of their associated resting-state functional architecture. 3T multiband resting state fMRI scans from 497 participants (60-82 years old) of the cross-sectional Whitehall II Imaging sub-study were analysed using an unbiased, data-driven network-modelling technique (FSLNets). Factor analysis was performed on the cognitive battery; the Hopkins Verbal Learning test and Rey-Osterreith Complex Figure test factors were used to assess verbal and visuospatial memory respectively. We present a map of the macroscopic functional connectome for the Whitehall II Imaging sub-study, comprising 58 functionally distinct nodes clustered into five major resting-state networks. Within this map we identified distinct functional connections associated with verbal and visuospatial memory. Functional anticorrelation between the hippocampal formation and the frontal pole was significantly associated with better verbal memory in an age-dependent manner. In contrast, hippocampus-motor and parietal-motor functional connections were associated with visuospatial memory independently of age. These relationships were not driven by grey matter volume and were unique to the respective memory domain. Our findings provide new insights into current models of brain-behaviour interactions, and suggest that while both episodic and visuospatial memory engage MTL nodes of the default mode network, the two memory domains differ in terms of the associated functional connections between the MTL and other resting-state brain networks.

Biphasic regulation of the transcription factor ABORTED MICROSPORES (AMS) is essential for tapetum and pollen development in Arabidopsis.

Viable pollen is essential for plant reproduction and crop yield. Its production requires coordinated expression at specific stages during anther development, involving early meiosis-associated events and late pollen wall formation. The ABORTED MICROSPORES (AMS) transcription factor is a master regulator of sporopollenin biosynthesis, secretion and pollen wall formation in Arabidopsis. Here we show that it has complex regulation and additional essential roles earlier in pollen formation. An inducible-AMS reporter was created for functional rescue, protein expression pattern analysis, and to distinguish between direct and indirect targets. Mathematical modelling was used to create regulatory networks based on wild-type RNA and protein expression. Dual activity of AMS was defined by biphasic protein expression in anther tapetal cells, with an initial peak around pollen meiosis and then later during pollen wall development. Direct AMS-regulated targets exhibit temporal regulation, indicating that additional factors are associated with their regulation. We demonstrate that AMS biphasic expression is essential for pollen development, and defines distinct functional activities during early and late pollen development. Mathematical modelling suggests that AMS may competitively form a protein complex with other tapetum-expressed transcription factors, and that biphasic regulation is due to repression of upstream regulators and promotion of AMS protein degradation.

Systems genetics: From GWAS to disease pathways.

Most common diseases are complex, involving multiple genetic and environmental factors and their interactions. In the past decade, genome-wide association studies (GWAS) have successfully identified thousands of genetic variants underlying susceptibility to complex diseases. However, the results from these studies often do not provide evidence on how the variants affect downstream pathways and lead to the disease. Therefore, in the post-GWAS era the greatest challenge lies in combining GWAS findings with additional molecular data to functionally characterize the associations. The advances in various ~omics techniques have made it possible to investigate the effect of risk variants on intermediate molecular levels, such as gene expression, methylation, protein abundance or metabolite levels. As disease aetiology is complex, no single molecular analysis is expected to fully unravel the disease mechanism. Multiple molecular levels can interact and also show plasticity in different physiological conditions, cell types and disease stages. There is therefore a great need for new integrative approaches that can combine data from different molecular levels and can help construct the causal inference from genotype to phenotype. Systems genetics is such an approach; it is used to study genetic effects within the larger scope of systems biology by integrating genotype information with various ~omics datasets as well as with environmental and physiological variables. In this review, we describe this approach and discuss how it can help us unravel the molecular mechanisms through which genetic variation causes disease. This article is part of a Special Issue entitled: From Genome to Function.

Nodal approach reveals differential impact of lateralized focal epilepsies on hub reorganization.

The impact of the hemisphere affected by impairment in models of network disease is not fully understood. Among such models, focal epilepsies are characterised by recurrent seizures generated in epileptogenic areas also responsible for wider network dysfunction between seizures. Previous work focusing on functional connectivity within circumscribed networks suggests a divergence of network integrity and compensatory capacity between epilepsies as a function of the laterality of seizure onset. We evaluated the ability of complex network theory to reveal changes in focal epilepsy in global and nodal parameters using graph theoretical analysis of functional connectivity data obtained with resting-state fMRI. Graphs of functional connectivity networks were derived from 19 right and 13 left focal epilepsy patients and 15 controls. Topological metrics (degree, local efficiency, global efficiency and modularity) were computed for a whole-brain, atlas-defined network. We also calculated a hub disruption index for each graph metric, measuring the capacity of the brain network to demonstrate increased connectivity in some nodes for decreased connectivity in others. Our data demonstrate that the patient group as a whole is characterised by network-wide pattern of reorganization, even while global parameters fail to distinguish between groups. Furthermore, multiple metrics indicate that epilepsies with differently lateralized epileptic networks are asymmetric in their burden on functional brain networks; with left epilepsy patients being characterised by reduced efficiency and modularity, while in right epilepsy patients we provide the first evidence that functional brain networks are characterised by enhanced connectivity and efficiency at some nodes whereas reduced in others.

Raccoon contact networks predict seasonal susceptibility to rabies outbreaks and limitations of vaccination.

Infectious disease transmission often depends on the contact structure of host populations. Although it is often challenging to capture the contact structure in wild animals, new technology has enabled biologists to obtain detailed temporal information on wildlife social contacts. In this study, we investigated the effects of raccoon contact patterns on rabies spread using network modelling. Raccoons (Procyon lotor) play an important role in the maintenance of rabies in the United States. It is crucial to understand how contact patterns influence the spread of rabies in raccoon populations in order to design effective control measures and to prevent transmission to human populations and other animals. We constructed a dynamic system of contact networks based on empirical data from proximity logging collars on a wild suburban raccoon population and then simulated rabies spread across these networks. Our contact networks incorporated the number and duration of raccoon interactions. We included differences in contacts according to sex and season, and both short-term acquaintances and long-term associations. Raccoons may display different behaviours when infectious, including aggression (furious behaviour) and impaired mobility (dumb behaviour); the network model was used to assess the impact of potential behavioural changes in rabid raccoons. We also tested the effectiveness of different vaccination coverage levels. Our results demonstrate that when rabies enters a suburban raccoon population, the likelihood of a disease outbreak affecting the majority of the population is high. Both the magnitude of rabies outbreaks and the speed of rabies spread depend strongly on the time of year that rabies is introduced into the population. When there is a combination of dumb and furious behaviours in the rabid raccoon population, there are similar outbreak sizes and speed of spread to when there are no behavioural changes due to rabies infection. By incorporating detailed data describing the variation in raccoon contact rates into a network modelling approach, we were able to show that suburban raccoon populations are highly susceptible to rabies outbreaks, that the risk of large outbreaks varies seasonally and that current vaccination target levels may be inadequate to prevent the spread of rabies within these populations. Our findings provide new insights into rabies dynamics in raccoon populations and have important implications for disease control.

Creating a Sustainability Scorecard as a predictive tool for measuring the complex social, economic and environmental impacts of industries, a case study: assessing the viability and sustainability of the dairy industry.

Sustainability is a key driver for decisions in the management and future development of industries. The World Commission on Environment and Development (WCED, 1987) outlined imperatives which need to be met for environmental, economic and social sustainability. Development of strategies for measuring and improving sustainability in and across these domains, however, has been hindered by intense debate between advocates for one approach fearing that efforts by those who advocate for another could have unintended adverse impacts. Studies attempting to compare the sustainability performance of countries and industries have also found ratings of performance quite variable depending on the sustainability indices used. Quantifying and comparing the sustainability of industries across the triple bottom line of economy, environment and social impact continues to be problematic. Using the Australian dairy industry as a case study, a Sustainability Scorecard, developed as a Bayesian network model, is proposed as an adaptable tool to enable informed assessment, dialogue and negotiation of strategies at a global level as well as being suitable for developing local solutions.

Exploring gene regulation and biological processes in insects: Insights from omics data using gene regulatory network models.

Gene regulatory network (GRN) comprises complicated yet intertwined gene-regulator relationships. Understanding the GRN dynamics will unravel the complexity behind the observed gene expressions. Insect gene regulation is often complicated due to their complex life cycles and diverse ecological adaptations. The main interest of this review is to have an update on the current mathematical modelling methods of GRNs to explain insect science. Several popular GRN architecture models are discussed, together with examples of applications in insect science. In the last part of this review, each model is compared from different aspects, including network scalability, computation complexity, robustness to noise and biological relevancy.

Application of Multivariate Regression and Artificial Neural Network Modelling for Prediction of Physicochemical Properties of Grape-Skin Compost.

The reusability of by-products in the food industry is consistent with sustainable and greener production; therefore, the aim of this paper was to evaluate the applicability of multiple linear regression (MLR), piecewise linear regression (PLR) and artificial neural network models (ANN) to the prediction of grape-skin compost's physicochemical properties (moisture, dry matter, organic matter, ash content, carbon content, nitrogen content, C/N ratio, total colour change of compost samples, pH, conductivity, total dissolved solids and total colour change of compost extract samples) during in-vessel composting based on the initial composting conditions (air-flow rate, moisture content and day of sampling). Based on the coefficient of determination for prediction, the adjusted coefficient of determination for calibration, the root-mean-square error of prediction (RMSEP), the standard error of prediction (SEP), the ratio of prediction to deviation (RPD) and the ratio of the error range (RER), it can be concluded that all developed MLR and PLR models are acceptable for process screening. Furthermore, the ANN model developed for predicting moisture and dry-matter content can be used for quality control (RER >11). The obtained results show the great potential of multivariate modelling for analysis of the physicochemical properties of compost during composting, confirming the high applicability of modelling in greener production processes.

On the characteristics of natural hydraulic dampers: An image-based approach to study the fluid flow behaviour inside the human meniscal tissue.

The meniscal tissue is a layered material with varying properties influenced by collagen content and arrangement. Understanding the relationship between structure and properties is crucial for disease management, treatment development, and biomaterial design. The internal layer of the meniscus is softer and more deformable than the outer layers, thanks to interconnected collagen channels that guide fluid flow. To investigate these relationships, we propose an integrated approach that combines Computational Fluid Dynamics (CFD) with Image Analysis (CFD-IA). We analyze fluid flow in the internal architecture of the human meniscus across a range of inlet velocities (0.1 mm/s to 1.6 m/s) using high-resolution 3D micro-computed tomography scans. Statistical correlations are observed between architectural parameters (tortuosity, connectivity, porosity, pore size) and fluid flow parameters (Re number distribution, permeability). Some channels exhibit Re values of 1400 at an inlet velocity of 1.6 m/s, and a transition from Darcy's regime to a non-Darcian regime occurs around an inlet velocity of 0.02 m/s. Location-dependent permeability ranges from 20-32 Darcy. Regression modelling reveals a strong correlation between fluid velocity and tortuosity at high inlet velocities, as well as with channel diameter at low inlet velocities. At higher inlet velocities, flow paths deviate more from the preferential direction, resulting in a decrease in the concentration parameter by an average of 0.4. This research provides valuable insights into the fluid flow behaviour within the meniscus and its structural influences. 3D models and image stack are available to download at https://doi.org/10.5281/zenodo.10401592. STATEMENT OF SIGNIFICANCE: The meniscus is a highly porous soft tissue with remarkable properties of load transfer and energy absorption. We give insight on the mechanism of energy absorption from high resolution uCT scans, never presented before, and a new method which combine CFD and image. The structure is similar to a sandwich structure with a stiff outside layer and a soft internal layer made of collagen channels oriented in a preferential direction guiding the fluid flow, enabling it to accommodate deformation and dissipate energy, making it a potentially optimized damping system. We investigate architectural/ fluid flow parameters- fluid regimes relationship, which is of interest of the readers working on designing suitable biomimetic systems that can be adopted for replacement.