ASU Scholarship Showcase

It is known that in classical fluids turbulence typically occurs at high Reynolds numbers. But can turbulence occur at low Reynolds numbers? Here we investigate the transition to turbulence in the classic Taylor-Couette system in which the rotating fluids are manufactured ferrofluids with magnetized nanoparticles embedded in liquid carriers. We find that, in the presence of a magnetic field transverse to the symmetry axis of the system, turbulence can occur at Reynolds numbers that are at least one order of magnitude smaller than those in conventional fluids. This is established by extensive computational ferrohydrodynamics through a detailed investigation of transitions in the flow structure, and characterization of behaviors of physical quantities such as the energy, the wave number, and the angular momentum through the bifurcations. A finding is that, as the magnetic field is increased, onset of turbulence can be determined accurately and reliably. Our results imply that experimental investigation of turbulence may be feasible by using ferrofluids. Our study of transition to and evolution of turbulence in the Taylor-Couette ferrofluidic flow system provides insights into the challenging problem of turbulence control.

Plasmodium vivax is the most prevalent malarial species in South America and exerts a substantial burden on the populations it affects. The control and eventual elimination of P. vivax are global health priorities. Genomic research contributes to this objective by improving our understanding of the biology of P. vivax and through the development of new genetic markers that can be used to monitor efforts to reduce malaria transmission. Here we analyze whole-genome data from eight field samples from a region in Cordóba, Colombia where malaria is endemic. We find considerable genetic diversity within this population, a result that contrasts with earlier studies suggesting that P. vivax had limited diversity in the Americas. We also identify a selective sweep around a substitution known to confer resistance to sulphadoxine-pyrimethamine (SP). This is the first observation of a selective sweep for SP resistance in this species. These results indicate that P. vivax has been exposed to SP pressure even when the drug is not in use as a first line treatment for patients afflicted by this parasite. We identify multiple non-synonymous substitutions in three other genes known to be involved with drug resistance in Plasmodium species. Finally, we found extensive microsatellite polymorphisms. Using this information we developed 18 polymorphic and easy to score microsatellite loci that can be used in epidemiological investigations in South America.

Inbreeding in hermaphroditic plants can occur through two different mechanisms: biparental inbreeding, when a plant mates with a related individual, or self-fertilization, when a plant mates with itself. To avoid inbreeding, many hermaphroditic plants have evolved self-incompatibility (SI) systems which prevent or limit self-fertilization. One particular SI system—homomorphic SI—can also reduce biparental inbreeding. Homomorphic SI is found in many angiosperm species, and it is often assumed that the additional benefit of reduced biparental inbreeding may be a factor in the success of this SI system. To test this assumption, we developed a spatially-explicit, individual-based simulation of plant populations that displayed three different types of homomorphic SI. We measured the total level of inbreeding avoidance by comparing each population to a self-compatible population (NSI), and we measured biparental inbreeding avoidance by comparing to a population of self-incompatible plants that were free to mate with any other individual (PSI).

Because biparental inbreeding is more common when offspring dispersal is limited, we examined the levels of biparental inbreeding over a range of dispersal distances. We also tested whether the introduction of inbreeding depression affected the level of biparental inbreeding avoidance. We found that there was a statistically significant decrease in autozygosity in each of the homomorphic SI populations compared to the PSI population and, as expected, this was more pronounced when seed and pollen dispersal was limited. However, levels of homozygosity and inbreeding depression were not reduced. At low dispersal, homomorphic SI populations also suffered reduced female fecundity and had smaller census population sizes. Overall, our simulations showed that the homomorphic SI systems had little impact on the amount of biparental inbreeding in the population especially when compared to the overall reduction in inbreeding compared to the NSI population. With further study, this observation may have important consequences for research into the origin and evolution of homomorphic self-incompatibility systems.

Under models of isolation-by-distance, population structure is determined by the probability of identity-by-descent between pairs of genes according to the geographic distance between them. Well established analytical results indicate that the relationship between geographical and genetic distance depends mostly on the neighborhood size of the population which represents a standardized measure of gene flow. To test this prediction, we model local dispersal of haploid individuals on a two-dimensional landscape using seven dispersal kernels: Rayleigh, exponential, half-normal, triangular, gamma, Lomax and Pareto. When neighborhood size is held constant, the distributions produce similar patterns of isolation-by-distance, confirming predictions. Considering this, we propose that the triangular distribution is the appropriate null distribution for isolation-by-distance studies. Under the triangular distribution, dispersal is uniform over the neighborhood area which suggests that the common description of neighborhood size as a measure of an effective, local panmictic population is valid for popular families of dispersal distributions. We further show how to draw random variables from the triangular distribution efficiently and argue that it should be utilized in other studies in which computational efficiency is important.

We investigate fundamental nonlinear dynamics of ferrofluidic Taylor-Couette flow - flow confined be-tween two concentric independently rotating cylinders - consider small aspect ratio by solving the ferro-hydrodynamical equations, carrying out systematic bifurcation analysis. Without magnetic field, we find steady flow patterns, previously observed with a simple fluid, such as those containing normal one- or two vortex cells, as well as anomalous one-cell and twin-cell flow states. However, when a symmetry-breaking transverse magnetic field is present, all flow states exhibit stimulated, finite two-fold mode. Various bifurcations between steady and unsteady states can occur, corresponding to the transitions between the two-cell and one-cell states. While unsteady, axially oscillating flow states can arise, we also detect the emergence of new unsteady flow states. In particular, we uncover two new states: one contains only the azimuthally oscillating solution in the configuration of the twin-cell flow state, and an-other a rotating flow state. Topologically, these flow states are a limit cycle and a quasiperiodic solution on a two-torus, respectively. Emergence of new flow states in addition to observed ones with classical fluid, indicates that richer but potentially more controllable dynamics in ferrofluidic flows, as such flow states depend on the external magnetic field.

We investigate the dynamics of ferrofluidic wavy vortex flows in the counter-rotating Taylor-Couette system, with a focus on wavy flows with a mixture of the dominant azimuthal modes. Without external magnetic field flows are stable and pro-grade with respect to the rotation of the inner cylinder. More complex behaviors can arise when an axial or a transverse magnetic field is applied. Depending on the direction and strength of the field, multi-stable wavy states and bifurcations can occur. We uncover the phenomenon of flow pattern reversal as the strength of the magnetic field is increased through a critical value. In between the regimes of pro-grade and retrograde flow rotations, standing waves with zero angular velocities can emerge. A striking finding is that, under a transverse magnetic field, a second reversal in the flow pattern direction can occur, where the flow pattern evolves into pro-grade rotation again from a retrograde state. Flow reversal is relevant to intriguing phenomena in nature such as geomagnetic reversal. Our results suggest that, in ferrofluids, flow pattern reversal can be induced by varying a magnetic field in a controlled manner, which can be realized in laboratory experiments with potential applications in the development of modern fluid devices.

Introduction: Abundance of immune cells has been shown to have prognostic and predictive significance in many tumor types. Beyond abundance, the spatial organization of immune cells in relation to cancer cells may also have significant functional and clinical implications. However there is a lack of systematic methods to quantify spatial associations between immune and cancer cells.

Methods: We applied ecological measures of species interactions to digital pathology images for investigating the spatial associations of immune and cancer cells in breast cancer. We used the Morisita-Horn similarity index, an ecological measure of community structure and predator–prey interactions, to quantify the extent to which cancer cells and immune cells colocalize in whole-tumor histology sections. We related this index to disease-specific survival of 486 women with breast cancer and validated our findings in a set of 516 patients from different hospitals.

Results: Colocalization of immune cells with cancer cells was significantly associated with a disease-specific survival benefit for all breast cancers combined. In HER2-positive subtypes, the prognostic value of immune-cancer cell colocalization was highly significant and exceeded those of known clinical variables. Furthermore, colocalization was a significant predictive factor for long-term outcome following chemotherapy and radiotherapy in HER2 and Luminal A subtypes, independent of and stronger than all known clinical variables.

Conclusions: Our study demonstrates how ecological methods applied to the tumor microenvironment using routine histology can provide reproducible, quantitative biomarkers for identifying high-risk breast cancer patients. We found that the clinical value of immune-cancer interaction patterns is highly subtype-specific but substantial and independent to known clinicopathologic variables that mostly focused on cancer itself. Our approach can be developed into computer-assisted prediction based on histology samples that are already routinely collected.

In this study, we present a novel methodology to infer indel parameters from multiple sequence alignments (MSAs) based on simulations. Our algorithm searches for the set of evolutionary parameters describing indel dynamics which best fits a given input MSA. In each step of the search, we use parametric bootstraps and the Mahalanobis distance to estimate how well a proposed set of parameters fits input data. Using simulations, we demonstrate that our methodology can accurately infer the indel parameters for a large variety of plausible settings. Moreover, using our methodology, we show that indel parameters substantially vary between three genomic data sets: Mammals, bacteria, and retroviruses. Finally, we demonstrate how our methodology can be used to simulate MSAs based on indel parameters inferred from real data sets.

In a meta-analysis published by myself and co-authors, we report differences in the life history risk factors for estrogen receptor negative (ER−) and estrogen receptor positive (ER+) breast cancers. Our meta-analysis did not find the association of ER− breast cancer risk with fast life history characteristics that Hidaka and Boddy suggest in their response to our article. There are a number of possible explanations for the differences between their conclusions and the conclusions we drew from our meta-analysis, including limitations of our meta-analysis and methodological challenges in measuring and categorizing estrogen receptor status. These challenges, along with the association of ER+ breast cancer with slow life history characteristics, may make it challenging to find a clear signal of ER− breast cancer with fast life history characteristics, even if that relationship does exist. The contradictory results regarding breast cancer risk and life history characteristics illustrate a more general challenge in evolutionary medicine: often different sub-theories in evolutionary biology make contradictory predictions about disease risk. In this case, life history models predict that breast cancer risk should increase with faster life history characteristics, while the evolutionary mismatch hypothesis predicts that breast cancer risk should increase with delayed reproduction. Whether life history tradeoffs contribute to ER− breast cancer is still an open question, but current models and several lines of evidence suggest that it is a possibility.

It has long been accepted that modern reproductive patterns are likely contributors to breast cancer susceptibility because of their influence on hormones such as estrogen and the importance of these hormones in breast cancer. We conducted a meta-analysis to assess whether this ‘evolutionary mismatch hypothesis’ can explain susceptibility to both estrogen receptor positive (ER-positive) and estrogen receptor negative (ER-negative) cancer. Our meta-analysis includes a total of 33 studies and examines parity, age of first birth and age of menarche broken down by estrogen receptor status. We found that modern reproductive patterns are more closely linked to ER-positive than ER-negative breast cancer. Thus, the evolutionary mismatch hypothesis for breast cancer can account for ER-positive breast cancer susceptibility but not ER-negative breast cancer.