Emergence of sensory selection mechanisms in Artificial Life simulations

Background The evolutionary advantages of selective attention are unclear. Since the study of selective attention began, it has been suggested that the nervous system only processes the most relevant stimuli because of its limited capacity [1]. An alternative proposal is that action planning requires the inhibition of irrelevant stimuli, which forces the nervous system to limit its processing [2]. An evolutionary approach might provide additional clues to clarify the role of selective attention. Methods We developed Artificial Life simulations wherein animals were repeatedly presented two objects, "left" and "right", each of which could be "food" or "non-food." The animals' neural networks (multilayer perceptrons) had two input nodes, one for each object, and two output nodes to determine if the animal ate each of the objects. The neural networks also had a variable number of hidden nodes, which determined whether or not it had enough capacity to process both stimuli (Table 1). The evolutionary relevance of the left and the right food objects could also vary depending on how much the animal's fitness was increased when ingesting them (Table 1). We compared sensory processing in animals with or without limited capacity, which evolved in simulations in which the objects had the same or different relevances. Table 1. Nine sets of simulations were performed, varying the values of food objects and the number of hidden nodes in the neural networks. The values of left and right food were swapped during the second half of the simulations. Non-food objects were always worth -3. The evolution of neural networks was simulated by a simple genetic algorithm. Fitness was a function of the number of food and non-food objects each animal ate and the chromosomes determined the node biases and synaptic weights. During each simulation, 10 populations of 20 individuals each evolved in parallel for 20,000 generations, then the relevance of food objects was swapped and the simulation was run again for another 20,000 generations. The neural networks were evaluated by their ability to identify the two objects correctly. The detectability (d') for the left and the right objects was calculated using Signal Detection Theory [3]. Results and conclusion When both stimuli were equally relevant, networks with two hidden nodes only processed one stimulus and ignored the other. With four or eight hidden nodes, they could correctly identify both stimuli. When the stimuli had different relevances, the d' for the most relevant stimulus was higher than the d' for the least relevant stimulus, even when the networks had four or eight hidden nodes. We conclude that selection mechanisms arose in our simulations depending not only on the size of the neuron networks but also on the stimuli's relevance for action.

Vascular wall stem cells. Selection and conditioning of progenitors useful for cell therapy. A pathological case study

The arterial wall contains MSCs with mesengenic and angiogenic abilities. These multipotent precursors have been isolated from variously-sized human adult segments, belying the notion that vessel wall is a relatively quiescent tissue. Recently, our group identified in normal human arteries a vasculogenic niche and subsequently isolated and characterized resident MSCs (VW-MSCs) with angiogenic ability and multilineage potential. To prove that VW-MSCs are involved in normal and pathological vascular remodeling, we used a long-term organ culture system; this method was of critical importance to follow spontaneous 3-D vascular remodeling without any influence of blood cells. Next we tried to identify and localize in situ the VW-MSCs and to understand their role in the vascular remodeling in failed arterial homografts. Subsequently, we isolated this cell population and tested in vitro their multilineage differentiation potential through immunohistochemical, immunofluorescence, RT-PCR and ultrastructural analysis. From 25-30cm2 of each vascular wall homograft sample, we isolated a cell population with MSCs properties; these cells expressed MSC lineage molecules (CD90, CD44, CD105, CD29, CD73), stemness (Notch-1, Oct-4, Sca-1, Stro-1) and pericyte markers (NG2) whilst were negative for hematopoietic and endothelial markers (CD34, CD133, CD45, KDR, CD146, CD31 and vWF). MSCs derived from failed homografts (H-MSCs) exhibited adipogenic, osteogenic and chondrogenic potential but scarce propensity to angiogenic and leiomyogenic differentiation. The present study demonstrates that failed homografts contain MSCs with morphological, phenotypic and functional MSCs properties; H-MSCs are long-lived in culture, highly proliferating and endowed with prompt ability to differentiate into adipocytes, osteocytes and chondrocytes; compared with VW-MSCs from normal arteries, H-MSCs show a failure in angiogenic and leiomyogenic differentiation. A switch in MSCs plasticity could be the basis of pathological remodeling and contribute to aneurysmal failure of arterial homografts. The study of VW-MSCs in a pathological setting indicate that additional mechanisms are involved in vascular diseases; their knowledge will be useful for opening new therapeutic options in cardiovascular diseases.

Evaluation of Antitumoral activity of bone targeted drugs/conventional chemotherapies and identification of biomarkers for the selection of patients with breast cancer for the bone targeted therapy in adjuvant setting

Bone metastases are responsible for different clinical complications defined as skeletal-related events (SREs) such as pathologic fractures, spinal cord compression, hypercalcaemia, bone marrow infiltration and severe bone pain requiring palliative radiotherapy. The general aim of these three years research period was to improve the management of patients with bone metastases through two different approaches of translational research. Firstly in vitro preclinical tests were conducted on breast cancer cells and on indirect co-colture of cancer cells and osteoclasts to evaluate bone targeted therapy singly and in combination with conventional chemotherapy. The study suggests that zoledronic acid has an antitumor activity in breast cancer cell lines. Its mechanism of action involves the decrease of RAS and RHO, as in osteoclasts. Repeated treatment enhances antitumor activity compared to non-repeated treatment. Furthermore the combination Zoledronic Acid + Cisplatin induced a high antitumoral activity in the two triple-negative lines MDA-MB-231 and BRC-230. The p21, pMAPK and m-TOR pathways were regulated by this combined treatment, particularly at lower Cisplatin doses. A co-colture system to test the activity of bone-targeted molecules on monocytes-breast conditioned by breast cancer cells was also developed. Another important criticism of the treatment of breast cancer patients, is the selection of patients who will benefit of bone targeted therapy in the adjuvant setting. A retrospective case-control study on breast cancer patients to find new predictive markers of bone metastases in the primary tumors was performed. Eight markers were evaluated and TFF1 and CXCR4 were found to discriminate between patients with relapse to bone respect to patients with no evidence of disease. In particular TFF1 was the most accurate marker reaching a sensitivity of 63% and a specificity of 79%. This marker could be a useful tool for clinicians to select patients who could benefit for bone targeted therapy in adjuvant setting.

Alternative reproductive tactics, sexual selection and the effects of inbreeding in the ant Hypoponera opacior

Inbreeding can lead to a fitness reduction due to the unmasking of deleterious recessive alleles and the loss of heterosis. Therefore, most sexually reproducing organisms avoid inbreeding, often by disperal. Besides the avoidance of inbreeding, dispersal lowers intraspecific competition on a local scale and leads to a spreading of genotypes into new habitats. In social insects, winged reproductives disperse and mate during nuptial flights. Therafter, queens independently found a new colony. However, some species also produce wingless sexuals as an alternative reproductive tactic. Wingless sexuals mate within or close to their colony and queens either stay in the nest or they found a new colony by budding. During this dependent colony foundation, wingless queens are accompanied by a fraction of nestmate workers. The production of wingless reproductives therefore circumvents the risks associated with dispersal and independent colony foundation. However, the absence of dispersal can lead to inbreeding and local competition.rnIn my PhD-project, I investigated the mating biology of Hypoponera opacior, an ant that produces winged and wingless reproductives in a population in Arizona. Besides the investigation of the annual reproductive cycle, I particularly focused on the consequences of wingless reproduction. An analysis of sex ratios in wingless sexuals should reveal the relative importance of local resource competition among queens (that mainly compete for the help of workers) and local mate competition among males. Further, sexual selection was expected to act on wingless males that were previously found to mate with and mate-guard pupal queens in response to local mate competition. We studied whether males are able to adapt their mating behaviour to the current competitive situation in the nest and which traits are under selection in this mating situation. Last, we investigated the extent and effects of inbreeding. As the species appeared to produce non-dispersive males and queens quite frequently, we assumed to find no or only weak negative effects of inbreeding and potentially mechanisms that moderate inbreeding levels despite frequent nest-matings.rnWe found that winged and wingless males and queens are produced during two separate seasons of the year. Winged sexuals emerge in early summer and conduct nuptial flights in July, when climate conditions due to frequent rainfalls lower the risks of dispersal and independent colony foundation. In fall, wingless sexuals are produced that reproduce within the colonies leading to an expansion on the local scale. The absence of dispersal during this second reproductive season resulted in a local genetic population viscosity and high levels of inbreeding within the colonies. Male-biased sex ratios in fall indicated a greater importance of local resource competition among queens than local mate competition among males. Males were observed to adjust mate-guarding durations to the competitive situation (i.e. the number of competing males and pupae) in the nest, an adaptation that helps maximising their reproductive success. Further, sexual selection was found to act on the timing of emergence as well as on body size in these males, i.e. earlier emerging and larger males show a higher mating success. Genetic analyses revealed that wingless males do not actively avoid inbreeding by choosing less related queens as mating partners. Further, we detected diploid males, a male type that is produced instead of diploid females if close relatives mate. In contrast to many other Hymenopteran species, diploid males were here viable and able to sire sterile triploid offspring. They did not differ in lifespan, body size and mating success from “normal” haploid males. Hence, diploid male production in H. opacior is less costly than in other social Hymenopteran species. No evidence of inbreeding depression was found on the colony level but more inbred colonies invested more resources into the production of sexuals. This effect was more pronounced in the dispersive summer generation. The increased investment in outbreeding sexuals can be regarded as an active strategy to moderate the extent and effects of inbreeding. rnIn summary, my thesis describes an ant species that has evolved alternative reproductive tactics as an adaptation to seasonal environmental variations. Hereby, the species is able to maintain its adaptive mating system without suffering from negative effects due to the absence of dispersal flights in fall.rn

Negative laparotomy in trauma: are we getting better?

One of the trauma surgeons' daily challenges is the balancing act between negative laparotomy and missed abdominal injury. We opted to characterize the indications that prompted a negative trauma exploratory laparotomy and the rate of missed abdominal injuries in an effort to optimize patient selection for laparotomy. At the Los Angeles County + University of Southern California Medical Center, negative laparotomies and missed injuries are consecutively captured and reviewed at the weekly mortality + morbidity (MM) conferences. All written reports of the MM meetings from January 2003 to December 2008 were reviewed to identify all patients who underwent a negative laparotomy or a laparotomy as a result of an initially missed abdominal injury. Over the 6-year study period, a total of 1871 laparotomies were performed, of which 73 (3.9%) were negative. The rate of missed injuries requiring subsequent laparotomy was 1.3 per cent (25 of 1871). The negative laparotomy rate and the rate of missed injuries did not vary significantly during the study period (2.8 to 4.7%, P = 0.875, and 0.7 to 2.9%, P = 0.689). Penetrating mechanisms accounted for the majority of negative laparotomies (58.9%). The primary indication for negative laparotomy was peritonitis (54.8%) followed by hypotension (28.8%) and suspicious computed tomographic scan findings (27.4%). The complication rate after negative laparotomy was 14.5 per cent, and of these, 10.1 per cent were directly related to the procedure. A low but steady rate of negative laparotomies and missed abdominal injuries after trauma remains. Negative laparotomies and missed abdominal injuries when they occur are still associated with significant complication rates and a prolonged length of stay.

Expression and diversity of Echinococcus multilocularis AgB genes in secondarily infected mice: evaluating the influence of T-cell immune selection on antigenic variation

The T-cell-mediated immune response exhibits a crucial function in the control of the intrahepatic proliferation of Echinococcus multilocularis larvae in mice and humans, both being natural intermediate hosts of the parasite. Antigen B (AgB), a metabolized Echinococcus spp. lipoprotein, contributes to the modulation of the T-cell immune response, and distinct sites of the corresponding AgB1, AgB3 and AgB4 genes were shown to be under positive selection pressure. Since several AgB gene variants are present in a single Echinococcus metacestode, we used secondary E. multilocularis infections in BALB/c and in athymic nude mice (devoid of T-cell responses) to analyze the effect of the cellular immune response on the expression and diversity of EmAgB1-EmAgB4 genes. We demonstrated hereby that EmAgB transcripts were less abundant in nude mice during the early phase of infection (at one month post-infection), and that EmAgB2 is simultaneously down-regulated when compared to the other three genes. A negative relationship exists between the level of transcription and diversity of EmAgB genes. Moreover, no excess of non-synonymous substitutions was found among the distinct EmAgB alleles from a single host. Together, these results pointed to the effect of purifying selection, which seemed to eliminate the detrimental AgB variants generated during the development of the metacestode within the peritoneal cavity of its intermediate host.

Isolated autosomal dominant growth hormone deficiency: stimulating mutant GH-1 gene expression drives GH-1 splice-site selection, cell proliferation, and apoptosis

The majority of mutations that cause isolated GH deficiency type II (IGHD II) affect splicing of GH-1 transcripts and produce a dominant-negative GH isoform lacking exon 3 resulting in a 17.5-kDa isoform, which further leads to disruption of the GH secretory pathway. A clinical variability in the severity of the IGHD II phenotype depending on the GH-1 gene alteration has been reported, and in vitro and transgenic animal data suggest that the onset and severity of the phenotype relates to the proportion of 17.5-kDa produced. The removal of GH in IGHD creates a positive feedback loop driving more GH expression, which may itself increase 17.5-kDa isoform productions from alternate splice sites in the mutated GH-1 allele. In this study, we aimed to test this idea by comparing the impact of stimulated expression by glucocorticoids on the production of different GH isoforms from wild-type (wt) and mutant GH-1 genes, relying on the glucocorticoid regulatory element within intron 1 in the GH-1 gene. AtT-20 cells were transfected with wt-GH or mutated GH-1 variants (5'IVS-3 + 2-bp T->C; 5'IVS-3 + 6 bp T->C; ISEm1: IVS-3 + 28 G->A) known to cause clinical IGHD II of varying severity. Cells were stimulated with 1 and 10 mum dexamethasone (DEX) for 24 h, after which the relative amounts of GH-1 splice variants were determined by semiquantitative and quantitative (TaqMan) RT-PCR. In the absence of DEX, only around 1% wt-GH-1 transcripts were the 17.5-kDa isoform, whereas the three mutant GH-1 variants produced 29, 39, and 78% of the 17.5-kDa isoform. DEX stimulated total GH-1 gene transcription from all constructs. Notably, however, DEX increased the amount of 17.5-kDa GH isoform relative to the 22- and 20-kDa isoforms produced from the mutated GH-1 variants, but not from wt-GH-1. This DEX-induced enhancement of 17.5-kDa GH isoform production, up to 100% in the most severe case, was completely blocked by the addition of RU486. In other studies, we measured cell proliferation rates, annexin V staining, and DNA fragmentation in cells transfected with the same GH-1 constructs. The results showed that that the 5'IVS-3 + 2-bp GH-1 gene mutation had a more severe impact on those measures than the splice site mutations within 5'IVS-3 + 6 bp or ISE +28, in line with the clinical severity observed with these mutations. Our findings that the proportion of 17.5-kDa produced from mutant GH-1 alleles increases with increased drive for gene expression may help to explain the variable onset progression, and severity observed in IGHD II.

Evaluation of the colorimetric VITEK 2 card for identification of gram-negative nonfermentative rods: comparison to 16S rRNA gene sequencing

Ninety strains of a collection of well-identified clinical isolates of gram-negative nonfermentative rods collected over a period of 5 years were evaluated using the new colorimetric VITEK 2 card. The VITEK 2 colorimetric system identified 53 (59%) of the isolates to the species level and 9 (10%) to the genus level; 28 (31%) isolates were misidentified. An algorithm combining the colorimetric VITEK 2 card and 16S rRNA gene sequencing for adequate identification of gram-negative nonfermentative rods was developed. According to this algorithm, any identification by the colorimetric VITEK 2 card other than Achromobacter xylosoxidans, Acinetobacter sp., Burkholderia cepacia complex, Pseudomonas aeruginosa, and Stenotrophomonas maltophilia should be subjected to 16S rRNA gene sequencing when accurate identification of nonfermentative rods is of concern.

An integral conservative gridding--algorithm using Hermitian curve interpolation

The problem of re-sampling spatially distributed data organized into regular or irregular grids to finer or coarser resolution is a common task in data processing. This procedure is known as 'gridding' or 're-binning'. Depending on the quantity the data represents, the gridding-algorithm has to meet different requirements. For example, histogrammed physical quantities such as mass or energy have to be re-binned in order to conserve the overall integral. Moreover, if the quantity is positive definite, negative sampling values should be avoided. The gridding process requires a re-distribution of the original data set to a user-requested grid according to a distribution function. The distribution function can be determined on the basis of the given data by interpolation methods. In general, accurate interpolation with respect to multiple boundary conditions of heavily fluctuating data requires polynomial interpolation functions of second or even higher order. However, this may result in unrealistic deviations (overshoots or undershoots) of the interpolation function from the data. Accordingly, the re-sampled data may overestimate or underestimate the given data by a significant amount. The gridding-algorithm presented in this work was developed in order to overcome these problems. Instead of a straightforward interpolation of the given data using high-order polynomials, a parametrized Hermitian interpolation curve was used to approximate the integrated data set. A single parameter is determined by which the user can control the behavior of the interpolation function, i.e. the amount of overshoot and undershoot. Furthermore, it is shown how the algorithm can be extended to multidimensional grids. The algorithm was compared to commonly used gridding-algorithms using linear and cubic interpolation functions. It is shown that such interpolation functions may overestimate or underestimate the source data by about 10-20%, while the new algorithm can be tuned to significantly reduce these interpolation errors. The accuracy of the new algorithm was tested on a series of x-ray CT-images (head and neck, lung, pelvis). The new algorithm significantly improves the accuracy of the sampled images in terms of the mean square error and a quality index introduced by Wang and Bovik (2002 IEEE Signal Process. Lett. 9 81-4).

Multiparameter flow cytometry characterization of MHC class I negative mouse bone marrow cells

BACKGROUND: MHC-I down-regulation was described in foetal liver progenitors, and two different subsets of adult bone marrow derived stem cells. These cells, namely, MHC-I-/Thy1+ bone marrow derived liver stem cells (BMDLSC) and the multipotent adult progenitors (MAPC) differentiated into functioning hepatocytes. The aim of this paper was to characterize the MHC-I negative bone marrow compartment as it pertains to BMDLSC and MAPC. MATERIAL/METHODS: We performed multiparameter flow-cytometry analyses of the MHC-I negative compartment using hematopoietic (CD45, Ter119), and stem cell markers (Thy1.2, c-Kit, IL-3R, CD34) in adult mice. RESULTS: When analysing CD45 and Ter119 expression, the MHC-I negative bone marrow compartment divides into four sub-populations: 1. CD45-/Ter119+: 86.0+/-4.4%; 2. CD45+/Ter119+: 0.2+/-0.1%; 3. CD45+/Ter119-: 11.6+/-3.0%; 4. CD45-/Ter119-: 2.0+/-2.1%. Stem cells markers were only expressed on MHC-I negative/ CD45+/Ter119- cells. In vivo, MAPC (Ter119-/CD45- cells) are composed of MHC-I negative (24%) and MHC-I positive cells and do not express any of the stem cell markers tested. CONCLUSIONS: In conclusion, mouse BMDLSC and MAPC are two distinct stem cell populations. Down-regulation of MHC-I was the only common characteristic found between BMDLSC and MAPC suggesting that selection of MHC-I negative cells might represent an efficient strategy to enrich for bone marrow stem cells with liver developmental potential.

Clinical and molecular features of methicillin-resistant, coagulase-negative staphylococci of pets and horses

OBJECTIVES To determine the antibiotic resistance and fingerprint profiles of methicillin-resistant coagulase-negative staphylococci (MRCoNS) from animal infections among different practices and examine the history of antibiotic treatment. METHODS Isolates were identified by mass spectrometry and tested for antimicrobial resistance by broth dilution, microarrays and sequence analysis of the topoisomerases. Diversity was assessed by PFGE, icaA PCR and staphylococcal cassette chromosome mec (SCCmec), arginine catabolic mobile element (ACME) and multilocus sequence typing. Clinical records were examined retrospectively. RESULTS MRCoNS were identified as Staphylococcus epidermidis (n=20), Staphylococcus haemolyticus (n=17), Staphylococcus hominis (n=3), Staphylococcus capitis (n=1), Staphylococcus cohnii (n=1) and Staphylococcus warneri (n=1). PFGE identified one clonal lineage in S. hominis isolates and several in S. haemolyticus and S. epidermidis. Fourteen sequence types were identified in S. epidermidis, with sequence type 2 (ST2) and ST5 being predominant. Ten isolates contained SCCmec IV, seven contained SCCmec V and the others were non-typeable. ACMEs were detected in 11 S. epidermidis isolates. One S. hominis and 10 S. epidermidis isolates were icaA positive. In addition to mecA-mediated β-lactam resistance, the most frequent resistance was to gentamicin/kanamycin [aac(6')-Ie-aph(2')-Ia, aph(3')-III] (n=34), macrolides/lincosamides [erm(C), erm(A), msr, lnu(A)] (n=31), tetracycline [tet(K)] (n=22), streptomycin [str, ant(6)-Ia] (n=20), trimethoprim [dfr(A), dfr(G)] (n=17), sulfamethoxazole (n = 34) and fluoroquinolones [amino acid substitutions in GyrA and GrlA] (n=30). Clinical data suggest selection through multiple antibiotic courses and emphasize the importance of accurate diagnosis and antibiograms. CONCLUSIONS MRCoNS from animal infection sites are genetically heterogeneous multidrug-resistant strains that represent a new challenge in the prevention and therapy of infections in veterinary clinics.

Extended-spectrum cephalosporin-resistant Gram-negative organisms in livestock: an emerging problem for human health?

Escherichia coli, Salmonella spp. and Acinetobacter spp. are important human pathogens. Serious infections due to these organisms are usually treated with extended-spectrum cephalosporins (ESCs). However, in the past two decades we have faced a rapid increasing of infections and colonization caused by ESC-resistant (ESC-R) isolates due to production of extended-spectrum-β-lactamases (ESBLs), plasmid-mediated AmpCs (pAmpCs) and/or carbapenemase enzymes. This situation limits drastically our therapeutic armamentarium and puts under peril the human health. Animals are considered as potential reservoirs of multidrug-resistant (MDR) Gram-negative organisms. The massive and indiscriminate use of antibiotics in veterinary medicine has contributed to the selection of ESC-R E. coli, ESC-R Salmonella spp. and, to less extent, MDR Acinetobacter spp. among animals, food, and environment. This complex scenario is responsible for the expansion of these MDR organisms which may have life-threatening clinical significance. Nowadays, the prevalence of food-producing animals carrying ESC-R E. coli and ESC-R Salmonella (especially those producing CTX-M-type ESBLs and the CMY-2 pAmpC) has reached worryingly high values. More recently, the appearance of carbapenem-resistant isolates (i.e., VIM-1-producing Enterobacteriaceae and NDM-1 or OXA-23-producing Acinetobacter spp.) in livestock has even drawn greater concerns. In this review, we describe the aspects related to the spread of the above MDR organisms among pigs, cattle, and poultry, focusing on epidemiology, molecular mechanisms of resistance, impact of antibiotic use, and strategies to contain the overall problem. The link and the impact of ESC-R organisms of livestock origin for the human scenario are also discussed.

Testing an alternative search algorithm for compound identification with the 'Wiley Registry of Tandem Mass Spectral Data. MSforID'

A tandem mass spectral database system consists of a library of reference spectra and a search program. State-of-the-art search programs show a high tolerance for variability in compound-specific fragmentation patterns produced by collision-induced decomposition and enable sensitive and specific 'identity search'. In this communication, performance characteristics of two search algorithms combined with the 'Wiley Registry of Tandem Mass Spectral Data, MSforID' (Wiley Registry MSMS, John Wiley and Sons, Hoboken, NJ, USA) were evaluated. The search algorithms tested were the MSMS search algorithm implemented in the NIST MS Search program 2.0g (NIST, Gaithersburg, MD, USA) and the MSforID algorithm (John Wiley and Sons, Hoboken, NJ, USA). Sample spectra were acquired on different instruments and, thus, covered a broad range of possible experimental conditions or were generated in silico. For each algorithm, more than 30,000 matches were performed. Statistical evaluation of the library search results revealed that principally both search algorithms can be combined with the Wiley Registry MSMS to create a reliable identification tool. It appears, however, that a higher degree of spectral similarity is necessary to obtain a correct match with the NIST MS Search program. This characteristic of the NIST MS Search program has a positive effect on specificity as it helps to avoid false positive matches (type I errors), but reduces sensitivity. Thus, particularly with sample spectra acquired on instruments differing in their Setup from tandem-in-space type fragmentation, a comparably higher number of false negative matches (type II errors) were observed by searching the Wiley Registry MSMS.

Genomic signatures of relaxed disruptive selection associated with speciation reversal in whitefish

Background: Speciation reversal: the erosion of species differentiation via an increase in introgressive hybridization due to the weakening of previously divergent selection regimes, is thought to be an important, yet poorly understood, driver of biodiversity loss. Our study system, the Alpine whitefish (Coregonus spp.) species complex is a classic example of a recent postglacial adaptive radiation: forming an array of endemic lake flocks, with the independent origination of similar ecotypes among flocks. However, many of the lakes of the Alpine radiation have been seriously impacted by anthropogenic nutrient enrichment, resulting in a collapse in neutral genetic and phenotypic differentiation within the most polluted lakes. Here we investigate the effects of eutrophication on the selective forces that have shaped this radiation, using population genomics. We studied eight sympatric species assemblages belonging to five independent parallel adaptive radiations, and one species pair in secondary contact. We used AFLP markers, and applied FST outlier (BAYESCAN, DFDIST) and logistic regression analyses (MATSAM), to identify candidate regions for disruptive selection in the genome and their associations with adaptive traits within each lake flock. The number of outlier and adaptive trait associated loci identified per lake were then regressed against two variables (historical phosphorus concentration and contemporary oxygen concentration) representing the strength of eutrophication. Results: Whilst we identify disruptive selection candidate regions in all lake flocks, we find similar trends, across analysis methods, towards fewer disruptive selection candidate regions and fewer adaptive trait/candidate loci associations in the more polluted lakes. Conclusions: Weakened disruptive selection and a concomitant breakdown in reproductive isolating mechanisms in more polluted lakes has lead to increased gene flow between coexisting Alpine whitefish species. We hypothesize that the resulting higher rates of interspecific recombination reduce either the number or extent of genomic islands of divergence surrounding loci evolving under disruptive natural selection. This produces the negative trend seen in the number of selection candidate loci recovered during genome scans of whitefish species flocks, with increasing levels of anthropogenic eutrophication: as the likelihood decreases that AFLP restriction sites will fall within regions of heightened genomic divergence and therefore be classified as FST outlier loci. This study explores for the first time the potential effects of human-mediated relaxation of disruptive selection on heterogeneous genomic divergence between coexisting species.

Frequency-Dependent Selection Predicts Patterns of Radiations and Biodiversity

Most empirical studies support a decline in speciation rates through time, although evidence for constant speciation rates also exists. Declining rates have been explained by invoking pre-existing niches, whereas constant rates have been attributed to non-adaptive processes such as sexual selection and mutation. Trends in speciation rate and the processes underlying it remain unclear, representing a critical information gap in understanding patterns of global diversity. Here we show that the temporal trend in the speciation rate can also be explained by frequency-dependent selection. We construct a frequency-dependent and DNA sequence-based model of speciation. We compare our model to empirical diversity patterns observed for cichlid fish and Darwin's finches, two classic systems for which speciation rates and richness data exist. Negative frequency-dependent selection predicts well both the declining speciation rate found in cichlid fish and explains their species richness. For groups like the Darwin's finches, in which speciation rates are constant and diversity is lower, speciation rate is better explained by a model without frequency-dependent selection. Our analysis shows that differences in diversity may be driven by incipient species abundance with frequency-dependent selection. Our results demonstrate that genetic-distance-based speciation and frequency-dependent selection are sufficient to explain the high diversity observed in natural systems and, importantly, predict decay through time in speciation rate in the absence of pre-existing niches.