Wybutosine hypomodification of tRNAphe activates HERVK and impairs neuronal differentiation.

We previously reported that loss of function of TYW1 led to cerebral palsy with severe intellectual disability through reduced neural proliferation. However, whether TYW1 loss affects neural differentiation is unknown. In this study, we first demonstrated that TYW1 loss blocked the formation of OHyW in tRNAphe and therefore affected the translation efficiency of UUU codon. Using the brain organoid model, we showed impaired neuron differentiation when TYW1 was depleted. Interestingly, retrotransposons were differentially regulated in TYW1-/- hESCs (human embryonic stem cells). In particular, one kind of human-specific endogenous retrovirus-K (HERVK/HML2), whose reactivation impaired human neurodevelopment, was significantly up-regulated in TYW1-/- hESCs. Consistently, a UUU codon-enriched protein, SMARCAD1, which was a key factor in controlling endogenous retroviruses, was reduced. Taken together, TYW1 loss leads to up-regulation of HERVK in hESCs by down-regulated SMARCAD1, thus impairing neuron differentiation.

A Multi-Information Spreading Model for One-Time Retweet Information in Complex Networks.

In the realm of online social networks, the spreading of information is influenced by a complex interplay of factors. To explore the dynamics of one-time retweet information spreading, we propose a Susceptible-Infected-Completed (SIC) multi-information spreading model. This model captures how multiple pieces of information interact in online social networks by introducing inhibiting and enhancement factors. The SIC model considers the completed state, where nodes cease to spread a particular piece of information after transmitting it. It also takes into account the impact of past and present information received from neighboring nodes, dynamically calculating the probability of nodes spreading each piece of information at any given moment. To analyze the dynamics of multiple information pieces in various scenarios, such as mutual enhancement, partial competition, complete competition, and coexistence of competition and enhancement, we conduct experiments on BA scale-free networks and the Twitter network. Our findings reveal that competing information decreases the likelihood of its spread while cooperating information amplifies the spreading of mutually beneficial content. Furthermore, the strength of the enhancement factor between different information pieces determines their spread when competition and cooperation coexist. These insights offer a fresh perspective for understanding the patterns of information propagation in multiple contexts.

Technical strategy for monozygotic twin discrimination by single-nucleotide variants.

Monozygotic (MZ) twins are theoretically genetically identical. Although they are revealed to accumulate mutations after the zygote splits, discriminating between twin genomes remains a formidable challenge in the field of forensic genetics. Single-nucleotide variants (SNVs) are responsible for a substantial portion of genetic variation, thus potentially serving as promising biomarkers for the identification of MZ twins. In this study, we sequenced the whole genome of a pair of female MZ twins when they were 27 and 33 years old to approximately 30 × coverage using peripheral blood on an Illumina NovaSeq 6000 Sequencing System. Potentially discordant SNVs supported by whole-genome sequencing were validated extensively by amplicon-based targeted deep sequencing and Sanger sequencing. In total, we found nine bona fide post-twinning SNVs, all of which were identified in the younger genomes and found in the older genomes. None of the SNVs occurred within coding exons, three of which were observed in introns, supported by whole-exome sequencing results. A double-blind test was employed, and the reliability of MZ twin discrimination by discordant SNVs was endorsed. All SNVs were successfully detected when input DNA amounts decreased to 0.25 ng, and reliable detection was limited to seven SNVs below 0.075 ng input. This comprehensive analysis confirms that SNVs could serve as cost-effective biomarkers for MZ twin discrimination.

Altered intestinal microbiota enhances adenoid hypertrophy by disrupting the immune balance.

Introduction: Adenoid hypertrophy (AH) is a common upper respiratory disorder in children. Disturbances of gut microbiota have been implicated in AH. However, the interplay of alteration of gut microbiome and enlarged adenoids remains elusive. Methods: 119 AH children and 100 healthy controls were recruited, and microbiome profiling of fecal samples in participants was performed using 16S rRNA gene sequencing. Fecal microbiome transplantation (FMT) was conducted to verify the effects of gut microbiota on immune response in mice. Results: In AH individuals, only a slight decrease of diversity in bacterial community was found, while significant changes of microbial composition were observed between these two groups. Compared with HCs, decreased abundances of Akkermansia, Oscillospiraceae and Eubacterium coprostanoligenes genera and increased abundances of Bacteroides, Faecalibacterium, Ruminococcus gnavus genera were revealed in AH patients. The abundance of Bacteroides remained stable with age in AH children. Notably, a microbial marker panel of 8 OTUs were identified, which discriminated AH from HC individuals with an area under the curve (AUC) of 0.9851 in the discovery set, and verified in the geographically different validation set, achieving an AUC of 0.9782. Furthermore, transfer of mice with fecal microbiota from AH patients dramatically reduced the proportion of Treg subsets within peripheral blood and nasal-associated lymphoid tissue (NALT) and promoted the expansion of Th2 cells in NALT. Conclusion: These findings highlight the effect of the altered gut microbiota in the AH pathogenesis.

Recent Advances in Bacterial Persistence Mechanisms.

The recurrence of bacterial infectious diseases is closely associated with bacterial persisters. This subpopulation of bacteria can escape antibiotic treatment by entering a metabolic status of low activity through various mechanisms, for example, biofilm, toxin-antitoxin modules, the stringent response, and the SOS response. Correspondingly, multiple new treatments are being developed. However, due to their spontaneous low abundance in populations and the lack of research on in vivo interactions between persisters and the host's immune system, microfluidics, high-throughput sequencing, and microscopy techniques are combined innovatively to explore the mechanisms of persister formation and maintenance at the single-cell level. Here, we outline the main mechanisms of persister formation, and describe the cutting-edge technology for further research. Despite the significant progress regarding study techniques, some challenges remain to be tackled.

Integrative analysis of single-cell embryo data reveals transcriptome signatures for the human pre-implantation inner cell mass.

As the source of embryonic stem cells (ESCs), inner cell mass (ICM) can form all tissues of the embryo proper, however, its role in early human lineage specification remains controversial. Although a stepwise differentiation model has been proposed suggesting the existence of ICM as a distinct developmental stage, the underlying molecular mechanism remains unclear. In the present study, we perform an integrated analysis on the public human preimplantation embryonic single-cell transcriptomic data and apply a trajectory inference algorithm to measure the cell plasticity. In our results, ICM population can be clearly discriminated on the dimension-reduced graph and confirmed by compelling evidences, thus validating the two-step hypothesis of lineage commitment. According to the branch probabilities and differentiation potential, we determine the precise time points for two lineage segregations. Further analysis on gene expression dynamics and regulatory network indicates that transcription factors including GSC, PRDM1, and SPIC may underlie the decisions of ICM fate. In addition, new human ICM marker genes, such as EPHA4 and CCR8 are discovered and validated by immunofluorescence. Given the potential clinical applications of ESCs, our analysis provides a further understanding of human ICM cells and facilitates the exploration of more unique characteristics in early human development.

Solute Transport across the Lymphatic Vasculature in a Soft Skin Tissue.

Convective transport of drug solutes in biological tissues is regulated by the interstitial fluid pressure, which plays a crucial role in drug absorption into the lymphatic system through the subcutaneous (SC) injection. In this paper, an approximate continuum poroelasticity model is developed to simulate the pressure evolution in the soft porous tissue during an SC injection. This poroelastic model mimics the deformation of the tissue by introducing the time variation of the interstitial fluid pressure. The advantage of this method lies in its computational time efficiency and simplicity, and it can accurately model the relaxation of pressure. The interstitial fluid pressure obtained using the proposed model is validated against both the analytical and the numerical solution of the poroelastic tissue model. The decreasing elasticity elongates the relaxation time of pressure, and the sensitivity of pressure relaxation to elasticity decreases with the hydraulic permeability, while the increasing porosity and permeability due to deformation alleviate the high pressure. An improved Kedem-Katchalsky model is developed to study solute transport across the lymphatic vessel network, including convection and diffusion in the multi-layered poroelastic tissue with a hybrid discrete-continuum vessel network embedded inside. At last, the effect of different structures of the lymphatic vessel network, such as fractal trees and Voronoi structure, on the lymphatic uptake is investigated. In this paper, we provide a novel and time-efficient computational model for solute transport across the lymphatic vasculature connecting the microscopic properties of the lymphatic vessel membrane to the macroscopic drug absorption.

Clinical Application of a Multiplex Droplet Digital PCR in the Rapid Diagnosis of Children with Suspected Bloodstream Infections.

Droplet digital PCR (ddPCR) recently has been shown to be a potential diagnostic tool for adults with bloodstream infections (BSIs); however, its application in children remains obscure. In this study, 76 blood samples of children with suspected BSIs were synchronously detected by traditional blood cultures (BCs) and ddPCRs. Our team validated the diagnostic performance of ddPCR including sensitivity, specificity, and positive and negative predictive values. The 76 pediatric patients from the hematology department (67.1%), the pediatric intensive care unit (PICU, 27.6%), and other departments (5.2%) were enrolled. The positive rate of ddPCR results was 47.9%, whereas that for BC was 6.6%. In addition, the time consumption of ddPCR was shorter, only for 4.7 ± 0.9 h, in comparison with the detection timing of BC (76.7 ± 10.4 h, p < 0.01). The levels of agreement and disagreement between BC and ddPCR were 96.1% and 4.2%, and the negative agreement reached 95.6%. The sensitivity of ddPCR was 100%, with corresponding specificities ranging from 95.3 to 100.0%. In addition, a total of nine viruses were identified by ddPCR. In China, the multiplexed ddPCR first could be a tool for the rapid and accurate diagnosis of children with suspected BSIs and can be an early indicator of the possibility of viraemia in children with immunosuppression.

Epigenetic regulation of human-specific gene expression in the prefrontal cortex.

BACKGROUND: Changes in gene expression levels during brain development are thought to have played an important role in the evolution of human cognition. With the advent of high-throughput sequencing technologies, changes in brain developmental expression patterns, as well as human-specific brain gene expression, have been characterized. However, interpreting the origin of evolutionarily advanced cognition in human brains requires a deeper understanding of the regulation of gene expression, including the epigenomic context, along the primate genome. Here, we used chromatin immunoprecipitation sequencing (ChIP-seq) to measure the genome-wide profiles of histone H3 lysine 4 trimethylation (H3K4me3) and histone H3 lysine 27 acetylation (H3K27ac), both of which are associated with transcriptional activation in the prefrontal cortex of humans, chimpanzees, and rhesus macaques. RESULTS: We found a discrete functional association, in which H3K4me3HP gain was significantly associated with myelination assembly and signaling transmission, while H3K4me3HP loss played a vital role in synaptic activity. Moreover, H3K27acHP gain was enriched in interneuron and oligodendrocyte markers, and H3K27acHP loss was enriched in CA1 pyramidal neuron markers. Using strand-specific RNA sequencing (ssRNA-seq), we first demonstrated that approximately 7 and 2% of human-specific expressed genes were epigenetically marked by H3K4me3HP and H3K27acHP, respectively, providing robust support for causal involvement of histones in gene expression. We also revealed the co-activation role of epigenetic modification and transcription factors in human-specific transcriptome evolution. Mechanistically, histone-modifying enzymes at least partially contribute to an epigenetic disturbance among primates, especially for the H3K27ac epigenomic marker. In line with this, peaks enriched in the macaque lineage were found to be driven by upregulated acetyl enzymes. CONCLUSIONS: Our results comprehensively elucidated a causal species-specific gene-histone-enzyme landscape in the prefrontal cortex and highlighted the regulatory interaction that drove transcriptional activation.

Recursive traffic percolation on urban transportation systems.

This paper proposes a recursive traffic percolation framework to capture the dynamics of cascading failures and analyze potential overloaded bottlenecks. In particular, compared to current work, the influence of external flow is considered, providing a new perspective for the study of regional commuting. Finally, we present an empirical study to verify the accuracy and effectiveness of our framework. Further analysis indicates that external flows from different regions affect the network. Our work requires only primary data and verifies the improvement of the functional network.

Association between Intestinal Colonization and Extraintestinal Infection with Carbapenem-Resistant Klebsiella pneumoniae in Children.

Carbapenem-resistant Klebsiella pneumoniae (CRKP) has become a critical public health threat. However, the association between intestinal colonization and parenteral infection among pediatric patients has not been elucidated. We collected 8 fecal CRKP strains and 10 corresponding CRKP strains responsible for extraintestinal infection from eight patients who did not manifest infection upon admission to the hospital. Paired isolates showed identical resistance to antimicrobials and identical virulence in vitro and in vivo. wzi capsule typing, multilocus sequence typing, and whole-genome sequencing (WGS) indicated high similarity between paired colonizing and infecting isolates. Mutations between colonizing and infecting isolate pairs found by WGS had a distinctive molecular signature of a high proportion of complex structural variants. The mutated genes were involved in pathways associated with infection-related physiological and pathogenic functions, including antibiotic resistance, virulence, and response to the extracellular environment. The latter is important for bacterial infection of environmental niches. Various mutations related to antibiotic resistance, virulence, and colonization that were not associated with any particular mutational hot spot correlated with an increased risk of extraintestinal infection. Notably, novel subclone carbapenem-resistant hypervirulent K. pneumoniae (CR-hvKP) KL19-ST15 exhibited hypervirulence in experimental assays that reflected the severe clinical symptoms of two patients infected with the clonal strains. Taken together, our findings indicate the association between CRKP intestinal colonization and extraintestinal infection, suggesting that active screening for colonization on admission could decrease infection risk in children. IMPORTANCE Carbapenem-resistant Klebsiella pneumoniae (CRKP) causes an increasing number of nosocomial infections, which can be life-threatening, as carbapenems are last-resort antibiotics. K. pneumoniae is part of the healthy human microbiome, and this provides a potential advantage for infection. This study demonstrated that CRKP intestinal colonization is strongly linked to extraintestinal infection, based on the evidence given by whole-genome sequencing data and phenotypic assays of antimicrobial resistance and virulence. Apart from these findings, our in-depth analysis of point mutations and chromosome structural variants in patient-specific infecting isolates compared with colonizing isolates may contribute insights into bacterial adaptation underlying CRKP infection. In addition, a novel subclone of carbapenem-resistant hypervirulent K. pneumoniae (CR-hvKP) was observed in the study. This finding highlights the importance of CRKP active surveillance among children, targeting in particular the novel high-risk CR-hvKP clone.

Lymphatic uptake of biotherapeutics through a 3D hybrid discrete-continuum vessel network in the skin tissue.

Subcutaneous administration is a common approach for the delivery of biotherapeutics, which is achieved mainly through the absorption across lymphatic vessels. In this paper, the drug transport and lymphatic uptake through a three-dimensional hybrid discrete-continuum vessel network in the skin tissue are investigated through high-fidelity numerical simulations. We find that the local lymphatic uptake through the explicit vessels significantly affects macroscopic drug absorption. The diffusion of drug solute through the explicit vessel network affects the lymphatic uptake after the injection. This effect, however, cannot be captured using previously developed continuum models. The lymphatic uptake is dominated by the convection due to lymphatic drainage driven by the pressure difference, which is rarely studied in experiments and simulations. Furthermore, the effects of injection volume and depth on the lymphatic uptake are investigated in a multi-layered domain. We find that the injection volume significantly affects the rate of lymphatic uptake through the heterogeneous vessel network, while the injection depth has little influence, which is consistent with the experimental results. At last, the binding and metabolism of drug molecules are studied to bridge the simulations to the drug clearance experients. We provide a new approach to study the diffusion and convection of drug molecules into the lymphatic system through the hybrid vessel network.

The gut microbiome and obstructive sleep apnea syndrome in children.

Obstructive sleep apnea syndrome (OSAS) in children has become a major public health problem that affects the physical and mental growth of children. OSAS can result in adverse outcomes during growth and development, inhibiting the normal development of the metabolic, cardiovascular, and immune systems. OSAS is characterized by partial or complete obstruction of the upper airway, and prolonged obstruction that causes intermittent hypoxia and sleep fragmentation in children. The human microbiota is a complex community that is in dynamic equilibrium in the human body. Intermittent hypoxia and sleep fragmentation induced by childhood OSAS alter the composition of the gut microbiome. At the same time, changes in the gut microbiome affect sleep patterns in children through immunomodulatory and metabolic mechanisms, and induce further comorbidities, such as obesity, hypertension, and cardiovascular disease. This article discusses recent progress in research into the mechanisms of OSAS-induced changes in the gut microbiota and its pathophysiology in children.

Nosocomial dissemination of hypervirulent Klebsiella pneumoniae with high-risk clones among children in Shanghai.

Objectives: Although hypervirulent Klebsiella pneumoniae (hvKp) is an increasing public health problem, there remains limited epidemiological information regarding hvKp infections in children. Here, we conducted a clinical, molecular and phenotypic surveillance of hvKp strains in a pediatric population. Methods: Non-repetitive K. pneumoniae (Kp) strains consecutively collected during 2019-2020 were screened for hypervirulence genes (prmpA, prmpA2, iucA, iroB, and peg344) using PCR. Positive strains were further characterized by four phenotypic assays (string test, serum killing assay, siderophore production, Galleria mellonella lethality assay), followed by murine sepsis model to determine virulence in vitro and in vivo. Also, capsular types, sequence types, plasmid replicon types, antimicrobial resistance determinants and susceptibility were analyzed. Results: A total of 352 isolates were collected, wherein 83 (23.6%) were hypervirulence genes-positive Kp (hgKp). A significant increase in KPC-2-producing KL47-ST11 among hgKp strains was observed, from 5.3% (1/19) in 2019 to 67.6% (25/37) in 2020 (P<.0001), suggesting the potential dissemination of the hybrid virulence and carbapenem-resistance encoding plasmid among children. Further, hgKp isolates were classified into hvKp (n = 27) and hgKp-low virulence (hgKp-Lv) (n = 56) based on virulence phenotypic assays. In hvKp, diverse genetic clones were observed and K1-ST23 or K2-ST25 strains with sensitivity to multiple antibiotics were prevalent (25.9%, 7/27). Compared with hgKp-Lv, hvKp infection had a higher propensity to involve severe pneumonia (22.2% vs. 12.5%) in elder children and significant higher mortality in mice (P = 0.0086). Additionally, either hvKp or hgKp-Lv infections were mostly healthcare-associated and hospital-acquired (74.1% vs. 91.9%). Conclusions: These data suggest that K1-ST23 and K2-ST25 are high-risk clones of hvKp, and the genetic convergence of virulence and carbapenem-resistance is increasing among children. Control measures are needed to prevent the dissemination in clinical settings.

Basins of attraction of chimera states on networks.

Networks of identical coupled oscillators display a remarkable spatiotemporal pattern, the chimera state, where coherent oscillations coexist with incoherent ones. In this paper we show quantitatively in terms of basin stability that stable and breathing chimera states in the original two coupled networks typically have very small basins of attraction. In fact, the original system is dominated by periodic and quasi-periodic chimera states, in strong contrast to the model after reduction, which can not be uncovered by the Ott-Antonsen ansatz. Moreover, we demonstrate that the curve of the basin stability behaves bimodally after the system being subjected to even large perturbations. Finally, we investigate the emergence of chimera states in brain network, through inducing perturbations by stimulating brain regions. The emerged chimera states are quantified by Kuramoto order parameter and chimera index, and results show a weak and negative correlation between these two metrics.

Inter-tissue convergence of gene expression during ageing suggests age-related loss of tissue and cellular identity.

Developmental trajectories of gene expression may reverse in their direction during ageing, a phenomenon previously linked to cellular identity loss. Our analysis of cerebral cortex, lung, liver, and muscle transcriptomes of 16 mice, covering development and ageing intervals, revealed widespread but tissue-specific ageing-associated expression reversals. Cumulatively, these reversals create a unique phenomenon: mammalian tissue transcriptomes diverge from each other during postnatal development, but during ageing, they tend to converge towards similar expression levels, a process we term Divergence followed by Convergence (DiCo). We found that DiCo was most prevalent among tissue-specific genes and associated with loss of tissue identity, which is confirmed using data from independent mouse and human datasets. Further, using publicly available single-cell transcriptome data, we showed that DiCo could be driven both by alterations in tissue cell-type composition and also by cell-autonomous expression changes within particular cell types.

Genetic Characteristics and Microbiological Profile of Hypermucoviscous Multidrug-Resistant Klebsiella variicola Coproducing IMP-4 and NDM-1 Carbapenemases.

We report here a hypermucoviscous, New Delhi metallo-ß-lactamase 1 (NDM-1) and imipenemase 4 (IMP-4) carbapenemases-coproducing Klebsiella variicola isolate obtained from a pediatric patient. This strain was resistant to carbapenems and most other ß-lactams. Although hypermucoviscous, this strain possessed attenuated virulence according to serum killing assay and Galleria mellonella infection model. Notably, two copies of blaNDM-1 were contained on two tandem ISCR1 elements and coexisted with blaIMP-4 in a novel hybrid multidrug resistance plasmid. This is the first description of the coexistence of blaNDM-1 and blaIMP-4 in a single plasmid of hypermucoviscous K. variicola. IMPORTANCE As an important member of the Klebsiella pneumoniae complex, Klebsiella variicola is poorly studied as an emerging human pathogen. We, for the first time, report a unique K. variicola isolated from a pediatric patient in China. This isolate exhibited hypermucoviscosity, a classic hypervirulence characteristic of K. pneumoniae, and contained multiple carbapenem-resistant genes, including blaIMP-1 and blaNDM-1. Interestingly, these antimicrobial resistance genes were located on a novel hybrid plasmid, and our results suggested that this plasmid might have been introduced from K. pneumoniae and undergone a series of integration and recombination evolutionary events. Overall, our study provides more insight into K. variicola and highlights its superior capability to acquire and maintain foreign resistance genes.

In-depth analysis reveals complex molecular aetiology in a cohort of idiopathic cerebral palsy.

Cerebral palsy is the most prevalent physical disability in children; however, its inherent molecular mechanisms remain unclear. In the present study, we performed in-depth clinical and molecular analysis on 120 idiopathic cerebral palsy families, and identified underlying detrimental genetic variants in 45% of these patients. In addition to germline variants, we found disease-related postzygotic mutations in ∼6.7% of cerebral palsy patients. We found that patients with more severe motor impairments or a comorbidity of intellectual disability had a significantly higher chance of harbouring disease-related variants. By a compilation of 114 known cerebral-palsy-related genes, we identified characteristic features in terms of inheritance and function, from which we proposed a dichotomous classification system according to the expression patterns of these genes and associated cognitive impairments. In two patients with both cerebral palsy and intellectual disability, we revealed that the defective TYW1, a tRNA hypermodification enzyme, caused primary microcephaly and problems in motion and cognition by hindering neuronal proliferation and migration. Furthermore, we developed an algorithm and demonstrated in mouse brains that this malfunctioning hypermodification specifically perturbed the translation of a subset of proteins involved in cell cycling. This finding provided a novel and interesting mechanism for congenital microcephaly. In another cerebral palsy patient with normal intelligence, we identified a mitochondrial enzyme GPAM, the hypomorphic form of which led to hypomyelination of the corticospinal tract in both human and mouse models. In addition, we confirmed that the aberrant Gpam in mice perturbed the lipid metabolism in astrocytes, resulting in suppressed astrocytic proliferation and a shortage of lipid contents supplied for oligodendrocytic myelination. Taken together, our findings elucidate novel aspects of the aetiology of cerebral palsy and provide insights for future therapeutic strategies.

Transport and Lymphatic Uptake of Biotherapeutics Through Subcutaneous Injection.

Drug transport and uptake in the subcutaneous tissue receives increasing attention in biomechanical and pharmaceutical researches, as subcutaneous administration becomes a common approach for the delivery of biotherapeutics, such as monoclonal antibodies. In this paper, high-fidelity numerical simulations are used to investigate the mechanisms governing drug transport and absorption in the subcutaneous tissue, which is expressed as a porous medium modeled by Darcy's law. The effects of tissue properties (permeability and porosity), the injection flow rate, and the vascular permeability of lymphatic vessels on the lymphatic uptake are studied. Additionally, an empirical formula for the lymphatic uptake during the injection is developed based on the numerical results. The roles of lymphatic drainage, blood perfusion, osmotic pressure, and the drug binding to the cells and the extracellular matrix in the lymphatic uptake are systematically studied. Furthermore, the drug distribution and absorption in a multi-layered porous medium are investigated to illustrate the effect of heterogeneity of permeability, as the permeability varies over a wide range in different layers of the tissue (such as dermis, subcutaneous tissue, muscle). While the interstitial pressure plays an essential role in the mechanisms regulating the absorption of free monoclonal antibodies, the binding and metabolism of drug proteins also affect the drug absorption by reducing the total free monoclonal antibodies.

Transport and lymphatic uptake of monoclonal antibodies after subcutaneous injection.

The subcutaneous injection has emerged to become a feasible self-administration practice for biotherapeutics due to the patient comfort and cost-effectiveness. However, the available knowledge about transport and absorption of these agents after subcutaneous injection is limited. Here, a mathematical framework to study the subcutaneous drug delivery of mAbs from injection to lymphatic uptake is presented. A three-dimensional poroelastic model is exploited to find the biomechanical response of the tissue by taking into account tissue deformation during the injection. The results show that including tissue deformability noticeably changes tissue poromechanical response due to the significant dependence of interstitial pressure on the tissue deformation. Moreover, the importance of the amount of lymph fluid at the injection site and the injection rate on the drug uptake to lymphatic capillaries is highlighted. Finally, variability of lymphatic uptake due to uncertainty in parameters including tissue poromechanical and lymphatic absorption parameters is evaluated. It is found that interstitial pressure due to injection is the major contributing factor in short-term lymphatic absorption, while the amount of lymph fluid at the site of injection determines the long-term absorption of the drug. Finally, it is shown that the lymphatic uptake results are consistent with experimental data available in the literature.