CopywriteR: DNA copy number detection from off-target sequence data.

Current methods for detection of copy number variants (CNV) and aberrations (CNA) from targeted sequencing data are based on the depth of coverage of captured exons. Accurate CNA determination is complicated by uneven genomic distribution and non-uniform capture efficiency of targeted exons. Here we present CopywriteR, which eludes these problems by exploiting 'off-target' sequence reads. CopywriteR allows for extracting uniformly distributed copy number information, can be used without reference, and can be applied to sequencing data obtained from various techniques including chromatin immunoprecipitation and target enrichment on small gene panels. CopywriteR outperforms existing methods and constitutes a widely applicable alternative to available tools.

A second generation radiation hybrid map to aid the assembly of the bovine genome sequence

BACKGROUND: Several approaches can be used to determine the order of loci on chromosomes and hence develop maps of the genome. However, all mapping approaches are prone to errors either arising from technical deficiencies or lack of statistical support to distinguish between alternative orders of loci. The accuracy of the genome maps could be improved, in principle, if information from different sources was combined to produce integrated maps. The publicly available bovine genomic sequence assembly with 6x coverage (Btau_2.0) is based on whole genome shotgun sequence data and limited mapping data however, it is recognised that this assembly is a draft that contains errors. Correcting the sequence assembly requires extensive additional mapping information to improve the reliability of the ordering of sequence scaffolds on chromosomes. The radiation hybrid (RH) map described here has been contributed to the international sequencing project to aid this process. RESULTS: An RH map for the 30 bovine chromosomes is presented. The map was built using the Roslin 3000-rad RH panel (BovGen RH map) and contains 3966 markers including 2473 new loci in addition to 262 amplified fragment-length polymorphisms (AFLP) and 1231 markers previously published with the first generation RH map. Sequences of the mapped loci were aligned with published bovine genome maps to identify inconsistencies. In addition to differences in the order of loci, several cases were observed where the chromosomal assignment of loci differed between maps. All the chromosome maps were aligned with the current 6x bovine assembly (Btau_2.0) and 2898 loci were unambiguously located in the bovine sequence. The order of loci on the RH map for BTA 5, 7, 16, 22, 25 and 29 differed substantially from the assembled bovine sequence. From the 2898 loci unambiguously identified in the bovine sequence assembly, 131 mapped to different chromosomes in the BovGen RH map. CONCLUSION: Alignment of the BovGen RH map with other published RH and genetic maps showed higher consistency in marker order and chromosome assignment than with the current 6x sequence assembly. This suggests that the bovine sequence assembly could be significantly improved by incorporating additional independent mapping information.

Complete Genome Sequence of Bovine Pestivirus Strain PG-2, a Second Member of the Tentative Pestivirus Species Giraffe

We report the complete genome sequence of bovine pestivirus strain PG-2. The sequence data from this virus showed that PG-2 is closely related to the giraffe pestivirus strain H138. PG-2 and H138 belong to one pestivirus species that should be considered an approved member of the genus Pestivirus.

Imputation of sequence level genotypes in the Franches-Montagnes horse breed

BACKGROUND A cost-effective strategy to increase the density of available markers within a population is to sequence a small proportion of the population and impute whole-genome sequence data for the remaining population. Increased densities of typed markers are advantageous for genome-wide association studies (GWAS) and genomic predictions. METHODS We obtained genotypes for 54 602 SNPs (single nucleotide polymorphisms) in 1077 Franches-Montagnes (FM) horses and Illumina paired-end whole-genome sequencing data for 30 FM horses and 14 Warmblood horses. After variant calling, the sequence-derived SNP genotypes (~13 million SNPs) were used for genotype imputation with the software programs Beagle, Impute2 and FImpute. RESULTS The mean imputation accuracy of FM horses using Impute2 was 92.0%. Imputation accuracy using Beagle and FImpute was 74.3% and 77.2%, respectively. In addition, for Impute2 we determined the imputation accuracy of all individual horses in the validation population, which ranged from 85.7% to 99.8%. The subsequent inclusion of Warmblood sequence data further increased the correlation between true and imputed genotypes for most horses, especially for horses with a high level of admixture. The final imputation accuracy of the horses ranged from 91.2% to 99.5%. CONCLUSIONS Using Impute2, the imputation accuracy was higher than 91% for all horses in the validation population, which indicates that direct imputation of 50k SNP-chip data to sequence level genotypes is feasible in the FM population. The individual imputation accuracy depended mainly on the applied software and the level of admixture.

Like will to like: abundances of closely related species can predict susceptibility to intestinal colonization by pathogenic and commensal bacteria

The intestinal ecosystem is formed by a complex, yet highly characteristic microbial community. The parameters defining whether this community permits invasion of a new bacterial species are unclear. In particular, inhibition of enteropathogen infection by the gut microbiota ( = colonization resistance) is poorly understood. To analyze the mechanisms of microbiota-mediated protection from Salmonella enterica induced enterocolitis, we used a mouse infection model and large scale high-throughput pyrosequencing. In contrast to conventional mice (CON), mice with a gut microbiota of low complexity (LCM) were highly susceptible to S. enterica induced colonization and enterocolitis. Colonization resistance was partially restored in LCM-animals by co-housing with conventional mice for 21 days (LCM(con21)). 16S rRNA sequence analysis comparing LCM, LCM(con21) and CON gut microbiota revealed that gut microbiota complexity increased upon conventionalization and correlated with increased resistance to S. enterica infection. Comparative microbiota analysis of mice with varying degrees of colonization resistance allowed us to identify intestinal ecosystem characteristics associated with susceptibility to S. enterica infection. Moreover, this system enabled us to gain further insights into the general principles of gut ecosystem invasion by non-pathogenic, commensal bacteria. Mice harboring high commensal E. coli densities were more susceptible to S. enterica induced gut inflammation. Similarly, mice with high titers of Lactobacilli were more efficiently colonized by a commensal Lactobacillus reuteri(RR) strain after oral inoculation. Upon examination of 16S rRNA sequence data from 9 CON mice we found that closely related phylotypes generally display significantly correlated abundances (co-occurrence), more so than distantly related phylotypes. Thus, in essence, the presence of closely related species can increase the chance of invasion of newly incoming species into the gut ecosystem. We provide evidence that this principle might be of general validity for invasion of bacteria in preformed gut ecosystems. This might be of relevance for human enteropathogen infections as well as therapeutic use of probiotic commensal bacteria.

Molecular epidemiology reveals long-term changes in HIV type 1 subtype B transmission in Switzerland

Sequence data from resistance testing offer unique opportunities to characterize the structure of human immunodeficiency virus (HIV) infection epidemics.

Rapid parallel adaptive radiations from a single hybridogenic ancestral population

The Alpine lake whitefish (Coregonus lavaretus) species complex is a classic example of a recent radiation, associated with colonization of the Alpine lakes following the glacial retreat (less than 15 kyr BP). They have formed a unique array of endemic lake flocks, each with one to six described sympatric species differing in morphology, diet and reproductive ecology. Here, we present a genomic investigation of the relationships between and within the lake flocks. Comparing the signal between over 1000 AFLP loci and mitochondrial control region sequence data, we use phylogenetic tree-based and population genetic methods to reconstruct the phylogenetic history of the group and to delineate the principal centres of genetic diversity within the radiation. We find significant cytonuclear discordance showing that the genomically monophyletic Alpine whitefish clade arose from a hybrid swarm of at least two glacial refugial lineages. Within this radiation, we find seven extant genetic clusters centred on seven lake systems. Most interestingly, we find evidence of sympatric speciation within and parallel evolution of equivalent phenotypes among these lake systems. However, we also find the genetic signature of human-mediated gene flow and diversity loss within many lakes, highlighting the fragility of recent radiations.

ABCtoolbox: a versatile toolkit for approximate Bayesian computations

Background The estimation of demographic parameters from genetic data often requires the computation of likelihoods. However, the likelihood function is computationally intractable for many realistic evolutionary models, and the use of Bayesian inference has therefore been limited to very simple models. The situation changed recently with the advent of Approximate Bayesian Computation (ABC) algorithms allowing one to obtain parameter posterior distributions based on simulations not requiring likelihood computations. Results Here we present ABCtoolbox, a series of open source programs to perform Approximate Bayesian Computations (ABC). It implements various ABC algorithms including rejection sampling, MCMC without likelihood, a Particle-based sampler and ABC-GLM. ABCtoolbox is bundled with, but not limited to, a program that allows parameter inference in a population genetics context and the simultaneous use of different types of markers with different ploidy levels. In addition, ABCtoolbox can also interact with most simulation and summary statistics computation programs. The usability of the ABCtoolbox is demonstrated by inferring the evolutionary history of two evolutionary lineages of Microtus arvalis. Using nuclear microsatellites and mitochondrial sequence data in the same estimation procedure enabled us to infer sex-specific population sizes and migration rates and to find that males show smaller population sizes but much higher levels of migration than females. Conclusion ABCtoolbox allows a user to perform all the necessary steps of a full ABC analysis, from parameter sampling from prior distributions, data simulations, computation of summary statistics, estimation of posterior distributions, model choice, validation of the estimation procedure, and visualization of the results.

Serodiagnosis of Echinococcus spp. infection: explorative selection of diagnostic antigens by peptide microarray

BACKGROUND: Production of native antigens for serodiagnosis of helminthic infections is laborious and hampered by batch-to-batch variation. For serodiagnosis of echinococcosis, especially cystic disease, most screening tests rely on crude or purified Echinococcus granulosus hydatid cyst fluid. To resolve limitations associated with native antigens in serological tests, the use of standardized and highly pure antigens produced by chemical synthesis offers considerable advantages, provided appropriate diagnostic sensitivity and specificity is achieved. METHODOLOGY/PRINCIPAL FINDINGS: Making use of the growing collection of genomic and proteomic data, we applied a set of bioinformatic selection criteria to a collection of protein sequences including conceptually translated nucleotide sequence data of two related tapeworms, Echinococcus multilocularis and Echinococcus granulosus. Our approach targeted alpha-helical coiled-coils and intrinsically unstructured regions of parasite proteins potentially exposed to the host immune system. From 6 proteins of E. multilocularis and 5 proteins of E. granulosus, 45 peptides between 24 and 30 amino acids in length were designed. These peptides were chemically synthesized, spotted on microarrays and screened for reactivity with sera from infected humans. Peptides reacting above the cut-off were validated in enzyme-linked immunosorbent assays (ELISA). Peptides identified failed to differentiate between E. multilocularis and E. granulosus infection. The peptide performing best reached 57% sensitivity and 94% specificity. This candidate derived from Echinococcus multilocularis antigen B8/1 and showed strong reactivity to sera from patients infected either with E. multilocularis or E. granulosus. CONCLUSIONS/SIGNIFICANCE: This study provides proof of principle for the discovery of diagnostically relevant peptides by bioinformatic selection complemented with screening on a high-throughput microarray platform. Our data showed that a single peptide cannot provide sufficient diagnostic sensitivity whereas pooling several peptide antigens improved sensitivity; thus combinations of several peptides may lead the way to new diagnostic tests that replace, or at least complement conventional immunodiagnosis of echinococcosis. Our strategy could prove useful for diagnostic developments in other pathogens.

Differences in Soil Fungal Communities between European Beech (Fagus sylvatica L.) Dominated Forests Are Related to Soil and Understory Vegetation

Fungi are important members of soil microbial communities with a crucial role in biogeochemical processes. Although soil fungi are known to be highly diverse, little is known about factors influencing variations in their diversity and community structure among forests dominated by the same tree species but spread over different regions and under different managements. We analyzed the soil fungal diversity and community composition of managed and unmanaged European beech dominated forests located in three German regions, the Schwäbische Alb in Southwestern, the Hainich-Dün in Central and the Schorfheide Chorin in the Northeastern Germany, using internal transcribed spacer (ITS) rDNA pyrotag sequencing. Multiple sequence quality filtering followed by sequence data normalization revealed 1655 fungal operational taxonomic units. Further analysis based on 722 abundant fungal OTUs revealed the phylum Basidiomycota to be dominant (54%) and its community to comprise 71.4% of ectomycorrhizal taxa. Fungal community structure differed significantly (p≤0.001) among the three regions and was characterized by non-random fungal OTUs co-occurrence. Soil parameters, herbaceous understory vegetation, and litter cover affected fungal community structure. However, within each study region we found no difference in fungal community structure between management types. Our results also showed region specific significant correlation patterns between the dominant ectomycorrhizal fungal genera. This suggests that soil fungal communities are region-specific but nevertheless composed of functionally diverse and complementary taxa.

Alboluxin, a snake C-type lectin from Trimeresurus albolabris venom is a potent platelet agonist acting via GPIb and GPVI

Alboluxin, a potent platelet activator, was purified from Trimeresurus albolabris venom with a mass of 120 kDa non-reduced and, after reduction, subunits of 17 and 24 kDa. Alboluxin induced a tyrosine phosphorylation profile in platelets that resembles those produced by collagen and convulxin, involving the time dependent tyrosine phosphorylation of Fc receptor gamma chain (Fc gamma), phospholipase Cgamma2 (PLCgamma2), LAT and p72SYK. Antibodies against both GPIb and GPVI inhibited platelet aggregation induced by alboluxin, whereas antibodies against alpha2beta1 had no effect. Inhibition of alphaIIb beta3 reduced the aggregation response to alboluxin, as well as tyrosine phosphorylation of platelet proteins, showing that activation of alphaIIb beta3 and binding of fibrinogen are involved in alboluxin-induced platelet aggregation and it is not simply agglutination. N-terminal sequence data from the beta-subunit of alboluxin indicates that it belongs to the snake C-type lectin family. The C-type lectin subunits are larger than usual possibly due to post-translational modifications such as glycosylation. Alboluxin is a hexameric (alphabeta)3 snake C-type lectin which activates platelets via both GPIb and GPVI.

Emended description of Actinobacillus capsulatus Arseculeratne 1962, 38AL

The taxonomic position of Actinobacillus capsulatus, a member of the family Pasteurellaceae found in rabbits, hares and hamsters, has been challenged. 16S rRNA gene (rrs) sequence data show the species to be heterogeneous. Using a polyphasic approach, 23 strains that were identified previously as belonging, or closely related, to A. capsulatus were analysed. Eighty characters were included in the phenotypic analysis. Phylogenetic analysis was done based on rrs, rpoB, infB and recN sequences. In addition, the recN sequence similarities were used to calculate the whole-genome sequence relatedness of all strains investigated as well as that with other members of the family Pasteurellaceae. The phenotypic analysis allowed identification of five groups. The major group of 17 strains could be classified as A. capsulatus. Two hamster isolates were closely related to A. capsulatus but differed in a few characters. Single isolates from a rabbit and snowshoe-hare were phenotypically related to Actinobacillus suis. One rabbit isolate was related to the genus Mannheimia, while another isolate could not be classified phenotypically with known taxa. The phylogenetic analysis confirmed the phenotypic grouping. In contrast to the rrs-based tree, the A. capsulatus strains clustered unambiguously with the type species and related species of the genus Actinobacillus in the rpoB-, infB- and recN-based trees. Genome similarity comparison using recN finally confirmed the high genomic relationship of the A. capsulatus strains with the type species and related species of the genus Actinobacillus and allowed a clear assignment of the other unrelated strains to the phenotypic and phylogenetic clusters outlined. The present findings allow the description of A. capsulatus to be emended and separate it more clearly from other species, both phenotypically and genotypically. The type strain of A. capsulatus is CCUG 12396(T) (=Frederiksen 243(T)=ATCC 51571(T)=NCTC 11408(T)=CIP 103283(T)).

Novel mutation in OTC gene causes neonatal death in twin brothers

Ornithine transcarbamylase (OTC) deficiency is the most common inborn error of the urea cycle. OTC locus is located in the short arm of X-chromosome. Authors report a case of a woman who gave birth to monozygotic male twins who later died because of severe neonatal-onset hyperammonaemic encephalopathy caused by a novel mutation of OTC gene. Post-mortem liver biopsy was taken from the second twin; afterwards, blood was drawn from the mother for examination. DNA sequence data showed that the mother was a carrier of the same novel mutation that was previously detected in the case of her son. In OTC deficiency, detection of female carriers is important for genetic counselling and eventual prenatal diagnosis.

A Nonsense Mutation in the IKBKG Gene in Mares with Incontinentia Pigmenti.

Ectodermal dysplasias (EDs) are a large and heterogeneous group of hereditary disorders characterized by abnormalities in structures of ectodermal origin. Incontinentia pigmenti (IP) is an ED characterized by skin lesions evolving over time, as well as dental, nail, and ocular abnormalities. Due to X-linked dominant inheritance IP symptoms can only be seen in female individuals while affected males die during development in utero. We observed a family of horses, in which several mares developed signs of a skin disorder reminiscent of human IP. Cutaneous manifestations in affected horses included the development of pruritic, exudative lesions soon after birth. These developed into wart-like lesions and areas of alopecia with occasional wooly hair re-growth. Affected horses also had streaks of darker and lighter coat coloration from birth. The observation that only females were affected together with a high number of spontaneous abortions suggested an X-linked dominant mechanism of transmission. Using next generation sequencing we sequenced the whole genome of one affected mare. We analyzed the sequence data for non-synonymous variants in candidate genes and found a heterozygous nonsense variant in the X-chromosomal IKBKG gene (c.184C>T; p.Arg62*). Mutations in IKBKG were previously reported to cause IP in humans and the homologous p.Arg62* variant has already been observed in a human IP patient. The comparative data thus strongly suggest that this is also the causative variant for the observed IP in horses. To our knowledge this is the first large animal model for IP.