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.     

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.     

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.     

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.     

MiR-27b promotes sheep skeletal muscle satellite cell proliferation by targeting <Emphasis Type="Italic">myostatin</Emphasis> gene

To investigate the role of miR-27b in sheep skeletal muscle development, here we first cloned the sequence of sheep pre-miR-27b, then further investigated its expression pattern in sheep skeletal muscle in vivo, the relationship of miR-27b expression and sheep skeletal muscle satellite cell proliferation and differentiation in vitro, and then finally confirmed its target gene during this development process. MiR-27b sequence, especially its mature sequence, was conservative among different species. MiR-27b highly expressed in sheep skeletal muscle than other tissues. In skeletal muscle of Suffolk and Bashbay sheep, miR-27b was upregulated during foetal period and downregulated during postnatal period significantly (\(P{<}0.01\)), but it still kept a relatively higher expression level in skeletal muscle of postnatal Suffolk sheep than Bashbay. There is a potential target site of miR-27b on \(3^\prime \)-UTR of sheep myostatin (MSTN) mRNA, and the double luciferase reporter assay proved that miR-27b could successfully bind on this site. When sheep satellite cells were in the proliferation status, miR-27b was upregulated and MSTN was downregulated significantly (\(P{<}0.01\)). When miR-27b mimics was transfected into sheep satellite cells, the cell proliferation was promoted and the protein level of MSTN was significantly downregulated (\(P{<}0.01\)). Moreover, miR-27b regulated its target gene MSTN by translation repression at an early step, and followed by inducing mRNA degradation in sheep satellite cells. Based on these results, we confirm that miR-27b could promote sheep skeletal muscle satellite cell proliferation by targeting MSTN and suppressing its expression.     

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.     

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.     

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.     

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.     

Caspase-1-Mediated Pyroptosis of the Predominance for Driving CD4$$^{}$$ T Cells Death: A Nonlocal Spatial Mathematical Model

Caspase-1-mediated pyroptosis is the predominance for driving CD4\(^{+}\) T cells death. Dying infected CD4\(^{+}\) T cells can release inflammatory signals which attract more uninfected CD4\(^{+}\) T cells to die. This paper is devoted to developing a diffusive mathematical model which can make useful contributions to understanding caspase-1-mediated pyroptosis by inflammatory cytokines IL-1\(\beta \) released from infected cells in the within-host environment. The well-posedness of solutions, basic reproduction number, threshold dynamics are investigated for spatially heterogeneous infection. Travelling wave solutions for spatially homogeneous infection are studied. Numerical computations reveal that the spatially heterogeneous infection can make \(\mathscr {R}_0>1\), that is, it can induce the persistence of virus compared to the spatially homogeneous infection. We also find that the random movements of virus have no effect on basic reproduction number for the spatially homogeneous model, while it may result in less infection risk for the spatially heterogeneous model, under some suitable parameters. Further, the death of infected CD4\(^{+}\) cells which are caused by pyroptosis can make \(\mathscr {R}_0<1\), that is, it can induce the extinction of virus, regardless of whether or not the parameters are spatially dependent.     

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.     

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.     

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.     

A Malaria Transmission Model with Temperature-Dependent Incubation Period

Malaria is an infectious disease caused by Plasmodium parasites and is transmitted among humans by female Anopheles mosquitoes. Climate factors have significant impact on both mosquito life cycle and parasite development. To consider the temperature sensitivity of the extrinsic incubation period (EIP) of malaria parasites, we formulate a delay differential equations model with a periodic time delay. We derive the basic reproduction ratio \(R_0\) and establish a threshold type result on the global dynamics in terms of \(R_0\), that is, the unique disease-free periodic solution is globally asymptotically stable if \(R_0<1\); and the model system admits a unique positive periodic solution which is globally asymptotically stable if \(R_0>1\). Numerically, we parameterize the model with data from Maputo Province, Mozambique, and simulate the long-term behavior of solutions. The simulation result is consistent with the obtained analytic result. In addition, we find that using the time-averaged EIP may underestimate the basic reproduction ratio.     

Ultrasound palpation for fast in-situ quantification of articular cartilage stiffness,thickness and relaxation capacity

Most current cartilage testing devices require the preparation of excised samples and therefore do not allow intra-operative application for diagnostic purposes. The gold standard during open or arthroscopic surgery is still the subjective perception of manual palpation. This work presents a new diagnostic method of ultrasound palpation (USP) to acquire applied stress and strain data during manual palpation of articular cartilage. With the proposed method, we obtain cartilage thickness and stiffness. Moreover, repeated palpations allow the quantification of relaxation effects. USP measurements on elastomer phantoms demonstrated very good repeatability for both, stage-guided (97.2%) and handheld (96.0%) applications. The USP measurements were compared with conventional indentation experiments and revealed very good agreement on elastomer phantoms (\(r = 0.98\)) and good agreement on porcine cartilage samples (\(r = 0.76\)). Artificially degenerated cartilage samples showed reduced stiffness, weak capacity to relax after palpation and an increase of stiffness of approximately 50% with each single palpation. Intact cartilage was measured by USP directly at the patella (in situ) and after excision and removal of the subchondral bone (ex situ), leading to stiffness values of \(12.1\pm 5.5\) and \(8.5\pm 5.9\,\hbox {MPa}\) (\(p<0.05\)), respectively. The results demonstrate the potential of the USP system for cartilage testing, its sensitivity to degenerative changes and as a method for quantifying relaxation processes by means of repeated palpations. Furthermore, the differences in the results of in-situ and ex-situ measurements are of general interest, since such comparison has not been reported previously. We point out the limited comparability of ex-situ cartilage with its in-situ biomechanical behavior.     

Contribution of left ventricular residual stress by myocytes and collagen: existence of inter-constituent mechanical interaction

We quantify the contribution of myocytes, collagen fibers and their interactions to the residual stress field found in the left ventricle (LV) using both experimental and theoretical methods. Ring tissue samples extracted from normal rat, male and female, LV were treated with collagenase and decellularization to isolate myocytes and collagen fibers, respectively. Opening angle tests were then performed on these samples as well as intact tissue samples containing both constituents that served as control. Our results show that the collagen fibers are the main contributor to the residual stress fields found in the LV. Specifically, opening angle measured in collagen-only samples (106.45\(^\circ \) ± 23.02\(^\circ \)) and myocytes-only samples (21.00\(^\circ \) ± 4.37\(^\circ \)) was significantly higher and lower than that of the control (57.88\(^\circ \) ± 12.29\(^\circ \)), respectively. A constrained mixture (CM) modeling framework was then used to infer these experimental results. We show that the framework cannot reproduce the opening angle found in the intact tissue with measurements made on the collagen-only and myocytes-only samples. Given that the CM framework assumes that each constituent contributes to the overall mechanics simply by their mere presence, this result suggests the existence of some myocyte–collagen mechanical interaction that cannot be ignored in the LV. We then propose an extended CM formulation that takes into account of the inter-constituent mechanical interaction in which constituents are deformed additionally when they are physically combined into a mixture. We show that the intact tissue opening angle can be recovered in this framework.     

A Mathematical Model Supports a Key Role for Ae4 (Slc4a9) in Salivary Gland Secretion

We develop a mathematical model of a salivary gland acinar cell with the objective of investigating the role of two \(\mathrm{Cl}^-/\mathrm{HCO}_3^-\) exchangers from the solute carrier family 4 (Slc4), Ae2 (Slc4a2) and Ae4 (Slc4a9), in fluid secretion. Water transport in this type of cell is predominantly driven by \(\mathrm{Cl}^-\) movement. Here, a basolateral \(\mathrm{Na}^+/ \mathrm{K}^+\) adenosine triphosphatase pump (NaK-ATPase) and a \(\mathrm{Na}^+\)–\(\mathrm{K}^+\)–\(2 \mathrm{Cl}^-\) cotransporter (Nkcc1) are primarily responsible for concentrating the intracellular space with \(\mathrm{Cl}^-\) well above its equilibrium potential. Gustatory and olfactory stimuli induce the release of \(\mathrm{Ca}^{2+}\) ions from the internal stores of acinar cells, which triggers saliva secretion. \(\mathrm{Ca}^{2+}\)-dependent \(\mathrm{Cl}^-\) and \(\mathrm{K}^+\) channels promote ion secretion into the luminal space, thus creating an osmotic gradient that promotes water movement in the secretory direction. The current model for saliva secretion proposes that \(\mathrm{Cl}^-/ \mathrm{HCO}_3^-\) anion exchangers (Ae), coupled with a basolateral \(\mathrm{Na}^+/\hbox {proton}\) (\(\hbox {H}^+\)) (Nhe1) antiporter, regulate intracellular pH and act as a secondary \(\mathrm{Cl}^-\) uptake mechanism (Nauntofte in Am J Physiol Gastrointest Liver Physiol 263(6):G823–G837, 1992; Melvin et al. in Annu Rev Physiol 67:445–469, 2005.  https://doi.org/10.1146/annurev.physiol.67.041703.084745). Recent studies demonstrated that Ae4 deficient mice exhibit an approximate \(30\%\) decrease in gland salivation (Peña-Münzenmayer et al. in J Biol Chem 290(17):10677–10688, 2015). Surprisingly, the same study revealed that absence of Ae2 does not impair salivation, as previously suggested. These results seem to indicate that the Ae4 may be responsible for the majority of the secondary \(\mathrm{Cl}^-\) uptake and thus a key mechanism for saliva secretion. Here, by using ‘in-silico’ Ae2 and Ae4 knockout simulations, we produced mathematical support for such controversial findings. Our results suggest that the exchanger’s cotransport of monovalent cations is likely to be important in establishing the osmotic gradient necessary for optimal transepithelial fluid movement.     

Competition of Multiple Species in Advective Environments

We study the effect of changes in flow speed on competition of an arbitrary number of species living in advective environments, such as streams and rivers. We begin with a spatial Lotka–Volterra model which is described by n reaction–diffusion–advection equations with Danckwerts boundary conditions. Using the dominant eigenvalue \(\lambda \le 0\) of the diffusion–advection operator subject to boundary conditions, we reduce the model to a system of ordinary differential equations. We impose a “transitive arrangement” of the competitors in terms of their interspecific coefficients and growth rates, which means that in the absence of advection, we have the following situation: for all \(1\le i<j\le n\), species i out-competes species j, while species j has higher intrinsic growth rate than species i. Changing advection speed in the original spatial model corresponds to changing the value of \(\lambda \) in the spatially implicit model. Considering the cases of the odd and even n separately, we obtain explicit intervals of the values of \(\lambda \) that allow all n species to be present in the habitat (coexistence interval). Stability of this equilibrium is shown for \(n\le 4\).