Solving the Problem of Comparing Whole Bacterial Genomes across Different Sequencing Platforms

Whole genome sequencing (WGS) shows great potential for real-time monitoring and identification of infectious disease outbreaks. However, rapid and reliable comparison of data generated in multiple laboratories and using multiple technologies is essential. So far studies have focused on using one technology because each technology has a systematic bias making integration of data generated from different platforms difficult. We developed two different procedures for identifying variable sites and inferring phylogenies in WGS data across multiple platforms. The methods were evaluated on three bacterial data sets and sequenced on three different platforms (Illumina, 454, Ion Torrent). We show that the methods are able to overcome the systematic biases caused by the sequencers and infer the expected phylogenies. It is concluded that the cause of the success of these new procedures is due to a validation of all informative sites that are included in the analysis. The procedures are available as web tools.     

ConPADE: Genome Assembly Ploidy Estimation from Next-Generation Sequencing Data

As a result of improvements in genome assembly algorithms and the ever decreasing costs of high-throughput sequencing technologies, new high quality draft genome sequences are published at a striking pace. With well-established methodologies, larger and more complex genomes are being tackled, including polyploid plant genomes. Given the similarity between multiple copies of a basic genome in polyploid individuals, assembly of such data usually results in collapsed contigs that represent a variable number of homoeologous genomic regions. Unfortunately, such collapse is often not ideal, as keeping contigs separate can lead both to improved assembly and also insights about how haplotypes influence phenotype. Here, we describe a first step in avoiding inappropriate collapse during assembly. In particular, we describe ConPADE (Contig Ploidy and Allele Dosage Estimation), a probabilistic method that estimates the ploidy of any given contig/scaffold based on its allele proportions. In the process, we report findings regarding errors in sequencing. The method can be used for whole genome shotgun (WGS) sequencing data. We also show applicability of the method for variant calling and allele dosage estimation. Results for simulated and real datasets are discussed and provide evidence that ConPADE performs well as long as enough sequencing coverage is available, or the true contig ploidy is low. We show that ConPADE may also be used for related applications, such as the identification of duplicated genes in fragmented assemblies, although refinements are needed.     

SomatiCA: Identifying,Characterizing and Quantifying Somatic Copy Number Aberrations from Cancer Genome Sequencing Data

Whole genome sequencing of matched tumor-normal sample pairs is becoming routine in cancer research. However, analysis of somatic copy-number changes from sequencing data is still challenging because of insufficient sequencing coverage, unknown tumor sample purity and subclonal heterogeneity. Here we describe a computational framework, named SomatiCA, which explicitly accounts for tumor purity and subclonality in the analysis of somatic copy-number profiles. Taking read depths (RD) and lesser allele frequencies (LAF) as input, SomatiCA will output 1) admixture rate for each tumor sample, 2) somatic allelic copy-number for each genomic segment, 3) fraction of tumor cells with subclonal change in each somatic copy number aberration (SCNA), and 4) a list of substantial genomic aberration events including gain, loss and LOH. SomatiCA is available as a Bioconductor R package at http://www.bioconductor.org/packages/2.13/bioc/html/SomatiCA.html.     

Genotyping by Genome Reducing and Sequencing for Outbred Animals

Next-generation sequencing (NGS) approaches are widely used in genome-wide genetic marker discovery and genotyping. However, current NGS approaches are not easy to apply to general outbred populations (human and some major farm animals) for SNP identification because of the high level of heterogeneity and phase ambiguity in the haplotype. Here, we reported a new method for SNP genotyping, called genotyping by genome reducing and sequencing (GGRS) to genotype outbred species. Through an improved procedure for library preparation and a marker discovery and genotyping pipeline, the GGRS approach can genotype outbred species cost-effectively and high-reproducibly. We also evaluated the efficiency and accuracy of our approach for high-density SNP discovery and genotyping in a large genome pig species (2.8 Gb), for which more than 70,000 single nucleotide polymorphisms (SNPs) can be identified for an expenditure of only $80 (USD)/sample.     

Population Genetic Inference from Personal Genome Data: Impact of Ancestry and Admixture on Human Genomic Variation

Full sequencing of individual human genomes has greatly expanded our understanding of human genetic variation and population history. Here, we present a systematic analysis of 50 human genomes from 11 diverse global populations sequenced at high coverage. Our sample includes 12 individuals who have admixed ancestry and who have varying degrees of recent (within the last 500 years) African, Native American, and European ancestry. We found over 21 million single-nucleotide variants that contribute to a 1.75-fold range in nucleotide heterozygosity across diverse human genomes. This heterozygosity ranged from a high of one heterozygous site per kilobase in west African genomes to a low of 0.57 heterozygous sites per kilobase in segments inferred to have diploid Native American ancestry from the genomes of Mexican and Puerto Rican individuals. We show evidence of all three continental ancestries in the genomes of Mexican, Puerto Rican, and African American populations, and the genome-wide statistics are highly consistent across individuals from a population once ancestry proportions have been accounted for. Using a generalized linear model, we identified subtle variations across populations in the proportion of neutral versus deleterious variation and found that genome-wide statistics vary in admixed populations even once ancestry proportions have been factored in. We further infer that multiple periods of gene flow shaped the diversity of admixed populations in the Americas—70% of the European ancestry in today’s African Americans dates back to European gene flow happening only 7–8 generations ago.     

不同地区HPV16 E7基因的克隆及序列差异分析

采用PCR扩增、pGEM T载体克隆和核苷酸序列分析的方法对一例武汉地区及两例五峰县高发区宫颈癌患者体内HPV16型的E7基因编码区进行序列分析并与野生型 (德国标准株 )及已发表的HPV16湖北株 (HPVHB)进行了比较。结果发现武汉地区HPV16型E7基因仅第 5 4位出现一个同义突变 ,而高发区HPV16型E7基因存在差异 ,第 77位氨基酸由精氨酸 (Arg)变为半胱氨酸 (Cys) ,第 96位由谷氨酰氨酸 (Gln)变为精氨酸 (Arg) ,E7蛋白的二级结构及亲、疏水性也相应改变 ,与野生型有较大差异     

Mutation Frequencies in HIV-1 Genome in Regions Containing Efficient RNAi Targets As Calculated from Ultra-Deep Sequencing Data

HIV-1 is one of the most variable viruses. The development of gene therapy technology using RNAi for AIDS/HIV-1 treatment is a potential alternative for traditional anti-retroviral therapy. Anti-HIV-1 siRNA should aim to exploit the most conserved viral targets. Using the deep sequencing of potential RNAi targets in 100-nt HIV-1 genome fragments from the clinical HIV-1 subtype A isolates in Russia, we found that the frequencies of all possible transversions and transitions in certain RNAi targets are 3–38 times lower than in adjacent sequences. Therefore, these targets are conserved. We propose the development of these RNAi targets for AIDS/HIV-1 treatment. Deep sequencing also enables the detection of the characteristic mutational bias of RT during the replication of viral RNA.     

不同地区HPV16E7基因的克隆及序列差异分析

采用PCR扩增、pGEM-T载体克隆和核苷酸序列分析的方法对一 例武汉地区及两例五峰县高发区宫颈癌患者体内HPV16型的E7基因编码区进行序列分析并与 野生型(德国标准株)及已发表的HPV16湖北株(HPVHB)进行了比较.结果发现武汉地区HPV16 型E7基因仅第54位出现一个同义突变,而高发区HPV16型E7基因存在差异,第77位氨基酸由 精氨酸(Arg)变为半胱氨酸(Cys),第96位由谷氨酰氨酸(Gln)变为精氨酸(Arg),E7蛋白的二级 结构及亲、疏水性也相应改变,与野生型有较大差异.     

Phylogenetics and Differentiation of Salmonella Newport Lineages by Whole Genome Sequencing

Salmonella Newport has ranked in the top three Salmonella serotypes associated with foodborne outbreaks from 1995 to 2011 in the United States. In the current study, we selected 26 S. Newport strains isolated from diverse sources and geographic locations and then conducted 454 shotgun pyrosequencing procedures to obtain 16–24 × coverage of high quality draft genomes for each strain. Comparative genomic analysis of 28 S. Newport strains (including 2 reference genomes) and 15 outgroup genomes identified more than 140,000 informative SNPs. A resulting phylogenetic tree consisted of four sublineages and indicated that S. Newport had a clear geographic structure. Strains from Asia were divergent from those from the Americas. Our findings demonstrated that analysis using whole genome sequencing data resulted in a more accurate picture of phylogeny compared to that using single genes or small sets of genes. We selected loci around the mutS gene of S. Newport to differentiate distinct lineages, including those between invH and mutS genes at the 3′ end of Salmonella Pathogenicity Island 1 (SPI-1), ste fimbrial operon, and Clustered, Regularly Interspaced, Short Palindromic Repeats (CRISPR) associated-proteins (cas). These genes in the outgroup genomes held high similarity with either S. Newport Lineage II or III at the same loci. S. Newport Lineages II and III have different evolutionary histories in this region and our data demonstrated genetic flow and homologous recombination events around mutS. The findings suggested that S. Newport Lineages II and III diverged early in the serotype evolution and have evolved largely independently. Moreover, we identified genes that could delineate sublineages within the phylogenetic tree and that could be used as potential biomarkers for trace-back investigations during outbreaks. Thus, whole genome sequencing data enabled us to better understand the genetic background of pathogenicity and evolutionary history of S. Newport and also provided additional markers for epidemiological response.     

中丹科学家完成古人类全基因组测序

由中国深圳华大基因研究院和丹麦哥本哈根大学联合创建的中丹基因组联合中心以Saqqaq古人为样本,完成了世界首例古人类全基因组的深度序列测定和解读.这一历史性成果以封面故事发表在2010年2月11日出版的科学期刊上.     

1株梅花鹿源狂犬病街毒株全基因组测序与分析

对1株梅花鹿源性狂犬病街毒株(DRV)进行全基因组克隆,对全长cDNA进行测序分析.RT-PCR扩增克隆覆盖全基因组9个重叠基因片段,基因组3'和5'末端采取3,-RACE和5'-RACE方法,9个重叠基因片段序列拼接得到DRV全基因组cDNA序列,共11 863个核苷酸.DRV毒株全基因组构成与其他狂犬病毒基因组构成相似,由5个编码区组成,基因起始位点和终止位点高度保守,在核蛋白和糖蛋白的重要抗原位点有个别氨基酸发生变异,对已完成全基因组测序的几个基因1型毒株分别进行了N、P、M、G、L基因核苷酸及氨基酸的同源性比较.与其他具有代表性的毒株进行N基因序列比较建立的系统进化树表明,DRV毒株属于基因1型,与中国人用疫苗株3aG同源性最高为94%,与分类位置未确定的北高加索毒株(WCBV)的同源性最低为71%.本研究结果可为狂犬病毒各项分子生物学研究提供理论参考.     

Complex Genotype Mixtures Analyzed by Deep Sequencing in Two Different Regions of Hepatitis B Virus

This study assesses the presence and outcome of genotype mixtures in the polymerase/surface and X/preCore regions of the HBV genome in patients with chronic hepatitis B virus (HBV) infection. Thirty samples from ten chronic hepatitis B patients were included. The polymerase/surface and X/preCore regions were analyzed by deep sequencing (UDPS) in the first available sample at diagnosis, a pre-treatment sample, and a sample while under treatment. HBV genotype was determined by phylogenesis. Quasispecies complexity was evaluated by mutation frequency and nucleotide diversity. The polymerase/surface and X/preCore regions were validated for genotyping from 113 GenBank reference sequences. UDPS yielded a median of 10,960 sequences per sample (IQR 16,645) in the polymerase/surface region and 11_,595 sequences per sample (IQR 14,682) in X/preCore. Genotype mixtures were more common in X/preCore (90%) than in polymerase/surface (30%) (p<0.001). On X/preCore genotyping, all samples were genotype A, whereas polymerase/surface yielded genotypes A (80%), D (16.7%), and F (3.3%) (p = 0.036). Genotype changes in polymerase/surface were observed in four patients during natural quasispecies dynamics and in two patients during treatment. There were no genotype changes in X/preCore. Quasispecies complexity was higher in X/preCore than in polymerase/surface (p = 0.004).The results provide evidence of genotype mixtures and differential genotype proportions in the polymerase/surface and X/preCore regions. The genotype dynamics in HBV infection and the different patterns of quasispecies complexity in the HBV genome suggest a new paradigm for HBV genotype classification.     

一株新的沙门菌烈性噬菌体IME-SAL1的分离、全基因组测序和比较基因组分析

目的:分离鉴定一株沙门菌的烈性噬菌体,观察其形态大小,完成全基因组测序,分析其基因组结构和进化关系,为治疗沙门菌感染提供新的策略和实验依据。方法:以沙门菌SAL95作为指示菌从医院废水中分离噬菌体,分离到的噬菌体经浓缩和纯化后采用透射电镜观察其形态大小,提取噬菌体的基因核酸并完成全基因组高通量测序,分析其全基因组的结构特征,通过比较基因组分析研究其进化关系。结果:从解放军307医院未经消毒处理的废水中分离到一株烈性沙门菌噬菌体。电镜观察显示,该噬菌体头部呈立体对称,有一不收缩的长尾。其基因组全长113 183 bp,比较基因组分析确定该噬菌体为一株新的沙门菌噬菌体,命名为IME-SAL1。结论:从医院废水中分离到一株烈性沙门菌噬菌体IME-SAL1,研究了该噬菌体的分类、基因组结构、进化关系,可为其实际应用提供参考。