Rapid and reliable genotyping of polymorphic loci modifying correct splicing of CFTR pre-mRNA using mass spectrometry

We describe a fast and unambiguous method for haplotyping the (TG)mTn repeat in IVS8 and determining three other single nucleotide polymorphisms (SNPs) in exons 10, 14a and 24 in the cystic fibrosis transmembrane conductance regulator (CFTR) gene affecting correct splicing of the CFTR pre-mRNA using primer extension and mass spectrometry. The diagnostic products are generated by primer extension (PEX) reactions, which require a single detection primer complementary to a region downstream of a target strand's variable site. On addition of a polymerase and an appropriate mixture of dNTP's and 2', 3'-dideoxynucleotide triphosphates (ddNTP's), the primer is extended through the mutation region until the first ddNTP is incorporated and the mass of the extension products determines the composition of the variable site. Analysis of patient DNA assigned the correct and unambiguous haplotype for the (TG)mTn repeat in intron 8 of the CFTR gene. Additional crucial SNPs influencing correct splicing in exon 10, 14 and 24 can easily be detected by biplexing the assay to genotype allelic variants important for correct splicing of the CFTR pre-mRNA. Different PEX reactions with subsequent mass spectrometry generate sufficient data, to enable unambiguous and easy haplotyping of the (TG)mTn repeat in the CFTR gene. The method can be easily extended to the inclusion of additional SNPs of interest by biplexing some of the PEX reactions. All experimental steps required for PEX are amenable to the high degree of automation desirable for a high-throughput diagnostic setting, facilitating the work of clinicians involved in the diagnosis of non-classic cystic fibrosis.

ASSOCIATON TESTS THAT ACCOMMODATE GENOTYPING ERRORS

High-throughput SNP arrays provide estimates of genotypes for up to one million loci, often used in genome-wide association studies. While these estimates are typically very accurate, genotyping errors do occur, which can influence in particular the most extreme test statistics and p-values. Estimates for the genotype uncertainties are also available, although typically ignored. In this manuscript, we develop a framework to incorporate these genotype uncertainties in case-control studies for any genetic model. We verify that using the assumption of a “local alternative” in the score test is very reasonable for effect sizes typically seen in SNP association studies, and show that the power of the score test is simply a function of the correlation of the genotype probabilities with the true genotypes. We demonstrate that the power to detect a true association can be substantially increased for difficult to call genotypes, resulting in improved inference in association studies.

Development of a genotyping microarray for Usher syndrome

BACKGROUND: Usher syndrome, a combination of retinitis pigmentosa (RP) and sensorineural hearing loss with or without vestibular dysfunction, displays a high degree of clinical and genetic heterogeneity. Three clinical subtypes can be distinguished, based on the age of onset and severity of the hearing impairment, and the presence or absence of vestibular abnormalities. Thus far, eight genes have been implicated in the syndrome, together comprising 347 protein-coding exons. METHODS: To improve DNA diagnostics for patients with Usher syndrome, we developed a genotyping microarray based on the arrayed primer extension (APEX) method. Allele-specific oligonucleotides corresponding to all 298 Usher syndrome-associated sequence variants known to date, 76 of which are novel, were arrayed. RESULTS: Approximately half of these variants were validated using original patient DNAs, which yielded an accuracy of >98%. The efficiency of the Usher genotyping microarray was tested using DNAs from 370 unrelated European and American patients with Usher syndrome. Sequence variants were identified in 64/140 (46%) patients with Usher syndrome type I, 45/189 (24%) patients with Usher syndrome type II, 6/21 (29%) patients with Usher syndrome type III and 6/20 (30%) patients with atypical Usher syndrome. The chip also identified two novel sequence variants, c.400C>T (p.R134X) in PCDH15 and c.1606T>C (p.C536S) in USH2A. CONCLUSION: The Usher genotyping microarray is a versatile and affordable screening tool for Usher syndrome. Its efficiency will improve with the addition of novel sequence variants with minimal extra costs, making it a very useful first-pass screening tool.

USING THE R PACKAGE crlmm FOR GENOTYPING AND COPY NUMBER ESTIMATION

Genotyping platforms such as Affymetrix can be used to assess genotype-phenotype as well as copy number-phenotype associations at millions of markers. While genotyping algorithms are largely concordant when assessed on HapMap samples, tools to assess copy number changes are more variable and often discordant. One explanation for the discordance is that copy number estimates are susceptible to systematic differences between groups of samples that were processed at different times or by different labs. Analysis algorithms that do not adjust for batch effects are prone to spurious measures of association. The R package crlmm implements a multilevel model that adjusts for batch effects and provides allele-specific estimates of copy number. This paper illustrates a workflow for the estimation of allele-specific copy number, develops markerand study-level summaries of batch effects, and demonstrates how the marker-level estimates can be integrated with complimentary Bioconductor software for inferring regions of copy number gain or loss. All analyses are performed in the statistical environment R. A compendium for reproducing the analysis is available from the author’s website (http://www.biostat.jhsph.edu/~rscharpf/crlmmCompendium/index.html).

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.

Genotyping of Francisella tularensis strains by pulsed-field gel electrophoresis, amplified fragment length polymorphism fingerprinting, and 16S rRNA gene sequencing

We evaluated three molecular methods for identification of Francisella strains: pulsed-field gel electrophoresis (PFGE), amplified fragment length polymorphism (AFLP) analysis, and 16S rRNA gene sequencing. The analysis was performed with 54 Francisella tularensis subsp. holarctica, 5 F. tularensis subsp. tularensis, 2 F. tularensis subsp. novicida, and 1 F. philomiragia strains. On the basis of the combination of results obtained by PFGE with the restriction enzymes XhoI and BamHI, PFGE revealed seven pulsotypes, which allowed us to discriminate the strains to the subspecies level and which even allowed us to discriminate among some isolates of F. tularensis subsp. holarctica. The AFLP analysis technique produced some degree of discrimination among F. tularensis subsp. holarctica strains (one primary cluster with three major subclusters and minor variations within subclusters) when EcoRI-C and MseI-A, EcoRI-T and MseI-T, EcoRI-A and MseI-C, and EcoRI-0 and MseI-CA were used as primers. The degree of similarity among the strains was about 94%. The percent similarities of the AFLP profiles of this subspecies compared to those of F. tularensis subsp. tularensis, F. tularensis subsp. novicida, and F. philomiragia were less than 90%, about 72%, and less than 24%, respectively, thus permitting easy differentiation of this subspecies. 16S rRNA gene sequencing revealed 100% similarity for all F. tularensis subsp. holarctica isolates compared in this study. These results suggest that although limited genetic heterogeneity among F. tularensis subsp. holarctica isolates was observed, PFGE and AFLP analysis appear to be promising tools for the diagnosis of infections caused by different subspecies of F. tularensis and suitable techniques for the differentiation of individual strains.

Occurrence, quantification, and genotyping of Mycoplasma conjunctivae in wild Caprinae with and without infectious keratoconjunctivitis

Mycoplasma conjunctivae, the causative agent of infectious keratoconjunctivitis (IKC), was recently detected in asymptomatic Alpine ibex (Capra ibex ibex). This suggested that an external source of infection may not be required for an IKC outbreak in wildlife but might be initiated by healthy carriers, which contradicted previous serologic investigations in chamois. Our aims were to 1) assess the prevalence of M. conjunctivae among asymptomatic ibex and Alpine chamois (Rupicapra rupicapra rupicapra) and its frequency in IKC-affected animals, 2) determine mycoplasma loads in different disease stages, and 3) characterize the M. conjunctivae strains involved. Eye swabs from 654 asymptomatic and 204 symptomatic animals were collected in diverse Swiss regions between 2008 and 2010, and tested by TaqMan real-time PCR. Data analysis was performed considering various patterns of IKC occurrence in the respective sampling regions. Strains from 24 animals were compared by cluster analysis. Prevalence of M. conjunctivae was 5.6% (95% confidence interval [CI]: 3.7-8.1%) in asymptomatic ibex and 5.8% (CI: 3.0-9.9%) in asymptomatic chamois, with significant differences between years and regions in both species. Detection frequency in symptomatic animals was significantly higher during IKC outbreaks than in nonepidemic situations (i.e., regular but low incidence or sporadic occurrence). Mycoplasma load was significantly lower in eyes from healthy carriers and animals with mild signs than from animals with moderate and severe signs. Although some strains were found in both asymptomatic and diseased animals of the same species, others apparently differed in their pathogenic potential depending on the infected species. Overall, we found a widespread occurrence of M. conjunctivae in wild Caprinae with and without IKC signs. Our results confirm the central role of M. conjunctivae in outbreaks but suggest that other infectious agents may be involved in IKC cases in nonepidemic situations. Additionally, presence and severity of signs are related to the quantity of M. conjunctivae in the eyes rather than to the strain. We propose that individual or environmental factors influence the clinical expression of the disease and that persistence of M. conjunctivae in populations of wild Caprinae cannot be excluded.

A multiplex PCR for the simultaneous detection and genotyping of the Echinococcus granulosus complex.

Echinococcus granulosus is characterized by high intra-specific variability (genotypes G1-G10) and according to the new molecular phylogeny of the genus Echinococcus, the E. granulosus complex has been divided into E. granulosus sensu stricto (G1-G3), E. equinus (G4), E. ortleppi (G5), and E. canadensis (G6-G10). The molecular characterization of E. granulosus isolates is fundamental to understand the spatio-temporal epidemiology of this complex in many endemic areas with the simultaneous occurrence of different Echinococcus species and genotypes. To simplify the genotyping of the E. granulosus complex we developed a single-tube multiplex PCR (mPCR) allowing three levels of discrimination: (i) Echinococcus genus, (ii) E. granulosus complex in common, and (iii) the specific genotype within the E. granulosus complex. The methodology was established with known DNA samples of the different strains/genotypes, confirmed on 42 already genotyped samples (Spain: 22 and Bulgaria: 20) and then successfully applied on 153 unknown samples (Tunisia: 114, Algeria: 26 and Argentina: 13). The sensitivity threshold of the mPCR was found to be 5 ng Echinoccoccus DNA in a mixture of up to 1 µg of foreign DNA and the specificity was 100% when template DNA from closely related members of the genus Taenia was used. Additionally to DNA samples, the mPCR can be carried out directly on boiled hydatid fluid or on alkaline-lysed frozen or fixed protoscoleces, thus avoiding classical DNA extractions. However, when using Echinococcus eggs obtained from fecal samples of infected dogs, the sensitivity of the mPCR was low (<40%). Thus, except for copro analysis, the mPCR described here has a high potential for a worldwide application in large-scale molecular epidemiological studies on the Echinococcus genus.

Clinical outcomes and molecular genotyping of Staphylococcus aureus isolated from milk samples of dairy primiparous Mediterranean buffaloes (Bubalus bubalis)

Staphylococcus aureus is one of the most important pathogens causing mastitis in dairy cows and in Mediterranean buffaloes. Genotype B (GTB) is contagious in dairy cows and may occur in up to 87% of cows of a dairy herd. It was the aim of this study to evaluate genotypes present, clinical outcomes, and prevalence of Staph. aureus in milk samples of primiparous Mediterranean dairy buffaloes. Two hundred composite milk samples originating from 40 primiparous buffaloes were collected from May to June 2012, at d 10, 30, 60, 90, and 150 d in milk (DIM) to perform somatic cell counts and bacteriological cultures. Daily milk yields were recorded. Before parturition until 40 to 50 DIM, all primiparous animals were housed separated from the pluriparous animals. Milking was performed in the same milking parlor, but the primiparous animals were milked first. After 50 DIM, the primiparous were mixed with the pluriparous animals, including the milking procedure. Individual quarter samples were collected from each animal, and aliquots of 1 mL were mixed and used for molecular identification and genotyping of Staph. aureus. The identification of Staph. aureus was performed verifying the presence of nuc gene by nuc gene PCR. All the nuc-positive isolates were subjected to genotype analysis by means of PCR amplification of the 16S-23S rRNA intergenic spacer region and analyzed by a miniaturized electrophoresis system. Of all 200 composite samples, 41 (20.5%) were positive for Staph. aureus, and no genotype other than GTB was identified. The prevalence of samples positive for Staph. aureus was 0% at 10 DIM and increased to a maximum of 22/40 (55%) at 90 DIM. During the period of interest, 14 buffaloes tested positive for Staph. aureus once, 6 were positive twice, and 5 were positive 3 times, whereas 15 animals were negative at every sampling. At 90 and 150 DIM, 7 (17.5%) and 3 buffaloes (7.5%), respectively, showed clinical mastitis (CM), and only 1 (2.5%) showed CM at both samplings. At 60, 90, and 150 DIM, 1 buffalo was found with subclinical mastitis at each sampling. At 30, 60, 90, and 150 DIM, 2.5 (1/40), 22.5 (9/40), 35 (14/40), and 10% (4/40) were considered affected by intramammary infection, respectively. Buffaloes with CM caused by Staph. aureus had statistically significantly higher mean somatic cell count values (6.06 ± 0.29, Log10 cells/mL ± standard deviation) and statistically significantly lower mean daily milk yields (7.15 ± 1.49, liters/animal per day) than healthy animals (4.69 ± 0.23 and 13.87 ± 2.64, respectively), buffaloes with IMI (4.82 ± 0.23 and 11.16 ± 1.80, respectively), or with subclinical mastitis (5.47 ± 0.10 and 10.33 ± 0.68, respectively). Based on our knowledge, this is the first time that Staph. aureus GTB has been identified in milk samples of dairy Mediterranean buffaloes.

Standard genotyping overestimates transmission of Mycobacterium tuberculosis among immigrants in a low incidence country.

Immigrants from high tuberculosis (TB) incidence regions are a risk group for TB in low-incidence countries such as Switzerland. In a previous analysis of a nationwide collection of 520 Mycobacterium tuberculosis isolates from 2000-2008, we identified 35 clusters comprising 90 patients based on standard genotyping (24-loci MIRU-VNTR and spoligotyping). Here, we used whole genome sequencing (WGS) to revisit these transmission clusters. Genome-based transmission clusters were defined as isolate pairs separated by ≤12 single nucleotide polymorphisms (SNPs). WGS confirmed 17/35 (49%) MIRU-VNTR clusters; the other 18 clusters contained pairs separated by >12 SNPs. Most transmission clusters (3/4) of Swiss-born patients were confirmed by WGS, as opposed to 25% (4/16) of clusters involving only foreign-born patients. The overall clustering proportion using standard genotyping was 17% (90 patients, 95% confidence interval [CI]: 14-21%), but only 8% (43 patients, 95% CI: 6-11%) using WGS. The clustering proportion was 17% (67/401, 95% CI: 13-21%) using standard genotyping and 7% (26/401, 95% CI: 4-9%) using WGS among foreign-born patients, and 19% (23/119, 95% CI: 13-28%) and 14% (17/119, 95% CI: 9-22%), respectively, among Swiss-born patients. Using weighted logistic regression, we found weak evidence for an association between birth origin and transmission (aOR 2.2, 95% CI: 0.9-5.5, comparing Swiss-born patients to others). In conclusion, standard genotyping overestimated recent TB transmission in Switzerland when compared to WGS, particularly among immigrants from high TB incidence regions, where genetically closely related strains often predominate. We recommend the use of WGS to identify transmission clusters in low TB incidence settings.

Whole-genome genotyping of grape using a panel of microsatellite

The use of microsatellite markers in large-scale genetic studies is limited by its low throughput and high cost and labor requirements. Here, we provide a panel of 45 multiplex PCRs for fast and cost-efficient genome-wide fluorescence-based microsatellite analysis in grapevine. The developed multiplex PCRs panel (with up to 15-plex) enables the scoring of 270 loci covering all the grapevine genome (9 to 20 loci/chromosome) using only 45 PCRs and sequencer runs. The 45 multiplex PCRs were validated using a diverse grapevine collection of 207 accessions, selected to represent most of the cultivated Vitis vinifera genetic diversity. Particular attention was paid to quality control throughout the whole process (assay replication, null allele detection, ease of scoring). Genetic diversity summary statistics and features of electrophoretic profiles for each studied marker are provided, as are the genotypes of 25 common cultivars that could be used as references in other studies.

Diversity Arrays: a solid state technology for sequence information independent genotyping

Here we present the successful application of the microarray technology platform to the analysis of DNA polymorphisms. Using the rice genome as a model, we demonstrate the potential of a high-throughput genome analysis method called Diversity Array Technology, DArT‘. In the format presented here the technology is assaying for the presence (or amount) of a specific DNA fragment in a representation derived from the total genomic DNA of an organism or a population of organisms. Two different approaches are presented: the first involves contrasting two representations on a single array while the second involves contrasting a representation with a reference DNA fragment common to all elements of the array. The Diversity Panels created using this method allow genetic fingerprinting of any organism or group of organisms belonging to the gene pool from which the panel was developed. Diversity Arrays enable rapid and economical application of a highly parallel, solid-state genotyping technology to any genome or complex genomic mixtures.

PCR candidate region mismatch scanning: adaptation to quantitative, high-throughput genotyping

Linkage and association analyses were performed to identify loci affecting disease susceptibility by scoring previously characterized sequence variations such as microsatellites and single nucleotide polymorphisms. Lack of markers in regions of interest, as well as difficulty in adapting various methods to high-throughput settings, often limits the effectiveness of the analyses. We have adapted the Escherichia coli mismatch detection system, employing the factors MutS, MutL and MutH, for use in PCR-based, automated, high-throughput genotyping and mutation detection of genomic DNA. Optimal sensitivity and signal-to-noise ratios were obtained in a straightforward fashion because the detection reaction proved to be principally dependent upon monovalent cation concentration and MutL concentration. Quantitative relationships of the optimal values of these parameters with length of the DNA test fragment were demonstrated, in support of the translocation model for the mechanism of action of these enzymes, rather than the molecular switch model. Thus, rapid, sequence-independent optimization was possible for each new genomic target region. Other factors potentially limiting the flexibility of mismatch scanning, such as positioning of dam recognition sites within the target fragment, have also been investigated. We developed several strategies, which can be easily adapted to automation, for limiting the analysis to intersample heteroduplexes. Thus, the principal barriers to the use of this methodology, which we have designated PCR candidate region mismatch scanning, in cost-effective, high-throughput settings have been removed.

DNA sequencing and genotyping by transcriptional synthesis of chain-terminated RNA ladders and MALDI-TOF mass spectrometry

Sets of RNA ladders can be synthesized by transcription of a bacteriophage-encoded RNA polymerase using 3′-deoxynucleotides as chain terminators. These ladders can be used for sequencing of DNA. Using a nicked form of phage SP6 RNA polymerase in this study substantially enhanced yields of transcriptional sequencing ladders. Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) of chain-terminated RNA ladders allowed DNA sequence determination of up to 56 nt. It is also demonstrated that A→G and C→T variations in heterozygous and homozygous samples can be unambiguously identified by the mass spectrometric analysis. As a step towards single-tube sequencing reactions, α-thiotriphosphate nucleotide analogs were used to overcome problems caused by chain terminator-independent, premature termination and by the small mass difference between natural pyrimidine nucleotides.