Hybridization between distant lineages increases adaptive variation during a biological invasion: stickleback in Switzerland

The three-spined stickleback is a widespread Holarctic species complex that radiated from the sea into freshwaters after the retreat of the Pleistocene ice sheets. In Switzerland, sticklebacks were absent with the exception of the far northwest, but different introduced populations have expanded to occupy a wide range of habitats since the late 19th century. A well-studied adaptive phenotypic trait in sticklebacks is the number of lateral plates. With few exceptions, freshwater and marine populations in Europe are fixed for either the low plated phenotype or the fully plated phenotype, respectively. Switzerland, in contrast, harbours in close proximity the full range of phenotypic variation known from across the continent. We addressed the phylogeographic origins of Swiss sticklebacks using mitochondrial partial cytochrome b and control region sequences. We found only five different haplotypes but these originated from three distinct European regions, fixed for different plate phenotypes. These lineages occur largely in isolation at opposite ends of Switzerland, but co-occur in a large central part. Across the country, we found a strong correlation between a microsatellite linked to the high plate ectodysplasin allele and the mitochondrial haplotype from a region where the fully plated phenotype is fixed. Phylogenomic and population genomic analysis of 481 polymorphic amplified fragment length polymorphism loci indicate genetic admixture in the central part of the country. The same part of the country also carries elevated within-population phenotypic variation. We conclude that during the recent invasive range expansion of sticklebacks in Switzerland, adaptive and neutral between-population genetic variation was converted into within-population variation, raising the possibility that hybridization between colonizing lineages contributed to the ecological success of sticklebacks in Switzerland.

Regulatory sequences of the porcine THBD gene facilitate endothelial-specific expression of bioactive human thrombomodulin in single- and multitransgenic pigs.

BACKGROUND Among other mismatches between human and pig, incompatibilities in the blood coagulation systems hamper the xenotransplantation of vascularized organs. The provision of the porcine endothelium with human thrombomodulin (hTM) is hypothesized to overcome the impaired activation of protein C by a heterodimer consisting of human thrombin and porcine TM. METHODS We evaluated regulatory regions of the THBD gene, optimized vectors for transgene expression, and generated hTM expressing pigs by somatic cell nuclear transfer. Genetically modified pigs were characterized at the molecular, cellular, histological, and physiological levels. RESULTS A 7.6-kb fragment containing the entire upstream region of the porcine THBD gene was found to drive a high expression in a porcine endothelial cell line and was therefore used to control hTM expression in transgenic pigs. The abundance of hTM was restricted to the endothelium, according to the predicted pattern, and the transgene expression of hTM was stably inherited to the offspring. When endothelial cells from pigs carrying the hTM transgene--either alone or in combination with an aGalTKO and a transgene encoding the human CD46-were tested in a coagulation assay with human whole blood, the clotting time was increased three- to four-fold (P<0.001) compared to wild-type and aGalTKO/CD46 transgenic endothelial cells. This, for the first time, demonstrated the anticoagulant properties of hTM on porcine endothelial cells in a human whole blood assay. CONCLUSIONS The biological efficacy of hTM suggests that the (multi-)transgenic donor pigs described here have the potential to overcome coagulation incompatibilities in pig-to-primate xenotransplantation.

Gene sequences of the lichen Cetraria aculeata

Lichens are symbioses between fungi (mycobionts) and photoautotrophic green algae or cyanobacteria (photobionts). Many lichens occupy large distributional ranges covering several climatic zones. So far, little is known about the large-scale phylogeography of lichen photobionts and their role in shaping the distributional ranges of lichens. We studied south polar, temperate and north polar populations of the widely distributed fruticose lichen Cetraria aculeata. Based on the DNA sequences from three loci for each symbiont, we compared the genetic structure of mycobionts and photobionts. Phylogenetic reconstructions and Bayesian clustering methods divided the mycobiont and photobiont data sets into three groups. An AMOVA shows that the genetic variance of the photobiont is best explained by differentiation between temperate and polar regions and that of the mycobiont by an interaction of climatic and geographical factors. By partialling out the relative contribution of climate, geography and codispersal, we found that the most relevant factors shaping the genetic structure of the photobiont are climate and a history of codispersal. Mycobionts in the temperate region are consistently associated with a specific photobiont lineage. We therefore conclude that a photobiont switch in the past enabled C. aculeata to colonize temperate as well as polar habitats. Rare photobiont switches may increase the geographical range and ecological niche of lichen mycobionts by associating them with locally adapted photobionts in climatically different regions and, together with isolation by distance, may lead to genetic isolation between populations and thus drive the evolution of lichens.

DNA alignment and Bayesian trees of partial sequences of COI, COIII and the haemocyanin functional unit f-g of 28 octopod species

Background: Octopods have successfully colonised the world's oceans from the tropics to the poles. Yet, successful persistence in these habitats has required adaptations of their advanced physiological apparatus to compensate impaired oxygen supply. Their oxygen transporter haemocyanin plays a major role in cold tolerance and accordingly has undergone functional modifications to sustain oxygen release at sub-zero temperatures. However, it remains unknown how molecular properties evolved to explain the observed functional adaptations. We thus aimed to assess whether natural selection affected molecular and structural properties of haemocyanin that explains temperature adaptation in octopods. Results: Analysis of 239 partial sequences of the haemocyanin functional units (FU) f and g of 28 octopod species of polar, temperate, subtropical and tropical origin revealed natural selection was acting primarily on charge properties of surface residues. Polar octopods contained haemocyanins with higher net surface charge due to decreased glutamic acid content and higher numbers of basic amino acids. Within the analysed partial sequences, positive selection was present at site 2545, positioned between the active copper binding centre and the FU g surface. At this site, methionine was the dominant amino acid in polar octopods and leucine was dominant in tropical octopods. Sites directly involved in oxygen binding or quaternary interactions were highly conserved within the analysed sequence. Conclusions: This study has provided the first insight into molecular and structural mechanisms that have enabled octopods to sustain oxygen supply from polar to tropical conditions. Our findings imply modulation of oxygen binding via charge-charge interaction at the protein surface, which stabilize quaternary interactions among functional units to reduce detrimental effects of high pH on venous oxygen release. Of the observed partial haemocyanin sequence, residue 2545 formed a close link between the FU g surface and the active centre, suggesting a role as allosteric binding site. The prevalence of methionine at this site in polar octopods, implies regulation of oxygen affinity via increased sensitivity to allosteric metal binding. High sequence conservation of sites directly involved in oxygen binding indicates that functional modifications of octopod haemocyanin rather occur via more subtle mechanisms, as observed in this study.

Metadata und statistic analysis of archael and bacterial sequences originating from sediments of the Håkon Mosby mud volcano (Z1-Z3, all habitats)

DNA extraction was carried out as described on the MICROBIS project pages (http://icomm.mbl.edu/microbis ) using a commercially available extraction kit. We amplified the hypervariable regions V4-V6 of archaeal and bacterial 16S rRNA genes using PCR and several sets of forward and reverse primers (http://vamps.mbl.edu/resources/primers.php). Massively parallel tag sequencing of the PCR products was carried out on a 454 Life Sciences GS FLX sequencer at Marine Biological Laboratory, Woods Hole, MA, following the same experimental conditions for all samples. Sequence reads were submitted to a rigorous quality control procedure based on mothur v30 (doi:10.1128/AEM.01541-09) including denoising of the flow grams using an algorithm based on PyroNoise (doi:10.1038/nmeth.1361), removal of PCR errors and a chimera check using uchime (doi:10.1093/bioinformatics/btr381). The reads were taxonomically assigned according to the SILVA taxonomy (SSURef v119, 07-2014; doi:10.1093/nar/gks1219) implemented in mothur and clustered at 98% ribosomal RNA gene V4-V6 sequence identity. V4-V6 amplicon sequence abundance tables were standardized to account for unequal sampling effort using 1000 (Archaea) and 2300 (Bacteria) randomly chosen sequences without replacement using mothur and then used to calculate inverse Simpson diversity indices and Chao1 richness (doi:10.2307/4615964). Bray-Curtis dissimilarities (doi:10.2307/1942268) between all samples were calculated and used for 2-dimensional non metric multidimensional scaling (NMDS) ordinations with 20 random starts (doi:10.1007/BF02289694). Stress values below 0.2 indicated that the multidimensional dataset was well represented by the 2D ordination. NMDS ordinations were compared and tested using Procrustes correlation analysis (doi:10.1007/BF02291478). All analyses were carried out with the R statistical environment and the packages vegan (available at: http://cran.r-project.org/package=vegan), labdsv (available at: http://cran.r-project.org/package=labdsv), as well as with custom R scripts. Operational taxonomic units at 98% sequence identity (OTU0.03) that occurred only once in the whole dataset were termed absolute single sequence OTUs (SSOabs; doi:10.1038/ismej.2011.132). OTU0.03 sequences that occurred only once in at least one sample, but may occur more often in other samples were termed relative single sequence OTUs (SSOrel). SSOrel are particularly interesting for community ecology, since they comprise rare organisms that might become abundant when conditions change.16S rRNA amplicons and metagenomic reads have been stored in the sequence read archive under SRA project accession number SRP042162.

Availability of sequences of aquatic and terrestrial taxa in genetic databases (GenBank and BOLD - the Barcode of Life Database) in 2010, 2012 and 2016

Correct species identifications are of tremendous importance for invasion ecology, as mistakes could lead to misdirecting limited resources against harmless species or inaction against problematic ones. DNA barcoding is becoming a promising and reliable tool for species identifications, however the efficacy of such molecular taxonomy depends on gene region(s) that provide a unique sequence to differentiate among species and on availability of reference sequences in existing genetic databases. Here, we assembled a list of aquatic and terrestrial non-indigenous species (NIS) and checked two leading genetic databases for corresponding sequences of six genome regions used for DNA barcoding. The genetic databases were checked in 2010, 2012, and 2016. All four aquatic kingdoms (Animalia, Chromista, Plantae and Protozoa) were initially equally represented in the genetic databases, with 64, 65, 69, and 61% of NIS included, respectively. Sequences for terrestrial NIS were present at rates of 58 and 78% for Animalia and Plantae, respectively. Six years later, the number of sequences for aquatic NIS increased to 75, 75, 74, and 63% respectively, while those for terrestrial NIS increased to 74 and 88% respectively. Genetic databases are marginally better populated with sequences of terrestrial NIS of plants compared to aquatic NIS and terrestrial NIS of animals. The rate at which sequences are added to databases is not equal among taxa. Though some groups of NIS are not detectable at all based on available data - mostly aquatic ones - encouragingly, current availability of sequences of taxa with environmental and/or economic impact is relatively good and continues to increase with time.

Scientific drilling projects in ancient lakes: Integrating geological and biological histories

Sedimentary sequences in ancient or long-lived lakes can reach several thousands of meters in thickness and often provide an unrivalled perspective of the lake's regional climatic, environmental, and biological history. Over the last few years, deep-drilling projects in ancient lakes became increasingly multi- and interdisciplinary, as, among others, seismological, sedimentological, biogeochemical, climatic, environmental, paleontological, and evolutionary information can be obtained from sediment cores. However, these multi- and interdisciplinary projects pose several challenges. The scientists involved typically approach problems from different scientific perspectives and backgrounds, and setting up the program requires clear communication and the alignment of interests. One of the most challenging tasks, besides the actual drilling operation, is to link diverse datasets with varying resolution, data quality, and age uncertainties to answer interdisciplinary questions synthetically and coherently. These problems are especially relevant when secondary data, i.e., datasets obtained independently of the drilling operation, are incorporated in analyses. Nonetheless, the inclusion of secondary information, such as isotopic data from fossils found in outcrops or genetic data from extant species, may help to achieve synthetic answers. Recent technological and methodological advances in paleolimnology are likely to increase the possibilities of integrating secondary information. Some of the new approaches have started to revolutionize scientific drilling in ancient lakes, but at the same time, they also add a new layer of complexity to the generation and analysis of sediment-core data. The enhanced opportunities presented by new scientific approaches to study the paleolimnological history of these lakes, therefore, come at the expense of higher logistic, communication, and analytical efforts. Here we review types of data that can be obtained in ancient lake drilling projects and the analytical approaches that can be applied to empirically and statistically link diverse datasets to create an integrative perspective on geological and biological data. In doing so, we highlight strengths and potential weaknesses of new methods and analyses, and provide recommendations for future interdisciplinary deep-drilling projects.

Thermal adaptation analyzed by comparison of protein sequences from mesophilic and extremely thermophilic Methanococcus species

The genome sequence of the extremely thermophilic archaeon Methanococcus jannaschii provides a wealth of data on proteins from a thermophile. In this paper, sequences of 115 proteins from M. jannaschii are compared with their homologs from mesophilic Methanococcus species. Although the growth temperatures of the mesophiles are about 50°C below that of M. jannaschii, their genomic G+C contents are nearly identical. The properties most correlated with the proteins of the thermophile include higher residue volume, higher residue hydrophobicity, more charged amino acids (especially Glu, Arg, and Lys), and fewer uncharged polar residues (Ser, Thr, Asn, and Gln). These are recurring themes, with all trends applying to 83–92% of the proteins for which complete sequences were available. Nearly all of the amino acid replacements most significantly correlated with the temperature change are the same relatively conservative changes observed in all proteins, but in the case of the mesophile/thermophile comparison there is a directional bias. We identify 26 specific pairs of amino acids with a statistically significant (P < 0.01) preferred direction of replacement.

Effective inhibition of influenza virus production in cultured cells by external guide sequences and ribonuclease P

The polymerase (PB2) and nucleocapsid (NP) genes encoded by the genome of influenza virus are essential for replication of the virus. When synthetic genes that express RNAs for external guide sequences targeted to the mRNAs of the PB2 and NP genes are stably incorporated into mouse cells in tissue culture, infection of these cells with influenza virus is nonproductive. Endogenous RNase P cleaves the targeted influenza virus mRNAs when they are in a complex with the external guide sequences. Targeting two different mRNAs simultaneously inhibits viral particle production more efficiently than does targeting only one mRNA.

The ability of HIV type 1 to use CCR-3 as a coreceptor is controlled by envelope V1/V2 sequences acting in conjunction with a CCR-5 tropic V3 loop

Although infection by primary HIV type 1 (HIV-1) isolates normally requires the functional interaction of the viral envelope protein with both CD4 and the CCR-5 coreceptor, a subset of such isolates also are able to use the distinct CCR-3 receptor. By analyzing the ability of a series of wild-type and chimeric HIV-1 envelope proteins to mediate CCR-3-dependent infection, we have determined that CCR-3 tropism maps to the V1 and V2 variable region of envelope. Although substitution of the V1/V2 region of a CCR-3 tropic envelope into the context of a CCR-5 tropic envelope is both necessary and sufficient to confer CCR-3 tropism, this same substitution has no phenotypic effect when inserted into a CXCR-4 tropic HIV-1 envelope context. However, this latter chimera acquires both CCR-3 and CCR-5 tropism when a CCR-5 tropic V3 loop sequence also is introduced. These data demonstrate that the V1/2 region of envelope can, like the V3 loop region, encode a particular coreceptor requirement and suggest that a functional envelope:CCR-3 interaction may depend on the cooperative interaction of CCR-3 with both the V1/V2 and the V3 region of envelope.

Direct isolation of human transcribed sequences from yeast artificial chromosomes through the application of RNA fingerprinting

The identification of cDNA clones from genomic regions known to contain human genes is usually the rate-limiting factor in positional cloning strategies. We demonstrate here that human genes present on yeast artificial chromosomes (YACs) are transcribed in yeast host cells. We have used the arbitrarily primed RNA (RAP) fingerprinting method to identify human-specific, transcribed sequences from YACs located in the 13q12 chromosome region. By comparing the RAP fingerprints generated using defined, arbitrary primers from various fragmented YACs, megaYACs, and host yeast, we were able to identify and map 20 products transcribed from the human YAC inserts. This method, therefore, permits the simultaneous isolation and mapping of novel expressed sequences directly from whole YACs.

Two sterol regulatory element-like sequences mediate up-regulation of caveolin gene transcription in response to low density lipoprotein free cholesterol

Caveolae form the terminus for a major pathway of intracellular free cholesterol (FC) transport. Caveolin mRNA levels in confluent human skin fibroblasts were up-regulated following increased uptake of low density lipoprotein (LDL) FC. The increase induced by FC was not associated with detectable change in mRNA stability, indicating that caveolin mRNA levels were mediated at the level of gene transcription. A total of 924 bp of 5′ flanking region of the caveolin gene were cloned and sequenced. The promoter sequence included three G+C-rich potential sterol regulatory elements (SREs), a CAAT sequence and a Sp1 consensus sequence. Deletional mutagenesis of individual SRE-like sequences indicated that of these two (at −646 and −395 bp) were essential for the increased transcription rates mediated by LDL-FC, whereas the third was inconsequential. Gel shift analysis of protein binding from nuclear extracts to these caveolin promoter DNA sequences, together with DNase I footprinting, confirmed nucleoprotein binding to the SRE-like elements as part of the transcriptional response to LDL-FC. A supershift obtained with antibody to SRE-binding protein 1 (SPEBP-1) indicated that this protein binds at −395 bp. There was no reaction at −395 bp with anti-Sp1 antibody nor with either antibody at −646 bp. The cysteine protease inhibitor N-acetyl-leu-leu-norleucinal (ALLN), which inhibits SREBP catabolism, superinhibited caveolin mRNA levels regardless of LDL-FC. This finding suggests that SREBP inhibits caveolin gene transcription in contrast to its stimulating effect on other promoters. The findings of this study are consistent with the postulated role for caveolin as a regulator of cellular FC homeostasis in quiescent peripheral cells, and the coordinate regulation by SREBP of FC influx and efflux.

A peptide derived from an extracellular domain selectively inhibits receptor internalization: Target sequences on insulin and insulin-like growth factor 1 receptors

Certain peptides derived from the α1 domain of the major histocompatibility class I antigen complex (MHC-I) inhibit receptor internalization, increasing the steady-state number of active receptors on the cell surface and thereby enhancing the sensitivity to hormones and other agonists. These peptides self-assemble, and they also bind to MHC-I at the same site from which they are derived, suggesting that they could bind to receptor sites with significant sequence similarity. Receptors affected by MHC-I peptides do, indeed, have such sequence similarity, as illustrated here by insulin receptor (IR) and insulin-like growth factor-1 receptor. A synthetic peptide with sequence identical to a certain extracellular receptor domain binds to that receptor in a ligand-dependent manner and inhibits receptor internalization. Moreover, each such peptide is selective for its cognate receptor. An antibody to the IR peptide not only binds to IR and competes with the peptide but also inhibits insulin-dependent internalization of IR. These observations, and binding studies with deletion mutants of IR, indicate that the sequence QILKELEESSF encoded by exon 10 plays a key role in IR internalization. Our results illustrate a principle for identifying receptor-specific sites of importance for receptor internalization, and for enhancing sensitivity to hormones and other agonists.

A structural census of the current population of protein sequences

We examine the occurrence of the ≈300 known protein folds in different groups of organisms. To do this, we characterize a large fraction of the currently known protein sequences (≈140,000) in structural terms, by matching them to known structures via sequence comparison (or by secondary-structure class prediction for those without structural homologues). Overall, we find that an appreciable fraction of the known folds are present in each of the major groups of organisms (e.g., bacteria and eukaryotes share 156 of 275 folds), and most of the common folds are associated with many families of nonhomologous sequences (i.e., >10 sequence families for each common fold). However, different groups of organisms have characteristically distinct distributions of folds. So, for instance, some of the most common folds in vertebrates, such as globins or zinc fingers, are rare or absent in bacteria. Many of these differences in fold usage are biologically reasonable, such as the folds of metabolic enzymes being common in bacteria and those associated with extracellular transport and communication being common in animals. They also have important implications for database-based methods for fold recognition, suggesting that an unknown sequence from a plant is more likely to have a certain fold (e.g., a TIM barrel) than an unknown sequence from an animal.

Evolution of gilled mushrooms and puffballs inferred from ribosomal DNA sequences

Homobasidiomycete fungi display many complex fruiting body morphologies, including mushrooms and puffballs, but their anatomical simplicity has confounded efforts to understand the evolution of these forms. We performed a comprehensive phylogenetic analysis of homobasidiomycetes, using sequences from nuclear and mitochondrial ribosomal DNA, with an emphasis on understanding evolutionary relationships of gilled mushrooms and puffballs. Parsimony-based optimization of character states on our phylogenetic trees suggested that strikingly similar gilled mushrooms evolved at least six times, from morphologically diverse precursors. Approximately 87% of gilled mushrooms are in a single lineage, which we call the “euagarics.” Recently discovered 90 million-year-old fossil mushrooms are probably euagarics, suggesting that (i) the origin of this clade must have occurred no later than the mid-Cretaceous and (ii) the gilled mushroom morphology has been maintained in certain lineages for tens of millions of years. Puffballs and other forms with enclosed spore-bearing structures (Gasteromycetes) evolved at least four times. Derivation of Gasteromycetes from forms with exposed spore-bearing structures (Hymenomycetes) is correlated with repeated loss of forcible spore discharge (ballistospory). Diverse fruiting body forms and spore dispersal mechanisms have evolved among Gasteromycetes. Nevertheless, it appears that Hymenomycetes have never been secondarily derived from Gasteromycetes, which suggests that the loss of ballistospory has constrained evolution in these lineages.