Study of biological communities subject to imperfect detection: bias and precision of community <Emphasis Type="Italic">N</Emphasis>-mixture abundance models in small-sample situations

Community N-mixture abundance models for replicated counts provide a powerful and novel framework for drawing inferences related to species abundance within communities subject to imperfect detection. To assess the performance of these models, and to compare them to related community occupancy models in situations with marginal information, we used simulation to examine the effects of mean abundance \((\bar{\lambda }\): 0.1, 0.5, 1, 5), detection probability \((\bar{p}\): 0.1, 0.2, 0.5), and number of sampling sites (n _site : 10, 20, 40) and visits (n _visit : 2, 3, 4) on the bias and precision of species-level parameters (mean abundance and covariate effect) and a community-level parameter (species richness). Bias and imprecision of estimates decreased when any of the four variables \((\bar{\lambda }\), \(\bar{p}\), n _site , n _visit ) increased. Detection probability \(\bar{p}\) was most important for the estimates of mean abundance, while \(\bar{\lambda }\) was most influential for covariate effect and species richness estimates. For all parameters, increasing n _site was more beneficial than increasing n _visit . Minimal conditions for obtaining adequate performance of community abundance models were n _site  ≥ 20, \(\bar{p}\) ≥ 0.2, and \(\bar{\lambda }\) ≥ 0.5. At lower abundance, the performance of community abundance and community occupancy models as species richness estimators were comparable. We then used additive partitioning analysis to reveal that raw species counts can overestimate β diversity both of species richness and the Shannon index, while community abundance models yielded better estimates. Community N-mixture abundance models thus have great potential for use with community ecology or conservation applications provided that replicated counts are available.     

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.     

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.     

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\).     

Genomewide identification of PPR gene family and prediction analysis on restorer gene in <Emphasis Type="BoldItalic">Gossypium</Emphasis>

Pentatricopeptide repeat (PPR) gene family plays an essential role in the regulation of plant growth and organelle gene expression. Some PPR genes are related to fertility restoration in plant, but there is no detailed information in Gossypium. In the present study, we identified 482 and 433 PPR homologues in Gossypium raimondii (\(\hbox {D}_{5}\)) and G. arboreum (\(\hbox {A}_{2}\)) genomes, respectively. Most PPR homologues showed an even distribution on the whole chromosomes. Given an evolutionary analysis to PPR genes from G. raimondii (\(\hbox {D}_{5}\)), G. arboreum (\(\hbox {A}_{2}\)) and G. hirsutum genomes, eight PPR genes were clustered together with restoring genes of other species. Most cotton PPR genes were qualified with no intron, high proportion of \(\upalpha \)-helix and classical tertiary structure of PPR protein. Based on bioinformatics analyses, eight PPR genes were targeted in mitochondrion, encoding typical P subfamily protein with protein binding activity and organelle RNA metabolism in function. Further verified by RNA-seq and quantitative real-time PCR (qRT-PCR) analyses, two PPR candidate genes, Gorai.005G0470 (\(\hbox {D}_{5}\)) and Cotton_A_08373 (\(\hbox {A}_{2}\)), were upregulated in fertile line than sterile line. These results reveal new insights into PPR gene evolution in Gossypium.     

Similarity of decay-associated spectra for tryptophan fluorescence of proteins with different structures

Tryptophan fluorescence lifetimes were analyzed for three proteins: human serum albumin, bovine serum albumin, and bacterial luciferase, which contain one, two, and seven tryptophan residues, respectively. For all of the proteins, the fluorescence decays were fitted by three lifetimes: τ₁ = 6–7 ns, τ₂ = 2.0–2.3 ns, and τ₃ ≤ 0.1 ns (the native state), and τ₁ = 4.4–4.6 ns, τ₂ = 1.7–1.8 ns, and τ₃ ≤ 0.1 ns (the denatured state). Corresponding decay-associated spectra had similar peak wavelengths and spectrum half-widths both in the native state (\(\lambda _{\max }^{{\tau _1}} = 324nm\), \(\lambda _{\max }^{{\tau _2}} = 328nm\), and \(\lambda _{\max }^{{\tau _3}} = 315nm\)), and in the denatured state (\(\lambda _{\max }^{{\tau _1}} = 350nm\), \(\lambda _{\max }^{{\tau _2}} = 343nm\), and \(\lambda _{\max }^{{\tau _3}} = 317nm\)). The differences in the steady-state spectra of the studied proteins were accounted for the individual ratio of the lifetime component contributions. The lifetime components were compared with a classification of tryptophan residues in the structure of these proteins within the discrete states model.     

A novel real-time driving fatigue detection system based on wireless dry EEG

Development of techniques for detection of mental fatigue has varied applications in areas where sustaining attention is of critical importance like security and transportation. The objective of this study is to develop a novel real-time driving fatigue detection methodology based on dry Electroencephalographic (EEG) signals. The study has employed two methods in the online detection of mental fatigue: power spectrum density (PSD) and sample entropy (SE). The wavelet packets transform (WPT) method was utilized to obtain the \(\theta \) (4–7 Hz), \(\alpha \) (8–12 Hz) and \(\beta \) (13–30 Hz) bands frequency components for calculating corresponding PSD of the selected channels. In order to improve the fatigue detection performance, the system was individually calibrated for each subject in terms of fatigue-sensitive channels selection. Two fatigue-related indexes: (\(\theta +\alpha \))/\(\beta \) and \(\theta \)/\(\beta \) were computed and then fused into an integrated metric to predict the degree of driving fatigue. In the case of SE extraction, the mean of SE averaged across two EEG channels (‘O1h’ and ‘O2h’) was used for fatigue detection. Ten healthy subjects participated in our study and each of them performed two sessions of simulated driving. In each session, subjects were required to drive simulated car for 90 min without any break. The results demonstrate that our proposed methods are effective for fatigue detection. The prediction of fatigue is consistent with the observation of reaction time that was recorded during simulated driving, which is considered as an objective behavioral measure.     

Complexity and algorithms for copy-number evolution problems

Background

Cancer is an evolutionary process characterized by the accumulation of somatic mutations in a population of cells that form a tumor. One frequent type of mutations is copy number aberrations, which alter the number of copies of genomic regions. The number of copies of each position along a chromosome constitutes the chromosome’s copy-number profile. Understanding how such profiles evolve in cancer can assist in both diagnosis and prognosis.

Results

We model the evolution of a tumor by segmental deletions and amplifications, and gauge distance from profile \(\mathbf {a}\) to \(\mathbf {b}\) by the minimum number of events needed to transform \(\mathbf {a}\) into \(\mathbf {b}\). Given two profiles, our first problem aims to find a parental profile that minimizes the sum of distances to its children. Given k profiles, the second, more general problem, seeks a phylogenetic tree, whose k leaves are labeled by the k given profiles and whose internal vertices are labeled by ancestral profiles such that the sum of edge distances is minimum.

Conclusions

For the former problem we give a pseudo-polynomial dynamic programming algorithm that is linear in the profile length, and an integer linear program formulation. For the latter problem we show it is NP-hard and give an integer linear program formulation that scales to practical problem instance sizes. We assess the efficiency and quality of our algorithms on simulated instances.

Availability

https://github.com/raphael-group/CNT-ILP

     

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.     

Isolation and characterization of microsatellite markers in the spiny lobster, <Emphasis Type="Italic">Panulirus echinatus</Emphasis> Smith, 1869 (Decapoda: Palinuridae) by Illumina MiSeq sequencing

Okra’s (Abelmoschus esculentus (L.) Moench) commercial cultivation is threatened in the tropics due to high incidence of yellow vein mosaic virus (YVMV) disease. Okra geneticists across the world tried to understand the inheritance pattern of YVMV disease tolerance without much success. Therefore, the inheritance pattern of YVMV disease in okra was revisited by employing six generations (\(\hbox {P}_{1}\), \(\hbox {P}_{2}\), \(\hbox {F}_{1}\), \(\hbox {F}_{2}\), \(\hbox {BC}_{1}\) and \(\hbox {BC}_{2}\)) of four selected crosses (one tolerant \(\times \) tolerant, two tolerant \(\times \) susceptible and one susceptible \(\times \) susceptible) using two tolerant (BCO-1 and Lal Bhendi) and two susceptible (Japanese Jhar Bhendi and PAN 2127) genotypes. Qualitative genetic analysis was done on the basis of segregation pattern of tolerant and susceptible plants in \(\hbox {F}_{2}\) and backcross generations of all the four crosses. It revealed that a single dominant gene along with some minor factors governed the disease tolerant trait in both the tolerant parents used. However, it was observed that genes governing disease tolerance identified in both the tolerant variety used was different. It could be concluded that the gene governing YVMV disease tolerance in okra was genotype specific. Further, duplicate gene action as evident from an approximate ratio of 15 : 1 (tolerant : susceptible) in the \(\hbox {F}_{2}\) population in the cross of two tolerant varieties gave a scope of increasing the tolerance level of the hybrid plants when both the tolerant genes are brought together. However, generation mean analysis revealed involvement of both additive and nonadditive effects in the inheritance of disease tolerance. Thus, the present study confirms that a complicated genetic inheritance pattern is involved in the disease tolerance against YVMV trait. The major tolerance genes could be transferred to other okra varieties, but the tolerance breaking virus strains might not allow them to achieve tolerance in stable condition. Therefore, accumulation of additional genes may be needed for a sustainable tolerance phenotype in okra.     

Deletion / insertion polymorphisms of the prion protein gene (<Emphasis Type="BoldItalic">PRNP</Emphasis>) in gayal (<Emphasis Type="BoldItalic">Bos frontalis</Emphasis>)

Resistance to fatal disease bovine spongiform encephalopathy (BSE), due to misfolded prion protein in cattle, is associated with a 23-bp indel polymorphism in the putative promoter and a 12-bp indel in intron 1 of the PRNP gene. Gayal (Bos frontalis) is an important semiwild bovid species and of great conservation concern, but till today these indel polymorphisms have not been evaluated in gayals. Therefore, we collected 225 samples of gayals and evaluated the genetic indel polymorphism in the two regions of this PRNP gene. The results revealed high allelic frequencies of insertions at these indel sites: 0.909 and 0.667 for, respectively, the 23 bp and 12 bp indels, both also with significant genotype frequencies (\(\chi ^{2}\): 9.81; 23 bp and \(\chi ^{2}\): 43.56; 12 bp). At the same time, the haplotype data showed indel polymorphisms with extremely low deletion (0.01) in both regions of the PRNP gene. We compared these data with those reported for healthy and BSE-affected cattle (Bos taurus) breeds from two European countries, Germany and Switzerland, and significant difference (\(P\,{<}\,0.001\)) was observed between BSE-affected as well as the healthy cattle. Further, our data were also extensively compared with previous reports on BSE and highly significant (\(P\,{<}\,0.001\)) outcomes were observed. This result suggested negligible genetic susceptibility to BSE in gayals. To the best of our knowledge, this study is the first comprehensive deciphering information about the PRNP indel polymorphisms of 23 bp and 12 bp in gayals, a semiwild species of China.     

A viscoplastic model for the active component in cardiac muscle

The HMK model (Hunter et al. in Prog Biophys Mol Biol 69:289–331, 1998) proposes mechanobiological equations for the influence of intracellular calcium concentration \(\hbox {Ca}_\mathrm{i}\) on the evolution of bound calcium concentration \(\hbox {Ca}_\mathrm{b}\) and the tropomyosin kinetics parameter z, which model processes in the active component of the tension in cardiac muscle. The inelastic response due to actin-myosin crossbridge kinetics is modeled in the HMK model with a function Q that depends on the history of the rate of total stretch of the muscle fiber. Here, an alternative model is proposed which models the active component of the muscle fiber as a viscoplastic material. In particular, an evolution equation is proposed for the elastic stretch \(\lambda _\mathrm{a} \) in the active component. Specific forms of the constitutive equations are proposed and used to match experimental data. The proposed viscoplastic formulation allows for separate modeling of three processes: the high rate deactivation of crossbridges causing rapid reduction in active tension; the high but lower rate reactivation of crossbridges causing recovery of active tension; and the low rate relaxation effects characterizing the Hill model of muscles.     

A new indicator-based many-objective ant colony optimizer for continuous search spaces

In this paper, we propose a novel multi-objective ant colony optimizer (called iMOACO\(_{\mathbb {R}}\)) for continuous search spaces, which is based on ACO\(_{\mathbb {R}}\) and the R2 performance indicator. iMOACO\(_{\mathbb {R}}\) is the first multi-objective ant colony optimizer (MOACO) specifically designed to tackle continuous many-objective optimization problems (i.e., multi-objective optimization problems having four or more objectives). Our proposed iMOACO\(_{\mathbb {R}}\) is compared to three state-of-the-art multi-objective evolutionary algorithms (NSGA-III, MOEA/D and SMS-EMOA) and a MOACO algorithm called MOACO\(_{\mathbb {R}}\) using standard test problems and performance indicators taken from the specialized literature. Our experimental results indicate that iMOACO\(_{\mathbb {R}}\) is very competitive with respect to NSGA-III and MOEA/D and it is able to outperform SMS-EMOA and MOACO\(_{\mathbb {R}}\) in most of the test problems adopted.     

Optimal foraging behavior with an explicit consideration of within-individual behavioral variation: an example of predation

Animal behavior is flexible, and the same individual can exhibit variable expressions under the equivalent ecological situations (i.e., within-individual behavioral variation). This study examines the evolution of within-individual behavioral variation using an individual-based model. A common predation scenario is considered where a predator spends a period h to handle and consume a captured prey. The model assumes the handling time of the predator to be a random variable. The average and within-individual variance of handling time are described by \(\mu _h\) and \(\sigma _h^2\), respectively, where each individual has its own unique \(\mu _h\) and \(\sigma _h^2\). Using a genetic algorithm, the evolution of \(\sigma _h^2\) is traced. The results show that natural selection acts on both \(\mu _h\) and \(\sigma _h^2\), and the optimal behavioral variation depends on the density of prey. In particular, individuals with high behavioral variance \(\sigma _h^2\) are more likely selected when prey density is low. Individual based modeling can be a useful tool for studying the ultimate significance of within-individual behavioral variation and generating empirically testable predictions. The mechanisms of the evolution of within-individual behavioral variation and their ecological implications are discussed.     

Sorting signed circular permutations by super short operations

Background

One way to estimate the evolutionary distance between two given genomes is to determine the minimum number of large-scale mutations, or genome rearrangements, that are necessary to transform one into the other. In this context, genomes can be represented as ordered sequences of genes, each gene being represented by a signed integer. If no gene is repeated, genomes are thus modeled as signed permutations of the form \(\pi =(\pi _1 \pi _2 \ldots \pi _n)\), and in that case we can consider without loss of generality that one of them is the identity permutation \(\iota _n =(1 2 \ldots n)\), and that we just need to sort the other (i.e., transform it into \(\iota _n\)). The most studied genome rearrangement events are reversals, where a segment of the genome is reversed and reincorporated at the same location; and transpositions, where two consecutive segments are exchanged. Many variants, e.g., combining different types of (possibly constrained) rearrangements, have been proposed in the literature. One of them considers that the number of genes involved, in a reversal or a transposition, is never greater than two, which is known as the problem of sorting by super short operations (or SSOs).

Results and conclusions

All problems considering SSOs in permutations have been shown to be in \(\mathsf {P}\), except for one, namely sorting signed circular permutations by super short reversals and super short transpositions. Here we fill this gap by introducing a new graph structure called cyclic permutation graph and providing a series of intermediate results, which allows us to design a polynomial algorithm for sorting signed circular permutations by super short reversals and super short transpositions.

     

Resolving the viscoelasticity and anisotropy dependence of the mechanical properties of skin from a porcine model

The mechanical response of skin to external loads is influenced by anisotropy and viscoelasticity of the tissue, but the underlying mechanisms remain unclear. Here, we report a study of the main effects of tissue orientation (TO, which is linked to anisotropy) and strain rate (SR, a measure of viscoelasticity), as well as the interaction effects between the two factors, on the tensile properties of skin from a porcine model. Tensile testing to rupture of porcine skin tissue was conducted to evaluate the sensitivity of the tissue modulus of elasticity (E) and fracture-related properties, namely maximum stress \((\sigma _{U})\) and strain \((\varepsilon _{U})\) at \(\sigma _{U}\), to varying SR and TO. Specimens were excised from the abdominal skin in two orientations, namely parallel (P) and right angle (R) to the torso midline. Each TO was investigated at three SR levels, namely 0.007–0.015 \(\hbox {s}^{-1}\) (low), 0.040 \(\hbox {s}^{-1}\) (mid) and 0.065 \(\hbox {s}^{-1}\) (high). Two-factor analysis of variance revealed that the respective parameters responded differently to varying SR and TO. Significant changes in the \(\sigma _{U}\) were observed with different TOs but not with SR. The \(\varepsilon _{U}\) decreased significantly with increasing SR, but no significant variation was observed for different TOs. Significant changes in E were observed with different TOs; E increased significantly with increasing SR. More importantly, the respective mechanical parameters were not significantly influenced by interactions between SR and TO. These findings suggest that the trends associated with the changes in the skin mechanical properties may be attributed partly to differences in the anisotropy and viscoelasticity but not through any interaction between viscoelasticity and anisotropy.     

Invasions Slow Down or Collapse in the Presence of Reactive Boundaries

Motivated by the propagation of thin bacterial films around planar obstacles, this paper considers the dynamics of travelling wave solutions to the Fisher–KPP equation \(u_t = u(1-u) + u_{xx} + u_{yy}\) in a planar strip \(-\infty< x < \infty \), \(0 \le y \le L\). We examine the propagation of fronts in the presence of a mixed boundary condition (also referred to as a ‘partially absorbing’ or ‘reactive’ boundary) \(u_y = \alpha u\), with \(\alpha >0\), at \(y=0\). The presence of boundary conditions of this kind leads to the development of front solutions that propagate in x but contain transverse structure in y. Motivated by the observation that the speed of propagation in the Fisher–KPP equation is determined (for exponentially decaying initial conditions) by the behaviour at the leading edge, we analyse the linearised Fisher–KPP equation in order to estimate the speed of the stable travelling front, a function of the width L and the imposed boundary conditions. For wide strips the speed estimate based on the linearised equation agrees well with the results of numerical simulations. For narrow channels numerical simulations indicate that the stable front propagates more slowly, and for sufficiently small L or sufficiently large \(\alpha \) the front speed falls to zero and the front collapses. The reason for the collapse is the non-existence, far behind the front, of a stable positive equilibrium solution u(x, y). While existence of these equilibrium states can be demonstrated via phase plane arguments, the investigation of stability is similar to calculations of critical patch sizes carried out in similar ecological models.