Assessing the Potential of Multiple Imputation in Sequence Analysis

Sequence analysis and optimal matching are useful heuristic tools for the descriptive analysis of heterogeneous individual pathways such as educational careers, job sequences or patterns of family formation. However, to date it remains unclear how to handle the inevitable problems caused by missing values with regard to such analysis. Multiple Imputation (MI) offers a possible solution for this problem but it has not been tested in the context of sequence analysis. Against this background, we contribute to the literature by assessing the potential of MI in the context of sequence analyses using an empirical example. Methodologically, we draw upon the work of Brendan Halpin and extend it to additional types of missing value patterns. Our empirical case is a sequence analysis of panel data with substantial attrition that examines the typical patterns and the persistence of sex segregation in school-to-work transitions in Switzerland. The preliminary results indicate that MI is a valuable methodology for handling missing values due to panel mortality in the context of sequence analysis. MI is especially useful in facilitating a sound interpretation of the resulting sequence types.

Novel approach for genotyping varicella-zoster virus strains from Germany

In this study, we present a novel genotyping scheme to classify German wild-type varicella-zoster virus (VZV) strains and to differentiate them from the Oka vaccine strain (genotype B). This approach is based on analysis of four loci in open reading frames (ORFs) 51 to 58, encompassing a total length of 1,990 bp. The new genotyping scheme produced identical clusters in phylogenetic analyses compared to full-genome sequences from well-characterized VZV strains. Based on genotype A, D, B, and C reference strains, a dichotomous identification key (DIK) was developed and applied for VZV strains obtained from vesicle fluid and liquor samples originating from 42 patients suffering from varicella or zoster between 2003 and 2006. Sequencing of regions in ORFs 51, 52, 53, 56, 57, and 58 identified 18 single-nucleotide polymorphisms (SNPs), including two novel ones, SNP 89727 and SNP 92792 in ORF51 and ORF52, respectively. The DIK as well as phylogenetic analysis by Bayesian inference showed that 14 VZV strains belonged to genotype A, and 28 VZV strains were classified as genotype D. Neither Japanese (vaccine)-like B strains nor recombinant-like C strains were found within the samples from Germany. The novel genotyping scheme and the DIK were demonstrated to be practical and simple and allow the highly efficient replication of phylogenetic patterns in VZV initially derived from full-genome DNA sequence analyses. Therefore, this approach may allow us to draw a more comprehensive picture of wild-type VZV strains circulating in Germany and Central Europe by high-throughput procedures in the future.

Experimental adaptation of wild-type canine distemper virus (CDV) to the human entry receptor CD150.

Canine distemper virus (CDV), a close relative of measles virus (MV), is widespread and well known for its broad host range. When the goal of measles eradication may be achieved, and when measles vaccination will be stopped, CDV might eventually cross the species barrier to humans and emerge as a new human pathogen. In order to get an impression how fast such alterations may occur, we characterized required adaptive mutations to the human entry receptors CD150 (SLAM) and nectin-4 as first step to infect human target cells. Recombinant wild-type CDV-A75/17(red) adapted quickly to growth in human H358 epithelial cells expressing human nectin-4. Sequencing of the viral attachment proteins (hemagglutinin, H, and fusion protein, F) genes revealed that no adaptive alteration was required to utilize human nectin-4. In contrast, the virus replicated only to low titres (10(2) pfu/ml) in Vero cells expressing human CD150 (Vero-hSLAM). After three passages using these cells virus was adapted to human CD150 and replicated to high titres (10(5) pfu/ml). Sequence analyses revealed that only one amino acid exchange in the H-protein at position 540 Asp→Gly (D540G) was required for functional adaptation to human CD150. Structural modelling suggests that the adaptive mutation D540G in H reflects the sequence alteration from canine to human CD150 at position 70 and 71 from Pro to Leu (P70L) and Gly to Glu (G71E), and compensates for the gain of a negative charge in the human CD150 molecule. Using this model system our data indicate that only a minimal alteration, in this case one adaptive mutation, is required for adaptation of CDV to the human entry receptors, and help to understand the molecular basis why this adaptive mutation occurs.

HLA-B⁎27 subtype specificity determines targeting and viral evolution of a hepatitis C virus-specific CD8+ T cell epitope

Background & Aims: HLA-B⁄27 is associated with spontaneous HCV genotype 1 clearance. HLA-B⁄27-restricted CD8+ T cells target three NS5B epitopes. Two of these epitopes are dominantly targeted in the majority of HLA-B⁄27+ patients. In chronic infection, viral escape occurs consistently in these two epitopes. The third epitope (NS5B2820) was dominantly targeted in an acutely infected patient. This was in contrast, however, to the lack of recognition and viral escape in the large majority of HLA-B⁄27+ patients. Here, we set out to determine the host factors contributing to selective targeting of this epitope. Methods: Four-digit HLA class I typing and viral sequence analyses were performed in 78 HLA-B⁄27+ patients with chronic HCV genotype 1 infection. CD8+ T cell analyses were performed in a subset of patients. In addition, HLA/peptide affinity was compared for HLA-B⁄27:02 and 05. Results: The NS5B2820 epitope is only restricted by the HLA-B⁄27 subtype HLA-B⁄27:02 (that is frequent in Mediterranean populations), but not by the prototype HLA-B⁄27 subtype B⁄27:05. Indeed, the epitope is very dominant in HLA-B⁄27:02+ patients and is associated with viral escape mutations at the anchor position for HLA-binding in 12 out of 13 HLA-B⁄27:02+ chronically infected patients. Conclusions: The NS5B2820 epitope is immunodominant in the context of HLA-B⁄27:02, but is not restricted by other HLA-B⁄27 subtypes. This finding suggests an important role of HLA subtypes in the restriction of HCV-specific CD8+ responses. With minor HLA subtypes covering up to 39% of specific populations, these findings may have important implications for the selection of epitopes for global vaccines.

Echinococcus P29 antigen: molecular characterization and implication on post-surgery follow-up of CE patients infected with different species of the Echinococcus granulosus complex.

The protein P29 is a potential serological marker for post-treatment monitoring of cystic echinococcosis (CE) especially in young patients. We now have demonstrated that P29 is encoded in the Echinococcus genus by a single gene consisting of 7 exons spanning 1.2 kb of DNA. Variability of the p29 gene at inter- and intra-species level was assessed with 50 cDNA and 280 genomic DNA clones isolated from different E. granulosus s.l. isolates (E. granulosus sensu stricto (G1), E. equinus (G4), E. ortleppi (G5), E. canadensis (G6), E. canadensis (G7) and E. canadensis (G10)) as well as four E. multilocularis isolates. Scarce interspecies polymorphism at the p29 locus was observed and affected predominantly E. granulosus s.s. (G1), where we identified two alleles (A1 and A2) coding for identical P29 proteins and yielding in three genotypes (A1/A1, A2/A2 and A1/A2). Genotypic frequencies expected under Hardy-Weinberg equilibrium revealed a high rate of heterozygosity (47%) that strongly supports the hypothesis that E. granulosus s.s. (G1) is predominantly outbreeding. Comparative sequence analyses of the complete p29 gene showed that phylogenetic relationships within the genus Echinococcus were in agreement with those of previous nuclear gene studies. At the protein level, the deduced P29 amino acid (AA) sequences exhibited a high level of conservation, ranging from 97.9% AA sequence identity among the whole E. granulosus s.l. group to 99.58% identity among E. multilocularis isolates. We showed that P29 proteins of these two species differ by three AA substitutions without implication for antigenicity. In Western-blot analyses, serum antibodies from a human CE patient infected with E. canadensis (G6) strongly reacted with recombinant P29 from E. granulosus s.s. (G1) (recEg(G1)P29). In the same line, human anti-Eg(G1)P29 antibodies bound to recEcnd(G6)P29. Thus, minor AA sequence variations appear not to impair the prognostic serological use of P29.

A 1.8-kb insertion in the 3'-UTR of RXFP2 is associated with polledness in sheep

Sheep breeds show a broad spectrum of different horn phenotypes. In most modern production breeds, sheep are polled (absence of horns), whereas horns occur mainly in indigenous breeds. Previous studies mapped the responsible locus to the region of the RXFP2 gene on ovine chromosome 10. A 4-kb region of the 3'-end of RXFP2 was amplified in horned and polled animals from seven Swiss sheep breeds. Sequence analysis identified a 1833-bp genomic insertion located in the 3'-UTR region of RXFP2 present in polled animals only. An efficient PCR-based genotyping method to determine the polled genotype of individual sheep is presented. Comparative sequence analyses revealed evidence that the polled-associated insertion adds a potential antisense RNA sequence of EEF1A1 to the 3'-end of RXFP2 transcripts.

Polydnavirus DNA of the braconid wasp Chelonus inanitus is integrated in the wasp's genome and excised only in later pupal and adult stages of the female

Many endoparasitic wasps inject, along with the egg, polydnavirus into their insect hosts, the virus being a prerequisite for successful parasitoid development. The genome of polydnaviruses consists of multiple circular dsDNA molecules of variable size. We show for a 12 kbp segment of the braconid Chelonus inanitus (CiV12) that it is integrated into the wasp genome. This is the first direct demonstration of integration for a bracovirus. PCR data indicated that the integrated form of CiV12 was present in all male and female stages investigated while the excised circular virus DNA only appeared in females after a specific stage in pupal-adult development. The data also indicated that after excision of virus DNA the genomic DNA was rejoined. This has not yet been reported for any polydnavirus. Sequence analyses in the junction regions revealed the presence of an imperfect consensus sequence of 15 nucleotides in CiV12, in each terminus of the integrated virus DNA and in the rejoined genomic DNA. Within these repeats two sequence types (ATA, TAC) were observed in the various virus clones and in the clones encompassing the rejoined genomic DNA; they corresponded to the sequence type in the right and left junction, respectively. To explain this, we propose a model of virus DNA replication in which the genomic DNA is folded to juxtapose the direct repeat of the left with that of the right junction; recombination at specific sites would then yield the two types of virus and rejoined genomic DNA.

Geochemical and Microbiological Studies of Nitrous Oxide Variations within the New NEEM Greenland Ice Core During the Last Glacial Period

Deep polar ice cores provide atmospheric records of nitrous oxide (N₂O) and other trace gases reflecting climate history along with a parallel archive of microbial cells transported with mineral dust, marine and volcanic aerosols from around the globe. Our interdisciplinary study of 32 samples from different depths of the recently drilled NEEM Greenland ice core addressed the question whether the identified microorganisms were capable of post-depositional biological production of N₂O in situ. We used high-resolution geochemical and microbiological approaches to examine the N₂O concentrations, the quantitative distributions of dust, Ca⁺², NH₄⁺ and NO₃⁻ ¡ons related to N cycle pathways, the microbial abundance and diversity at specific NEEM core depths from 1758 m to 1867.8 m. Results showed varying concentrations of N₂O (220 –271.5 ppb). Microbial abundance fluctuated between 3.3 x 10⁴ and 3.3 x 10⁶ cells mL⁻¹ in direct correlation with dust and Ca²⁺ concentrations with higher cell numbers deposited during colder periods. The average values of NH₄⁺ and NO₃⁻ indicated that substrates were available for the microorganisms capable of utilizing them. PCR amplification of selected functional genes involved in bacterial and archaeal nitrification and denitrification was not successful. Sanger and Illumina MiSeq sequence analyses of SSU rRNA genes showed variable representation of Alpha-, Beta- and Gammaproteobacteria, Firmicutes, Actinobacteria, chloroplasts and fungi. The metabolic potential of the dominant genera of Proteobacteria and Firmicutes as possible N₂O producers suggested that denitrification activity may have led to in-situ production and accumulation of N₂O.

The venom composition of the parasitic wasp Chelonus inanitus resolved by combined expressed sequence tags analysis and proteomic approach

Background Parasitic wasps constitute one of the largest group of venomous animals. Although some physiological effects of their venoms are well documented, relatively little is known at the molecular level on the protein composition of these secretions. To identify the majority of the venom proteins of the endoparasitoid wasp Chelonus inanitus (Hymenoptera: Braconidae), we have randomly sequenced 2111 expressed sequence tags (ESTs) from a cDNA library of venom gland. In parallel, proteins from pure venom were separated by gel electrophoresis and individually submitted to a nano-LC-MS/MS analysis allowing comparison of peptides and ESTs sequences. Results About 60% of sequenced ESTs encoded proteins whose presence in venom was attested by mass spectrometry. Most of the remaining ESTs corresponded to gene products likely involved in the transcriptional and translational machinery of venom gland cells. In addition, a small number of transcripts were found to encode proteins that share sequence similarity with well-known venom constituents of social hymenopteran species, such as hyaluronidase-like proteins and an Allergen-5 protein. An overall number of 29 venom proteins could be identified through the combination of ESTs sequencing and proteomic analyses. The most highly redundant set of ESTs encoded a protein that shared sequence similarity with a venom protein of unknown function potentially specific of the Chelonus lineage. Venom components specific to C. inanitus included a C-type lectin domain containing protein, a chemosensory protein-like protein, a protein related to yellow-e3 and ten new proteins which shared no significant sequence similarity with known sequences. In addition, several venom proteins potentially able to interact with chitin were also identified including a chitinase, an imaginal disc growth factor-like protein and two putative mucin-like peritrophins. Conclusions The use of the combined approaches has allowed to discriminate between cellular and truly venom proteins. The venom of C. inanitus appears as a mixture of conserved venom components and of potentially lineage-specific proteins. These new molecular data enrich our knowledge on parasitoid venoms and more generally, might contribute to a better understanding of the evolution and functional diversity of venom proteins within Hymenoptera.

Sequence capture and next-generation resequencing of multiple tagged nucleic acid samples for mutation screening of urea cycle disorders

Molecular genetic testing is commonly used to confirm clinical diagnoses of inherited urea cycle disorders (UCDs); however, conventional mutation screenings encompassing only the coding regions of genes may not detect disease-causing mutations occurring in regulatory elements and introns. Microarray-based target enrichment and next-generation sequencing now allow more-comprehensive genetic screening. We applied this approach to UCDs and combined it with the use of DNA bar codes for more cost-effective, parallel analyses of multiple samples.

Genotype and antibiotic resistance analyses of Campylobacter isolates from ceca and carcasses of slaughtered broiler flocks

To obtain genetic information about Campylobacter jejuni and Campylobacter coli from broilers and carcasses at slaughterhouses, we analyzed and compared 340 isolates that were collected in 2008 from the cecum right after slaughter or from the neck skin after processing. We performed rpoB sequence-based identification, multilocus sequence typing (MLST), and flaB sequence-based typing; we additionally analyzed mutations within the 23S rRNA and gyrA genes that confer resistance to macrolide and quinolone antibiotics, respectively. The rpoB-based identification resulted in a distribution of 72.0% C. jejuni and 28.0% C. coli. The MLST analysis revealed that there were 59 known sequence types (STs) and 6 newly defined STs. Most of the STs were grouped into 4 clonal complexes (CC) that are typical for poultry (CC21, CC45, CC257, and CC828), and these represented 61.8% of all of the investigated isolates. The analysis of 95 isolates from the cecum and from the corresponding carcass neck skin covered 44 different STs, and 54.7% of the pairs had matching genotypes. The data indicate that cross-contamination from various sources during slaughter may occur, although the majority of Campylobacter contamination on carcasses appeared to originate from the slaughtered flock itself. Mutations in the 23S rRNA gene were found in 3.1% of C. coli isolates, although no mutations were found in C. jejuni isolates. Mutations in the gyrA gene were observed in 18.9% of C. jejuni and 26.8% of C. coli isolates, which included two C. coli strains that carried mutations conferring resistance to both classes of antibiotics. A relationship between specific genotypes and antibiotic resistance/susceptibility was observed.

Sequence variations in the von Hippel-Lindau tumor suppressor gene in patients with intracranial aneurysms

The rupture of intracranial aneurysms leads to subarachnoid hemorrhage, which is often associated with poor outcome. Preventive treatment of unruptured intracranial aneurysms is possible and recommended. However, the lack of candidate genes precludes identifying patients at risk by genetic analyses. We observed intracranial aneurysms in 2 patients with von Hippel-Lindau (VHL) disease and the known disease-causing mutation c.292T > C (p.Tyr98His) in the VHL tumor suppressor gene. This study investigates whether the VHL gene is a possible candidate gene for aneurysm formation.

The mutation causing the black-and-tan pigmentation phenotype of Mangalitza pigs maps to the porcine ASIP locus but does not affect its coding sequence

The gene for agouti signaling protein (ASIP) is centrally involved in the expression of coat color traits in animals. The Mangalitza pig breed is characterized by a black-and-tan phenotype with black dorsal pigmentation and yellow or white ventral pigmentation. We investigated a Mangalitza x Piétrain cross and observed a coat color segregation pattern in the F2 generation that can be explained by virtue of two alleles at the MC1R locus and two alleles at the ASIP locus. Complete linkage of the black-and-tan phenotype to microsatellite alleles at the ASIP locus on SSC 17q21 was observed. Corroborated by the knowledge of similar mouse coat color mutants, it seems therefore conceivable that the black-and-tan pigmentation of Mangalitza pigs is caused by an ASIP allele a(t), which is recessive to the wild-type allele A. Toward positional cloning of the a(t) mutation, a 200-kb genomic BAC/PAC contig of this chromosomal region has been constructed and subsequently sequenced. Full-length ASIP cDNAs obtained by RACE differed in their 5' untranslated regions, whereas they shared a common open reading frame. Comparative sequencing of all ASIP exons and ASIP cDNAs between Mangalitza and Piétrain pigs did not reveal any differences associated with the coat color phenotype. Relative qRT-PCR analyses showed different dorsoventral skin expression intensities of the five ASIP transcripts in black-and-tan Mangalitza. The a(t) mutation is therefore probably a regulatory ASIP mutation that alters its dorsoventral expression pattern.

Identification and in silico analyses of novel TGFBR1 and TGFBR2 mutations in Marfan syndrome-related disorders

Very recently, heterozygous mutations in the genes encoding transforming growth factor beta receptors I (TGFBR1) and II (TGFBR2) have been reported in Loeys-Dietz aortic aneurysm syndrome (LDS). In addition, dominant TGFBR2 mutations have been identified in Marfan syndrome type 2 (MFS2) and familial thoracic aortic aneurysms and dissections (TAAD). In the past, mutations of these genes were associated with atherosclerosis and several human cancers. Here, we report a total of nine novel and one known heterozygous sequence variants in the TGFBR1 and TGFBR2 genes in nine of 70 unrelated individuals with MFS-like phenotypes who previously tested negative for mutations in the gene encoding the extracellular matrix protein fibrillin-1 (FBN1). To assess the pathogenic impact of these sequence variants, in silico analyses were performed by the PolyPhen, SIFT, and Fold-X algorithms and by means of a 3D homology model of the TGFBR2 kinase domain. Our results showed that in all but one of the patients the pathogenic effect of at least one sequence variant is highly probable (c.722C > T, c.799A > C, and c.1460G > A in TGFBR1 and c.773T > G, c.1106G > T, c.1159G > A, c.1181G > A, and c.1561T > C in TGFBR2). These deleterious alleles occurred de novo or segregated with the disease in the families, indicating a causative association between the sequence variants and clinical phenotypes. Since TGFBR2 mutations found in patients with MFS-related disorders cannot be distinguished from heterozygous TGFBR2 mutations reported in tumor samples, we emphasize the importance of segregation analysis in affected families. In order to be able to find the mutation that is indeed responsible for a MFS-related phenotype, we also propose that genetic testing for sequence alterations in TGFBR1 and TGFBR2 should be complemented by mutation screening of the FBN1 gene.