Evolutionary insight from whole-genome sequencing of experimentally evolved microbes

Experimental evolution (EE) combined with whole‐genome sequencing (WGS) has become a compelling approach to study the fundamental mechanisms and processes that drive evolution. Most EE‐WGS studies published to date have used microbes, owing to their ease of propagation and manipulation in the laboratory and relatively small genome sizes. These experiments are particularly suited to answer long‐standing questions such as: How many mutations underlie adaptive evolution, and how are they distributed across the genome and through time? Are there general rules or principles governing which genes contribute to adaptation, and are certain kinds of genes more likely to be targets than others? How common is epistasis among adaptive mutations, and what does this reveal about the variety of genetic routes to adaptation? How common is parallel evolution, where the same mutations evolve repeatedly and independently in response to similar selective pressures? Here, we summarize the significant findings of this body of work, identify important emerging trends and propose promising directions for future research. We also outline an example of a computational pipeline for use in EE‐WGS studies, based on freely available bioinformatics tools.     

Habitat selection as a mechanism of resource partitioning in two cryptic bat species Pipistrellus pipistrellus and Pipistrellus pygmaeus

Ecomorphological studies of bat communities often reveal the spatial and temporal coexistence of morphologically similar species, leading to suggestions that these communities are structured by non-deterministic processes. However, the diversification of echolocation call structure in bats allows for considerable morphological similarity while still permitting niche differentiation based on specialisation for prey type and habitat structure. The recent separation of a common Palaearctic bat, the pipistrelle, into Pipistrellus pipistrellus and P. pygmaeus, which are sympatrically distributed throughout their range, raises the question as to whether these two morphologically similar species partition resources in time and space.
To test the hypothesis that the coexistence of these cryptic species is facilitated by differential habitat use, 14 P . pipistrellus , and 12 P. pygmaeus were radio-tracked from adjacent maternity roosts, in northeast Scotland, from May to September 2002/2003. The two species showed distinct habitat partitioning with P. pygmaeus foraging predominantly in riparian woodland and over water, and P. pipistrellus foraging along woodland edges and short isolated tree lines. Inter-specific overlap in habitat use was low and consequently foraging ranges were segregated spatially.
The degree of habitat partitioning revealed in these species, which show considerable overlap in echolocation call parameters and functional morphology, suggests that morphological features, whilst useful in separating chiropteran species into coarse-grained foraging guilds, may not predict fine-grained ecological segregation.     

Looking for therapeutic antibodies in next-generation sequencing repositories

ABSTRACT

Recently it has become possible to query the great diversity of natural antibody repertoires using next-generation sequencing (NGS). These methods are capable of producing millions of sequences in a single experiment. Here we compare clinical-stage therapeutic antibodies to the ~1b sequences from 60 independent sequencing studies in the Observed Antibody Space database, which includes antibody sequences from NGS analysis of immunoglobulin gene repertoires. Of 242 post-Phase 1 antibodies, we found 16 with sequence identity matches of 95% or better for both heavy and light chains. There are also 54 perfect matches to therapeutic CDR-H3 regions in the NGS outputs, suggesting a nontrivial amount of convergence between naturally observed sequences and those developed artificially. This has potential implications for both the legal protection of commercial antibodies and the discovery of antibody therapeutics.     

Massive parallel sequencing in animal genetics: wherefroms and wheretos

Next generation sequencing (NGS) has revolutionized genomics research, making it difficult to overstate its impact on studies of Biology. NGS will immediately allow researchers working in non‐mainstream species to obtain complete genomes together with a comprehensive catalogue of variants. In addition, RNA‐seq will be a decisive way to annotate genes that cannot be predicted purely by computational or comparative approaches. Future applications include whole genome sequence association studies, as opposed to classical SNP‐based association, and implementing this new source of information into breeding programmes. For these purposes, one of the main advantages of sequencing vs. genotyping is the possibility of identifying copy number variants. Currently, experimental design is a topic of utmost interest, and here we discuss some of the options available, including pools and reduced representation libraries. Although bioinformatics is still an important bottleneck, this limitation is only transient and should not deter animal geneticists from embracing these technologies.     

Milk microbiome signatures of subclinical mastitis-affected cattle analysed by shotgun sequencing

Aims: Metagenomic analysis of milk samples collected from Kankrej, Gir (Bos indicus) and crossbred (Bos taurus × B. indicus) cattle harbouring subclinical mastitis was carried out by next‐generation sequencing 454 GS‐FLX technology to elucidate the microbial community structure of cattle milk. Methods and Results: Milk samples from Kankrej, Gir and crossbred cattle were subjected to metagenomic profiling by pyrosequencing. The Metagenomic analysis produced 63·07, 11·09 and 7·87 million base pairs (Mb) of sequence data, assembled in 264 798, 56 114 and 36 762 sequences with an average read length of 238, 197 and 214 nucleotides in Kankrej, Gir and crossbred cattle, respectively. Phylogenetic and metabolic profiles by the web‐based tool MG‐RAST revealed that the members of Enterobacteriales were predominant in mastitic milk followed by Pseudomonadales, Bacillales and Lactobacillales. Around 56 different species with varying abundance were detected in the subclinically infected milk. Escherichia coli was found to be the most predominant species in Kankrej and Gir cattle followed by Pseudomonas aeruginosa, Pseudomonas mendocina, Shigella flexneri and Bacillus cereus. In crossbred cattle, Staphylococcus aureus followed by Klebsiella pneumoniae, Staphylococcus epidermidis and E. coli were detected in descending order. Metabolic profiling indicated fluoroquinolones, methicillin, copper, cobalt–zinc–cadmium as the groups of antibiotics and toxic compounds to which the organisms showed resistance. Sequences indicating potential of organisms exhibiting multidrug resistance against antibiotics and resistance to toxic compounds were also present. Interestingly, presence of bacteriophages against Staph. aureus, E. coli, Enterobacter and Yersinia species was also observed. Conclusions: The analysis identified potential infectious organisms in mastitis, resistance of organisms to antibiotics and chemical compounds and the natural resistance potential of dairy cows. Significance and Impact of the Study: The findings of this study may help in formulating strategies for the prevention and treatment of mastitis in dairy animals and consequently in reducing economic losses incurred because of it.     

Rise of the machines--recommendations for ecologists when using next generation sequencing for microsatellite development

Next generation sequencing is revolutionizing molecular ecology by simplifying the development of molecular genetic markers, including microsatellites. Here, we summarize the results of the large-scale development of microsatellites for 54 nonmodel species using next generation sequencing and show that there are clear differences amongst plants, invertebrates and vertebrates for the number and proportion of motif types recovered that are able to be utilized as markers. We highlight that the heterogeneity within each group is very large. Despite this variation, we provide an indication of what number of sequences and consequent proportion of a 454 run are required for the development of 40 designable, unique microsatellite loci for a typical molecular ecological study. Finally, to address the challenges of choosing loci from the vast array of microsatellite loci typically available from partial genome runs (average for this study, 2341 loci), we provide a microsatellite development flowchart as a procedural guide for application once the results of a partial genome run are obtained.     

Partial short-read sequencing of a highly inbred Iberian pig and genomics inference thereof

Despite dramatic reduction in sequencing costs with the advent of next generation sequencing technologies, obtaining a complete mammalian genome sequence at sufficient depth is still costly. An alternative is partial sequencing. Here, we have sequenced a reduced representation library of an Iberian sow from the Guadyerbas strain, a highly inbred strain that has been used in numerous QTL studies because of its extreme phenotypic characteristics. Using the Illumina Genome Analyzer II (San Diego, CA, USA), we resequenced ～ 1% of the genome with average 4 × depth, identifying 68,778 polymorphisms. Of these, 55,457 were putative fixed differences with respect to the assembly, based on the genome of a Duroc pig, and 13,321 were heterozygous positions within Guadyerbas. Despite being highly inbred, the estimate of heterozygosity within Guadyerbas was ～ 0.78 kb(-1) in autosomes, after correcting for low depth. Nucleotide variability was consistently higher at the telomeric regions than on the rest of the chromosome, likely a result of increased recombination rates. Further, variability was 50% lower in the X-chromosome than in autosomes, which may be explained by a recent bottleneck or by selection. We divided the whole genome in 500 kb windows and we analyzed overrepresented gene ontology terms in regions of low and high variability. Multi organism process, pigmentation and cell killing were overrepresented in high variability regions and metabolic process ontology, within low variability regions. Further, a genome wide Hudson-Kreitman-Aguadé test was carried out per window; overall, variability was in agreement with neutral expectations.     

jMHC: software assistant for multilocus genotyping of gene families using next-generation amplicon sequencing

Genotyping of multilocus gene families, such as the major histocompatibility complex (MHC), may be challenging because of problems with assigning alleles to loci and copy number variation among individuals. Simultaneous amplification and genotyping of multiple loci may be necessary, and in such cases, next-generation deep amplicon sequencing offers a great promise as a genotyping method of choice. Here, we describe jMHC, a computer program developed for analysing and assisting in the visualization of deep amplicon sequencing data. Software operates on FASTA files; therefore, output from any sequencing technology may be used. jMHC was designed specifically for MHC studies but it may be useful for analysing amplicons derived from other multigene families or for genotyping other polymorphic systems. The program is written in Java with user-friendly graphical interface (GUI) and can be run on Microsoft Windows, Linux OS and Mac OS.     

新一代测序技术在植物转录组研究中的应用

随着DNA测序技术的发展,新一代测序技术以其高通量、低成本的特点,成为越来越多的生物学研究者在开展工作时的首选。在所有的新一代测序技术中,454测序系统是最早实现商业化且发展相对成熟的一种,目前被广泛的应用于各个领域的生物学研究中。文章以454测序系统为例,综述了新一代测序系统的原理、优缺点,及其在植物转录组研究中的应用,并对其在植物研究领域中可能的发展应用方向进行了展望。     

新一代测序技术应用于microRNA检测

MicroRNAs(miRNAs)是一类在进化上高度保守的非编码小分子单链RNA(~22nt),在基因转录后调控中发挥至关重要的作用。越来越多的证据表明,miRNAs参与很多重要的生理和病理过程,例如发育、器官形成、调亡、细胞增殖、肿瘤发生等。近年来飞速发展的新一代测序技术在miRNA检测方面具有重要的应用。文章简要介绍了新一代测序技术3大平台的基本步骤和原理,测序数据的生物信息学分析方法以及新一代测序技术在miRNA方向的主要应用。相比于传统的miRNA检测方法,新一代测序技术具有通量高、对遗传物质检测完全且准确度高,可重复性好等优点,在探索新miRNA、miRNA互补链、miRNA编辑、miRNA异构体检测以及miRNA靶基因检测等方面具有巨大优势。随着新一代测序技术的不断发展,测序成本不断降低,在未来几年,新一代测序技术的使用率或将大大增加。新一代测序技术的不断应用将进一步促进人类对于miRNA在各种生理病理过程中的功能和调控的认识。     

Diel variation in feeding and vertical distribution of ten co-occurring fish species: consequences for resource partitioning

The diel variations in feeding behaviour and vertical distribution were determined for ten species of a tropical fish community in a shallow SE Sri Lankan reservoir. The fish community consisted of two introduced exotic tilapias and eight indigenous riverine species including five cyprinids, one clarid, one hemiramphid and one gobid. Multi-mesh gillnets dividing the water column in six depth strata were set with a six-hour interval throughout four 24-hour periods. The study of the diel variation in feeding behaviour was based on analysis of catch data, gut contents and gut fullness per setting time. The diel variation in vertical distribution was determined per species for each of the six depth-strata. Feeding behaviour and vertical distribution differed significantly among species and for most species throughout the 24-hour period. Feeding behaviour and vertical distribution were correlated. The status of whether or not a fish is feeding determines if the trophic and spatial resource dimensions are interrelated. Resource partitioning along the spatial dimension may reduce exploitative competition at the time of feeding, whereas it may limit interference competition or predation during non-feeding periods. This revised version was published online in July 2006 with corrections to the Cover Date.     

High-throughput microsatellite isolation through 454 GS-FLX Titanium pyrosequencing of enriched DNA libraries

Microsatellites (or SSRs: simple sequence repeats) are among the most frequently used DNA markers in many areas of research. The use of microsatellite markers is limited by the difficulties involved in their de novo isolation from species for which no genomic resources are available. We describe here a high-throughput method for isolating microsatellite markers based on coupling multiplex microsatellite enrichment and next-generation sequencing on 454 GS-FLX Titanium platforms. The procedure was calibrated on a model species (Apis mellifera) and validated on 13 other species from various taxonomic groups (animals, plants and fungi), including taxa for which severe difficulties were previously encountered using traditional methods. We obtained from 11,497 to 34,483 sequences depending on the species and the number of detected microsatellite loci ranged from 199 to 5791. We thus demonstrated that this procedure can be readily and successfully applied to a large variety of taxonomic groups, at much lower cost than would have been possible with traditional protocols. This method is expected to speed up the acquisition of high-quality genetic markers for nonmodel organisms.     

Development and characterization of polymorphic microsatellite markers for Neolitsea sericea using Illumina paired‐end draft sequencing data

Simple sequence repeat (SSR) markers were developed and characterized for Neolitsea sericea (Bl.) Koidz. (Lauraceae). Out of 196 designed primer pairs, a total of 144 pairs showed amplification, of which 44 had clear and stable chromatograms. Polymorphism of these 44 loci was tested using 32 individuals sampled from a single population of N. sericea. The number of alleles and the polymorphism information content varied from 3 to 12 and 0.271 to 0.853, respectively. A significant departure from the Hardy‐Weinberg equilibrium was observed in one of the 44 loci. These SSR markers are useful for population genetic studies and parentage analysis in N. sericea, which is one of the most common evergreen species in coastal Pinus thunbergii forests in central‐western Japan.     

Isolation and characterization of microsatellite markers for Viola mirabilis (Violaceae) using a next‐generation sequencing platform

We isolated and characterized microsatellite loci in Viola mirabilis (Violaceae), an endangered species from South Korea. Twenty‐three polymorphic microsatellite loci were developed and tested in Korean, Chinese and Japanese populations. The number of alleles per locus varied from two to eight. The observed and expected heterozygosities within the three populations were 0.000–0.625 and 0.469–0.695, respectively. A total of six loci in the Korean population, one locus in the Chinese population and seven loci in the Japanese population deviated from Hardy–Weinberg equilibrium. We expect that these newly developed microsatellite markers will contribute to understanding the phylogeography and population genetics of V. mirabilis, which will aid in developing conservation strategies for this species.     

Identification of high utility SNPs for population assignment and traceability purposes in the pig using high-throughput sequencing

The objectives of this study were to develop breed-specific single nucleotide polymorphisms (SNPs) in five pig breeds sequenced with Illumina's Genome Analyzer and to investigate their usefulness for breed assignment purposes. DNA pools were prepared for Duroc, Landrace, Large White, Pietrain and Wild Boar. The total number of animals used for sequencing was 153. SNP discovery was performed by aligning the filtered reads against Build 7 of the pig genome. A total of 313,964 high confidence SNPs were identified and analysed for the presence of breed-specific SNPs (defined in this context as SNPs for which one of the alleles was detected in only one breed). There were 29,146 putative breed-specific SNPs identified, of which 4441 were included in the PorcineSNP60 beadchip. Upon re-examining the genotypes obtained using the beadchip, 193 SNPs were confirmed as being breed specific. These 193 SNPs were subsequently used to assign an additional 490 individuals from the same breeds, using the sequenced individuals as reference populations. In total, four breed assignment tests were performed. Results showed that for all methods tested 99% of the animals were correctly assigned, with an average probability of assignment of at least 99.2%, indicating the high utility of breed-specific markers for breed assignment and traceability. This study provides a blueprint for the way next-generation sequencing technologies can be used for the identification of breed-specific SNPs, as well as evidence that these SNPs may be a powerful tool for breed assignment and traceability of animal products to their breeds of origin.