Reproducibility and Consistency of Proteomic Experiments on Natural Populations of a Non-Model Aquatic Insect

Population proteomics has a great potential to address evolutionary and ecological questions, but its use in wild populations of non-model organisms is hampered by uncontrolled sources of variation. Here we compare the response to temperature extremes of two geographically distant populations of a diving beetle species (Agabus ramblae) using 2-D DIGE. After one week of acclimation in the laboratory under standard conditions, a third of the specimens of each population were placed at either 4 or 27°C for 12 h, with another third left as a control. We then compared the protein expression level of three replicated samples of 2–3 specimens for each treatment. Within each population, variation between replicated samples of the same treatment was always lower than variation between treatments, except for some control samples that retained a wider range of expression levels. The two populations had a similar response, without significant differences in the number of protein spots over- or under-expressed in the pairwise comparisons between treatments. We identified exemplary proteins among those differently expressed between treatments, which proved to be proteins known to be related to thermal response or stress. Overall, our results indicate that specimens collected in the wild are suitable for proteomic analyses, as the additional sources of variation were not enough to mask the consistency and reproducibility of the response to the temperature treatments.     

Establishment of a coastal fish in the Azores: recent colonisation or sudden expansion of an ancient relict population?

The processes and timescales associated with ocean-wide changes in the distribution of marine species have intrigued biologists since Darwin''s earliest insights into biogeography. The Azores, a mid-Atlantic volcanic archipelago located >1000 km off the European continental shelf, offers ideal opportunities to investigate phylogeographic colonisation scenarios. The benthopelagic sparid fish known as the common two-banded seabream (Diplodus vulgaris) is now relatively common along the coastline of the Azores archipelago, but was virtually absent before the 1990s. We employed a multiple genetic marker approach to test whether the successful establishment of the Azorean population derives from a recent colonisation from western continental/island populations or from the demographic explosion of an ancient relict population. Results from nuclear and mtDNA sequences show that all Atlantic and Mediterranean populations belong to the same phylogroup, though microsatellite data indicate significant genetic divergence between the Azorean sample and all other locations, as well as among Macaronesian, western Iberian and Mediterranean regions. The results from Approximate Bayesian Computation indicate that D. vulgaris has likely inhabited the Azores for ∼40 (95% confidence interval (CI): 5.5–83.6) to 52 (95% CI: 6.32–89.0) generations, corresponding to roughly 80–150 years, suggesting near-contemporary colonisation, followed by a more recent demographic expansion that could have been facilitated by changing climate conditions. Moreover, the lack of previous records of this species over the past century, together with the absence of lineage separation and the presence of relatively few private alleles, do not exclude the possibility of an even more recent colonisation event.     

New insight into the SSC8 genetic determination of fatty acid composition in pigs

Background

Fat content and fatty acid composition in swine are becoming increasingly studied because of their effect on sensory and nutritional quality of meat. A QTL (quantitative trait locus) for fatty acid composition in backfat was previously detected on porcine chromosome 8 (SSC8) in an Iberian x Landrace F₂ intercross. More recently, a genome-wide association study detected the same genomic region for muscle fatty acid composition in an Iberian x Landrace backcross population. ELOVL6, a strong positional candidate gene for this QTL, contains a polymorphism in its promoter region (ELOVL6:c.-533C < T), which is associated with percentage of palmitic and palmitoleic acids in muscle and adipose tissues. Here, a combination of single-marker association and the haplotype-based approach was used to analyze backfat fatty acid composition in 470 animals of an Iberian x Landrace F₂ intercross genotyped with 144 SNPs (single nucleotide polymorphisms) distributed along SSC8.

Results

Two trait-associated SNP regions were identified at 93 Mb and 119 Mb on SSC8. The strongest statistical signals of both regions were observed for palmitoleic acid (C16:1(n-7)) content and C18:0/C16:0 and C18:1(n-7)/C16:1(n-7) elongation ratios. MAML3 and SETD7 are positional candidate genes in the 93 Mb region and two novel microsatellites in MAML3 and nine SNPs in SETD7 were identified. No significant association for the MAML3 microsatellite genotypes was detected. The SETD7:c.700G > T SNP, although statistically significant, was not the strongest signal in this region. In addition, the expression of MAML3 and SETD7 in liver and adipose tissue varied among animals, but no association was detected with the polymorphisms in these genes. In the 119 Mb region, the ELOVL6:c.-533C > T polymorphism showed a strong association with percentage of palmitic and palmitoleic fatty acids and elongation ratios in backfat.

Conclusions

Our results suggest that the polymorphisms studied in MAML3 and SETD7 are not the causal mutations for the QTL in the 93 Mb region. However, the results for ELOVL6 support the hypothesis that the ELOVL6:c.-533C > T polymorphism has a pleiotropic effect on backfat and intramuscular fatty acid composition and that it has a role in the determination of the QTL in the 119 Mb region.     

Deux nouveaux représentants du genre Ascidiella(Clavatoraceae,Charophyta)

Two new species of the genus AscidiellaGrambast are described: A. inflata et A. irregularis. They come from the Lower Cretaceous of the Iberian peninsula.     

Epigenetic regulation of the ELOVL6 gene is associated with a major QTL effect on fatty acid composition in pigs

Background

In previous studies on an Iberian x Landrace cross, we have provided evidence that supported the porcine ELOVL6 gene as the major causative gene of the QTL on pig chromosome 8 for palmitic and palmitoleic acid contents in muscle and backfat. The single nucleotide polymorphism (SNP) ELOVL6:c.-533C > T located in the promoter region of ELOVL6 was found to be highly associated with ELOVL6 expression and, accordingly, with the percentages of palmitic and palmitoleic acids in longissimus dorsi and adipose tissue. The main goal of the current work was to further study the role of ELOVL6 on these traits by analyzing the regulation of the expression of ELOVL6 and the implication of ELOVL6 polymorphisms on meat quality traits in pigs.

Results

High-throughput sequencing of BAC clones that contain the porcine ELOVL6 gene coupled to RNAseq data re-analysis showed that two isoforms of this gene are expressed in liver and adipose tissue and that they differ in number of exons and 3’UTR length. Although several SNPs in the 3’UTR of ELOVL6 were associated with palmitic and palmitoleic acid contents, this association was lower than that previously observed with SNP ELOVL6:c.-533C > T. This SNP is in full linkage disequilibrium with SNP ELOVL6:c.-394G > A that was identified in the binding site for estrogen receptor alpha (ERα). Interestingly, the ELOVL6:c.-394G allele is associated with an increase in methylation levels of the ELOVL6 promoter and with a decrease of ELOVL6 expression. Therefore, ERα is clearly a good candidate to explain the regulation of ELOVL6 expression through dynamic epigenetic changes in the binding site of known regulators of ELOVL6 gene, such as SREBF1 and SP1.

Conclusions

Our results strongly suggest the ELOVL6:c.-394G > A polymorphism as the causal mutation for the QTL on pig chromosome 8 that affects fatty acid composition in pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0111-y) contains supplementary material, which is available to authorized users.     

SARS-CoV-2 infection is associated with a pro-thrombotic platelet phenotype

Novel coronavirus disease 2019 (COVID-19) is associated with a hypercoagulable state, characterized by abnormal coagulation parameters and by increased incidence of cardiovascular complications. With this study, we aimed to investigate the activation state and the expression of transmembrane proteins in platelets of hospitalized COVID-19 patients. We investigated transmembrane proteins expression with a customized mass cytometry panel of 21 antibodies. Platelets of 8 hospitalized COVID-19 patients not requiring intensive care support and without pre-existing conditions were compared to platelets of healthy controls (11 donors) with and without in vitro stimulation with thrombin receptor-activating peptide (TRAP). Mass cytometry of non-stimulated platelets detected an increased surface expression of activation markers P-Selectin (0.67 vs. 1.87 median signal intensity for controls vs. patients, p = 0.0015) and LAMP-3 (CD63, 0.37 vs. 0.81, p = 0.0004), the GPIIb/IIIa complex (4.58 vs. 5.03, p < 0.0001) and other adhesion molecules involved in platelet activation and platelet–leukocyte interactions. Upon TRAP stimulation, mass cytometry detected a higher expression of P-selectin in COVID-19 samples compared to controls (p < 0.0001). However, we observed a significantly reduced capacity of COVID-19 platelets to increase the expression of activation markers LAMP-3 and P-Selectin upon stimulation with TRAP. We detected a hyperactivated phenotype in platelets during SARS-CoV-2 infection, consisting of highly expressed platelet activation markers, which might contribute to the hypercoagulopathy observed in COVID-19. In addition, several transmembrane proteins were more highly expressed compared to healthy controls. These findings support research projects investigating antithrombotic and antiplatelet treatment regimes in COVID-19 patients, and provide new insights on the phenotypical platelet expression during SARS-CoV-2 infection.Subject terms: Mechanisms of disease, Viral infection     

The cell cycle- and insulin-signaling-inhibiting miRNA expression pattern of very small embryonic-like stem cells contributes to their quiescent state

Murine Oct4⁺, very small embryonic-like stem cells (VSELs), are a quiescent stem cell population that requires a supportive co-culture layer to proliferate and/or to differentiate in vitro. Gene expression studies have revealed that the quiescence of these cells is due to changes in expression of parentally imprinted genes, including genes involved in cell cycle regulation and insulin and insulin-like growth factor signaling (IIS). To investigate the role of microRNAs (miRNAs) in VSEL quiescence, we performed miRNA studies in highly purified VSELs and observed a unique miRNA expression pattern in these cells. Specifically, we observed significant differences in the expression of certain miRNA species (relative to a reference cell population), including (i) miRNA-25_1 and miRNA-19 b, whose downregulation has the effect of upregulating cell cycle checkpoint genes and (ii) miRNA-675-3 p and miRNA-675-5 p, miRNA-292-5 p, miRNA-184, and miRNA-125 b, whose upregulation attenuates IIS. These observations are important for understanding the biology of these cells and for developing efficient ex vivo expansion strategies for VSELs isolated from adult tissues.     

Synovial membrane protein expression differs between juvenile idiopathic arthritis subtypes in early disease

Introduction

Juvenile idiopathic arthritis (JIA) is the most common rheumatological disease of childhood with a prevalence of around 1 in 1,000. Without appropriate treatment it can have devastating consequences including permanent disability from joint destruction and growth deformities. Disease aetiology remains unknown. Investigation of disease pathology at the level of the synovial membrane is required if we want to begin to understand the disease at the molecular and biochemical level. The synovial membrane proteome from early disease-stage, treatment naive JIA patients was compared between polyarticular and oligoarticular subgroups.

Methods

Protein was extracted from 15 newly diagnosed, treatment naive JIA synovial membrane biopsies and separated by two dimensional fluorescent difference in-gel electrophoresis. Proteins displaying a two-fold or greater change in expression levels between the two subgroups were identified by matrix assisted laser desorption ionization-time of flight mass spectrometry with expression further verified by Western blotting and immunohistochemistry.

Results

Analysis of variance analysis (P ≤ 0.05) revealed 25 protein spots with a two-fold or greater difference in expression levels between polyarticular and oligoarticular patients. Hierarchical cluster analysis with Pearson ranked correlation revealed two distinctive clusters of proteins. Some of the proteins that were differentially expressed included: integrin alpha 2b (P = 0.04); fibrinogen D fragment (P = 0.005); collagen type VI (P = 0.03); fibrinogen gamma chain (P = 0.05) and peroxiredoxin 2 (P = 0.02). The identified proteins are involved in a number of different processes including platelet activation and the coagulation system.

Conclusions

The data indicate distinct synovial membrane proteome profiles between JIA subgroups at an early stage in the disease process. The identified proteins also provide insight into differentially perturbed pathways which could influence pathological events at the joint level.     

Waves of dispersal in island-hopping Chondrina species (Gastropoda, Pulmonata, Chondrinidae)

The aim of this study was to unravel the historical biogeography of the speciose land snail genus Chondrina. To this end phylogenetic hypotheses were tested using mitochondrial DNA sequence data.Mitochondrial DNA sequences of the Cytochrome Oxidase subunit I region were obtained for 89 individuals, representing just over 70% of the extant Chondrina species. The extent of molecular genetic diversity and phylogeographical patterns were investigated by using neighbour joining, parsimony and bayesian methods for phylogeny reconstructions. The resulting data were used to infer historical biogeographical patterns for the genus Chondrina.The three phylogenetic methods yielded congruent topologies for the phylogeny reconstruction. Six clades were identified, each of which with at least one taxon that is known from the Iberian peninsula. The most parsimonious scenario indicates at least three waves of dispersal out of the Iberian peninsula into the North and East of Europe and Northern Africa.The phylogenetic relationships combined with the distributional patterns of the various species, indicate that only vicariance events cannot explain the actual situation. Apparently, separate waves of dispersal and subsequent speciation occurred, each time starting from the southwestern part of the present generic range. Until recently, this was obscured by repetitive cases of parallel or convergent evolution in shell characters, as became evident with the use of molecular methods.     

Incipient radiation within the dominant Hawaiian tree Metrosideros polymorpha

Although trees comprise a primary component of terrestrial species richness, the drivers and temporal scale of divergence in trees remain poorly understood. We examined the landscape-dominant tree, Metrosideros polymorpha, for variation at nine microsatellite loci across 23 populations on young Hawai''i Island, sampling each of the island''s five varieties throughout its full geographic range. For four varieties, principal coordinate analysis revealed strong clustering of populations by variety across the 10 430 km² island, indicating partitioning of the species into multiple evolutionarily significant units. The single island-endemic form, riparian var. newellii, showed especially strong differentiation from other varieties despite occurring in sympatry with other varieties and likely evolved from a bog form on the oldest volcano, Kohala, within the past 500 000 years. Along with comparable riparian forms on other Pacific Islands, var. newellii appears to represent parallel incipient ecological speciation within Metrosideros. Greater genetic distance among the more common varieties on the oldest volcano and an inverse relationship between allelic diversity and substrate age appear consistent with colonization of Hawai''i Island by older, partially diverged varieties followed by increased hybridization among varieties on younger volcanoes. This study demonstrates that broad population-level sampling is required to uncover patterns of diversification within a ubiquitous and long-lived tree species. Hawaiian Metrosideros appears to be a case of incipient radiation in trees and thus should be useful for studies of divergence and the evolution of reproductive isolating barriers at the early stages of speciation.     

Inferring historical survivals of climate relicts: the effects of climate changes,geography, and population-specific factors on herbaceous hydrangeas

Climate relicts hold considerable importance because they have resulted from numerous historical changes. However, there are major interspecific variations among the ways by which they survived climate changes. Therefore, investigating the factors and timing that affected population demographics can expand our understanding of how climate relicts responded to historical environmental changes. Here, we examined herbaceous hydrangeas of genus Deinanthe in East Asia, which show limited distributions and a remarkable disjunction between Japan and central China. Chloroplast genome and restriction site-associated DNA sequencing revealed that speciation event occurred in the late Miocene (ca. 7–9 Mya) in response to global climate change. Two lineages apparently remained not branched until the middle Quaternary, and afterwards started to diverge to regional population groups. The narrow endemic species in central China showed lower genetic diversity (He = 0.082), as its population size rapidly decreased during the Holocene due to isolation in montane refugia. Insular populations in the three Japanese islands (He = 0.137–0.160) showed a genetic structure that was inconsistent with sea barriers, indicating that it was shaped in the glacial period when its range retreated to coastal refugia on the exposed sea floor. Demographic modelling by stairway-plot analysis reconstructed variable responses of Japanese populations: some experienced glacial bottlenecks in refugial isolation, while post-glacial range expansion seemingly exerted founder effects on other populations. Overall, this study demonstrated the involvement of not just one, but multiple factors, such as the interplay between climate changes, geography, and other population-specific factors, that determine the demographics of climate relicts.Subject terms: Population genetics, Plant genetics     

Historical biogeography of the genus Rhadinaea (Squamata: Dipsadinae)

Multiple geological and climatic events have created geographical or ecological barriers associated with speciation events, playing a role in biological diversification in North and Central America. Here, we evaluate the influence of the Neogene and Quaternary geological events, as well as the climatic changes in the diversification of the colubrid snake genus Rhadinaea using molecular dating and ancestral area reconstruction. A multilocus sequence dataset was generated for 37 individuals of Rhadinaea from most of the biogeographical provinces where the genus is distributed, representing 19 of the 21 currently recognized species, and two undescribed species. Our analyses show that the majority of the Rhadinaea species nest in two main clades, herein identified as “Eastern” and “Southern”. These clades probably diverged from each other in the early Miocene, and their divergence was followed by 11 divergences during the middle to late Miocene, three divergences during the Pliocene, and six divergences in the Pleistocene. The ancestral distribution of Rhadinaea was reconstructed across the Sierra Madre del Sur. Our phylogenetic analyses do not support the monophyly of Rhadinaea. The Miocene and Pliocene geomorphology, perhaps in conjunction with climate change, appears to have triggered the diversification of the genus, while the climatic changes during the Miocene probably induced the diversification of Rhadinaea in the Sierra Madre del Sur. Our analysis suggests that the uplifting of the Trans‐Mexican Volcanic Belt and Chiapan–Guatemalan highlands in this same period resulted in northward and southward colonization events. This was followed by more recent, independent colonization events in the Pliocene and Pleistocene involving the Balsas Basin, Chihuahuan Desert, Pacific Coast, Sierra Madre Occidental, Sierra Madre Oriental, Sierra Madre del Sur, Trans‐Mexican Volcanic Belt, and Veracruz provinces, probably driven by the climatic fluctuations of the time.     

Genomic Signatures for Species-Specific Adaptation in Lake Victoria Cichlids Derived from Large-Scale Standing Genetic Variation

The cichlids of Lake Victoria are a textbook example of adaptive radiation, as >500 endemic species arose in just 14,600 years. The degree of genetic differentiation among species is very low due to the short period of time after the radiation, which allows us to ascertain highly differentiated genes that are strong candidates for driving speciation and adaptation. Previous studies have revealed the critical contribution of vision to speciation by showing the existence of highly differentiated alleles in the visual opsin gene among species with different habitat depths. In contrast, the processes of species-specific adaptation to different ecological backgrounds remain to be investigated. Here, we used genome-wide comparative analyses of three species of Lake Victoria cichlids that inhabit different environments—Haplochromis chilotes, H. sauvagei, and Lithochromis rufus—to elucidate the processes of adaptation by estimating population history and by searching for candidate genes that contribute to adaptation. The patterns of changes in population size were quite distinct among the species according to their habitats. We identified many novel adaptive candidate genes, some of which had surprisingly long divergent haplotypes between species, thus showing the footprint of selective sweep events. Molecular phylogenetic analyses revealed that a large fraction of the allelic diversity among Lake Victoria cichlids was derived from standing genetic variation that originated before the adaptive radiation. Our analyses uncovered the processes of species-specific adaptation of Lake Victoria cichlids and the complexity of the genomic substrate that facilitated this adaptation.