Analysis of two susceptibility SNPs in HLA region and evidence of interaction between rs6457617 in <Emphasis Type="BoldItalic">HLA-DQB1</Emphasis> and <Emphasis Type="BoldItalic">HLA-DRB1</Emphasis>*04 locus on Tunisian rheumatoid arthritis

Previous genomewide association studies (GWAS) and meta-analyses have enumerated several genes/loci in major histocompatibility complex region, which are consistently associated with rheumatoid arthritis (RA) in different ethnic populations. Given the genetic heterogeneity of the disease, it is necessary to replicate these susceptibility loci in other populations. In this case, we investigate the analysis of two SNPs, rs13192471 and rs6457617, from the human leukocyte antigen (HLA) region with the risk of RA in Tunisian population. These SNPs were previously identified to have a strong RA association signal in several GWAS studies. A case–control sample composed of 142 RA patients and 123 healthy controls was analysed. Genotyping of rs13192471 and rs6457617 was carried out using real-time PCR methods by TaqMan allelic discrimination assay. A trend of significant association was found in rs6457617 TT genotype with susceptibility to RA (\(P = 0.04\), \(p_{c} = 0.08\), \(\hbox {OR} = 1.73\)). Moreover, using multivariable analysis, the combination of rs6457617*TT–HLA-DRB1*\(04^{+}\) increased risk of RA (\(\hbox {OR} = 2.38\)), which suggest a gene–gene interaction event between rs6457617 located within the HLA-DQB1 and HLA-DRB1. Additionally, haplotypic analysis highlighted a significant association of rs6457617*T–HLA-DRB1*\(04^{+}\) haplotype with susceptibility to RA (\(P = 0.018\), \(p_{c} = 0.036\), \(\hbox {OR} = 1.72\)). An evidence of association was shown subsequently in \(\hbox {antiCCP}^{+}\) subgroup with rs6457617 both in T allele and TT genotype (\(P = 0.01\), \(p_{c} = 0.03\), \(\hbox {OR} = 1.66\) and \(P = 0.008\), \(p_{c} = 0.024\), \(\hbox {OR} = 1.28\), respectively). However, no association was shown for rs13192471 polymorphism with susceptibility and severity to RA. This study suggests the involvement of rs6457617 locus as risk variant for susceptibility/severity to RA in Tunisian population. Secondly, it highlights the gene–gene interaction between HLA-DQB1 and HLA-DRB1.     

Cell Volume Regulation in the Proximal Tubule of Rat Kidney

We developed a dynamic model of a rat proximal convoluted tubule cell in order to investigate cell volume regulation mechanisms in this nephron segment. We examined whether regulatory volume decrease (RVD), which follows exposure to a hyposmotic peritubular solution, can be achieved solely via stimulation of basolateral K\(^+\) and \(\hbox {Cl}^-\) channels and \(\hbox {Na}^+\)–\(\hbox {HCO}_3^-\) cotransporters. We also determined whether regulatory volume increase (RVI), which follows exposure to a hyperosmotic peritubular solution under certain conditions, may be accomplished by activating basolateral \(\hbox {Na}^+\)/H\(^+\) exchangers. Model predictions were in good agreement with experimental observations in mouse proximal tubule cells assuming that a 10% increase in cell volume induces a fourfold increase in the expression of basolateral K\(^+\) and \(\hbox {Cl}^-\) channels and \(\hbox {Na}^+\)–\(\hbox {HCO}_3^-\) cotransporters. Our results also suggest that in response to a hyposmotic challenge and subsequent cell swelling, \(\hbox {Na}^+\)–\(\hbox {HCO}^-_3\) cotransporters are more efficient than basolateral K\(^+\) and \(\hbox {Cl}^-\) channels at lowering intracellular osmolality and reducing cell volume. Moreover, both RVD and RVI are predicted to stabilize net transcellular \(\hbox {Na}^+\) reabsorption, that is, to limit the net \(\hbox {Na}^+\) flux decrease during a hyposmotic challenge or the net \(\hbox {Na}^+\) flux increase during a hyperosmotic challenge.     

Genetic variability and structure of an isolated population of <Emphasis Type="Italic">Ambystoma altamirani</Emphasis>, a mole salamander that lives in the mountains of one of the largest urban areas in the world

Amphibians are globally threatened by habitat loss and fragmentation; species within the order Ambystoma are not the exception, as there are 18 species of mole salamanders in México, of which 16 are endemic and all species are under some national or international status of protection. The mole salamander, Ambystoma altamirani is a microendemic species, which is distributed in central México, within the trans-Mexican volcanic belt, and is one of the most threatened species due to habitat destruction and the introduction of exotic species. Nine microsatellite markers were used to determine the genetic structure, genetic variability, effective population size, presence of bottlenecks and inbreeding coefficient of one population of A. altamirani to generate information which might help to protect and conserve this threatened species. We found two genetic subpopulations with significant level of genetic structure (\(F_{\mathrm{ST}}= 0.005\)) and high levels of genetic variability (\(H_{\mathrm{o}}= 0.883\); \(H_{\mathrm{e}}= 0.621\)); we also found a small population size (\(N_{\mathrm{e}} = 8.8\)), the presence of historical (\(M =\) 0.486) and recent bottlenecks under IAM and TPM models, with a low, but significant coefficient of inbreeding (\(F_{\mathrm{IS}} = -\)0.451). This information will help us to raise conservation strategies of this microendemic mole salamander species.     

Oestrogen regulates the expression of cathepsin E-A-like gene through ER$$\upbeta $$ in liver of chicken ( Gallus gallus)

The cathepsin E-A-like, also known as ‘similar to nothepsin’, is a new member of the aspartic protease family, which may take part in processing of egg yolk macromolecules, due to it was identified in the chicken egg-yolk. Previously, studies have suggested that the expression of cathepsin E-A-like increased gradually during sexual maturation of pullets, but the exact regulation mechanism is poorly understood. In this study, to gain insight into the function and regulation mechanism of the gene in egg-laying hen, we cloned the cathepsin E-A-like gene and evaluated its evolutionary origin by using both phylogenetic and syntenic methods. The mode of the gene expression regulation was analysed through stimulating juvenile hens with \(17\upbeta \)-estradiol and chicken embryo hepatocytes with \(17\upbeta \)-estradiol combined with oestrogen receptor antagonists including MPP, ICI 182,780 and tamoxifen. Our results showed that cathepsin E-A-like was an orthologoues gene with nothepsin, which is present in birds but not in mammals. The expression of cathepsin E-A-like significantly increased in a dose-dependent manner after the juvenile hens were treated with \(17\upbeta \)-estradiol (\(P~<~0.05\)). Compared with the \(17\upbeta \)-estradiol treatment group, the expression of cathepsin E-A-like was not significantly changed when the hepatocytes were treated with \(17\upbeta \)-estradiol combined with MPP (\(P~<~0.05\)). In contrast, compared with the \(17\upbeta \)-estradiol combined with MPP treatment group, the expression of cathepsin E-A-like was significantly downregulated when the hepatocytes were treated with \(17\upbeta \)-estradiol combined with tamoxifen or ICI 182,780 (\(P~<~0.05\)). These results demonstrated that cathepsin E-A-like shared the same evolutionary origin with nothepsin. The expression of cathepsin E-A-like was regulated by oestrogen, and the regulative effect was predominantly mediated through ER-\(\upbeta \) in liver of chicken.     

Glutathione <Emphasis Type="BoldItalic">S</Emphasis>-transferase P1 gene polymorphisms and susceptibility to coronary artery disease in a subgroup of north Indian population

The present study aimed to investigate the association of \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) and \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) polymorphisms of GSTP1 with coronary artery disease (CAD) in a subgroup of north Indian population. In the present case–control study, CAD patients (\(n = 200\)) and age-matched, sex-matched and ethnicity-matched healthy controls (\(n = 200\)) were genotyped for polymorphisms in GSTP1 using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method. Genotype distribution of \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) and \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) polymorphisms of GSTP1 gene was significantly different between cases and controls (\(P = 0.005\) and 0.024, respectively). Binary logistic regression analysis showed significant association of A/G (odds ratio (OR): 1.6, 95% CI: 1.08–2.49, \(P = 0.020\)) and G/G (OR: 3.1, 95% CI: 1.41–6.71, P \(=\) 0.005) genotypes of GSTP1 \(\hbox {g}.313\hbox {A}{\!>\!}\hbox {G}\), and C/T (OR: 5.8, 95% CI: 1.26–26.34, \(P = 0.024\)) genotype of GSTP1 \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) with CAD. The A/G and G/G genotypes of \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) and C/T genotype of \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) conferred 6.5-fold increased risk for CAD (OR: 6.5, 95% CI: 1.37–31.27, \(P = 0.018\)). Moreover, the recessive model of GSTP1 \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) is the best fit inheritance model to predict the susceptible gene effect (OR: 2.3, 95% CI: 1.11–4.92, \(P = 0.020\)). In conclusion, statistically significant associations of GSTP1 \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) (A/G, G/G) and \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) (C/T) genotypes with CAD were observed.     

Trend analysis of variations in carbon stock using stock big data

Changes in land use affect the terrestrial carbon stock through changes in the land cover. Research on land use and analysis of variations in carbon stock have practical applications in the optimization of land use and the mitigation of climate change effects. This study was conducted in Baixiang and Julu counties in the Taihang Piedmont by employing the trend analysis method to characterize the variation in county land use and carbon stock. The findings show that in both counties, agricultural and unused land areas decreased while built-up land area increased, and the reduction in cropland was the main reason behind the agricultural land reduction. An inflection point appeared on the cropland curves of Julu, because the cropland area decreased by 1576.97 hm\(^{2}\) from 2004 to 2006. Cropland area in Baixiang decreased from 1996 to 1998 by a total of 129.89 hm\(^{2}\) and then remained relatively stable after 1998. The total carbon storage and variation in land use in the two counties displayed similar trends. Total carbon reserves in Julu increased by 2.76 \(\times \) 10\(^{4}\) tC (carbon equivalent), while those in Baixiang decreased by 0.63 \(\times \) 10\(^{4}\) tC. Carbon stock of built-up land in Julu and Baixiang increased by 2.44 \(\times \) 10\(^{4}\) and 1.22 \(\times \) 10\(^{4}\) tC, respectively.     

Joint non-uniform sampling of all incremented time delays for quicker acquisition in protein relaxation studies

NMR relaxometry plays crucial role in studies of protein dynamics. The measurement of longitudinal and transverse relaxation rates of \(^{15}\)N is the main source of information on backbone motions. However, even the most basic approach exploiting a series of \(^{15}\)N HSQC spectra can require several hours of measurement time. Standard non-uniform sampling (NUS), i.e. random under-sampling of indirect time domain, typically cannot reduce this by more than 2–4\(\times\) due to relatively low “compressibility” of these spectra. In this paper we propose an extension of NUS to relaxation delays. The two-dimensional space of \(t_1\)/\(t_{relax}\) is sampled in a way similar to NUS of \(t_1\)/\(t_2\) domain in 3D spectra. The signal is also processed in a way similar to that known from 3D NUS spectra i.e. using one of the most popular compressed sensing algorithms, iterative soft thresholding. The 2D Fourier transform matrix is replaced with mixed inverse Laplace-Fourier transform matrix. The peak positions in resulting 3D spectrum are characterized by two frequency coordinates and relaxation rate and thus no additional fitting of exponential curves is required. The method is tested on three globular proteins, providing satisfactory results in a time corresponding to acquisition of two conventional \(^{15}\)N HSQC spectra.     

Improved identifiability of myocardial material parameters by an energy-based cost function

Myocardial stiffness is a valuable clinical biomarker for the monitoring and stratification of heart failure (HF). Cardiac finite element models provide a biomechanical framework for the assessment of stiffness through the determination of the myocardial constitutive model parameters. The reported parameter intercorrelations in popular constitutive relations, however, obstruct the unique estimation of material parameters and limit the reliable translation of this stiffness metric to clinical practice. Focusing on the role of the cost function (CF) in parameter identifiability, we investigate the performance of a set of geometric indices (based on displacements, strains, cavity volume, wall thickness and apicobasal dimension of the ventricle) and a novel CF derived from energy conservation. Our results, with a commonly used transversely isotropic material model (proposed by Guccione et al.), demonstrate that a single geometry-based CF is unable to uniquely constrain the parameter space. The energy-based CF, conversely, isolates one of the parameters and in conjunction with one of the geometric metrics provides a unique estimation of the parameter set. This gives rise to a new methodology for estimating myocardial material parameters based on the combination of deformation and energetics analysis. The accuracy of the pipeline is demonstrated in silico, and its robustness in vivo, in a total of 8 clinical data sets (7 HF and one control). The mean identified parameters of the Guccione material law were \(C_1=3000\pm 1700\,\hbox {Pa}\) and \(\alpha =45\pm 25\) (\(b_f=25\pm 14\), \(b_{ft}=11\pm 6\), \(b_{t}=9\pm 5\)) for the HF cases and \(C_1=1700\,\hbox {Pa}\) and \(\alpha =15\) (\(b_f=8\), \(b_{ft}=4\), \(b_{t}=3\)) for the healthy case.     

Backbone resonance assignments for the SET domain of the human methyltransferase NSD2

Aberrant NSD2 methyltransferase activity is implicated as the oncogenic driver in multiple myeloma, suggesting opportunities for novel therapeutic intervention. The methyltransferase activity of NSD2 resides in its catalytic SET domain, which is conserved among most lysine methyltransferases. Here we report the backbone \(\hbox {H}^{\mathrm{N}}\), N, C\(^{\prime }\), \(\hbox {C}^\alpha\) and side-chain \(\hbox {C}^\beta\) assignments of a 25 kDa NSD2 SET domain construct, spanning residues 991–1203. A chemical shift analysis of C\(^{\prime }\), \(\hbox {C}^\alpha\) and \(\hbox {C}^\beta\) resonances predicts a secondary structural pattern that is in agreement with homology models.     

Reduced Winter Snowpack and Greater Soil Frost Reduce Live Root Biomass and Stimulate Radial Growth and Stem Respiration of Red Maple (<Emphasis Type="Italic">Acer rubrum</Emphasis>) Trees in a Mixed-Hardwood Forest

Northeastern U.S. forests are currently net carbon (C) sinks, but rates of C loss from these ecosystems may be altered by the projected reduction in snowpack and increased soil freezing over the next century. Soil freezing damages fine roots, which may reduce radial tree growth and stem respiration. We conducted a snow removal experiment at Harvard Forest, MA to quantify effects of a reduced winter snowpack and increased soil freezing on root biomass, stem radial growth and respiration in a mixed-hardwood forest. The proportion of live fine root biomass during spring (late-April) declined with increasing soil frost severity (P = 0.05). Basal area increment index was positively correlated with soil frost severity for Acer rubrum, but not Quercus rubra. Rates of stem respiration in the growing season correlated positively with soil frost duration in the previous winter, (\( R^{2}_{{{\text{LMM}}({\text{m}})}} \) = 0.15 and 0.24 for Q. rubra and A. rubrum, respectively). Losses of C from stem respiration were comparable to or greater than C storage from radial growth of Q. rubra and A. rubrum, respectively. Overall, our findings suggest that in mixed-hardwood forests (1) soil freezing has adverse effects on spring live root biomass, but at least in the short-term could stimulate aboveground processes such as stem respiration and radial growth for A. rubrum more than Q. rubra, (2) stem respiration is an important ecosystem C flux and (3) the increasing abundance of A. rubrum relative to Q. rubra may have important implications for C storage in tree stem biomass.     

Association between SSR markers and fibre traits in sea island cotton (<Emphasis Type="BoldItalic">Gossypium barbadense</Emphasis>) germplasm resources

Identification of molecular markers associated with fibre traits can accelerate cotton marker-assisted selection (MAS) programmes. In this study, Gossypium barbadense germplasm accessions with diverse origins (\(n = 123\)) were used to perform association analysis of fibre traits with 120 polymorphic simple sequence repeat (SSR) markers. In total, 120 polymorphic primer pairs amplified 258 loci with a mean of 2.15 loci per primer. Population structure analysis identified three main clusters for the accessions, which indicated agreement of genetic and predefined populations. Marker–trait associations (\(n= 58\)) were detected for 10 fibre traits with 26 SSR markers located on 15 chromosomes. The \(R^{2}\) (phenotypic variation explained) ranged from 3.19 to 15.21%. Two markers (NAU5465 and NAU3013) were found to be stably associated with boll number per plant (BNP) and fibre uniformity (UI), respectively. Four markers (BNL252, NAU3424, NAU3324 and CGR5202) associated with fibre quality traits preferentially clustered on the D8 chromosome, which was thus identified as an important candidate region for study molecular mechanisms underlying fibre quality and for use in breeding cotton cultivars for improving fibre quality. This study generated molecular data with a potential for better understanding of the genetic basis of the fibre traits and provided new markers for MAS in G. barbadense breeding programmes.     

Mathematical modelling of bone adaptation of the metacarpal subchondral bone in racehorses

In Thoroughbred racehorses, fractures of the distal limb are commonly catastrophic. Most of these fractures occur due to the accumulation of fatigue damage from repetitive loading, as evidenced by microdamage at the predilection sites for fracture. Adaptation of the bone in response to training loads is important for fatigue resistance. In order to better understand the mechanism of subchondral bone adaptation to its loading environment, we utilised a square root function defining the relationship between bone volume fraction \((f_{BM} )\) and specific surface \((S_v )\) of the subchondral bone of the lateral condyles of the third metacarpal bone (MCIII) of the racehorse, and using this equation, developed a mathematical model of subchondral bone that adapts to loading conditions observed in vivo. The model is expressed as an ordinary differential equation incorporating a formation rate that is dependent on strain energy density. The loading conditions applied to a selected subchondral region, i.e. volume of interest, were estimated based on joint contact forces sustained by racehorses in training. For each of the initial conditions of \(f_{BM} \) we found no difference between subsequent homoeostatic \(f_{BM} \) at any given loading condition, but the time to reach equilibrium differed by initial \(f_{BM} \) and loading condition. We found that the observed values for \(f_{BM} \) from the mathematical model output were a good approximation to the existing data for racehorses in training or at rest. This model provides the basis for understanding the effect of changes to training strategies that may reduce the risk of racehorse injury.     

Prediction of heterosis in rice based on divergence of morphological and molecular markers

Identifying the best performing hybrid without a field test was essential to save resources and time. In this study, the genetic divergence was estimated using morphological and expressed sequence tag (EST)-derived simple sequence repeats (SSR) markers. Cluster analysis showed that APMS6A and RPHR 1005 belong to groups I and II, respectively, and the hybrid combination recorded the highest mean grain yield of 32.25 g among generated 40 \(\hbox {F}_{1}\hbox {s}\) with standard heterosis of 8.4% over hybrid check, KRH2. The coefficient of marker polymorphism (CMP) value was calculated based on EST-SSR markers; it ranged from 0.40 to 0.80, and a higher CMP value of 0.80 was obtained for the parental combination APMS6A \(\times \) RPHR1005. We predicted heterosis for 40 \(\hbox {F}_{1}\hbox {s}\) based on correlation between CMP and standard heterosis in different traits with standard varietal and hybrid checks indicating positive correlation and significant value for grain yield per plant (\(r=0.58\)**), productivity per day (\(r=0.54\)**), productive tillers (\(r=0.34\)*) and panicle weight (\(r=0.42\)**). This study revealed that the relationship of molecular marker heterozygosity, along with the combining ability, high mean value of different traits, grouping of parental lines based on morphological and molecular characterization is helpful to identify heterotic patterns in rice.     

Association of <Emphasis Type="BoldItalic">lactase</Emphasis> 13910 C/T polymorphism with bone mineral density and fracture risk: a meta-analysis

A number of studies have investigated the association of lactase (LCT, C/T-13910) gene polymorphism with bone mineral density (BMD) and fracture risk, but previous results were inconclusive. In this study, a meta-analysis was performed to quantify the association of LCT (C/T-13910) polymorphism with BMD and fracture risk. Eligible publications were searched in the PubMed, Web of Science, Embase databases, Google Scholar, Yahoo and Baidu. Pooled weighed mean difference (WMD) or odds ratio (OR) with their 95% confidence interval (CI) were calculated using a fixed-effects or random-effects model. A total of nine articles with 8871 subjects were investigated in the present meta-analysis. Overall, the TT/TC genotypes of LCT 13910 C/T polymorphism showed significantly higher BMD than those with the CC genotype at femur neck (FN) (\(\hbox {WMD} = 0.011\,\hbox {g/cm}^{2}\), 95% CI \(=\) 0.004–0.018, \(P = 0.003\)). Besides, LCT 13910 C/T polymorphism may decrease the risk of any site fractures (for TT versus TC \(+\) CC, OR \(=\) 0.813, 95% CI \(=\) 0.704–0.938, \(P = 0.005\); for T allele versus C allele, OR \(=\) 0.885, 95% CI \(=\) 0.792–0.989, \(P = 0.032\)). However, there was no significant association of LCT 13910 C/T polymorphism with BMD at lumbar spine and risk of vertebral fractures under all genetic contrast models (all P values were \({>}0.05\)). The meta-analysis suggests that there are significant effects of LCT 13910 C/T polymorphism on BMD and fracture risk. Large-scale studies with different ethnic populations will be needed to further investigate the possible race-specific effect of LCT 13910 C/T polymorphism on BMD and fracture risk.     

Fine root dynamics and soil carbon accretion under thinned and un-thinned Cupressus lusitanica stands in, Southern Ethiopia

Background and aims

Forest management activities influences stand nutrient budgets, belowground carbon allocation and storage in the soil. A field experiment was carried out in Southern Ethiopia to investigate the effect of thinning on fine root dynamics and associated soil carbon accretion of 6-year old C. lusitanica stands.

Methods

Fine roots (≤2 mm in diameter) were sampled seasonally to a depth of 40 cm using sequential root coring method. Fine root biomass and necromass, vertical distribution, seasonal dynamics, annual turnover and soil carbon accretion were quantified.

Results

Fine root biomass and necromass showed vertical and temporal variations. More than 70 % of the fine root mass was concentrated in the top 20 cm soil depth. Fine root biomass showed significant seasonal variation with peaks at the end of the major rainy season and short rainy season. Thinning significantly increased fine root necromass, annual fine root production and turnover. Mean annual soil carbon accretion, through fine root necromass, in the thinned stand was 63 % higher than that in the un-thinned stand.

Conclusions

The temporal dynamics in fine roots is driven by the seasonality in precipitation. Thinning of C. lusitanica plantation would increase soil C accretion considerably through increased fine root necromass and turnover.     

Simulation of extracellular matrix remodeling by fibroblast cells in soft three-dimensional bioresorbable scaffolds

To culture functional soft tissues and organs in three-dimensional scaffolds, it is essential to elucidate the optimal scaffold mechanical properties. However, mechanoregulated soft tissue remodeling is not well understood. In this study, we hypothesized that individual cells are capable of remodeling extracellular matrix within a short proximity of themselves in order to match the stiffness of the broader surrounding matrix. This theory was implemented in a three-dimensional finite element model to simulate soft tissue remodeling of human fibroblast cells in two collagen–chitosan scaffolds with different mechanical properties. Simulation results closely matched with previously reported experimental data, showing that soft tissue growth in compliant (Scaf-A, 4.30 kPa) and stiff (Scaf-B, 17.03 kPa) scaffolds led to an almost eightfold difference in the resulting stiffnesses after 10 days (8.40 kPa for Scaf-A, 59.25 kPa for Scaf-B). Furthermore, varying the simulated rate for tissue remodeling by \(\pm \)50 % caused unequal changes in the resulting stiffness (+3.6 and \(-\)23 % for Scaf-A, +5 and \(-\)17 % for Scaf-B), and \(\pm \)50 % changes in the assumed upper limit on tissue stiffness only had larger effects on the stiff scaffold (+42 and \(-\)44 % for Scaf-B). These results reinforce the notion that soft tissue remodeling is governed by the stiffness of the surrounding matrix, until meeting an upper limit on tissue stiffness.     

Towards the prediction of flow-induced shear stress distributions experienced by breast cancer cells in the lymphatics

Tumour metastasis in the lymphatics is a crucial step in the progression of breast cancer. The dynamics by which breast cancer cells (BCCs) travel in the lymphatics remains poorly understood. The goal of this work is to develop a model capable of predicting the shear stresses metastasising BCCs experience using numerical and experimental techniques. This paper models the fluidic transport of large particles (\(\eta =d_{\mathrm{p}}/W=0.1-0.4\) where \(d_{\mathrm{p}}\) is the particle diameter and W is the channel width) subjected to lymphatic flow conditions (\({ Re}=0.04\)), in a \(100\times 100\,\upmu \hbox {m}\) microchannel. The feasibility of using the dynamic fluid body interaction (DFBI) method to predict particle motion was assessed, and particle tracking experiments were performed. The experiments found that particle translational velocity decreased from the undisturbed fluid velocity with increasing particle size (5–14% velocity lag for \(\eta =0.1-0.3\)). DFBI simulations were found to better predict particle behaviour than theoretical predictions; however, mesh restrictions in the near-wall region (\(0.2\,\mathrm{W}>y>0.8\,\mathrm{W}\)) result in computationally expensive models. The simulations were in good agreement with the experiments (\(<12\%\) difference) across the channel (\(0.2\,\mathrm{W}\le y\le 0.8\,\mathrm{W}\)), with differences up to 25% in the near-wall region. Particles experience a range of shear stresses (0.002–0.12 Pa) and spatial shear gradients (\(0.004-0.137\,\hbox {Pa}/\upmu \hbox {m}\)) depending on their size and radial position. The predicted shear gradients are far in excess of values associated with BCC apoptosis (\(0.004-0.023\,\hbox {Pa}/\upmu \hbox {m}\)). Increasing our understanding of the shear stress magnitudes and gradients experienced by BCCs could be leveraged to elucidate whether a particular BCC size or location exists that encourages metastasis within the lymphatics.     

The Role of <Emphasis Type="BoldItalic">Hyalomma Truncatum</Emphasis> on the Dynamics of Rift Valley Fever: Insights from a Mathematical Epidemic Model

To date, our knowledge of Rift Valley fever (RVF) disease spread and maintenance is still limited, as flooding, humid weather and presence of biting insects such as mosquitoes, have not completely explained RVF outbreaks. We propose a model that includes livestock, mosquitoes and ticks compartments structured according to their questing and feeding behaviour in order to study the possible role of ticks on the dynamics of RVF. To quantify disease transmission at the initial stage of the epidemic, we derive an explicit formula of the basic reproductive number, \(R_0\). Using the concept of Metzler matrix, we state necessary conditions for global asymptotic stability of the disease-free equilibrium. Results suggest that although host-ticks interactions may serve as disease reservoirs or disease amplifiers, the Aedes reproductive number should be kept under unity if disease post-epizootics activities are to be controlled. Results of both local and global sensitivity analysis of selected model parameters indicate that \(R_0\) is more sensitive to the ticks attachment and detachment rates, probability of transmission from ticks to host and from host to ticks, length of infection in livestock and ticks death rate. Furthermore, when comparing the mean value of \(R_0\) with that from previous studies which did not include ticks we found that our \(R_0\) is very much larger resulting in an increase in the exponential phase of an outbreak. These findings suggest that if ticks are capable of transmitting the virus, they may be contributing to disease outbreaks and endemicity.     

Critical zone properties control the fate of nitrogen during experimental rainfall in montane forests of the Colorado Front Range

Several decades of research in alpine ecosystems have demonstrated links among the critical zone, hydrologic response, and the fate of elevated atmospheric nitrogen (N) deposition. Less research has occurred in mid-elevation forests, which may be important for retaining atmospheric N deposition. To explore the fate of N in the montane zone, we conducted plot-scale experimental rainfall events across a north–south transect within a catchment of the Boulder Creek Critical Zone Observatory. Rainfall events mimicked relatively common storms (20–50% annual exceedance probability) and were labeled with ¹⁵N-nitrate (\( {\text{NO}}_{3}^{ - } \)) and lithium bromide tracers. For 4 weeks, we measured soil–water and leachate concentrations of Br^?, \( {}^{15}{\text{NO}}_{3}^{ - } , \) and \( {\text{NO}}_{3}^{ - } \) daily, followed by recoveries of ¹⁵N species in bulk soils and microbial biomass. Tracers moved immediately into the subsurface of north-facing slope plots, exhibiting breakthrough at 10 and 30 cm over 22 days. Conversely, little transport of Br^? or \( {}^{15}{\text{NO}}_{3}^{ - } \) occurred in south-facing slope plots; tracers remained in soil or were lost via pathways not measured. Hillslope position was a significant determinant of soil ¹⁵N-\( {\text{NO}}_{3}^{ - } \) recoveries, while soil depth and time were significant determinants of ¹⁵N recovery in microbial biomass. Overall, ¹⁵N recovery in microbial biomass and leachate was greater in upper north-facing slope plots than lower north-facing (toeslope) and both south-facing slope plots in August; by October, ¹⁵N recovery in microbial N biomass within south-facing slope plots had increased substantially. Our results point to the importance of soil properties in controlling the fate of N in mid-elevation forests during the summer season.