新城疫不同毒株交叉鸡胚中和指数及其与F和HN基因变异的相关性

选取13株国内2001～2004年分离的新城疫流行病毒(Newcastle disease virus,NDV),经蚀斑纯化,克隆其融合蛋白(F)和血凝素.神经氨酸酶(HN)基因,结合疫苗株La Sota、Clone30和国内标准强毒株F48E9等的基因序列,进行遗传变异分析.利用纯化的病毒制备特异阳性血清,进行鸡胚交叉中和试验,确定不同NDV毒株之间的抗原相关性,并与NDV不同毒株之间的HN和F基因核苷酸(氨基酸)同源性进行相关比较.结果表明:病毒中和指数与HN基因的核苷酸(氨基酸)同源性显著相关(P<0.01,r=-0.35),与F基因呈弱相关(P<0.05,r=0.20),而与F基因前374bp的核甘酸同源性不相关.这表明,NDV的分子变异已经对NDV的抗原性变异产生了影响,研制新型的疫苗成为必然.     

Understanding animal welfare

Background

Hepatitis E virus (HEV) is a zoonotic pathogen of which swine was reported as major reservoirs. HEV has been divided into 4 different genotypes according to phylogenetic analysis. Recent reports showed that genotype 4 HEV is freely transmitted between humans and swine in eastern China, including Shanghai area. This paper investigated the recent infection status of HEV among swine population of Shanghai area in China.

Methods

480 swine faecal specimens were collected from 23 farms which distribute all over Shanghai from September to November, 2007 and tested for the presence of HEV RNA by the polymerase chain reaction (PCR).

Results

Our results showed that 26.1% (6/23) of the swine farms were positive for HEV RNA and the positive rate of the six farms were ranged from 9.1% to 33.3%. The HEV RNA positive rate for total samples were 5% (24/480). The resulted positive band specific for HEV was sequenced and sequence analysis indicated that all of these isolates belonged to genotype 4 HEV. Phylogenetic analysis showed that the 24 isolates clustered into 4 distinct subgroups, sharing 83.3–89.7% inter-subgroup and 97–99% intra-subgroup identities. More over, isolates in three of the four subgroups closely clustered with previous identified strains, sharing up high to 97% identity with them.

Conclusion

These results suggested that there were 4 different subgenotypes of HEV prevalent in Shanghai, and some of them may not be indigenous to Shanghai but introduced from other geographic regions.     

Rapid virulence prediction and identification of Newcastle disease virus genotypes using third-generation sequencing

Background

Newcastle disease (ND) outbreaks are global challenges to the poultry industry. Effective management requires rapid identification and virulence prediction of the circulating Newcastle disease viruses (NDV), the causative agent of ND. However, these diagnostics are hindered by the genetic diversity and rapid evolution of NDVs.

Methods

An amplicon sequencing (AmpSeq) workflow for virulence and genotype prediction of NDV samples using a third-generation, real-time DNA sequencing platform is described here. 1D MinION sequencing of barcoded NDV amplicons was performed using 33 egg-grown isolates, (15 NDV genotypes), and 15 clinical swab samples collected from field outbreaks. Assembly-based data analysis was performed in a customized, Galaxy-based AmpSeq workflow. MinION-based results were compared to previously published sequences and to sequences obtained using a previously published Illumina MiSeq workflow.

Results

For all egg-grown isolates, NDV was detected and virulence and genotype were accurately predicted. For clinical samples, NDV was detected in ten of eleven NDV samples. Six of the clinical samples contained two mixed genotypes as determined by MiSeq, of which the MinION method detected both genotypes in four samples. Additionally, testing a dilution series of one NDV isolate resulted in NDV detection in a dilution as low as 10¹ 50% egg infectious dose per milliliter. This was accomplished in as little as 7 min of sequencing time, with a 98.37% sequence identity compared to the expected consensus obtained by MiSeq.

Conclusion

The depth of sequencing, fast sequencing capabilities, accuracy of the consensus sequences, and the low cost of multiplexing allowed for effective virulence prediction and genotype identification of NDVs currently circulating worldwide. The sensitivity of this protocol was preliminary tested using only one genotype. After more extensive evaluation of the sensitivity and specificity, this protocol will likely be applicable to the detection and characterization of NDV.

     

The partitioned LASSO-patternsearch algorithm with application to gene expression data

Background

The generation of multiple sequence alignments (MSAs) is a crucial step for many bioinformatic analyses. Thus improving MSA accuracy and identifying potential errors in MSAs is important for a wide range of post-genomic research. We present a novel method called MergeAlign which constructs consensus MSAs from multiple independent MSAs and assigns an alignment precision score to each column.

Results

Using conventional benchmark tests we demonstrate that on average MergeAlign MSAs are more accurate than MSAs generated using any single matrix of sequence substitution. We show that MergeAlign column scores are related to alignment precision and hence provide an ab initio method of estimating alignment precision in the absence of curated reference MSAs. Using two novel and independent alignment performance tests that utilise a large set of orthologous gene families we demonstrate that increasing MSA performance leads to an increase in the performance of downstream phylogenetic analyses.

Conclusion

Using multiple tests of alignment performance we demonstrate that this novel method has broad general application in biological research.     

HAlign-II: efficient ultra-large multiple sequence alignment and phylogenetic tree reconstruction with distributed and parallel computing

Background

Multiple sequence alignment (MSA) plays a key role in biological sequence analyses, especially in phylogenetic tree construction. Extreme increase in next-generation sequencing results in shortage of efficient ultra-large biological sequence alignment approaches for coping with different sequence types.

Methods

Distributed and parallel computing represents a crucial technique for accelerating ultra-large (e.g. files more than 1 GB) sequence analyses. Based on HAlign and Spark distributed computing system, we implement a highly cost-efficient and time-efficient HAlign-II tool to address ultra-large multiple biological sequence alignment and phylogenetic tree construction.

Results

The experiments in the DNA and protein large scale data sets, which are more than 1GB files, showed that HAlign II could save time and space. It outperformed the current software tools. HAlign-II can efficiently carry out MSA and construct phylogenetic trees with ultra-large numbers of biological sequences. HAlign-II shows extremely high memory efficiency and scales well with increases in computing resource.

Conclusions

THAlign-II provides a user-friendly web server based on our distributed computing infrastructure. HAlign-II with open-source codes and datasets was established at http://lab.malab.cn/soft/halign.

     

Hepatitis C virus genotype 3a with phylogenetically distinct origin is circulating in Pakistan

Background

Hepatitis C virus (HCV) is one of the leading causes of viral hepatitis worldwide and its genotype 3a is predominant in vast areas of Pakistan.

Findings

The present study reports the first full sequence of HCV 3a isolate PK-1 from Pakistan. This nucleotide sequence was compared with six other HCV genotype 3a full length sequences from different regions of the world by using statistical methods of phylogenetic analysis.

Conclusion

The nucleotide difference of these seven sequences shows that HCV genotype 3a of phylogenetically distinct origin is circulating in Pakistan.     

In vivo - in vitro toxicogenomic comparison of TCDD-elicited gene expression in Hepa1c1c7 mouse hepatoma cells and C57BL/6 hepatic tissue

Background

By comparing the quail genome with that of chicken, chromosome rearrangements that have occurred in these two galliform species over 35 million years of evolution can be detected. From a more practical point of view, the definition of conserved syntenies helps to predict the position of genes in quail, based on information taken from the chicken sequence, thus enhancing the utility of this species in biological studies through a better knowledge of its genome structure. A microsatellite and an Amplified Fragment Length Polymorphism (AFLP) genetic map were previously published for quail, as well as comparative cytogenetic data with chicken for macrochromosomes. Quail genomics will benefit from the extension and the integration of these maps.

Results

The integrated linkage map presented here is based on segregation analysis of both anonymous markers and functional gene loci in 1,050 quail from three independent F2 populations. Ninety-two loci are resolved into 14 autosomal linkage groups and a Z chromosome-specific linkage group, aligned with the quail AFLP map. The size of linkage groups ranges from 7.8 cM to 274.8 cM. The total map distance covers 904.3 cM with an average spacing of 9.7 cM between loci. The coverage is not complete, as macrochromosome CJA08, the gonosome CJAW and 23 microchromosomes have no marker assigned yet. Significant sequence identities of quail markers with chicken enabled the alignment of the quail linkage groups on the chicken genome sequence assembly. This, together with interspecific Fluorescence In Situ Hybridization (FISH), revealed very high similarities in marker order between the two species for the eight macrochromosomes and the 14 microchromosomes studied.

Conclusion

Integrating the two microsatellite and the AFLP quail genetic maps greatly enhances the quality of the resulting information and will thus facilitate the identification of Quantitative Trait Loci (QTL). The alignment with the chicken chromosomes confirms the high conservation of gene order that was expected between the two species for macrochromosomes. By extending the comparative study to the microchromosomes, we suggest that a wealth of information can be mined in chicken, to be used for genome analyses in quail.     

Monoclonal antibodies to hemagglutinin-neuraminidase and fusion glycoproteins of Newcastle disease virus: relationship between glycosylation and reactivity.

Eighteen hybridoma lines obtained by immunization of mice with Newcastle disease virus (NDV) lentogenic strain La Sota or velogenic strain Italien produced hemagglutinating monoclonal antibodies. The 18 monoclones were divided into four groups according to their reactivity toward native hemagglutinin neuraminidase protein (HN), nonglycosylated HN precursor, and heat-denatured HN blotted on nitrocellulose membranes. Only group II reagents were reactive toward their targets in all conditions tested. They were considered sequence-specific antibodies. Group I antibodies did not require glycosylation but lacked reactivity towards the denatured glycosylated antigen. Monoclonal antibodies from group III recognized only the native HN. Group IV was made up of a single monoclone that lacked reactivity with NDV Italien but recognized the La Sota strain in hemagglutination inhibition and enzyme-linked immunosorbent assays. Five hybridoma lines produced monoclonal antibodies which neutralized viral infectivity but failed to inhibit hemagglutination. One monoclonal antibody obtained after immunization of mice with NDV La Sota showed a low neutralization index versus NDV Italien. Four monoclonal antibodies derived from mice immunized with NDV Italien showed higher neutralization indices towards this strain. Neither the denatured F protein nor its nonglycosylated precursor was reacted against by the five monoclonal antibodies.     

Mapping the sequence mutations of the 2009 H1N1 influenza A virus neuraminidase relative to drug and antibody binding sites

 

In this work, we study the consequences of sequence variations of the "2009 H1N1" (swine or Mexican flu) influenza A virus strain neuraminidase for drug treatment and vaccination. We find that it is phylogenetically more closely related to European H1N1 swine flu and H5N1 avian flu rather than to the H1N1 counterparts in the Americas. Homology-based 3D structure modeling reveals that the novel mutations are preferentially located at the protein surface and do not interfere with the active site. The latter is the binding cavity for 3 currently used neuraminidase inhibitors: oseltamivir (Tamiflu^®), zanamivir (Relenza^®) and peramivir; thus, the drugs should remain effective for treatment. However, the antigenic regions of the neuraminidase relevant for vaccine development, serological typing and passive antibody treatment can differ from those of previous strains and already vary among patients.

Reviewers

This article was reviewed by Sandor Pongor and L. Aravind.     

Aberrant antigenic expression in extranodal NK/T-cell lymphoma: a multi-parameter study from Thailand

Background

Extraosseous plasmacytoma, so called extramedullary plasmacytoma (EMP) is relatively rare in China. The aim was investigate the clinicopathologic features of EMP and the role of Immunophenotype and genotype detection in diagnosis of EMP.

Methods

Thirty-two cases of EMP were investigated retrospectively by histopathology, immunophenotype, genotype and survival analysis.

Results

Clinically, the mean age of the patients was 53.4. Most of the patients received no treatment after the diagnosis was established, and the prognosis was relatively poor. Histologically, in 40% of the cases, the neoplastic cells were grade II or III. The neoplastic cells expressed one or more PC associated antigens. The immunophenotype of EMP and inflammation of sinonasal regions with numerous PC infiltrations were compared and showed some difference in expression of CD45, CD27, CD44v6 and Bcl-2 as well. Ig light chain restriction was detected in 87.5% of the cases.

Conclusions

we described 32 Chinese cases of EMP, compare with that reported in the literature, some differences are presented, including higher percentage of grade II and III cases, clinically inconsistent treatment and management as well as poor outcome of the disease.     

Generation by Reverse Genetics of an Effective,Stable, Live-Attenuated Newcastle Disease Virus Vaccine Based on a Currently Circulating,Highly Virulent Indonesian Strain

Newcastle disease virus (NDV) can cause severe disease in chickens. Although NDV vaccines exist, there are frequent reports of outbreaks in vaccinated chickens. During 2009–2010, despite intense vaccination, NDV caused major outbreaks among commercial poultry farms in Indonesia. These outbreaks raised concern regarding the protective immunity of current vaccines against circulating virulent strains in Indonesia. In this study, we investigated whether a recombinant attenuated Indonesian NDV strain could provide better protection against prevalent Indonesian viruses. A reverse genetics system for the highly virulent NDV strain Banjarmasin/010/10 (Ban/010) isolated in Indonesia in 2010 was constructed. The Ban/010 virus is classified in genotype VII of class II NDV, which is genetically distinct from the commercial vaccine strains B1 and LaSota, which belong to genotype II, and shares only 89 and 87% amino acid identity for the protective antigens F and HN, respectively. A mutant virus, named Ban/AF, was developed in which the virulent F protein cleavage site motif “RRQKR↓F” was modified to an avirulent motif “GRQGR↓L” by three amino acid substitutions (underlined). The Ban/AF vaccine virus did not produce syncytia or plaques in cell culture, even in the presence of added protease. Pathogenicity tests showed that Ban/AF was completely avirulent. Ban/AF replicated efficiently during 10 consecutive passages in chickens and remained genetically stable. Serological analysis showed that Ban/AF induced higher neutralization and hemagglutination inhibition antibody titers against the prevalent viruses than the commercial vaccines B1 or LaSota. Both Ban/AF and commercial vaccines provided protection against clinical disease and mortality after challenge with virulent NDV strain Ban/010 (genotype VII) or GB Texas (genotype II). However, Ban/AF significantly reduced challenge virus shedding from the vaccinated birds compared to B1 vaccine. These results suggest that Ban/AF can provide better protection than commercial vaccines and is a promising vaccine candidate against NDV strains circulating in Indonesia.     

High-throughput peptide quantification using mTRAQ reagent triplex

Background

Continuing research into the global multiple sequence alignment problem has resulted in more sophisticated and principled alignment methods. Unfortunately these new algorithms often require large amounts of time and memory to run, making it nearly impossible to run these algorithms on large datasets. As a solution, we present two general methods, Crumble and Prune, for breaking a phylogenetic alignment problem into smaller, more tractable sub-problems. We call Crumble and Prune meta-alignment methods because they use existing alignment algorithms and can be used with many current alignment programs. Crumble breaks long alignment problems into shorter sub-problems. Prune divides the phylogenetic tree into a collection of smaller trees to reduce the number of sequences in each alignment problem. These methods are orthogonal: they can be applied together to provide better scaling in terms of sequence length and in sequence depth. Both methods partition the problem such that many of the sub-problems can be solved independently. The results are then combined to form a solution to the full alignment problem.

Results

Crumble and Prune each provide a significant performance improvement with little loss of accuracy. In some cases, a gain in accuracy was observed. Crumble and Prune were tested on real and simulated data. Furthermore, we have implemented a system called Job-tree that allows hierarchical sub-problems to be solved in parallel on a compute cluster, significantly shortening the run-time.

Conclusions

These methods enabled us to solve gigabase alignment problems. These methods could enable a new generation of biologically realistic alignment algorithms to be applied to real world, large scale alignment problems.     

不同基因型新城疫病毒HN基因保守区域的克隆表达及免疫活性的比较

对于NDV不同基因型毒株的HN基因的氨基酸序列比较后发现,基因Ⅶ型的毒株在第65~75位的氨基酸序列较为保守,而其他各基因型在该区域则不尽相同。因此本文克隆了NDV基因Ⅱ、Ⅶ、Ⅸ型的代表毒株La Sota、GX-2、F48E9的含有该区域的序列,并进行蛋白表达,通过特异性抗血清对表达肽段进行抗原性测定。应用表达的蛋白免疫SPF鸡,用ELISA检测各组的抗体水平。结果表明3个毒株在反应原性上存在差异。应用NDV F48E9强毒株攻毒,结果显示各多肽蛋白的保护率不同,表明在该区域这3个毒株之间存在着免疫原性差异。     

Sargassum fusiforme fraction is a potent and specific inhibitor of HIV-1 fusion and reverse transcriptase

Background

Hepatitis C virus (HCV) circulates in an infected individual as a heterogeneous mixture of closely related viruses called quasispecies. The E1/E2 region of the HCV genome is hypervariable (HVR1) and is targeted by the humoral immune system. Hepatitis C virions are found in two forms: antibody associated or antibody free. The objective of this study was to investigate if separation of Hepatitis C virions into antibody enriched and antibody depleted fractions segregates quasispecies populations into distinctive swarms.

Results

A HCV genotype 4a specimen was fractionated into IgG-depleted and IgG-enriched fractions by use of Albumin/IgG depletion spin column. Clonal analysis of these two fractions was performed and then compared to an unfractionated sample. Following sequence analysis it was evident that the antibody depleted fraction was significantly more heterogeneous than the antibody enriched fraction, revealing a unique quasispecies profile. An in-frame 3 nt insertion was observed in 26% of clones in the unfractionated population and in 64% of clones in the IgG-depleted fraction. In addition, an in-frame 3 nt indel event was observed in 10% of clones in the unfractionated population and in 9% of clones in the IgG-depleted fraction. Neither of these latter events, which are rare occurrences in genotype 4a, was identified in the IgG-enriched fraction.

Conclusion

In conclusion, the homogeneity of the IgG-enriched species is postulated to represent a sequence that was strongly recognised by the humoral immune system at the time the sample was obtained. The heterogeneous nature of the IgG-depleted fraction is discussed in the context of humoral escape.     

Comparison of solution-based exome capture methods for next generation sequencing

Background

Techniques enabling targeted re-sequencing of the protein coding sequences of the human genome on next generation sequencing instruments are of great interest. We conducted a systematic comparison of the solution-based exome capture kits provided by Agilent and Roche NimbleGen. A control DNA sample was captured with all four capture methods and prepared for Illumina GAII sequencing. Sequence data from additional samples prepared with the same protocols were also used in the comparison.

Results

We developed a bioinformatics pipeline for quality control, short read alignment, variant identification and annotation of the sequence data. In our analysis, a larger percentage of the high quality reads from the NimbleGen captures than from the Agilent captures aligned to the capture target regions. High GC content of the target sequence was associated with poor capture success in all exome enrichment methods. Comparison of mean allele balances for heterozygous variants indicated a tendency to have more reference bases than variant bases in the heterozygous variant positions within the target regions in all methods. There was virtually no difference in the genotype concordance compared to genotypes derived from SNP arrays. A minimum of 11× coverage was required to make a heterozygote genotype call with 99% accuracy when compared to common SNPs on genome-wide association arrays.

Conclusions

Libraries captured with NimbleGen kits aligned more accurately to the target regions. The updated NimbleGen kit most efficiently covered the exome with a minimum coverage of 20×, yet none of the kits captured all the Consensus Coding Sequence annotated exons.     

Association studies using family pools of outcrossing crops based on allele-frequency estimates from DNA sequencing

Key message

We propose a method in which GBS data can be conveniently analyzed without calling genotypes.

Abstract

F2 families are frequently used in breeding of outcrossing species, for instance to obtain trait measurements on plots. We propose to perform association studies by obtaining a matching “family genotype” from sequencing a pooled sample of the family, and to directly use allele frequencies computed from sequence read-counts for mapping. We show that, under additivity assumptions, there is a linear relationship between the family phenotype and family allele frequency, and that a regression of family phenotype on family allele frequency will estimate twice the allele substitution effect at a locus. However, medium-to-low sequencing depth causes underestimation of the true allele substitution effect. An expression for this underestimation is derived for the case that parents are diploid, such that F2 families have up to four dosages of every allele. Using simulation studies, estimation of the allele effect from F2-family pools was verified and it was shown that the underestimation of the allele effect is correctly described. The optimal design for an association study when sequencing budget would be fixed is obtained using large sample size and lower sequence depth, and using higher SNP density (resulting in higher LD with causative mutations) and lower sequencing depth. Therefore, association studies using genotyping by sequencing are optimal and use low sequencing depth per sample. The developed framework for association studies using allele frequencies from sequencing can be modified for other types of family pools and is also directly applicable for association studies in polyploids.     

Genome skimming identifies polymorphism in tern populations and species

Background

Large-scale protein structure alignment, an indispensable tool to structural bioinformatics, poses a tremendous challenge on computational resources. To ensure structure alignment accuracy and efficiency, efforts have been made to parallelize traditional alignment algorithms in grid environments. However, these solutions are costly and of limited accessibility. Others trade alignment quality for speedup by using high-level characteristics of structure fragments for structure comparisons.

Findings

We present ppsAlign, a parallel protein structure Alignment framework designed and optimized to exploit the parallelism of Graphics Processing Units (GPUs). As a general-purpose GPU platform, ppsAlign could take many concurrent methods, such as TM-align and Fr-TM-align, into the parallelized algorithm design. We evaluated ppsAlign on an NVIDIA Tesla C2050 GPU card, and compared it with existing software solutions running on an AMD dual-core CPU. We observed a 36-fold speedup over TM-align, a 65-fold speedup over Fr-TM-align, and a 40-fold speedup over MAMMOTH.

Conclusions

ppsAlign is a high-performance protein structure alignment tool designed to tackle the computational complexity issues from protein structural data. The solution presented in this paper allows large-scale structure comparisons to be performed using massive parallel computing power of GPU.     

Immunoproteomic analysis of bacterial proteins of Actinobacillus pleuropneumoniae serotype 1

Background

Actinobacillus pleuropneumoniae (APP) is one of the most important swine pathogens worldwide. Identification and characterization of novel antigenic APP vaccine candidates are underway. In the present study, we use an immunoproteomic approach to identify APP protein antigens that may elicit an immune response in serotype 1 naturally infected swine and serotype 1 virulent strain S259-immunized rabbits.

Results

Proteins from total cell lysates of serotype 1 APP were separated by two-dimensional electrophoresis (2DE). Western blot analysis revealed 21 immunoreactive protein spots separated in the pH 4-7 range and 4 spots in the pH 7-11 range with the convalescent sera from swine; we found 5 immunoreactive protein spots that separated in the pH 4-7 range and 2 in the pH 7-11 range with hyperimmune sera from S259-immunized rabbits. The proteins included the known antigens ApxIIA, protective surface antigen D15, outer membrane proteins P5, subunit NqrA. The remaining antigens are being reported as immunoreactive proteins in APP for the first time, to our knowledge.

Conclusions

We identified a total of 42 immunoreactive proteins of the APP serotype 1 virulent strain S259 which represented 32 different proteins, including some novel immunoreactive factors which could be researched as vaccine candidates.