Applications of next generation sequencing in molecular ecology of non-model organisms

As most biologists are probably aware, technological advances in molecular biology during the last few years have opened up possibilities to rapidly generate large-scale sequencing data from non-model organisms at a reasonable cost. In an era when virtually any study organism can 'go genomic', it is worthwhile to review how this may impact molecular ecology. The first studies to put the next generation sequencing (NGS) to the test in ecologically well-characterized species without previous genome information were published in 2007 and the beginning of 2008. Since then several studies have followed in their footsteps, and a large number are undoubtedly under way. This review focuses on how NGS has been, and can be, applied to ecological, population genetic and conservation genetic studies of non-model species, in which there is no (or very limited) genomic resources. Our aim is to draw attention to the various possibilities that are opening up using the new technologies, but we also highlight some of the pitfalls and drawbacks with these methods. We will try to provide a snapshot of the current state of the art for this rapidly advancing and expanding field of research and give some likely directions for future developments.     

The utility of PacBio circular consensus sequencing for characterizing complex gene families in non-model organisms

Background

Molecular characterization of highly diverse gene families can be time consuming, expensive, and difficult, especially when considering the potential for relatively large numbers of paralogs and/or pseudogenes. Here we investigate the utility of Pacific Biosciences single molecule real-time (SMRT) circular consensus sequencing (CCS) as an alternative to traditional cloning and Sanger sequencing PCR amplicons for gene family characterization. We target vomeronasal gene receptors, one of the most diverse gene families in mammals, with the goal of better understanding intra-specific V1R diversity of the gray mouse lemur (Microcebus murinus). Our study compares intragenomic variation for two V1R subfamilies found in the mouse lemur. Specifically, we compare gene copy variation within and between two individuals of M. murinus as characterized by different methods for nucleotide sequencing. By including the same individual animal from which the M. murinus draft genome was derived, we are able to cross-validate gene copy estimates from Sanger sequencing versus CCS methods.

Results

We generated 34,088 high quality circular consensus sequences of two diverse V1R subfamilies (here referred to as V1RI and V1RIX) from two individuals of Microcebus murinus. Using a minimum threshold of 7× coverage, we recovered approximately 90% of V1RI sequences previously identified in the draft M. murinus genome (59% being identical at all nucleotide positions). When low coverage sequences were considered (i.e. < 7× coverage) 100% of V1RI sequences identified in the draft genome were recovered. At least 13 putatively novel V1R loci were also identified using CCS technology.

Conclusions

Recent upgrades to the Pacific Biosciences RS instrument have improved the CCS technology and offer an alternative to traditional sequencing approaches. Our results suggest that the Microcebus murinus V1R repertoire has been underestimated in the draft genome. In addition to providing an improved understanding of V1R diversity in the mouse lemur, this study demonstrates the utility of CCS technology for characterizing complex regions of the genome. We anticipate that long-read sequencing technologies such as PacBio SMRT will allow for the assembly of multigene family clusters and serve to more accurately characterize patterns of gene copy variation in large gene families, thus revealing novel micro-evolutionary patterns within non-model organisms.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-720) contains supplementary material, which is available to authorized users.     

Chemical genetics of zipper-interacting protein kinase reveal myosin light chain as a bona fide substrate in permeabilized arterial smooth muscle

Zipper-interacting protein kinase (ZIPK) has been implicated in Ca(2+)-independent smooth muscle contraction, although its specific role is unknown. The addition of ZIPK to demembranated rat caudal arterial strips induced an increase in force, which correlated with increases in LC(20) and MYPT1 phosphorylation. However, because of the number of kinases capable of phosphorylating LC(20) and MYPT1, it has proven difficult to identify the mechanism underlying ZIPK action. Therefore, we set out to identify bona fide ZIPK substrates using a chemical genetics method that takes advantage of ATP analogs with bulky substituents at the N(6) position and an engineered ZIPK capable of utilizing such substrates. (32)P-Labeled 6-phenyl-ATP and ZIPK-L93G mutant protein were added to permeabilized rat caudal arterial strips, and substrate proteins were detected by autoradiography following SDS-PAGE. Mass spectrometry identified LC(20) as a direct target of ZIPK in situ for the first time. Tissues were also exposed to 6-phenyl-ATP and ZIPK-L93G in the absence of endogenous ATP, and putative ZIPK substrates were identified by Western blotting. LC(20) was thereby confirmed as a direct target of ZIPK; however, no phosphorylation of MYPT1 was detected. We conclude that ZIPK is involved in the regulation of smooth muscle contraction through direct phosphorylation of LC(20).     

基于RNA-Seq转录组测序技术揭示肉兔生长发育调控的分子机制

为研究肉兔肌肉生长发育调控的分子机制,本试验以84日龄齐卡巨型白兔和齐兴肉兔为研究对象,屠宰后取背最长肌组织并提取总RNA,反转录建立c DNA文库,利用Illumina HiSeq 2500测序平台进行测序,并进行序列分析和注释。筛选与肉兔肌肉生长发育相关的信号通路及差异表达基因,最后进行荧光定量PCR验证。结果表明:齐卡巨型白兔和齐兴肉兔背最长肌中显著差异表达基因总数为833个,其中齐卡巨型白兔中显著上调基因325个,显著下调基因508个。KEGG功能注释发现差异基因显著富集到PI3K-AKT通路、癌症通路和黏着斑通路,最终筛选出4个可能与生长发育显著相关的差异表达基因,为进一步对肉兔进行遗传改良提供了理论基础。     

Mouse bestrophin-2 is a bona fide Cl(-) channel: identification of a residue important in anion binding and conduction

Bestrophins have recently been proposed to comprise a new family of Cl(-) channels. Our goal was to test whether mouse bestrophin-2 (mBest2) is a bona fide Cl(-) channel. We expressed mBest2 in three different mammalian cell lines. mBest2 was trafficked to the plasma membrane as shown by biotinylation and immunoprecipitation, and induced a Ca(2+)-activated Cl(-) current in all three cell lines (EC(50) for Ca(2+) = 230 nM). The permeability sequence was SCN(-): I(-): Br(-): Cl(-): F(-) (8.2: 1.9: 1.4: 1: 0.5). Although SCN(-) was highly permeant, its conductance was approximately 10% that of Cl(-) and SCN(-) blocked Cl(-) conductance (IC(50) = 12 mM). Therefore, SCN(-) entered the pore more easily than Cl(-), but bound more tightly than Cl(-). Mutations in S79 altered the relative permeability and conductance for SCN(-) as expected if S79 contributed to an anion binding site in the channel. P(SCN)/P(Cl) = 8.2 +/- 1.3 for wild-type and 3.9 +/- 0.4 for S79C. G(SCN)/G(Cl) = 0.14 +/- 0.03 for wild-type and 0.94 +/- 0.04 for S79C. In the S79 mutants, SCN(-) did not block Cl(-) conductance. This suggested that the S79C mutation altered the affinity of an anion binding site for SCN(-). Additional evidence that S79 was located in the conduction pathway was provided by the finding that modification of the sulfhydryl group in S79C with MTSET(+) or MTSES(-) increased conductance significantly. Because the effect of positively and negatively charged MTS reagents was similar, electrostatic interactions between the permeant anion and the channel at this residue were probably not critical in anion selectivity. These data provide strong evidence that mBest2 forms part of the novel Cl(-) conduction pathway in mBest2-transfected cells and that S79 plays an important role in anion binding in the pore of the channel.     

Functional relevance of Gedunin as a bona fide ligand of NADPH oxidase 5 and ROS scavenger: An in silico and in vitro assessment in a hyperglycemic RBC model

Clinical evidence suggests that type 2 diabetes therapy can greatly benefit from the suppression of reactive oxygen species generation and the activation or restoration of cellular antioxidant mechanisms. In human, NADPH oxidase (NOX) is the main producer of reactive oxygen species (ROS) that supress the activity of endogenous antioxidant enzymes. In the present study, the antioxidant potential of Gedunin was studied. In silico findings reveal its strong binding affinity with NOX5 C terminal HSP90 binding site that disrupts NOX5 stability and its ability to generate ROS, leading to restoration antioxidant enzymes activities. It was found that Gedunin suppressed hyperglycaemia induced oxidative stress in an in vitro RBC model and markedly reversed glucose induced changes including haemoglobin glycosylation and lipid peroxidation. A significant restoration of activities of cellular antioxidant enzymes; superoxide dismutase, catalase and glutathione peroxidase in the presence of Gedunin revealed its ability to reduce oxidative stress. These results substantiated Gedunin as a bona fide inhibitor of human NOX5 and a ROS scavenging antioxidant with promising therapeutic attributes including its natural origin and inhibition of multiple diabetic targets.     

The primordial differentiation of tumor-specific memory CD8+ T cells as bona fide responders to PD-1/PD-L1 blockade in draining lymph nodes

1. Download : Download high-res image (226KB)
2. Download : Download full-size image

     

Inferring the colonization of a mountain range--refugia vs. nunatak survival in high alpine ground beetles

It has long been debated whether high alpine specialists survived ice ages in situ on small ice-free islands of habitat, so-called nunataks, or whether glacial survival was restricted to larger massifs de refuge at the periphery. We evaluate these alternative hypotheses in a local radiation of high alpine carabid beetles (genus Trechus) in the Orobian Alps, Northern Italy. While summits along the northern ridge of this mountain range were surrounded by the icesheet as nunataks during the last glacial maximum, southern areas remained unglaciated. We analyse a total of 1366 bp of mitochondrial (Cox1 and Cox2) data sampled from 150 individuals from twelve populations and 530 bp of nuclear (PEPCK) sequence sampled for a subset of 30 individuals. Using Bayesian inference, we estimate ancestral location states in the gene trees, which in turn are used to infer the most likely order of recolonization under a model of sequential founder events from a massif de refuge from the mitochondrial data. We test for the paraphyly expected under this model and for reciprocal monophyly predicted by a contrasting model of prolonged persistence of nunatak populations. We find that (i) only three populations are incompatible with the paraphyly of the massif de refuge model, (ii) both mitochondrial and nuclear data support separate refugial origins for populations on the western and eastern ends of the northern ridge, and (iii) mitochondrial node ages suggest persistence on the northern ridge for part of the last ice age.     

Robustness of predictions of extremely thermally stable proteins in ancient organisms

A number of studies have addressed the environmental temperatures experienced by ancient life. Computational studies using a nonhomogeneous evolution model have estimated ancestral G + C contents of ribosomal RNAs and the amino acid compositions of ancestral proteins, generating hypotheses regarding the mesophilic last universal common ancestor. In contrast, our previous study computationally reconstructed ancestral amino acid sequences of nucleoside diphosphate kinases using a homogeneous model and then empirically resurrected the ancestral proteins. The thermal stabilities of these ancestral proteins were equivalent to or greater than those of extant homologous thermophilic proteins, supporting the thermophilic universal ancestor theory. In this study, we reinferred ancestral sequences using a dataset from which hyperthermophilic sequences were excluded. We also reinferred ancestral sequences using a nonhomogeneous evolution model. The newly reconstructed ancestral proteins are still thermally stable, further supporting the hypothesis that the ancient organisms contained thermally stable proteins and therefore that they were thermophilic.     

Biomarker/bioindicator response profiles of organisms can help differentiate between sources of anthropogenic stressors in aquatic ecosystems

Aquatic ecosystems can be chronically stressed by multiple environmental factors which originate from a variety of point and non-point sources. In addition, these stressors may vary both spatially and temporally, and, combined with synergestic and cumulative interactions of these stressors, complicate the interpretation and evaluation of stress responses in organisms. To help identify and differentiate between sources of anthropogenic stressors in aquatic systems, a diagnostic approach based on exposure-response profiles in sentinel organisms was developed from the known effects of various anthropogenic activities on biological systems. To generate these exposure-effects profiles, biomarkers of exposure were plotted against bioindicators of corresponding effects for several major anthropogenic activities including petrochemical, pulp and paper, domestic sewage, mining operations, land-development, and agricultural activities. Biomarkers of exposure to environmental stressors varied widely depending on the type of anthropogenic activity involved. Bioindicator effects, however, including histopathological lesions, bioenergetic status, growth, reproductive impairment, and community-level endpoints were similar among several of the major anthropogenic activities because responses at these higher levels are less specific to stressors than are biomarkers. This approach appears useful for helping to identify and diagnose sources of stress in environments impacted by multiple stressors. By identifying the types and sources of environmental stressors impacting key components of biological systems, aquatic ecosystems can be more effectively protected, regulated, and managed to help improve and maintain environmental quality and ecosystem fitness.     

Distribution and evolution of chitinase genes in Streptomyces species: involvement of gene-duplication and domain-deletion

Streptomyces coelicolor A3(2) possesses nine genes for family 18 chitinases and two for family 19, showing high multiplicity. By hybridization analyses, distribution of those chitinase genes was investigated in six other Streptomyces species covering the whole phylogenetic range based on 16S rDNA sequences. All strains showed high-multiplicity of chitinase genes, like S. coelicolor A3(2). The phylogeny and gene organization of the family 18 chitinase genes cloned from Streptomyces species so far were then analyzed to investigate the gene evolution. It was concluded that Streptomyces already possessed a variety of chitinase genes prior to branching into many species, and that the ancestral genes of chiA and chiB have been generated by gene-duplication. In the course of the analyses, evidence that the chi30 and chi40 genes of S. thermoviolaceus were derived from their corresponding original chitinase genes by losing gene parts for substrate-binding domains and fibronectin type III-like domains was obtained. It was thus shown that gene-duplication and domain-deletion were implicated in generating the high diversity and multiplicity of chitinase genes in Streptomyces species. This revised version was published online in June 2006 with corrections to the Cover Date.     

Role of a reaction coordinate in free energy calculations

Abstract

The free energy calculation method emerges as a viable technique for ‘in-silico’ calorimetry. Efficient sampling techniques and the good choice of a reaction path connecting the reactant and the product state enable accurate computations of the free energy differences. We argue that in many cases the thermodynamic integration technique has the lowest variance when the transformation between the reactant and the product state proceeds along the natural path of the studied chemical reaction. We provide examples of free energy calculations for the fragmentation of the charged clusters and the swapping reaction of oligomer formation in proteins that follow a tentative reaction mechanism.     

Application of three-dimensional molecular hydrophobicity potential to the analysis of spatial organization of membrane domains in proteins: I. Hydrophobic properties of transmembrane segments of Na+, K+-ATPase

A new computer-aided molecular modeling approach based on the concept of three-dimensional (3D) molecular hydrophobicity potential has been developed to calculate the spatial organization of intramembrane domains in proteins. The method has been tested by calculating the arrangement of membrane-spanning segments in the photoreaction center ofRhodopseudomonas viridis and comparing the results obtained with those derived from the X-ray data. We have applied this computational procedure to the analysis of interhelical packing in membrane moiety of Na⁺, K⁺-ATPase. The work consists of three parts. In Part I, 3D distributions of electrostatic and molecular hydrophobicity potentials on the surfaces of transmembrane helical peptides were computed and visualized. The hydrophobic and electrostatic properties of helices are discussed from the point of view of their possible arrangement within the protein molecule. Interlocation of helical segments connected with short extramembrane loops found by means of optimization of their hydrophobic/hydrophilic contacts is considered in Part II. The most probable 3D model of packing of helical peptides in the membrane domain of Na⁺, K⁺-ATPase is discussed in the final part of the work.