Variant Identification in Multi-Sample Pools by Illumina Genome Analyzer Sequencing

Multi-sample pooling and Illumina Genome Analyzer (GA) sequencing allows high throughput sequencing of multiple samples to determine population sequence variation. A preliminary experiment, using the RET proto-oncogene as a model, predicted ≤30 samples could be pooled to reliably detect singleton variants without requiring additional confirmation testing. This report used 30 and 50 sample pools to test the hypothesized pooling limit and also to test recent protocol improvements, Illumina GAIIx upgrades, and longer read chemistry. The SequalPrep^TM method was used to normalize amplicons before pooling. For comparison, a single ‘control’ sample was run in a different flow cell lane. Data was evaluated by variant read percentages and the subtractive correction method which utilizes the control sample. In total, 59 variants were detected within the pooled samples, which included all 47 known true variants. The 15 known singleton variants due to Sanger sequencing had an average of 1.62±0.26% variant reads for the 30 pool (expected 1.67% for a singleton variant [unique variant within the pool]) and 1.01±0.19% for the 50 pool (expected 1%). The 76 base read lengths had higher error rates than shorter read lengths (33 and 50 base reads), which eliminated the distinction of true singleton variants from background error. This report demonstrated pooling limits from 30 up to 50 samples (depending on error rates and coverage), for reliable singleton variant detection. The presented pooling protocols and analysis methods can be used for variant discovery in other genes, facilitating molecular diagnostic test design and interpretation.     

Rapid Development of Microsatellite Markers for Callosobruchus chinensis Using Illumina Paired-End Sequencing

Background

The adzuki bean weevil, Callosobruchus chinensis L., is one of the most destructive pests of stored legume seeds such as mungbean, cowpea, and adzuki bean, which usually cause considerable loss in the quantity and quality of stored seeds during transportation and storage. However, a lack of genetic information of this pest results in a series of genetic questions remain largely unknown, including population genetic structure, kinship, biotype abundance, and so on. Co-dominant microsatellite markers offer a great resolving power to determine these events. Here, we report rapid microsatellite isolation from C. chinensis via high-throughput sequencing.

Principal Findings

In this study, 94,560,852 quality-filtered and trimmed reads were obtained for the assembly of genome using Illumina paired-end sequencing technology. In total, the genome with total length of 497,124,785 bp, comprising 403,113 high quality contigs was generated with de novo assembly. More than 6800 SSR loci were detected and a suit of 6303 primer pair sequences were designed and 500 of them were randomly selected for validation. Of these, 196 pair of primers, i.e. 39.2%, produced reproducible amplicons that were polymorphic among 8 C. chinensis genotypes collected from different geographical regions. Twenty out of 196 polymorphic SSR markers were used to analyze the genetic diversity of 18 C. chinensis populations. The results showed the twenty SSR loci were highly polymorphic among these populations.

Conclusions

This study presents a first report of genome sequencing and de novo assembly for C. chinensis and demonstrates the feasibility of generating a large scale of sequence information and SSR loci isolation by Illumina paired-end sequencing. Our results provide a valuable resource for C. chinensis research. These novel markers are valuable for future genetic mapping, trait association, genetic structure and kinship among C. chinensis.     

Identification and Characteristics of Batrachuperus karlschmidti miRNA Using Illumina Deep Sequencing

Russian Journal of Bioorganic Chemistry - MicroRNAs (miRNAs) play critical roles in regulating many biological processes of eukaryotes. Batrachuperus karlschmidti, as an endangered amphibian...     

Identification of Optimum Sequencing Depth Especially for De Novo Genome Assembly of Small Genomes Using Next Generation Sequencing Data

Next Generation Sequencing (NGS) is a disruptive technology that has found widespread acceptance in the life sciences research community. The high throughput and low cost of sequencing has encouraged researchers to undertake ambitious genomic projects, especially in de novo genome sequencing. Currently, NGS systems generate sequence data as short reads and de novo genome assembly using these short reads is computationally very intensive. Due to lower cost of sequencing and higher throughput, NGS systems now provide the ability to sequence genomes at high depth. However, currently no report is available highlighting the impact of high sequence depth on genome assembly using real data sets and multiple assembly algorithms. Recently, some studies have evaluated the impact of sequence coverage, error rate and average read length on genome assembly using multiple assembly algorithms, however, these evaluations were performed using simulated datasets. One limitation of using simulated datasets is that variables such as error rates, read length and coverage which are known to impact genome assembly are carefully controlled. Hence, this study was undertaken to identify the minimum depth of sequencing required for de novo assembly for different sized genomes using graph based assembly algorithms and real datasets. Illumina reads for E.coli (4.6 MB) S.kudriavzevii (11.18 MB) and C.elegans (100 MB) were assembled using SOAPdenovo, Velvet, ABySS, Meraculous and IDBA-UD. Our analysis shows that 50X is the optimum read depth for assembling these genomes using all assemblers except Meraculous which requires 100X read depth. Moreover, our analysis shows that de novo assembly from 50X read data requires only 6–40 GB RAM depending on the genome size and assembly algorithm used. We believe that this information can be extremely valuable for researchers in designing experiments and multiplexing which will enable optimum utilization of sequencing as well as analysis resources.     

Oxford Nanopore MinION Sequencing and Genome Assembly

The revolution of genome sequencing is continuing after the successful second-generation sequencing (SGS) technology. The third-generation sequencing (TGS) technology, led by Pacific Biosciences (PacBio), is progressing rapidly, moving from a technology once only capable of providing data for small genome analysis, or for performing targeted screening, to one that pro-mises high quality de novo assembly and structural variation detection for human-sized genomes. In 2014, the MinION, the first commercial sequencer using nanopore technology, was released by Oxford Nanopore Technologies (ONT). MinION identifies DNA bases by measuring the changes in electrical conductivity generated as DNA strands pass through a biological pore. Its portability, affordability, and speed in data production makes it suitable for real-time applications, the release of the long read sequencer MinION has thus generated much excitement and interest in the geno-mics community. While de novo genome assemblies can be cheaply produced from SGS data, assem-bly continuity is often relatively poor, due to the limited ability of short reads to handle long repeats. Assembly quality can be greatly improved by using TGS long reads, since repetitive regions can be easily expanded into using longer sequencing lengths, despite having higher error rates at the base level. The potential of nanopore sequencing has been demonstrated by various studies in gen-ome surveillance at locations where rapid and reliable sequencing is needed, but where resources are limited.     

基于高通量测序的厚唇裸重唇鱼微卫星分子标记筛选

本研究对厚唇裸重唇鱼(Gymnodiptychus pachycheilus)基因组进行Illumina HiSeq高通量测序,共获得85 905条序列,包含567 200个微卫星位点,从中筛选出15个微卫星位点,采用雅砻江新龙种群和黄河渭河种群对其多态性进行了验证。新龙种群中,15个位点的平均等位基因数为6.9(3 ~ 13),观测杂合度(Ho)为0.712 4,多态信息含量(PIC)为0.630 3,有8个位点显著偏离哈迪-温伯格平衡,1对位点表现出连锁。渭河种群中,15个位点的平均等位基因数为8.4(4 ~ 16),观测杂合度(Ho)为0.719 5,多态信息含量(PIC)为0.680 7,有7个位点显著偏离了哈迪-温伯格平衡,4对位点表现出连锁。筛选获得的15个微卫星位点多态性较高,适合用于厚唇裸重唇鱼种群遗传学研究。     

微生物全基因组测序组装中的g印填补方法

目的：研究和开发高通量测序全基因组组装过程中的填补gap的方法。方法：研究组装软件的算法,使用Perl语言编写自动填补gap的程序,并建立全基因组组装的流程。结果：提出了填补gap的末端延伸法,并使用Perl语言进行了编程;在对立克次体高通量测序的组装过程中,这些方法能大大减少gap的数量。结论：本研究提出的末端延伸法能够高效填补全序列组装过程中出现的gap,具有很强的实用性。     

水稻基因组测序及基因功能的鉴定

水稻是重要的粮食作物。作为单子叶模式植物,水稻基因组的大规模测序具有巨大的理论价值和现实意义。目前已获得了籼稻“93—11”和粳稻“日本晴”高质量的基因组数据,这为在基因组水平上深入研究其生长、发育、抗病和高产等的遗传机理提供了便利,从而为进一步解决世界粮食危机提供了新的突破口和契机。随着水稻基因组计划的顺利结束,其研究重心也已由建立高分辨率的遗传、物理和转录图谱为主的结构基因组学转向基因功能的研究。结构基因组学研究获得的大量序列数据为揭示和开发功能基因开辟了广阔的前景。目前,利用图位克隆和电子克隆等方法已成功分离了多个水稻抗病、抗虫、抗逆境、抗倒伏、高产、优质等重要农艺性状相关的基因,对培育水稻新品种,促进农业的可持续发展意义重大。据估计,水稻至少拥有3．7万个非转座因子相关的蛋白编码基因。因此,完成全基因组序列测定后,重要基因功能的鉴定已成为当前基因组学研究的主要目标。反向遗传学、大规模基因功能表达谱分析和蛋白质组研究等策略已在研究水稻重要基因的功能方面发挥了重要作用。文章综述了水稻基因组测序及基因功能研究的现状,并就新基因发掘和基因功能注释的方法作了评述,期待为水稻遗传工程和育种实践提供参考。     

Genome Sequencing and Identification of Gene Function in Rice

Rice is known to be one of the most important crops for human consumption. As the model cereal crop, large-scale sequencing of rice genome must play quite important roles both in theoretical research and practical application in rice breeding, which announces the opening of another new way to resolve the world food crisis. At present, the emphasis of rice genome research has been transferred from structure genomics to functional analysis. The discovery of new genes and annotation of gene function was believed to be an important issue in functional genomics research. In this article, the sequencing and functional research of the rice genome were reviewed. These results may provide some useful clues for rice genetic engineering and breeding practices.     

真菌基因组de novo测序组装的方法与实践

真菌基因组较其他真核生物基因组结构简单,长度短,易于测序、组装与注释,因此真菌基因组是研究真核生物基因组的模型。为研究真菌基因组组装策略,本研究基于Illumina HiSeq测序平台对烟曲霉菌株An16007基因组测序,分别使用5种de novo组装软件ABySS、SOAP-denovo、Velvet、MaSuRCA和IDBA-UD组装基因组,然后通过Augustus软件进行基因预测,BUSCO软件评估组装结果。研究发现,5种组装软件对基因组组装结果不同,ABySS组装的基因组较其他4种组装软件具有更高的完整性和准确性,且预测的基因数量较高,因此,ABySS更适合本研究基因组的组装。本研究提供了真菌de novo测序、组装及组装质量评估的技术流程,为基因组<100 Mb的真菌或其他生物基因组的研究提供参考。     

Characterization of the Soybean Genome Using EST-derived Microsatellite Markers

We generated a high-density genetic linkage map of soybean usingexpressed sequence tag (EST)-derived microsatellite markers.A total of 6920 primer pairs (10.9%) were designed to amplifysimple sequence repeats (SSRs) from 63 676 publicly availablenon-redundant soybean ESTs. The polymorphism of two parent plants,the Japanese cultivar ‘Misuzudaizu’ and the Chineseline ‘Moshidou Gong 503’, were examined using 10%polyacrylamide gel electrophoresis. Primer pairs showing polymorphismwere then used for genotyping 94 recombinant inbred lines (RILs)derived from a cross between the parents. In addition to previouslyreported markers, 680 EST-derived microsatellite markers wereselected and subjected to linkage analysis. As a result, 935marker loci were mapped successfully onto 20 linkage groups,which totaled 2700.3 cM in length; 693 loci were detected usingthe 668 EST-derived microsatellite markers developed in thisstudy, the other 242 loci were detected with 105 RFLP markers,136 genome-derived microsatellite markers, and one phenotypicmarker. We examined allelic variation among 23 soybean cultivars/linesand a wild soybean line using 668 mapped EST-derived microsatellitemarkers (corresponding to 686 marker loci), in order to determinethe transferability of the markers among soybean germplasms.A limited degree of macrosynteny was observed at the segmentallevel between the genomes of soybean and the model legume Lotusjaponicus, which suggests that considerable genome shufflingoccurred after separation of the species and during establishmentof the paleopolyploid soybean genome.     

Species Identification and Profiling of Complex Microbial Communities Using Shotgun Illumina Sequencing of 16S rRNA Amplicon Sequences

The high throughput and cost-effectiveness afforded by short-read sequencing technologies, in principle, enable researchers to perform 16S rRNA profiling of complex microbial communities at unprecedented depth and resolution. Existing Illumina sequencing protocols are, however, limited by the fraction of the 16S rRNA gene that is interrogated and therefore limit the resolution and quality of the profiling. To address this, we present the design of a novel protocol for shotgun Illumina sequencing of the bacterial 16S rRNA gene, optimized to amplify more than 90% of sequences in the Greengenes database and with the ability to distinguish nearly twice as many species-level OTUs compared to existing protocols. Using several in silico and experimental datasets, we demonstrate that despite the presence of multiple variable and conserved regions, the resulting shotgun sequences can be used to accurately quantify the constituents of complex microbial communities. The reconstruction of a significant fraction of the 16S rRNA gene also enabled high precision (>90%) in species-level identification thereby opening up potential application of this approach for clinical microbial characterization.