Genomic signature of natural and anthropogenic stress in wild populations of the waterflea Daphnia magna: validation in space, time and experimental evolution

Natural populations are confronted with multiple selection pressures resulting in a mosaic of environmental stressors at the landscape level. Identifying the genetic underpinning of adaptation to these complex selection environments and assigning causes of natural selection within multidimensional selection regimes in the wild is challenging. The water flea Daphnia is a renowned ecological model system with its well-documented ecology, the possibility to analyse subfossil dormant egg banks and the short generation time allowing an experimental evolution approach. Capitalizing on the strengths of this model system, we here link candidate genome regions to three selection pressures, known to induce micro-evolutionary responses in Daphnia magna: fish predation, parasitism and land use. Using a genome scan approach in space, time and experimental evolution trials, we provide solid evidence of selection at the genome level under well-characterized environmental gradients in the wild and identify candidate genes linked to the three environmental stressors. Our study reveals differential selection at the genome level in Daphnia populations and provides evidence for repeatable patterns of local adaptation in a geographic mosaic of environmental stressors fuelled by standing genetic variation. Our results imply high evolutionary potential of local populations, which is relevant to understand the dynamics of trait changes in natural populations and their impact on community and ecosystem responses through eco-evolutionary feedbacks.     

Recent population decline and selection shape diversity of taxol-related genes

Taxanes are defensive metabolites produced by Taxus species (yews) and used in anticancer therapies. Despite their medical interest, patterns of natural diversity in taxane-related genes are unknown. We examined variation at five main genes of Taxus baccata in the Iberian Peninsula, a region where unique yew genetic resources are endangered. We looked at several gene features and applied complementary neutrality tests, including diversity/divergence tests, tests solely based on site frequency spectrum (SFS) and Zeng's compound tests. To account for specific demography, microsatellite data were used to infer historical changes in population size based on an Approximate Bayesian Computation (ABC) approach. Polymorphism-divergence tests pointed to positive selection for genes TBT and TAT and balancing selection for DBAT. In addition, neutrality tests based on SFS found that while a recent reduction in population size may explain most statistics' values, selection may still be in action in genes TBT and DBAT, at least in some populations. Molecular signatures on taxol genes suggest the action of frequent selective waves with different direction or intensity, possibly related to varying adaptive pressures produced by the host-enemy co-evolution on defence-related genes. Such natural selection processes may have produced taxane variants still undiscovered.     

How lizards survived blizzards: phylogeography of the Liolaemus lineomaculatus group (Liolaemidae) reveals multiple breaks and refugia in southern Patagonia and their concordance with other codistributed taxa

Patagonia was shaped by a complex geological history, including the Miocene uplift of the Andes, followed by volcanism, marine introgressions, and extreme climatic oscillations during Pliocene–Pleistocene glaciation–deglaciation cycles. The distributional patterns and phylogenetic relationships of southern patagonian animals and plants were affected in different ways, and those imprints are reflected in the seven phylogeographic breaks and eight refugia that have been previously proposed. In this study, we estimated time‐calibrated phylogenetic/phylogeographic patterns in lizards of the Liolaemus lineomaculatus group and relate them to historical Miocene‐to‐Pleistocene events of Patagonia and the previously proposed phylogeographic patterns. Individuals from 51 localities were sequenced for the mitochondrial marker (cyt‐b) and a subsample of individuals from each mitochondrial lineage was sequenced for one nuclear (LDA12D) and one slow evolving mitochondrial gene (12S). Our analyses revealed strong phylogeographic structure among lineages and, in most cases, no signal of demographic changes through time. The lineomaculatus group is composed of three strongly supported clades (lineomaculatus, hatcheri and kolengh + silvanae), and divergence estimates suggested their origins associated with the oldest known Patagonian glaciation (7–5 Ma); subsequent diversification within the lineomaculatus clade coincided with the large Pliocene glaciations (~3.5 Ma). The lineomaculatus clade includes nine strongly genetically and geographically structured lineages, five of which are interpreted as candidate species. Our findings suggest that some Liolaemus lineages have persisted in situ, each of them in a different refugium, through several glaciation–deglaciation cycles without demographic fluctuations. We also summarize and update qualitative evidence of some shared phylogeographic breaks and refugia among plants, rodents and lizards.     

Morphological evolution in land plants: new designs with old genes

The colonization and radiation of multicellular plants on land that started over 470 Ma was one of the defining events in the history of this planet. For the first time, large amounts of primary productivity occurred on the continental surface, paving the way for the evolution of complex terrestrial ecosystems and altering global biogeochemical cycles; increased weathering of continental silicates and organic carbon burial resulted in a 90 per cent reduction in atmospheric carbon dioxide levels. The evolution of plants on land was itself characterized by a series of radical transformations of their body plans that included the formation of three-dimensional tissues, de novo evolution of a multicellular diploid sporophyte generation, evolution of multicellular meristems, and the development of specialized tissues and organ systems such as vasculature, roots, leaves, seeds and flowers. In this review, we discuss the evolution of the genes and developmental mechanisms that drove the explosion of plant morphologies on land. Recent studies indicate that many of the gene families which control development in extant plants were already present in the earliest land plants. This suggests that the evolution of novel morphologies was to a large degree driven by the reassembly and reuse of pre-existing genetic mechanisms.     

QTL and quantitative genetic analysis of beak morphology reveals patterns of standing genetic variation in an Estrildid finch

The intra- and interspecific diversity of avian beak morphologies is one of the most compelling examples for the power of natural selection acting on a morphological trait. The development and diversification of the beak have also become a textbook example for evolutionary developmental biology, and variation in expression levels of several genes is known to causally affect beak shape. However, until now, no genomic polymorphisms have been identified, which are related to beak morphology in birds. QTL mapping does reveal the location of causal polymorphisms, albeit with poor spatial resolution. Here, we estimate heritability and genetic correlations for beak length, depth and width and perform a QTL linkage analysis for these traits based on 1404 informative single-nucleotide polymorphisms genotyped in a four-generation pedigree of 992 captive zebra finches (Taeniopygia guttata). Beak size, relative to body size, was sexually dimorphic (larger in males). Heritability estimates ranged from 0.47 for beak length to 0.74 for beak width. QTL mapping revealed four to five regions of significant or suggestive genome-wide linkage for each of the three beak dimensions (nine different regions in total). Eight out of 11 genes known to influence beak morphology are located in these nine peak regions. Five QTL do not cover known candidates demonstrating that yet unknown genes or regulatory elements may influence beak morphology in the zebra finch.     

Proteomic identification of predictive biomarkers of resistance to neoadjuvant chemotherapy in luminal breast cancer: a possible role for 14-3-3 theta/tau and tBID?

Introduction

Chemotherapy resistance is a major obstacle in effective neoadjuvant treatment for estrogen receptor-positive breast cancer. The ability to predict tumour response would allow chemotherapy administration to be directed towards only those patients who would benefit, thus maximising treatment efficiency. We aimed to identify putative protein biomarkers associated with chemotherapy resistance, using fresh tumour samples with antibody microarray analysis and then to perform pilot clinical validation experiments.

Materials and methods

Chemotherapy resistant and chemotherapy sensitive tumour samples were collected from breast cancer patients who had received anthracycline based neoadjuvant therapy consisting of epirubicin with cyclophosphamide followed by docetaxel. A total of 5 comparative proteomics experiments were performed using invasive ductal carcinomas which demonstrated estrogen receptor positivity (luminal subtype). Protein expression was compared between chemotherapy resistant and chemotherapy sensitive tumour samples using the Panorama XPRESS Profiler725 antibody microarray containing 725 antibodies from a wide variety of cell signalling and apoptosis pathways. A pilot series of archival samples was used for clinical validation of putative predictive biomarkers.

Results

AbMA analysis revealed 38 differentially expressed proteins which demonstrated at least 1.8 fold difference in expression in chemotherapy resistant tumours and 7 of these proteins (Zyxin, 14-3-3 theta/tau, tBID, Pinin, Bcl-xL, RIP and MyD88) were found in at least 2 experiments. Clinical validation in a pilot series of archival samples revealed 14-3-3 theta/tau and tBID to be significantly associated with chemotherapy resistance.

Conclusions

For the first time, antibody microarrays have been used to identify proteins associated with chemotherapy resistance using fresh breast cancer tissue. We propose a potential role for 14-3-3 theta/tau and tBID as predictive biomarkers of neoadjuvant chemotherapy resistance in breast cancer. Further validation in a larger sample series is now required.     

Comparative proteomic analysis of salt response proteins in seedling roots of two wheat varieties

A comparative proteomic analysis was made of salt response in seedling roots of wheat cultivars Jing-411 (salt tolerant) and Chinese Spring (salt sensitive) subjected to a range of salt stress concentrations (0.5%, 1.5% and 2.5%) for 2 days. One hundred and ninety eight differentially expressed protein spots (DEPs) were located with at least two-fold differences in abundance on 2-DE maps, of which 144 were identified by MALDI-TOF-TOF MS. These proteins were involved primarily in carbon metabolism (31.9%), detoxification and defense (12.5%), chaperones (5.6%) and signal transduction (4.9%). Comparative analysis showed that 41 DEPs were salt responsive with significant expression changes in both varieties under salt stress, and 99 (52 in Jing-411 and 47 in Chinese Spring) were variety specific. Only 15 and 9 DEPs in Jing-411 and Chinese Spring, respectively, were up-regulated in abundance under all three salt concentrations. All dynamics of the DEPs were analyzed across all treatments. Some salt responsive DEPs, such as guanine nucleotide-binding protein subunit beta-like protein, RuBisCO large subunit-binding protein subunit alpha and pathogenesis related protein 10, were up-regulated significantly in Jing-411 under all salt concentrations, whereas they were down-regulated in salinity-stressed Chinese Spring.     

Kidney proteome responses in the teleost fish Paralichthys olivaceus indicate a putative immune response against Streptococcus parauberis

The proteomic response to bacterial infection in a teleost fish (Paralichthys olivaceus) infected with Streptococcus parauberis was analyzed using label-free protein quantitation coupled with LC-MS(E) tandem mass spectrometry. A total of 82 proteins from whole kidney, a major lymphoid organ in this fish, were found to be differentially expressed between healthy and diseased fish analyzed 6, 24, 72 and 120 h post-infection. Among the differentially expressed proteins, those involved in mediating immune responses (e.g., heat shock proteins, cathepsins, goose-type lysozyme and complement components) were most significantly up-regulated by infection. In addition, cell division cycle 48 (CDC48) and calreticulin, which are associated with cellular recovery and glycoprotein synthesis, were up-regulated in the universal protein group, whereas the other proteins in that group were down-regulated. There was continuous activation of expression of immune-associated proteins during infection, but there was also loss of expression of proteins not involved in immune function. We expect that our findings regarding immune response at the protein level would offer new insight into the systemic response to bacterial infection of a major immune organ in teleost fish.     

Self/nonself recognition in Tuber melanosporum is not mediated by a heterokaryon incompatibility system

Vegetative incompatibility is a widespread phenomenon in filamentous ascomycetes, which limits formation of viable heterokaryons. Whether this phenomenon plays a role in maintaining the homokaryotic state of the hyphae during the vegetative growth of Tuber spp. Gene expression, polymorphism analysis as well as targeted in vitro experiments allowed us to test whether a heterokaryon incompatibility (HI) system operates in Tuber melanosporum. HI is controlled by different genetic systems, often involving HET domain genes and their partners whose interaction can trigger a cell death reaction. Putative homologues to HI-related genes previously characterized in Neurospora crassa and Podospora anserina were identified in the T. melanosporum genome. However, only two HET domain genes were found. In many other ascomycetes HET domains have been found within different genes including some members of the NWD (NACHT and WD-repeat associated domains) gene family of P. anserina. More than 50 NWD homologues were found in T. melanosporum but none of these contain a HET domain. All these T. melanosporum paralogs showed a conserved gene organization similar to the microexon genes only recently characterized in Schistosoma mansoni. Expression data of the annotated HI-like genes along with low allelic polymorphism suggest that they have cellular functions unrelated to HI. Moreover, morphological analyses did not provide evidence for HI reactions between pairs of genetically different T. melanosporum strains. Thus, the maintenance of the genetic integrity during the vegetative growth of this species likely depends on mechanisms that act before hyphal fusion.