鸭源鸡杆菌的分离鉴定及其全基因组序列分析

【背景】鸭源鸡杆菌作为一种条件致病菌能引起家禽卵巢炎、输卵管炎和腹膜炎等疾病,严重威胁养殖业的发展。【目的】四川某养鸡场送检了一批疑似感染鸭源鸡杆菌的病死鸡,为探究其感染机制与防治方法,对该菌进行分离鉴定及全基因组测序分析。【方法】从病料中分离并纯化细菌,再依次进行生化试验、16S rRNA基因序列分析和药敏试验,同时通过全基因组测序对其进行物种分型与毒力、耐药等基因功能注释及遗传进化分析。【结果】该分离菌被鉴定为鸭源鸡杆菌,菌株命名为TS0001,药敏试验显示其仅对硝基呋喃类和少数β-内酰胺类药物敏感,对部分β-内酰胺类、氯霉素、部分氨基糖苷类、大环内酯类、四环素类和磺胺类药物具有耐药性。全基因组序列长度为2 626 722 bp,蛋白质编码基因功能注释显示其有较强的自我加工修饰能力,全基因组注释到83个与毒力因子和耐药性相关的基因,包含4个前噬菌体区域。序列类型分析结果显示,该菌株为ST69型,而且管家基因联合建树表明其与墨西哥普通家鸡分离株7990一致性最高。【结论】本研究为鸭源鸡杆菌的感染机制与防治研究提供了参考,丰富了后续研究的分子生物学背景。     

Comparative genomics and selection analysis of Yeonsan Ogye black chicken with whole-genome sequencing

Yeonsan Ogye (OGYE; Gallus gallus domesticus) is a rare indigenous chicken breed that inhabits the Korean Peninsula. This breed has completely black coloring, including plumage, skin, eyes, beak, and internal organs. Despite these unique morphological characteristics, the population of OGYE has declined without in-depth research into their genome research. Therefore, this study aimed to compare the whole genome of OGYE to 12 other chicken populations, including ancestral breed, commercial breeds, Chinese indigenous breeds, and Korean native chickens. We focused on revealing the selection signature of OGYE, which has occurred through environmental pressures in the Korean Peninsula. Genome-wide selection analysis has identified local adaptation traits, such as egg development, that contribute to fetal viability and innate immune response to prevent viral and microbes infection in OGYE. In particular, SPP1 (Secreted Phosphoprotein 1), HSP90AA1 (Heat Shock Protein 90 Alpha Family Class A Member 1), and P2RX4 (Purinergic Receptor P2X 4) could have considerable involvement in egg development and RNASEL (Ribonuclease L), BRIP1 (BRCA1 Interacting Protein C-terminal Helicase 1), and TLR4 (Toll-Like Receptor 4) are crucial for the determination of the innate immune response. This study revealed the unique genetic diversity of OGYE at the genome-wide level. Furthermore, we emphasized the sustainable management of genetic resources and formulated breeding strategies for livestock on the Korean Peninsula.     

Integration of whole-genome DNA methylation data with RNA sequencing data to identify markers for bull fertility

Predicting bull fertility prior to breeding is a current challenge for the dairy industry. The use of molecular biomarkers has been previously assessed. However, the integration of this information has not been performed to extract biologically relevant markers. The goal of this study was to integrate DNA methylation data with previously published RNA-sequencing results in order to identify candidate markers for sire fertility. A total of 1765 differentially methylated cytosines were found between high- and low-fertility sires. Ten genes associated with 11 differentially methylated cytosines were found in a previous study of gene expression between high- and low-fertility sires. Additionally, two of these genes code for proteins found exclusively in bull seminal plasma. Collectively, our results reveal 10 genes that could be used in the future as a panel for predicting bull fertility.     

High throughput genomic sequencing of bioaerosols in broiler chicken production facilities

Chronic inhalation exposure to agricultural dust promotes the development of chronic respiratory diseases among poultry workers. Poultry dust is composed of dander, chicken feed, litter bedding and microbes. However, the microbial composition and abundance has not been fully elucidated. Genomic DNA was extracted from settled dust and personal inhalable dust collected while performing litter sampling or mortality collection tasks. DNA libraries were sequenced using a paired‐end sequencing‐by‐synthesis approach on an Illumina HiSeq 2500. Sequencing data showed that poultry dust is predominantly composed of bacteria (64–67%) with a small quantity of avian, human and feed DNA (< 2% of total reads). Staphylococcus sp. AL1, Salinicoccus carnicancri and Lactobacillus crispatus were the most abundant bacterial species in personal exposure samples of inhalable dust. Settled dust had a moderate relative abundance of these species as well as Staphylococcus lentus and Lactobacillus salivarius. There was a statistical difference between the microbial composition of aerosolized and settled dust. Unlike settled dust composition, aerosolized dust composition had little variance between samples. These data provide an extensive analysis of the microbial composition and relative abundance in personal inhalable poultry dust and settled poultry dust.     

Performance comparison of whole-genome sequencing platforms

Whole-genome sequencing is becoming commonplace, but the accuracy and completeness of variant calling by the most widely used platforms from Illumina and Complete Genomics have not been reported. Here we sequenced the genome of an individual with both technologies to a high average coverage of ～76×, and compared their performance with respect to sequence coverage and calling of single-nucleotide variants (SNVs), insertions and deletions (indels). Although 88.1% of the ～3.7 million unique SNVs were concordant between platforms, there were tens of thousands of platform-specific calls located in genes and other genomic regions. In contrast, 26.5% of indels were concordant between platforms. Target enrichment validated 92.7% of the concordant SNVs, whereas validation by genotyping array revealed a sensitivity of 99.3%. The validation experiments also suggested that >60% of the platform-specific variants were indeed present in the genome. Our results have important implications for understanding the accuracy and completeness of the genome sequencing platforms.     

Efficient identification of genomic insertions and flanking regions through whole-genome sequencing in three transgenic soybean events

Transgenic Research - Genomic insertions and flanking regions of transgenes in host genomes constitute a critical component of precise molecular characterization and event-specific detection, which...     

Identifying insertion mutations by whole-genome sequencing

Insertion mutagenesis via mobile genetic element is a common technique for the analysis of gene function in model organisms. Next-generation sequencing offers an attractive approach for localizing the site of insertion, but alignment-based mapping of mobile genetic elements is challenging. A computational method for identifying insertion sites is reported herein. The technique was validated by mapping transposons in both bacterial and nematode species. The approach should be extensible to other systems that employ mobile genetic elements to generate mutations.     

Bioinformatics for whole-genome shotgun sequencing of microbial communities

The application of whole-genome shotgun sequencing to microbial communities represents a major development in metagenomics, the study of uncultured microbes via the tools of modern genomic analysis. In the past year, whole-genome shotgun sequencing projects of prokaryotic communities from an acid mine biofilm, the Sargasso Sea, Minnesota farm soil, three deep-sea whale falls, and deep-sea sediments have been reported, adding to previously published work on viral communities from marine and fecal samples. The interpretation of this new kind of data poses a wide variety of exciting and difficult bioinformatics problems. The aim of this review is to introduce the bioinformatics community to this emerging field by surveying existing techniques and promising new approaches for several of the most interesting of these computational problems.     

Toward genomic prediction from whole-genome sequence data: impact of sequencing design on genotype imputation and accuracy of predictions

Genomic prediction from whole-genome sequence data is attractive, as the accuracy of genomic prediction is no longer bounded by extent of linkage disequilibrium between DNA markers and causal mutations affecting the trait, given the causal mutations are in the data set. A cost-effective strategy could be to sequence a small proportion of the population, and impute sequence data to the rest of the reference population. Here, we describe strategies for selecting individuals for sequencing, based on either pedigree relationships or haplotype diversity. Performance of these strategies (number of variants detected and accuracy of imputation) were evaluated in sequence data simulated through a real Belgian Blue cattle pedigree. A strategy (AHAP), which selected a subset of individuals for sequencing that maximized the number of unique haplotypes (from single-nucleotide polymorphism panel data) sequenced gave good performance across a range of variant minor allele frequencies. We then investigated the optimum number of individuals to sequence by fold coverage given a maximum total sequencing effort. At 600 total fold coverage (x 600), the optimum strategy was to sequence 75 individuals at eightfold coverage. Finally, we investigated the accuracy of genomic predictions that could be achieved. The advantage of using imputed sequence data compared with dense SNP array genotypes was highly dependent on the allele frequency spectrum of the causative mutations affecting the trait. When this followed a neutral distribution, the advantage of the imputed sequence data was small; however, when the causal mutations all had low minor allele frequencies, using the sequence data improved the accuracy of genomic prediction by up to 30%.     

Research ethics and the challenge of whole-genome sequencing

The recent completion of the first two individual whole-genome sequences is a research milestone. As personal genome research advances, investigators and international research bodies must ensure ethical research conduct. We identify three major ethical considerations that have been implicated in whole-genome research: the return of research results to participants; the obligations, if any, that are owed to participants' relatives; and the future use of samples and data taken for whole-genome sequencing. Although the issues are not new, we discuss their implications for personal genomics and provide recommendations for appropriate management in the context of research involving individual whole-genome sequencing.     

Caenorhabditis elegans mutant allele identification by whole-genome sequencing

Identification of the molecular lesion in Caenorhabditis elegans mutants isolated through forward genetic screens usually involves time-consuming genetic mapping. We used Illumina deep sequencing technology to sequence a complete, mutant C. elegans genome and thus pinpointed a single-nucleotide mutation in the genome that affects a neuronal cell fate decision. This constitutes a proof-of-principle for using whole-genome sequencing to analyze C. elegans mutants.     

Amplification-free whole-genome bisulfite sequencing by post-bisulfite adaptor tagging

DNA methylation plays a key role in epigenetic regulation of eukaryotic genomes. Hence the genome-wide distribution of 5-methylcytosine, or the methylome, has been attracting intense attention. In recent years, whole-genome bisulfite sequencing (WGBS) has enabled methylome analysis at single-base resolution. However, WGBS typically requires microgram quantities of DNA as well as global PCR amplification, thereby precluding its application to samples of limited amounts. This is presumably because bisulfite treatment of adaptor-tagged templates, which is inherent to current WGBS methods, leads to substantial DNA fragmentation. To circumvent the bisulfite-induced loss of intact sequencing templates, we conceived an alternative method termed Post-Bisulfite Adaptor Tagging (PBAT) wherein bisulfite treatment precedes adaptor tagging by two rounds of random primer extension. The PBAT method can generate a substantial number of unamplified reads from as little as subnanogram quantities of DNA. It requires only 100 ng of DNA for amplification-free WGBS of mammalian genomes. Thus, the PBAT method will enable various novel applications that would not otherwise be possible, thereby contributing to the rapidly growing field of epigenomics.     

Allele-specific copy-number discovery from whole-genome and whole-exome sequencing

Copy-number variants (CNVs) are a major form of genetic variation and a risk factor for various human diseases, so it is crucial to accurately detect and characterize them. It is conceivable that allele-specific reads from high-throughput sequencing data could be leveraged to both enhance CNV detection and produce allele-specific copy number (ASCN) calls. Although statistical methods have been developed to detect CNVs using whole-genome sequence (WGS) and/or whole-exome sequence (WES) data, information from allele-specific read counts has not yet been adequately exploited. In this paper, we develop an integrated method, called AS-GENSENG, which incorporates allele-specific read counts in CNV detection and estimates ASCN using either WGS or WES data. To evaluate the performance of AS-GENSENG, we conducted extensive simulations, generated empirical data using existing WGS and WES data sets and validated predicted CNVs using an independent methodology. We conclude that AS-GENSENG not only predicts accurate ASCN calls but also improves the accuracy of total copy number calls, owing to its unique ability to exploit information from both total and allele-specific read counts while accounting for various experimental biases in sequence data. Our novel, user-friendly and computationally efficient method and a complete analytic protocol is freely available at https://sourceforge.net/projects/asgenseng/.     

Development of genomic resources for Nothofagus species using next‐generation sequencing data

Using next‐generation sequencing, we developed the first whole‐genome resources for two hybridizing Nothofagus species of the Patagonian forests that crucially lack genomic data, despite their ecological and industrial value. A de novo assembly strategy combining base quality control and optimization of the putative chloroplast gene map yielded ~32 000 contigs from 43% of the reads produced. With 12.5% of assembled reads, we covered ~96% of the chloroplast genome and ~70% of the mitochondrial gene content, providing functional and structural annotations for 112 and 52 genes, respectively. Functional annotation was possible on 15% of the contigs, with ~1750 potentially novel nuclear genes identified for Nothofagus species. We estimated that the new resources (13.41 Mb in total) included ~4000 gene regions representing ~6.5% of the expected genic partition of the genome, the remaining contigs potentially being nongenic DNA. A high‐quality single nucleotide polymorphisms resource was developed by comparing various filtering methods, and preliminary results indicate a strong conservation of cpDNA genomes in contrast to numerous exclusive nuclear polymorphisms in both species. Finally, we characterized 2274 potential simple sequence repeat (SSR) loci, designed primers for 769 of them and validated nine of 29 loci in 42 individuals per species. Nothofagus obliqua had more alleles (4.89) on average than N. nervosa (2.89), 8 SSRs were efficient to discriminate species, and three were successfully transferred in three other Nothofagus species. These resources will greatly help for future inferences of demographic, adaptive and hybridizing events in Nothofagus species, and for conserving and managing natural populations.     

An automated system for genome analysis to support microbial whole-genome shotgun sequencing

We developed a semi-automated genome analysis system called GAMBLER in order to support the current whole-genome sequencing project focusing on alkaliphilic Bacillus halodurans C-125. GAMBLER was designed to reduce the human intervention required and to reduce the complications in annotating thousands of ORFs in the microbial genome. GAMBLER automates three major routines: analyzing assembly results provided by genome assembler software, assigning ORFs, and homology searching. GAMBLER is equipped with an interface for convenience of annotation. All processes and options are manipulatable through a WWW browser that enables scientists to share their genome analysis results without choosing computer platforms.     

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.     

Generation of EST and microarray resources for functional genomic studies on chicken intestinal health

Expressed sequenced tags (ESTs) and microarray resources have a great impact on the ability to study host response in mice and humans. Unfortunately, these resources are not yet available for domestic farm animals. The aim of this study was to provide genomic resources to study chicken intestinal health, in particular malabsorption syndrome (MAS), which affects mainly the intestine. Therefore a normalized and subtracted cDNA library containing more than 7000 clones was prepared. Randomly chosen clones were sequenced for control purposes. New ESTs were found and multiple ESTs not identified in the chicken intestine before were observed. The number of non-specific ESTs in this cDNA library was low. Based on this normalized and subtracted library a cDNA microarray was made. In a preliminary hybridization experiment with the microarray, genes were identified to be up- or downregulated in MAS infected chickens. This indicates that the generated resources are valuable tools to investigate chicken intestinal health by whole genome expression analysis approaches.     

Comparative whole-genome hybridization reveals genomic islands in Brucella species

Brucella species are responsible for brucellosis, a worldwide zoonotic disease causing abortion in domestic animals and Malta fever in humans. Based on host preference, the genus is divided into six species. Brucella abortus, B. melitensis, and B. suis are pathogenic to humans, whereas B. ovis and B. neotomae are nonpathogenic to humans and B. canis human infections are rare. Limited genome diversity exists among Brucella species. Comparison of Brucella species whole genomes is, therefore, likely to identify factors responsible for differences in host preference and virulence restriction. To facilitate such studies, we used the complete genome sequence of B. melitensis 16M, the species highly pathogenic to humans, to construct a genomic microarray. Hybridization of labeled genomic DNA from Brucella species to this microarray revealed a total of 217 open reading frames (ORFs) altered in five Brucella species analyzed. These ORFs are often found in clusters (islands) in the 16M genome. Examination of the genomic context of these islands suggests that many are horizontally acquired. Deletions of genetic content identified in Brucella species are conserved in multiple strains of the same species, and genomic islands missing in a given species are often restricted to that particular species. These findings suggest that, whereas the loss or gain of genetic material may be related to the host range and virulence restriction of certain Brucella species for humans, independent mechanisms involving gene inactivation or altered expression of virulence determinants may also contribute to these differences.