Aneuploidy in dizygotic twin sheep detected using genome-wide single nucleotide polymorphism data from two commonly used commercial vendors

Early detection of karyotype abnormalities, including aneuploidy, could aid producers in identifying animals which, for example, would not be suitable candidate parents. Genome-wide genetic marker data in the form of single nucleotide polymorphisms (SNPs) are now being routinely generated on animals. The objective of the present study was to describe the statistics that could be generated from the allele intensity values from such SNP data to diagnose karyotype abnormalities; of particular interest was whether detection of aneuploidy was possible with both commonly used genotyping platforms in agricultural species, namely the Applied Biosystems^TM Axiom^TM and the Illumina platform. The hypothesis was tested using a case study of a set of dizygotic X-chromosome monosomy 53,X sheep twins. Genome-wide SNP data were available from the Illumina platform (11 082 autosomal and 191 X-chromosome SNPs) on 1848 male and 8954 female sheep and available from the Axiom^TM platform (11 128 autosomal and 68 X-chromosome SNPs) on 383 female sheep. Genotype allele intensity values, either as their original raw values or transformed to logarithm intensity ratio (LRR), were used to accurately diagnose two dizygotic (i.e. fraternal) twin 53,X sheep, both of which received their single X chromosome from their sire. This is the first reported case of 53,X dizygotic twins in any species. Relative to the X-chromosome SNP genotype mean allele intensity values of normal females, the mean allele intensity value of SNP genotypes on the X chromosome of the two females monosomic for the X chromosome was 7.45 to 12.4 standard deviations less, and were easily detectable using either the Axiom^TM or Illumina genotype platform; the next lowest mean allele intensity value of a female was 4.71 or 3.3 standard deviations less than the population mean depending on the platform used. Both 53,X females could also be detected based on the genotype LRR although this was more easily detectable when comparing the mean LRR of the X chromosome of each female to the mean LRR of their respective autosomes. On autopsy, the ovaries of the two sheep were small for their age and evidence of prior ovulation was not appreciated. In both sheep, the density of primordial follicles in the ovarian cortex was lower than normally found in ovine ovaries and primary follicle development was not observed. Mammary gland development was very limited. Results substantiate previous studies in other species that aneuploidy can be readily detected using SNP genotype allele intensity values generally already available, and the approach proposed in the present study was agnostic to genotype platform.     

The new sequencer on the block: comparison of Life Technology’s Proton sequencer to an Illumina HiSeq for whole-exome sequencing

We assessed the performance of the new Life Technologies Proton sequencer by comparing whole-exome sequence data in a Centre d’Etude du Polymorphisme Humain trio (family 1463) to the Illumina HiSeq instrument. To simulate a typical user’s results, we utilized the standard capture, alignment and variant calling methods specific to each platform. We restricted data analysis to include the capture region common to both methods. The Proton produced high quality data at a comparable average depth and read length, and the Ion Reporter variant caller identified 96 % of single nucleotide polymorphisms (SNPs) detected by the HiSeq and GATK pipeline. However, only 40 % of small insertion and deletion variants (indels) were identified by both methods. Usage of the trio structure and segregation of platform-specific alleles supported this result. Further comparison of the trio data with Complete Genomics sequence data and Illumina SNP microarray genotypes documented high concordance and accurate SNP genotyping of both Proton and Illumina platforms. However, our study underscored the problem of accurate detection of indels for both the Proton and HiSeq platforms.     

An efficient quantitation method of next-generation sequencing libraries by using MiSeq sequencer

Library quantitation is a critical step to obtain high data output in Illumina HiSeq sequencers. Here, we introduce a library quantitation method that uses the Illumina MiSeq sequencer designated as quantitative MiSeq (qMiSeq). In this procedure, 96 dual-index libraries, including control samples, are denatured, pooled in equal volume, and sequenced by MiSeq. We found that relative concentration of each library can be determined based on the observed index ratio and can be used to determine HiSeq run condition for each library. Thus, qMiSeq provides an efficient way to quantitate a large number of libraries at a time.     

Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms

DNA sequencing continues to decrease in cost with the Illumina HiSeq2000 generating up to 600 Gb of paired-end 100 base reads in a ten-day run. Here we present a protocol for community amplicon sequencing on the HiSeq2000 and MiSeq Illumina platforms, and apply that protocol to sequence 24 microbial communities from host-associated and free-living environments. A critical question as more sequencing platforms become available is whether biological conclusions derived on one platform are consistent with what would be derived on a different platform. We show that the protocol developed for these instruments successfully recaptures known biological results, and additionally that biological conclusions are consistent across sequencing platforms (the HiSeq2000 versus the MiSeq) and across the sequenced regions of amplicons.     

A bioinformatics approach for determining sample identity from different lanes of high-throughput sequencing data

The ability to generate whole genome data is rapidly becoming commoditized. For example, a mammalian sized genome (～3Gb) can now be sequenced using approximately ten lanes on an Illumina HiSeq 2000. Since lanes from different runs are often combined, verifying that each lane in a genome's build is from the same sample is an important quality control. We sought to address this issue in a post hoc bioinformatic manner, instead of using upstream sample or "barcode" modifications. We rely on the inherent small differences between any two individuals to show that genotype concordance rates can be effectively used to test if any two lanes of HiSeq 2000 data are from the same sample. As proof of principle, we use recent data from three different human samples generated on this platform. We show that the distributions of concordance rates are non-overlapping when comparing lanes from the same sample versus lanes from different samples. Our method proves to be robust even when different numbers of reads are analyzed. Finally, we provide a straightforward method for determining the gender of any given sample. Our results suggest that examining the concordance of detected genotypes from lanes purported to be from the same sample is a relatively simple approach for confirming that combined lanes of data are of the same identity and quality.     

Comparison of Sequencing Platforms for Single Nucleotide Variant Calls in a Human Sample

Next-generation sequencings platforms coupled with advanced bioinformatic tools enable re-sequencing of the human genome at high-speed and large cost savings. We compare sequencing platforms from Roche/454(GS FLX), Illumina/HiSeq (HiSeq 2000), and Life Technologies/SOLiD (SOLiD 3 ECC) for their ability to identify single nucleotide substitutions in whole genome sequences from the same human sample. We report on significant GC-related bias observed in the data sequenced on Illumina and SOLiD platforms. The differences in the variant calls were investigated with regards to coverage, and sequencing error. Some of the variants called by only one or two of the platforms were experimentally tested using mass spectrometry; a method that is independent of DNA sequencing. We establish several causes why variants remained unreported, specific to each platform. We report the indel called using the three sequencing technologies and from the obtained results we conclude that sequencing human genomes with more than a single platform and multiple libraries is beneficial when high level of accuracy is required.     

Validation of multiple single nucleotide variation calls by additional exome analysis with a semiconductor sequencer to supplement data of whole-genome sequencing of a human population

Background

Validation of single nucleotide variations in whole-genome sequencing is critical for studying disease-related variations in large populations. A combination of different types of next-generation sequencers for analyzing individual genomes may be an efficient means of validating multiple single nucleotide variations calls simultaneously.

Results

Here, we analyzed 12 independent Japanese genomes using two next-generation sequencing platforms: the Illumina HiSeq 2500 platform for whole-genome sequencing (average depth 32.4×), and the Ion Proton semiconductor sequencer for whole exome sequencing (average depth 109×). Single nucleotide polymorphism (SNP) calls based on the Illumina Human Omni 2.5-8 SNP chip data were used as the reference. We compared the variant calls for the 12 samples, and found that the concordance between the two next-generation sequencing platforms varied between 83% and 97%.

Conclusions

Our results show the versatility and usefulness of the combination of exome sequencing with whole-genome sequencing in studies of human population genetics and demonstrate that combining data from multiple sequencing platforms is an efficient approach to validate and supplement SNP calls.

Electronic supplementary material

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

濒危珍稀植物半枫荷的转录组分析

为加强中国特有濒危植物半枫荷资源的保护与利用工作,采用高通量测序平台Illumina HiSeq 2500对其进行转录组测序,将得到的数据过滤后进行de novo组装并聚类去冗余,获得77 629个Unigenes,通过九大功能数据库比对、分析、注释,最终有45 293个Unigenes获得注释信息;其中在KOG按功能分为25个子类,获得25 253个功能注释信息;GO功能注释可分为细胞组分、生物学过程、分子功能3大类、分别可细分为22、26、17亚类（共65个亚类）;与KEGG数据库对比,共发现286条代谢通路,其中发现可能与半枫荷药用活性成分相关的次生代谢产物生物合成的177条途径;根据组装结果预测出88个基因家族共1 547个编码转录因子的Unigenes,发现控制药效合成的转录因子家族。另外,根据注释结果检测到12 579个SNP多态位点和预测出57 671个CDS位点。本研究首次对半枫荷转录组进行分析,为深入开展半枫荷分子生物学研究提供基础数据来源。     

基于比较转录组分析海洋弧菌X511的褐藻胶代谢途径

【目的】从北海涠洲岛海域腐烂的马尾藻中分离得到的海洋弧菌(Vibrio X511)具有较强的利用褐藻胶能力,本文利用转录组测序的方法以研究弧菌X511的褐藻胶代谢途径。【方法】采用Illumina Hi Seq2500测序平台对菌株在褐藻胶及葡萄糖培养下的转录组进行测序;比较和分析差异转录本,利用荧光定量PCR验证测序结果;采用GO(Gene Ontology)和KEGG(Kyoto Encyclopedia of Genes and Genomes)对差异转录本进行功能和Pathway注释。【结果】经比较发现,菌株在褐藻胶培养下相对于葡萄糖的培养共有2024个差异表达基因,其中1066个基因上调,958个基因下调;某些普遍存在于代谢途径中的基因在不同培养条件下也存在差异表达;海洋弧菌X511中涉及褐藻胶利用的所有基因以及合成乙醇的关键基因其转录量均有一定程度的上调;此外,通过分析发现该菌株具有独特的褐藻胶利用方式,其中的一个代谢过程尚未在弧菌中被报道。【结论】成功解析了海洋弧菌X511的褐藻胶代谢途径,丰富了生物方法降解褐藻胶的研究,为大型海藻生物质能源的研究提供有价值的数据支持。     

Multiplex sequencing of bacterial artificial chromosomes for assembling complex plant genomes

Hierarchical shotgun sequencing remains the method of choice for assembling high‐quality reference sequences of complex plant genomes. The efficient exploitation of current high‐throughput technologies and powerful computational facilities for large‐insert clone sequencing necessitates the sequencing and assembly of a large number of clones in parallel. We developed a multiplexed pipeline for shotgun sequencing and assembling individual bacterial artificial chromosomes (BACs) using the Illumina sequencing platform. We illustrate our approach by sequencing 668 barley BACs (Hordeum vulgare L.) in a single Illumina HiSeq 2000 lane. Using a newly designed parallelized computational pipeline, we obtained sequence assemblies of individual BACs that consist, on average, of eight sequence scaffolds and represent >98% of the genomic inserts. Our BAC assemblies are clearly superior to a whole‐genome shotgun assembly regarding contiguity, completeness and the representation of the gene space. Our methods may be employed to rapidly obtain high‐quality assemblies of a large number of clones to assemble map‐based reference sequences of plant and animal species with complex genomes by sequencing along a minimum tiling path.     

紫背天葵MBW相关调控因子转录组测序分析

为获取紫背天葵(Cynura bicolor)花青素合成相关转录调控因子MBW家族基因,采用二代高通量测序技术进行全转录功能基因组测序后组装,再通过Pfam、Swiss Prot和Nr数据库搜索,共获得138个MBW相关Unigene,分别有42个MYB、67个b HLH、15个b HLH-MYB和14个WD40,其中目前已报道与花青素合成代谢相关的MYB、b HLH、b HLH-MYB和WD40分别为11、33、6和3个。这为进一步研究紫背天葵花青素合成的调控机理和相关基因克隆等奠定基础。     

An Efficient Genotyping Method in Chicken Based on Genome Reducing and Sequencing

Single nucleotide polymorphisms (SNPs) are essential for identifying the genetic mechanisms of complex traits. In the present study, we applied genotyping by genome reducing and sequencing (GGRS) method to construct a 252-plex sequencing library for SNP discovery and genotyping in chicken. The library was successfully sequenced on an Illumina HiSeq 2500 sequencer with a paired-end pattern; approximately 400 million raw reads were generated, and an average of approximately 1.4 million good reads per sample were generated. A total of 91,767 SNPs were identified after strict filtering, and all of the 252 samples and all of the chromosomes were well represented. Compared with the Illumina 60K chicken SNP chip data, approximately 34,131 more SNPs were identified using GGRS, and a higher SNP density was found using GGRS, which could be beneficial for downstream analysis. Using the GGRS method, more than 3528 samples can be sequenced simultaneously, and the cost is reduced to $18 per sample. To the best of our knowledge, this study describes the first report of such highly multiplexed sequencing in chicken, indicating potential applications for genome-wide association and genomic selection in chicken.     

Reconstructing the plant mitochondrial genome for marker discovery: a case study using Pinus

Whole‐genome‐shotgun (WGS) sequencing of total genomic DNA was used to recover ~1 Mbp of novel mitochondrial (mtDNA) sequence from Pinus sylvestris (L.) and three members of the closely related Pinus mugo species complex. DNA was extracted from megagametophyte tissue from six mother trees from locations across Europe, and 100‐bp paired‐end sequencing was performed on the Illumina HiSeq platform. Candidate mtDNA sequences were identified by their size and coverage characteristics, and by comparison with published plant mitochondrial genomes. Novel variants were identified, and primers targeting these loci were trialled on a set of 28 individuals from across Europe. In total, 31 SNP loci were successfully resequenced, characterizing 15 unique haplotypes. This approach offers a cost‐effective means of developing marker resources for mitochondrial genomes in other plant species where reference sequences are unavailable.     

Selection of the sex‐linked inhibitor of apoptosis in mountain pine beetle (Dendroctonus ponderosae) driven by enhanced expression during early overwintering

The mountain pine beetle (Dendroctonus ponderosae) is an insect native to western North America; however, its geographical range has recently expanded north in BC and east into Alberta. To understand the population structure in the areas of expansion, 16 gene‐linked microsatellites were screened and compared to neutral microsatellites using outlier analyses of F_st and F_ct values. One sex‐linked gene, inhibitor of apoptosis (IAP), showed a strong signature of positive selection for neo‐X alleles and was analyzed for evidence of adaptive variation. Alleles of IAP were sequenced, and differences between the neo‐X and neo‐Y alleles were consistent with neutral evolution suggesting that the neo‐Y allele may not be under functional constraints. Neo‐Y alleles were amplified from gDNA, but not effectively from cDNA, suggesting that there was little IAP expression from neo‐Y alleles. There were no differences in overall IAP expression between males and females with the common northern neo‐X allele suggesting that the neo‐X allele in males compensates for the reduced expression of neo‐Y alleles. However, males lacking the most common northern neo‐X allele thought to be selected for in northern populations had reduced overall IAP expression in early October—at a time when beetles are preparing for overwintering. This suggests that the most common allele may have more rapid upregulation. The reduced function of neo‐Y alleles of IAP suggested by both sequence differences and lower levels of expression may foster a highly selective environment for neo‐X alleles such as the common northern allele with more efficient upregulation.     

Engineering glyceraldehyde-3-phosphate dehydrogenase for switching control of glycolysis in Escherichia coli

Glycolysis has evolved to be a highly robust mechanism for maintaining the cellular metabolism of living organisms. However, relevant modifications of glycolytic activity are required to intentionally modulate cellular phenotypes. Here, we designed a platform that allows switching control of glycolysis in Escherichia coli in response to an environmental signal, in this case, temperature. This system functions by regulating the expression of gapA, which encodes glyceraldehyde‐3‐phosphate dehydrogenase (GAPDH), one of the key glycolytic enzymes. Because a very low level of gapA expression is capable of maintaining cellular physiology, we also modified GAPDH through directed evolution to provide sensitive regulation of glycolytic activity. The switching control of glycolysis was successfully demonstrated by regulating the expression of engineered gapA through changes in temperature. This system offers potential control over the cell's central carbon‐metabolism switch, providing the ability to perform reprogrammed tasks with desired timing depending on environmental signals. Biotechnol. Bioeng. 2012; 109: 2612–2619. © 2012 Wiley Periodicals, Inc.