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.

8p23 beta-defensin copy number determination by single-locus pseudogene-based paralog ratio tests risk bias due to low-frequency sequence variations.

BACKGROUND The copy number variation (CNV) in beta-defensin genes (DEFB) on human chromosome 8p23 has been proposed to contribute to the phenotypic differences in inflammatory diseases. However, determination of exact DEFB CN is a major challenge in association studies. Quantitative real-time PCR (qPCR), paralog ratio tests (PRT) and multiplex ligation-dependent probe amplification (MLPA) have been extensively used to determine DEFB CN in different laboratories, but inter-method inconsistencies were observed frequently. In this study we asked which one is superior among the three methods for DEFB CN determination. RESULTS We developed a clustering approach for MLPA and PRT to statistically correlate data from a single experiment. Then we compared qPCR, a newly designed PRT and MLPA for DEFB CN determination in 285 DNA samples. We found MLPA had the best convergence and clustering results of the raw data and the highest call rate. In addition, the concordance rates between MLPA or PRT and qPCR (32.12% and 37.99%, respectively) were unacceptably low with underestimated CN by qPCR. Concordance rate between MLPA and PRT (90.52%) was high but PRT systematically underestimated CN by one in a subset of samples. In these samples a sequence variant which caused complete PCR dropout of the respective DEFB cluster copies was found in one primer binding site of one of the targeted paralogous pseudogenes. CONCLUSION MLPA is superior to PRT and even more to qPCR for DEFB CN determination. Although the applied PRT provides in most cases reliable results, such a test is particularly sensitive to low-frequency sequence variations preferably accumulating in loci like pseudogenes which are most likely not under selective pressure. In the light of the superior performance of multiplex assays, the drawbacks of such single PRTs could be overcome by combining more test markers.

Multiple regulatory variants located in cell type-specific enhancers within the PKP2 locus form major risk and protective haplotypes for canine atopic dermatitis in German shepherd dogs.

BACKGROUND Canine atopic dermatitis (CAD) is a chronic inflammatory skin disease triggered by allergic reactions involving IgE antibodies directed towards environmental allergens. We previously identified a ~1.5 Mb locus on canine chromosome 27 associated with CAD in German shepherd dogs (GSDs). Fine-mapping indicated association closest to the PKP2 gene encoding plakophilin 2. RESULTS Additional genotyping and association analyses in GSDs combined with control dogs from five breeds with low-risk for CAD revealed the top SNP 27:19,086,778 (p = 1.4 × 10(-7)) and a rare ~48 kb risk haplotype overlapping the PKP2 gene and shared only with other high-risk CAD breeds. We selected altogether nine SNPs (four top-associated in GSDs and five within the ~48 kb risk haplotype) that spanned ~280 kb forming one risk haplotype carried by 35 % of the GSD cases and 10 % of the GSD controls (OR = 5.1, p = 5.9 × 10(-5)), and another haplotype present in 85 % of the GSD cases and 98 % of the GSD controls and conferring a protective effect against CAD in GSDs (OR = 0.14, p = 0.0032). Eight of these SNPs were analyzed for transcriptional regulation using reporter assays where all tested regions exerted regulatory effects on transcription in epithelial and/or immune cell lines, and seven SNPs showed allelic differences. The DNA fragment with the top-associated SNP 27:19,086,778 displayed the highest activity in keratinocytes with 11-fold induction of transcription by the risk allele versus 8-fold by the control allele (pdifference = 0.003), and also mapped close (~3 kb) to an ENCODE skin-specific enhancer region. CONCLUSIONS Our experiments indicate that multiple CAD-associated genetic variants located in cell type-specific enhancers are involved in gene regulation in different cells and tissues. No single causative variant alone, but rather multiple variants combined in a risk haplotype likely contribute to an altered expression of the PKP2 gene, and possibly nearby genes, in immune and epithelial cells, and predispose GSDs to CAD.

Membrane disrupting antimicrobial peptide dendrimers with multiple amino termini

Antimicrobial peptide dendrimer H1 Leu8(Lys-Leu)4(Lys-Phe)2Lys-LysNH2 (Lys = branching lysine) was identified by screening a 6750-membered combinatorial library by the bead-diffusion assay. Sequence variations also revealed dendrimer bH1 Leu8(Dap-Leu)4(Dap-Phe)2Dap-LysNH2 (Dap = branching 2,3-diaminopropanoic acid) as a more potent analog. H1 and bH1 showed good antimicrobial activities mediated by membrane disruption (MIC = 2–4 μg mL−1 on Bacillus subtilis and Escherichia coli) but low hemolytic activity (MHC = 310 μg mL−1 respectively >2000 μg mL−1).

Rapid typing of Moraxella catarrhalis subpopulations based on outer membrane proteins using mass spectrometry

Moraxella catarrhalis is a major mucosal pathogen of the human respiratory tract both in children and in adults. Two subpopulations of this organism have been described that differ in 16S rRNA gene sequence and virulence traits. Three 16S rRNA types have been defined. 2-DE followed by protein identification by MS revealed significant differences in the outer membrane protein (OMP) patterns of each M. catarrhalis 16S rRNA type. Approximately 130 features were detected on the 2-DE map of each M. catarrhalis 16S rRNA type. However, only 50 features were expressed by all strains. Furthermore, direct profiling of isolated OMP using MALDI-TOF MS resulted in a characteristic spectral fingerprint for each 16S rRNA type. Fingerprints remained identical when intact cells instead of isolated OMP were analyzed. This finding suggests that the source of desorbed ions is the outer membrane. Based on the fingerprint we were able to assign 18 well-characterized clinical M. catarrhalis isolates to the correct subpopulation. Therefore, MALDI-TOF of intact M. catarrhalis provides a rapid and robust tool for M. catarrhalis strain typing that could be applied in epidemiological studies.

Diversity of Mycoplasma hyopneumoniae in pig farms revealed by direct molecular typing of clinical material

Mycoplasma hyopneumoniae is the etiological agent of enzootic pneumonia in swine. Various reports indicate that different strains are circulating in the swine population. We investigated the variety of M. hyopneumoniae strains by a newly developed genetic typing method based on the polyserine repeat motif of the LppS homolog P146. PCR amplification using M. hyopneumoniae specific, conserved primers flanking the region encoding the repeat motif, followed by sequencing and cluster analysis was carried out. The study included strains isolated from different geographic regions as well as lysates from lung swabs from a series of pig farms in Switzerland. High diversity of M. hyopneumoniae was observed but farms being in close geographic or operative contact generally seemed to be affected by the same strains. Moreover, analysis of multiple samples from single pig farms indicated that these harbored the same, farm-specific strain. The results indicate that multiple strains of M. hyopneumoniae are found in the swine population but that specific strains or clones are responsible for local outbreaks. The method presented is a highly reproducible epidemiologic tool allowing direct typing of M. hyopneumoniae from clinical material without prior isolation and cultivation of strains.

The locus for bovine dilated cardiomyopathy maps to chromosome 18

Bovine dilated cardiomyopathy (BDCMP) is a severe and terminal disease of the heart muscle observed in Holstein-Friesian cattle over the last 30 years. There is strong evidence for an autosomal recessive mode of inheritance for BDCMP. The objective of this study was to genetically map BDCMP, with the ultimate goal of identifying the causative mutation. A whole-genome scan using 199 microsatellite markers and one SNP revealed an assignment of BDCMP to BTA18. Fine-mapping on BTA18 refined the candidate region to the MSBDCMP06-BMS2785 interval. The interval containing the BDCMP locus was confirmed by multipoint linkage analysis using the software loki. The interval is about 6.7 Mb on the bovine genome sequence (Btau 3.1). The corresponding region of HSA19 is very gene-rich and contains roughly 200 genes. Although telomeric of the marker interval, TNNI3 is a possible positional and a functional candidate for BDCMP given its involvement in a human form of dilated cardiomyopathy. Sequence analysis of TNNI3 in cattle revealed no mutation in the coding sequence, but there was a G-to-A transition in intron 6 (AJ842179:c.378+315G>A). The analysis of this SNP using the study's BDCMP pedigree did not conclusively exclude TNNI3 as a candidate gene for BDCMP. Considering the high density of genes on the homologous region of HSA19, further refinement of the interval on BTA18 containing the BDCMP locus is needed.

Whole-body MRI of multiple myeloma: comparison of different MRI sequences in assessment of different growth patterns

PURPOSE: To determine sensitivity, specificity and inter-observer variability of different whole-body MRI (WB-MRI) sequences in patients with multiple myeloma (MM). METHODS AND MATERIALS: WB-MRI using a 1.5T MRI scanner was performed in 23 consecutive patients (13 males, 10 females; mean age 63+/-12 years) with histologically proven MM. All patients were clinically classified according to infiltration (low-grade, n=7; intermediate-grade, n=7; high-grade, n=9) and to the staging system of Durie and Salmon PLUS (stage I, n=12; stage II, n=4; stage III, n=7). The control group consisted of 36 individuals without malignancy (25 males, 11 females; mean age 57+/-13 years). Two observers independently evaluated the following WB-MRI sequences: T1w-TSE (T1), T2w-TIRM (T2), and the combination of both sequences, including a contrast-enhanced T1w-TSE with fat-saturation (T1+/-CE/T2). They had to determine growth patterns (focal and/or diffuse) and the MRI sequence that provided the highest confidence level in depicting the MM lesions. Results were calculated on a per-patient basis. RESULTS: Visual detection of MM was as follows: T1, 65% (sensitivity)/85% (specificity); T2, 76%/81%; T1+/-CE/T2, 67%/88%. Inter-observer variability was as follows: T1, 0.3; T2, 0.55; T1+/-CE/T2, 0.55. Sensitivity improved depending on infiltration grade (T1: 1=60%; 2=36%; 3=83%; T2: 1=70%; 2=71%; 3=89%; T1+/-CE/T2: 1=50%; 2=50%; 3=89%) and clinical stage (T1: 1=58%; 2=63%; 3=79%; T2: 1=58%; 2=88%; 3=100%; T1+/-CE/T2: 1=50%; 2=63%; 3=100%). T2w-TIRM sequences achieved the best reliability in depicting the MM lesions (65% in the mean of both readers). CONCLUSIONS: T2w-TIRM sequences achieved the highest level of sensitivity and best reliability, and thus might be valuable for initial assessment of MM. For an exact staging and grading the examination protocol should encompass unenhanced and enhanced T1w-MRI sequences, in addition to T2w-TIRM.

Functional properties of multiple isoforms of human divalent metal-ion transporter 1 (DMT1)

DMT1 (divalent metal-ion transporter 1) is a widely expressed metal-ion transporter that is vital for intestinal iron absorption and iron utilization by most cell types throughout the body, including erythroid precursors. Mutations in DMT1 cause severe microcytic anaemia in animal models. Four DMT1 isoforms that differ in their N- and C-termini arise from mRNA transcripts that vary both at their 5'-ends (starting in exon 1A or exon 1B) and at their 3'-ends giving rise to mRNAs containing (+) or lacking (-) the 3'-IRE (iron-responsive element) and resulting in altered C-terminal coding sequences. To determine whether these variations result in functional differences between isoforms, we explored the functional properties of each isoform using the voltage clamp and radiotracer assays in cRNA-injected Xenopus oocytes. 1A/IRE+-DMT1 mediated Fe2+-evoked currents that were saturable (K(0.5)(Fe) approximately 1-2 microM), temperature-dependent (Q10 approximately 2), H+-dependent (K(0.5)(H) approximately 1 muM) and voltage-dependent. 1A/IRE+-DMT1 exhibited the provisional substrate profile (ranked on currents) Cd2+, Co2+, Fe2+, Mn2+>Ni2+, V3+>>Pb2+. Zn2+ also evoked large currents; however, the zinc-evoked current was accounted for by H+ and Cl- conductances and was not associated with significant Zn2+ transport. 1B/IRE+-DMT1 exhibited the same substrate profile, Fe2+ affinity and dependence on the H+ electrochemical gradient. Each isoform mediated 55Fe2+ uptake and Fe2+-evoked currents at low extracellular pH. Whereas iron transport activity varied markedly between the four isoforms, the activity for each correlated with the density of anti-DMT1 immunostaining in the plasma membrane, and the turnover rate of the Fe2+ transport cycle did not differ between isoforms. Therefore all four isoforms of human DMT1 function as metal-ion transporters of equivalent efficiency. Our results reveal that the N- and C-terminal sequence variations among the DMT1 isoforms do not alter DMT1 functional properties. We therefore propose that these variations serve as tissue-specific signals or cues to direct DMT1 to the appropriate subcellular compartments (e.g. in erythroid cells) or the plasma membrane (e.g. in intestine).

Clustering subspecies of Aeromonas salmonicida using IS630 typing.

BACKGROUND The insertion element IS630 found in Aeromonas salmonicida belongs to the IS630-Tc1-mariner superfamily of transposons. It is present in multiple copies and represents approximately half of the IS present in the genome of A. salmonicida subsp. salmonicida A449. RESULTS By using High Copy Number IS630 Restriction Fragment Length Polymorphism (HCN-IS630-RFLP), strains of various subspecies of Aeromonas salmonicida showed conserved or clustering patterns, thus allowing their differentiation from each other. Fingerprints of A. salmonicida subsp. salmonicida showed the highest homogeneity while 'atypical' A. salmonicida strains were more heterogeneous. IS630 typing also differentiated A. salmonicida from other Aeromonas species. The copy number of IS630 in Aeromonas salmonicida ranges from 8 to 35 and is much lower in other Aeromonas species. CONCLUSIONS HCN-IS630-RFLP is a powerful tool for subtyping of A. salmonicida. The high stability of IS630 insertions in A. salmonicida subsp. salmonicida indicates that it might have played a role in pathoadaptation of A. salmonicida which has reached an optimal configuration in the highly virulent and specific fish pathogen A. salmonicida subsp. salmonicida.

A non-coding genomic duplication at the HMX1 locus is associated with crop ears in highland cattle.

Highland cattle with congenital crop ears have notches of variable size on the tips of both ears. In some cases, cartilage deformation can be seen and occasionally the external ears are shortened. We collected 40 cases and 80 controls across Switzerland. Pedigree data analysis confirmed a monogenic autosomal dominant mode of inheritance with variable expressivity. All affected animals could be traced back to a single common ancestor. A genome-wide association study was performed and the causative mutation was mapped to a 4 Mb interval on bovine chromosome 6. The H6 family homeobox 1 (HMX1) gene was selected as a positional and functional candidate gene. By whole genome re-sequencing of an affected Highland cattle, we detected 6 non-synonymous coding sequence variants and two variants in an ultra-conserved element at the HMX1 locus with respect to the reference genome. Of these 8 variants, only a non-coding 76 bp genomic duplication (g.106720058_106720133dup) located in the conserved region was perfectly associated with crop ears. The identified copy number variation probably results in HMX1 misregulation and possible gain-of-function. Our findings confirm the role of HMX1 during the development of the external ear. As it is sometimes difficult to phenotypically diagnose Highland cattle with slight ear notches, genetic testing can now be used to improve selection against this undesired trait.

DQB1 locus alone explains most of the risk and protection in narcolepsy with cataplexy in Europe

STUDY OBJECTIVE Prior research has identified five common genetic variants associated with narcolepsy with cataplexy in Caucasian patients. To replicate and/or extend these findings, we have tested HLA-DQB1, the previously identified 5 variants, and 10 other potential variants in a large European sample of narcolepsy with cataplexy subjects. DESIGN Retrospective case-control study. SETTING A recent study showed that over 76% of significant genome-wide association variants lie within DNase I hypersensitive sites (DHSs). From our previous GWAS, we identified 30 single nucleotide polymorphisms (SNPs) with P < 10(-4) mapping to DHSs. Ten SNPs tagging these sites, HLADQB1, and all previously reported SNPs significantly associated with narcolepsy were tested for replication. PATIENTS AND PARTICIPANTS For GWAS, 1,261 narcolepsy patients and 1,422 HLA-DQB1*06:02-matched controls were included. For HLA study, 1,218 patients and 3,541 controls were included. MEASUREMENTS AND RESULTS None of the top variants within DHSs were replicated. Out of the five previously reported SNPs, only rs2858884 within the HLA region (P < 2x10(-9)) and rs1154155 within the TRA locus (P < 2x10(-8)) replicated. DQB1 typing confirmed that DQB1*06:02 confers an extraordinary risk (odds ratio 251). Four protective alleles (DQB1*06:03, odds ratio 0.17, DQB1*05:01, odds ratio 0.56, DQB1*06:09 odds ratio 0.21, DQB1*02 odds ratio 0.76) were also identified. CONCLUSION An overwhelming portion of genetic risk for narcolepsy with cataplexy is found at DQB1 locus. Since DQB1*06:02 positive subjects are at 251-fold increase in risk for narcolepsy, and all recent cases of narcolepsy after H1N1 vaccination are positive for this allele, DQB1 genotyping may be relevant to public health policy.

Unravelling the multiple roles of FUS in RNA processing and on ALS pathogenesis

Two genes with related functions in RNA biogenesis were recently reported in patients with familial ALS: the FUS/TLS gene at the ALS6 locus and the TARDBP/TDP-43 gene at the ALS10 locus [1, 2]. FUS has been implicated to function in several steps of gene expression, including transcription regulation [3], RNA splicing [4, 5], mRNA transport in neurons [6] and, interestingly, in microRNA (miRNA) processing [7]. The goal of this project is to identify the molecular mechanisms leading to the development of FUS mutations-associated ALS. Specifically, we want to test the hypothesis that these FUS mutations misregulate miRNA levels that in turn affect the expression of genes critical for motor neuron survival. In addition we want to test whether misregulation of the miRNA profile is a common feature in ALS. We have performed immunoprecipitations from total extracts of 293T cells expressing FLAG-tagged FUS to characterize its interactome by mass spectrometry. This proteomic study not only revealed a strong interaction of FUS with splicing factors, but shows that FUS might be involved in many, quite different pathways. To map which parts of the FUS protein contribute to the interaction with splicing factors, we have performed a set of experiments with a series of missense and deletion mutants. With this approach, we will not only gain information on the binding partners of FUS along with a map of the required domains for the interactions, but it will also help to unravel whether certain ALS-associated FUS mutations lead to a loss or gain of function due to gain or loss of interactors. Additionally, we have performed quantitative interactomics using SILAC to identify interactome differences of ALS-associated FUS mutants. To this end we have performed immunoprecipitations of total extract from 293T cells, stably transduced with constructs expressing wild-type FUS-FLAG as well as three different ALS-associated mutants (G156E, R244C, P525L). First results indicate striking differences in the interactome with certain RNA binding proteins. We are now validating these candidates in order to reveal the importance of these differential interactions in the context of ALS.

Global phylogenomic analysis of nonencapsulated Streptococcus pneumoniae reveals a deep-branching classic lineage that is distinct from multiple sporadic lineages.

The surrounding capsule of Streptococcus pneumoniae has been identified as a major virulence factor and is targeted by pneumococcal conjugate vaccines (PCV). However, nonencapsulated Streptococcus pneumoniae (Non-Ec-Sp) have also been isolated globally, mainly in carriage studies. It is unknown if Non-Ec-Sp evolve sporadically, if they have high antibiotic non-susceptiblity rates and a unique, specific gene content. Here, whole genome sequencing of 131 Non-Ec-Sp isolates sourced from 17 different locations around the world was performed. Results revealed a deep-branching classic lineage that is distinct from multiple sporadic lineages. The sporadic lineages clustered with a previously sequenced, global collection of encapsulated S. pneumoniae (Ec-Sp) isolates while the classic lineage is comprised mainly of the frequently identified multi-locus sequences types ST344 (n=39) and ST448 (n=40). All ST344 and nine ST448 isolates had high non-susceptiblity rates to β-lactams and other antimicrobials. Analysis of the accessory genome reveals that the classic Non-Ec-Sp contained an increased number of mobile elements, than Ec-Sp and sporadic Non-Ec-Sp. Performing adherence assays to human epithelial cells for selected classic and sporadic Non-Ec-Sp revealed that the presence of a integrative conjugative element (ICE) results in increased adherence to human epithelial cells (P=0.005). In contrast, sporadic Non-Ec-Sp lacking the ICE had greater growth in vitro possibly resulting in improved fitness. In conclusion, Non-Ec-Sp isolates from the classic lineage have evolved separately. They have spread globally, are well adapted to nasopharyngeal carriage and are able to coexist with Ec-Sp. Due to continued use of pneumococcal conjugate vaccines, Non-Ec-Sp may become more prevalent.

An mRNA-derived ncRNA targets and regulates the ribosome

Small non-protein-coding RNA (ncRNA) molecules are key players in controlling gene expression at multiple steps in all domains of life. While the list of validated ncRNAs that regulate crucial cellular processes grows steadily (such as micro RNAs and small-interfering RNAs), not a single ncRNA has been identified that directly interacts and regulates the ribosome during protein biosynthesis (with the notable exceptions of 7SL RNA and tmRNA). This is unexpected, given the central position the ribosome plays during gene expression. To investigate whether such a class of regulatory ncRNAs does exist we performed genomic screens for small ribosome-associated RNAs in various model organisms of all three domains [1,2]. Here we show that an mRNA-derived 18 nucleotide long ncRNA is capable of down-regulating translation in Saccharomyces cerevisiae by directly targeting the ribosome [3]. This 18-mer ncRNA binds to polysomes upon salt stress and is crucial for efficient growth under hyperosmotic conditions. Although the 18-mer RNA originates from the TRM10 locus, which encodes a tRNA methyltransferase, genetic analyses revealed the 18-mer RNA nucleotide sequence, rather than the mRNA-encoded enzyme, as the translation regulator under these stress conditions. Our data reveal the ribosome as a target for small regulatory ncRNAs and unveil the existence of a novel mechanism of translation regulation. Analogous genomic screens in organisms spanning all three domains of life demonstrate the existence of thousands of ncRNA candidates putatively regulating the ribosome. We therefore anticipate that ribosome-bound ncRNAs are capable of fine tuning translation and might represent a so far largely unexplored class of regulatory ncRNAs.