Airfoil flow and noise computation using monotonically integrated large eddy simulation and acoustic analogy: Effect of the grid resolution.

A new scalable Monotonically Integrated Large Eddy Simulation (MILES) method based on the Compact Accurately Boundary-Adjusting high-REsolution Technique (CABARET) has been applied for the simulation of unsteady flow around NACA0012 airfoil at Re = 400,000 and M = 0.058. The flow solution is coupled with the Ffowcs Williams-Hawkings formulation for far-field noise prediction. The computational modeling results are presented for several computational grid resolutions: 8, 16, and 32 million grid cells and compared with the experimental data available.

Approximate Bayesian Computation for Smoothing

We consider a method for approximate inference in hidden Markov models (HMMs). The method circumvents the need to evaluate conditional densities of observations given the hidden states. It may be considered an instance of Approximate Bayesian Computation (ABC) and it involves the introduction of auxiliary variables valued in the same space as the observations. The quality of the approximation may be controlled to arbitrary precision through a parameter ε > 0. We provide theoretical results which quantify, in terms of ε, the ABC error in approximation of expectations of additive functionals with respect to the smoothing distributions. Under regularity assumptions, this error is, where n is the number of time steps over which smoothing is performed. For numerical implementation, we adopt the forward-only sequential Monte Carlo (SMC) scheme of [14] and quantify the combined error from the ABC and SMC approximations. This forms some of the first quantitative results for ABC methods which jointly treat the ABC and simulation errors, with a finite number of data and simulated samples. © Taylor & Francis Group, LLC.

Global analysis of dynamical decision-making models through local computation around the hidden saddle

Bistable dynamical switches are frequently encountered in mathematical modeling of biological systems because binary decisions are at the core of many cellular processes. Bistable switches present two stable steady-states, each of them corresponding to a distinct decision. In response to a transient signal, the system can flip back and forth between these two stable steady-states, switching between both decisions. Understanding which parameters and states affect this switch between stable states may shed light on the mechanisms underlying the decision-making process. Yet, answering such a question involves analyzing the global dynamical (i.e., transient) behavior of a nonlinear, possibly high dimensional model. In this paper, we show how a local analysis at a particular equilibrium point of bistable systems is highly relevant to understand the global properties of the switching system. The local analysis is performed at the saddle point, an often disregarded equilibrium point of bistable models but which is shown to be a key ruler of the decision-making process. Results are illustrated on three previously published models of biological switches: two models of apoptosis, the programmed cell death and one model of long-term potentiation, a phenomenon underlying synaptic plasticity. © 2012 Trotta et al.

Cubically convergent iterations for invariant subspace computation

We propose a Newton-like iteration that evolves on the set of fixed dimensional subspaces of ℝ n and converges locally cubically to the invariant subspaces of a symmetric matrix. This iteration is compared in terms of numerical cost and global behavior with three other methods that display the same property of cubic convergence. Moreover, we consider heuristics that greatly improve the global behavior of the iterations.

Riemannian geometry of Grassmann manifolds with a view on algorithmic computation

We give simple formulas for the canonical metric, gradient, Lie derivative, Riemannian connection, parallel translation, geodesics and distance on the Grassmann manifold of p-planes in ℝn. In these formulas, p-planes are represented as the column space of n × p matrices. The Newton method on abstract Riemannian manifolds proposed by Smith is made explicit on the Grassmann manifold. Two applications - computing an invariant subspace of a matrix and the mean of subspaces - are worked out.

A Newton algorithm for invariant subspace computation with large basins of attraction

We study the global behaviour of a Newton algorithm on the Grassmann manifold for invariant subspace computation. It is shown that the basins of attraction of the invariant subspaces may collapse in case of small eigenvalue gaps. A Levenberg-Marquardt-like modification of the algorithm with low numerical cost is proposed. A simple strategy for choosing the parameter is shown to dramatically enlarge the basins of attraction of the invariant subspaces while preserving the fast local convergence.

Evolutionary Conservation of Dazl Genomic Organization and its Continuous and Dynamic Distribution Throughout Germline Development in Gynogenetic Gibel Carp

To investigate germline development and germ cell specification, we identified a Dazl homolog (CagDazl) from gynogenetic gibel carp (Carassius auratus gibelio). Its cDNA sequence and BAC clone sequence analyses revealed the genomic organization conservation and conserved synteny of the Dazl family members and their neighborhood genes among vertebrates, especially in fish. Moreover, a polyclonal antibody specific to CagDazl was produced and used to examine its expression and distribution throughout germline development at protein level. Firstly, ovary-specific expression pattern of CagDazl was confirmed in adult tissues by RT-PCR and Western blot. In addition, in situ hybridization and immunofluorescence localization demonstrated its specific expression in germ cells, and both its transcript and protein were localized to germ plasm. Then, co-localization of CagDazl and mitochondrial cloud was found, confirming that CagDazl transcript and its protein are germ plasm component and move via METRO pathway during oogenesis. Furthermore, the CagDazl is abundant and continuous throughout germline development and germ cell specification including primordial germ cell (PGC) formation, oogonium differentiation, oocyte development, and embryogenesis, and the dynamic distribution occurs at different development stages. The data suggest that maternal CagDazl might play an important role in gibel carp PGC formation. Therefore, CagDazl is a useful and specific marker for tracing germ plasm and germ cell development in the gynogenetic gibel carp. In addition, in comparison with previous studies in sexual reproduction species, the continuous and dynamic distribution of CagDazl protein in the germ plasm throughout the life cycle seems to have significant implication in sex evolution of vertebrates. J. Exp. Zool. (Mol. Deu. Euol.) 312B:855-871, 2009. (C) 2009 Wiley-Liss, Inc.

Distinct Evolutionary Patterns Between Two Duplicated Color Vision Genes Within Cyprinid Fishes

We investigated the molecular evolution of duplicated color vision genes (LWS-1 and SWS2) within cyprinid fish, focusing on the most cavefish-rich genus-Sinocyclocheilus. Maximum likelihood-based codon substitution approaches were used to analyze the evolution of vision genes. We found that the duplicated color vision genes had unequal evolutionary rates, which may lead to a related function divergence. Divergence of LWS-1 was strongly influenced by positive selection causing an accelerated rate of substitution in the proportion of pocket-forming residues. The SWS2 pigment experienced divergent selection between lineages, and no positively selected site was found. A duplicate copy of LWS-1 of some cyprinine species had become a pseudogene, but all SWS2 sequences remained intact in the regions examined in the cyprinid fishes examined in this study. The pseudogenization events did not occur randomly in the two copies of LWS-1 within Sinocyclocheilus species. Some cave species of Sinocyclocheilus with numerous morphological specializations that seem to be highly adapted for caves, retain both intact copies of color vision genes in their genome. We found some novel amino acid substitutions at key sites, which might represent interesting target sites for future mutagenesis experiments. Our data add to the increasing evidence that duplicate genes experience lower selective constraints and in some cases positive selection following gene duplication. Some of these observations are unexpected and may provide insights into the effect of caves on the evolution of color vision genes in fishes.

Phylogenetic relationships of Danio within the order cypriniformes: A framework for comparative and evolutionary studies of a model species

The evolutionary relationships of species of Danio and the monophyly and phylogenetic placement of the genus within the family Cyprinidae and subfamily Rasborinae provide fundamentally important phyloinformatics necessary for direct evaluations of an array of pertinent questions in modern comparative biology. Although the genus Danio is not one of the most diverse within the family, Danio rerio is one of the most important model species in biology. Many investigations have used this species or presumed close relatives to address specific questions that have lasting impact on the hypothesis and theory of development in vertebrates. Largely lacking from this approach has been a holistic picture of the exact phylogenetic or evolutionary relationships of this species and its close relatives. One thing that has been learned over the previous century is that many organismal attributes (e.g., developmental pathways, ecologies, behaviors, speciation) are historically constrained and their origins and functions are best explained via a phylogenetic approach. Herein, we provide a molecular evaluation of the phylogenetic placement of the model species Danio rerio within the genus Danio and among hypothesized closely related species and genera. Our analysis is derived from data using two nuclear genes (RAG1, rhodopsin) and five mitochondrial genes (ND4, ND4L, ND5, COI, cyt b) evaluated using parsimony, maximum likelihood, and Bayesian analyses. The family Cyprinidae is resolved as monophyletic but the subfamily Rasborinae (priority over Danioinae) is an unnatural assemblage. Danio is identified as a monophyletic group sister to a clade inclusive of the genera Chela, Microrasbora, Devario, and Inlecypris, not Devario nor Esomus as hypothesized in previous studies. Danio rerio is sister to D. kyathit among the species of Danio evaluated in this analysis. Microrasbora and Rasbora are non-monophyletic assemblages; however, Boraras is monophyletic.

Identification of novel SINEs from Cyprinidae and their evolutionary significance

Short interspersed nuclear elements (SINEs) are widespread among eukaryotic genomes. They are repetitive DNA sequences that have been amplified by retrotransposition. In this study, a class of SINEs were isolated from the Opsariichthys bidens genome, and named Opsar. Sequence analysis confirmed that Opsar is a new class of typical SINEs derived from tRNA molecules. With the tRNA-derived region of Opsar and through BLASTN search, we further identified Zb-SINEs from the zebrafish genome, which includes two groups: Zb-SINE-A and Zb-SINE-B. The Zb-SINE-A group comprises subfamilies of -Al--A5, and the Zb-SINE-B group is a dimer of the tRNA(Ala)-derived region and shares a similar dimeric composition to Alu. Zb-SINEs are composed of three distinct regions: a 5 end tRNA-derived region, a tRNA-unrelated region and a 3 end AT-rich region. The flanking regions are AT rich. The average length of Zb-SINEs elements is about 340 6p. Zb-SINEs account for as much as 0.1% of the whole zebrafish genome. About 70% of the Zb-SINEs are on chromosomes 11, 18, and 19. These Zb-SINEs were characterized by PCR and dot hybridization. The distribution pattern of Zb-SINEs in genome strongly supports the master genes model. The tRNA-derived regions of Opsar and Zb-SINEs were compared with the tRNA(Ala) gene, and they showed 76% similarity, indicating that Opsar and Zb-SINEs originated from an inactive tRNA(Ala) sequence or a tRNA(Ala)-like sequence. In view of the evolutionary status of zebrafish in the Cyprinidae, we deduced that Zb-SINEs were a very old class of interspersed sequences.

A possible evolutionary lag in the relationship between freshwater mussels and European bitterling

The coevolutionary dynamics between European bitterling Rhodeus amarus and freshwater unionid mussels, which the former parasitize by laying eggs on their gills, were tested. In a series of experiments fish preferences and mussel responses were compared in parasites and hosts of recent (Europe) and ancient (Asia) sympatry. Rhodeus amarus readily oviposited on the gills of all mussel species tested. Fish that laid their eggs on the gills of Asian Anodonta woodiana, however, suffered a dramatic reduction in reproductive success compared to fish that oviposited on the gills of European mussels: Unio pictorum, Unio tumidus, Anodonta anatina and Anodonta cygnea. This difference was the result of egg ejection behaviour by mussels rather than the unsuitability of the internal gill environment for European bitterling embryo development. The ejection response of mussels with a long sympatry with European bitterling was considerably more pronounced than that of mussels with a substantially shorter sympatry. The data support a coevolutionary arms race between bitterling and mussels and point to an evolutionary lag in the relationship between R. amarus and its European mussel hosts. (c) 2007 The Authors. Journal compilation (c) 2007 The Fisheries Society of the British Isles.

The co-evolutionary relationship between bitterling fishes and freshwater mussels: insights from interspecific comparisons

Hypothesis: In parasites that use hosts for offspring development, adults may base oviposition decisions on a range of host traits related either to host quality or the co-evolutionary relationship between parasite and host. We examined whether host quality or co-evolutionary dynamics drive the use of hosts in the bitterling-mussel relationship. Organisms: Six species of bitterling fish (Acheilognathinae) and eight species of freshwater mussels (Unionidae, Corbiculidae) that are used by bitterling for oviposition. Site of experiments: Experimental tanks in Wuhan, China, at the site of the natural distribution of the studied species. Methods: Three experiments that controlled for host accessibility and interspecific interactions were conducted to identify host preferences among bitterling fishes and their mussel hosts. We started with a broad interspecific comparison. We then tested bitterling behavioural choices, their temporal stability, and mussel host ejection behaviour of the eggs of generalist and specialist bitterling species. Finally, we measured host mussel quality based on respiration rate and used published studies on mussel gill structure to infer mussel suitability as hosts for bitterling eggs. Results: We found significant interspecific differences among bitterling species in their use of mussel hosts. Bitterling species varied in their level of host specificity and identity of preferred hosts. Host preferences were flexible even among apparently specialized species and fishes switched their preferences adaptively when the quality of individuals of preferred host species declined. Mussels varied considerably in their response to oviposition through egg ejections. Host preference by a generalist bitterling species correlated positively with host quality measured as the efficiency of the mussel gills to extract oxygen from inhaled water. Host ability to eject bitterling eggs correlated positively with their relative respiration rate, probably due to a higher velocity of water circulating in the mussel gill chamber.

Comparison of evolutionary rates in the mitochondrial DNA cytochrome b gene and control region and their implications for phylogeny of the Cobitoidea (Teleostei : Cypriniformes)

It is widely accepted that mitochondrial DNA (mtDNA) control region evolves faster than protein encoding genes with few exceptions. In the present study, we sequenced the mitochondrial cytochrome b gene (cyt b) and control region (CR) and compared their rates in 93 specimens representing 67 species of loaches and some related taxa in the Cobitoidea (Order Cypriniformes). The results showed that sequence divergences of the CR were broadly higher than those of the cyt b (about 1.83 times). However, in considering only closely related species, CR sequence evolution was slower than that of cyt b gene (ratio of CR/cyt b is 0.78), a pattern that is found to be very common in Cypriniformes. Combined data of the cyt b and CR were used to estimate the phylogenetic relationship of the Cobitoidea by maximum parsimony, neighbor-joining, and Bayesian methods. With Cyprinus carpio and Danio rerio as outgroups, three analyses identified the same four lineages representing four subfamilies of loaches, with Botiinae on the basal-most clade. The phylogenctic relationship of the Cobitoidea was ((Catostomidae + Gyrinocheilidae) + (Botiinae + (Balitorinae + (Cobitinae + Nemacheilinae)))), which indicated that Sawada's Cobitidae (including Cobitinae and Botiinae) was not monophyletic. Our molecular phylogenetic analyses are in very close agreement with the phylogenetic results based on the morphological data proposed by Nalbant and Bianco, wherein these four subfamilies were elevated to the family level as Botiidae, Balitoridae, Cobitidae, and Nemacheilidae. (c) 2005 Elsevier Inc. All rights reserved.