Evolutionary dynamics of molecular markers during local adaptation: a case study in Drosophila subobscura

Background: Natural selection and genetic drift are major forces responsible for temporal genetic changes in populations. Furthermore, these evolutionary forces may interact with each other. Here we study the impact of an ongoing adaptive process at the molecular genetic level by analyzing the temporal genetic changes throughout 40 generations of adaptation to a common laboratory environment. Specifically, genetic variability, population differentiation and demographic structure were compared in two replicated groups of Drosophila subobscura populations recently sampled from different wild sources. Results: We found evidence for a decline in genetic variability through time, along with an increase in genetic differentiation between all populations studied. The observed decline in genetic variability was higher during the first 14 generations of laboratory adaptation. The two groups of replicated populations showed overall similarity in variability patterns. Our results also revealed changing demographic structure of the populations during laboratory evolution, with lower effective population sizes in the early phase of the adaptive process. One of the ten microsatellites analyzed showed a clearly distinct temporal pattern of allele frequency change, suggesting the occurrence of positive selection affecting the region around that particular locus. Conclusion: Genetic drift was responsible for most of the divergence and loss of variability between and within replicates, with most changes occurring during the first generations of laboratory adaptation. We also found evidence suggesting a selective sweep, despite the low number of molecular markers analyzed. Overall, there was a similarity of evolutionary dynamics at the molecular level in our laboratory populations, despite distinct genetic backgrounds and some differences in phenotypic evolution.

Odorant Receptor (Or) genes: polymorphism and divergence in the D. melanogaster and D. pseudoobscura Lineages

Background: In insects, like in most invertebrates, olfaction is the principal sensory modality, which provides animals with essential information for survival and reproduction. Odorant receptors are involved in this response, mediating interactions between an individual and its environment, as well as between individuals of the same or different species. The adaptive importance of odorant receptors renders them good candidates for having their variation shaped by natural selection. Methodology/Principal Findings: We analyzed nucleotide variation in a subset of eight Or genes located on the 3L chromosomal arm of Drosophila melanogaster in a derived population of this species and also in a population of Drosophila pseudoobscura. Some heterogeneity in the silent polymorphism to divergence ratio was detected in the D. melanogaster/D. simulans comparison, with a single gene (Or67b) contributing ~37% to the test statistic. However, no other signals of a very recent selective event were detected at this gene. In contrast, at the speciation timescale, the MK test uncovered the footprint of positive selection driving the evolution of two of the encoded proteins in both D. melanogaster ¿OR65c and OR67a ¿and D. pseudoobscura ¿OR65b1 and OR67c. Conclusions: The powerful polymorphism/divergence approach provided evidence for adaptive evolution at a rather high proportion of the Or genes studied after relatively recent speciation events. It did not provide, however, clear evidence for very recent selective events in either D. melanogaster or D. pseudoobscura.

Adaptive colour polymorphism of Acrida ungarica H. (Orthoptera: Acrididae) in a spatially heterogeneous environment

Intra-specific colour polymorphism provides a cryptic camouflage from predators in heterogeneous habitats. The orthoptera species, Acrida ungarica (Herbst, 1786) possess two well-distinguished colour morphs: brown and green and displays several disruptive colouration patterns within each morph to improve the crypsis. This study focused on how the features of the background environment relate to the proportion of the two morphs and to the intensity of disruptive colouration patterns in A. ungarica. As the two sexes are very distinct with respect to mass and length, we also distinctively tested the relationship for each sex. In accordance with the background matching hypothesis, we found that, for both sexes, the brown morph was in higher proportion at sites with a brown-dominant environment, and green morphs were in higher proportion in green-dominant environments. Globally, individuals in drier sites and in the drier year also had more intense disruptive colouration patterns, and brown morphs and females were also more striped. Colour patterns differed largely between populations and were significantly correlated with relevant environmental features. Even if A. ungarica is a polymorphic specialist, disruptive colouration still appears to provide strong benefits, particularly in some habitats. Moreover, because females are larger, they are less able to flee, which might explain the difference between sexes

Thermal evolution of gene expression profiles in Drosophila subobscura

Background: Despite its pervasiveness, the genetic basis of adaptation resulting in variation directly or indirectly related to temperature (climatic) gradients is poorly understood. By using 3-fold replicated laboratory thermal stocks covering much of the physiologically tolerable temperature range for the temperate (i.e., cold tolerant) species Drosophila subobscura we have assessed whole-genome transcriptional responses after three years of thermal adaptation, when the populations had already diverged for inversion frequencies, pre-adult life history components, and morphological traits. Total mRNA from each population was compared to a reference pool mRNA in a standard, highly replicated two-colour competitive hybridization experiment using cDNA microarrays.Results: A total of 306 (6.6%) cDNA clones were identified as 'differentially expressed' (following a false discovery rate correction) after contrasting the two furthest apart thermal selection regimes (i.e., 13°C vs . 22°C), also including four previously reported candidate genes for thermotolerance in Drosophila (Hsp26, Hsp68, Fst, and Treh). On the other hand, correlated patterns of gene expression were similar in cold- and warm-adapted populations. Analysis of functional categories defined by the Gene Ontology project point to an overrepresentation of genes involved in carbohydrate metabolism, nucleic acids metabolism and regulation of transcription among other categories. Although the location of differently expressed genes was approximately at random with respect to chromosomes, a physical mapping of 88 probes to the polytene chromosomes of D. subobscura has shown that a larger than expected number mapped inside inverted chromosomal segments.Conclusion: Our data suggest that a sizeable number of genes appear to be involved in thermal adaptation in Drosophila, with a substantial fraction implicated in metabolism. This apparently illustrates the formidable challenge to understanding the adaptive evolution of complex trait variation. Furthermore, some clustering of genes within inverted chromosomal sections was detected. Disentangling the effects of inversions will be obviously required in any future approach if we want to identify the relevant candidate genes.

Comparative genomic analysis of the odorant-binding protein family in 12 Drosophila genomes: purifying selection and birth-and-death evolution

Background: Chemoreception is a widespread mechanism that is involved in critical biologic processes, including individual and social behavior. The insect peripheral olfactory system comprises three major multigene families: the olfactory receptor (Or), the gustatory receptor (Gr), and the odorant-binding protein (OBP) families. Members of the latter family establish the first contact with the odorants, and thus constitute the first step in the chemosensory transduction pathway.Results: Comparative analysis of the OBP family in 12 Drosophila genomes allowed the identification of 595 genes that encode putative functional and nonfunctional members in extant species, with 43 gene gains and 28 gene losses (15 deletions and 13 pseudogenization events). The evolution of this family shows tandem gene duplication events, progressive divergence in DNA and amino acid sequence, and prevalence of pseudogenization events in external branches of the phylogenetic tree. We observed that the OBP arrangement in clusters is maintained across the Drosophila species and that purifying selection governs the evolution of the family; nevertheless, OBP genes differ in their functional constraints levels. Finally, we detect that the OBP repertoire evolves more rapidly in the specialist lineages of the Drosophila melanogaster group (D. sechellia and D. erecta) than in their closest generalists.Conclusion: Overall, the evolution of the OBP multigene family is consistent with the birth-and-death model. We also found that members of this family exhibit different functional constraints, which is indicative of some functional divergence, and that they might be involved in some of the specialization processes that occurred through the diversification of the Drosophila genus.

An adaptive multiscale method for density-driven instabilities

Accurate modeling of flow instabilities requires computational tools able to deal with several interacting scales, from the scale at which fingers are triggered up to the scale at which their effects need to be described. The Multiscale Finite Volume (MsFV) method offers a framework to couple fine-and coarse-scale features by solving a set of localized problems which are used both to define a coarse-scale problem and to reconstruct the fine-scale details of the flow. The MsFV method can be seen as an upscaling-downscaling technique, which is computationally more efficient than standard discretization schemes and more accurate than traditional upscaling techniques. We show that, although the method has proven accurate in modeling density-driven flow under stable conditions, the accuracy of the MsFV method deteriorates in case of unstable flow and an iterative scheme is required to control the localization error. To avoid large computational overhead due to the iterative scheme, we suggest several adaptive strategies both for flow and transport. In particular, the concentration gradient is used to identify a front region where instabilities are triggered and an accurate (iteratively improved) solution is required. Outside the front region the problem is upscaled and both flow and transport are solved only at the coarse scale. This adaptive strategy leads to very accurate solutions at roughly the same computational cost as the non-iterative MsFV method. In many circumstances, however, an accurate description of flow instabilities requires a refinement of the computational grid rather than a coarsening. For these problems, we propose a modified iterative MsFV, which can be used as downscaling method (DMsFV). Compared to other grid refinement techniques the DMsFV clearly separates the computational domain into refined and non-refined regions, which can be treated separately and matched later. This gives great flexibility to employ different physical descriptions in different regions, where different equations could be solved, offering an excellent framework to construct hybrid methods.

Evaluation of the RapidAir 457 Air Void Analyzer, March 2012

An adequate air void system is imperative to produce concrete with freeze-thaw durability in a wet freeze environment such as found in Iowa. Specifications rely on a percentage of air obtained in the plastic state by the pressure meter. Actual, in place air contents, of some concrete pavements in Iowa, have been found with reduced air content due to a number of factors such as excessive vibration and inadequate mixing. Determining hardened air void parameters is a time consuming process involving potential for human error. The RapidAir 457 air void analyzer is an automated device used to determine hardened air void parameters. The device is used in Europe and has been shown to quickly produce accurate and repeatable hardened air results. This research investigates how well the RapidAir 457 results correlate to plastic air content and the image analysis air technique. The repeatability and operator variation were also investigated, as well as, the impact of aggregate porosity and selection of threshold value on hardened air results.

Evaluating the Effectiveness of Red Light Running Camera Enforcement in Cedar Rapids and Developing Guidelines for Selection and Use of Red Light Running Countermeasures Final Report, November 2011

Red light running (RLR) is a problem in the US that has resulted in 165,000 injuries and 907 fatalities annually. In Iowa, RLR-related crashes make up 24.5 percent of all crashes and account for 31.7 percent of fatal and major injury crashes at signalized intersections. RLR crashes are a safety concern due to the increased likelihood of injury compared to other types of crashes. One tool used to combat red light running is automated enforcement in the form of RLR cameras. Automated enforcement, while effective, is often controversial. Cedar Rapids, Iowa installed RLR and speeding cameras at seven intersections across the city. The intersections were chosen based on crash rates and whether cameras could feasibly be placed at the intersection approaches. The cameras were placed starting in February 2010 with the last one becoming operational in December 2010. An analysis of the effect of the cameras on safety at these intersections was determined prudent in helping to justify the installation and effectiveness of the cameras. The objective of this research was to assess the safety effectiveness of the RLR program that has been implemented in Cedar Rapids. This was accomplished by analyzing data to determine changes in the following metrics:  Reductions in red light violation rates based on overall changes, time of day changes, and changes by lane  Effectiveness of the cameras over time  Time in which those running the red light enter the intersection  Changes in the average headway between vehicles entering the intersection

Low Water Stream Crossing: Design and Construction Recommendations; A Selection and Design Guide, February 2003

Most Iowa counties maintain low volume roads with at least one bridge or culvert that is structurally deficient or obsolete. In some counties the percentage of deficient drainage structures may be as high as 62%. Replacement with structures of similar size would require large capital expenditures that many counties cannot afford. Low water stream crossings (LWSCs) may be an acceptable lowcost alternative in some cases.

Automatic selection of tidal volume, respiratory frequency and minute ventilation in intubated ICU patients as start up procedure for closed-loop controlled ventilation.

OBJECTIVE: Before a patient can be connected to a mechanical ventilator, the controls of the apparatus need to be set up appropriately. Today, this is done by the intensive care professional. With the advent of closed loop controlled mechanical ventilation, methods will be needed to select appropriate start up settings automatically. The objective of our study was to test such a computerized method which could eventually be used as a start-up procedure (first 5-10 minutes of ventilation) for closed-loop controlled ventilation. DESIGN: Prospective Study. SETTINGS: ICU's in two adult and one children's hospital. PATIENTS: 25 critically ill adult patients (age > or = 15 y) and 17 critically ill children selected at random were studied. INTERVENTIONS: To stimulate 'initial connection', the patients were disconnected from their ventilator and transiently connected to a modified Hamilton AMADEUS ventilator for maximally one minute. During that time they were ventilated with a fixed and standardized breath pattern (Test Breaths) based on pressure controlled synchronized intermittent mandatory ventilation (PCSIMV). MEASUREMENTS AND MAIN RESULTS: Measurements of airway flow, airway pressure and instantaneous CO2 concentration using a mainstream CO2 analyzer were made at the mouth during application of the Test-Breaths. Test-Breaths were analyzed in terms of tidal volume, expiratory time constant and series dead space. Using this data an initial ventilation pattern consisting of respiratory frequency and tidal volume was calculated. This ventilation pattern was compared to the one measured prior to the onset of the study using a two-tailed paired t-test. Additionally, it was compared to a conventional method for setting up ventilators. The computer-proposed ventilation pattern did not differ significantly from the actual pattern (p > 0.05), while the conventional method did. However the scatter was large and in 6 cases deviations in the minute ventilation of more than 50% were observed. CONCLUSIONS: The analysis of standardized Test Breaths allows automatic determination of an initial ventilation pattern for intubated ICU patients. While this pattern does not seem to be superior to the one chosen by the conventional method, it is derived fully automatically and without need for manual patient data entry such as weight or height. This makes the method potentially useful as a start up procedure for closed-loop controlled ventilation.

Some spectral properties of the canonical solution operator to $\bar\partial$ on weighted Fock spaces

We characterize the Schatten class membership of the canonical solution operator to $\overline{\partial}$ acting on $L^2(e^{-2\phi})$, where $\phi$ is a subharmonic function with $\Delta\phi$ a doubling measure. The obtained characterization is in terms of $\Delta\phi$. As part of our approach, we study Hankel operators with anti-analytic symbols acting on the corresponding Fock space of entire functions in $L^2(e^{-2\phi})$

Life-history evolution and the polyphenic regulation of somatic maintenance and survival.

Here we discuss life-history evolution from the perspective of adaptive phenotypic plasticity, with a focus on polyphenisms for somatic maintenance and survival. Polyphenisms are adaptive discrete alternative phenotypes that develop in response to changes in the environment. We suggest that dauer larval diapause and its associated adult phenotypes in the nematode (Caenorhabditis elegans), reproductive dormancy in the fruit fly (Drosophila melanogaster) and other insects, and the worker castes of the honey bee (Apis mellifera) are examples of what may be viewed as the polyphenic regulation of somatic maintenance and survival. In these and other cases, the same genotype can--depending upon its environment--express either of two alternative sets of life-history phenotypes that differ markedly with respect to somatic maintenance, survival ability, and thus life span. This plastic modulation of somatic maintenance and survival has traditionally been underappreciated by researchers working on aging and life history. We review the current evidence for such adaptive life-history switches and their molecular regulation and suggest that they are caused by temporally and/or spatially varying, stressful environments that impose diversifying selection, thereby favoring the evolution of plasticity of somatic maintenance and survival under strong regulatory control. By considering somatic maintenance and survivorship from the perspective of adaptive life-history switches, we may gain novel insights into the mechanisms and evolution of aging.

The Effect of Inoculum Size on Selection of In Vitro Resistance to Vancomycin, Daptomycin, and Linezolid in Methicillin-Resistant Staphylococcus aureus.

OBJECTIVES: The inoculum effect (IE) is an increase in the minimum inhibitory concentration (MIC) at high bacterial densities. The effect of three inoculum sizes on the selection of resistance to vancomycin, daptomycin, and linezolid was investigated in methicillin-resistant Staphylococcus aureus (MRSA). METHODS: Low (10(4) CFU/ml), medium (10(6) CFU/ml), and high (10(8) CFU/ml) inocula of MRSA were exposed to twofold increasing concentrations of either drug during 15 days of cycling. MICs for low (MICL), medium (MICM), and high (MICH) inocula were determined daily. Conventional MICs were measured at days 1, 5, 10, and 15. Experiments were performed in triplicate. RESULTS: At the beginning of the experiment a small IE was observed for vancomycin (MICL=1 μg/ml, MICM=1-2 μg/ml, and MICH=2 μg/ml) and a significant IE for daptomycin (MICL=0.25 μg/ml, MICM=0.25-0.5 μg/ml, and MICH=2 μg/ml). Linezolid exhibited no IE at low and medium inocula (MICL=1 μg/ml and MICM=1-2 μg/ml), but with the high inoculum, concentrations up to 2,048 μg/ml did not fully inhibit visual growth. During cycling, increase of MIC was observed for all antibiotics. At day 15, MICL, MICM, and MICH of vancomycin were 2-4, 4-8, and 4-16 μg/ml and of daptomycin were 0.5-2, 8-128, and 64-256 μg/ml, respectively. MICL and MICM of linezolid were 1 and 2-4 μg/ml, respectively. Conventional MICs showed vancomycin and daptomycin selection of resistance since day 5 depending on the inocula. No selection of linezolid resistance was observed. CONCLUSIONS: Our results showed the importance of the inoculum size in the development of resistance. Measures aimed at lowering the inoculum at the site of infection should be used whenever possible in parallel to antimicrobial therapy.

Autoreactive T cells bypass negative selection and respond to self-antigen stimulation during infection.

Central and peripheral tolerance prevent autoimmunity by deleting the most aggressive CD8(+) T cells but they spare cells that react weakly to tissue-restricted antigen (TRA). To reveal the functional characteristics of these spared cells, we generated a transgenic mouse expressing the TCR of a TRA-specific T cell that had escaped negative selection. Interestingly, the isolated TCR matches the affinity/avidity threshold for negatively selecting T cells, and when developing transgenic cells are exposed to their TRA in the thymus, only a fraction of them are eliminated but significant numbers enter the periphery. In contrast to high avidity cells, low avidity T cells persist in the antigen-positive periphery with no signs of anergy, unresponsiveness, or prior activation. Upon activation during an infection they cause autoimmunity and form memory cells. Unexpectedly, peptide ligands that are weaker in stimulating the transgenic T cells than the thymic threshold ligand also induce profound activation in the periphery. Thus, the peripheral T cell activation threshold during an infection is below that of negative selection for TRA. These results demonstrate the existence of a level of self-reactivity to TRA to which the thymus confers no protection and illustrate that organ damage can occur without genetic predisposition to autoimmunity.