Multiplex strategy for multilocus sequence typing, fla typing, and genetic determination of antimicrobial resistance of Campylobacter jejuni and Campylobacter coli isolates collected in Switzerland

We present an optimized multilocus sequence typing (MLST) scheme with universal primer sets for amplifying and sequencing the seven target genes of Campylobacter jejuni and Campylobacter coli. Typing was expanded by sequence determination of the genes flaA and flaB using optimized primer sets. This approach is compatible with the MLST and flaA schemes used in the PubMLST database and results in an additional typing method using the flaB gene sequence. An identification module based on the 16S rRNA and rpoB genes was included, as well as the genetic determination of macrolide and quinolone resistances based on mutations in the 23S rRNA and gyrA genes. Experimental procedures were simplified by multiplex PCR of the 13 target genes. This comprehensive approach was evaluated with C. jejuni and C. coli isolates collected in Switzerland. MLST of 329 strains resulted in 72 sequence types (STs) among the 186 C. jejuni strains and 39 STs for the 143 C. coli isolates. Fourteen (19%) of the C. jejuni and 20 (51%) of the C. coli STs had not been found previously. In total, 35% of the C. coli strains collected in Switzerland contained mutations conferring antibiotic resistance only to quinolone, 15% contained mutations conferring resistance only to macrolides, and 6% contained mutations conferring resistance to both classes of antibiotics. In C. jejuni, these values were 31% and 0% for quinolone and macrolide resistance, respectively. The rpoB sequence allowed phylogenetic differentiation between C. coli and C. jejuni, which was not possible by 16S rRNA gene analysis. An online Integrated Database Network System (SmartGene, Zug, Switzerland)-based platform for MLST data analysis specific to Campylobacter was implemented. This Web-based platform allowed automated allele and ST designation, as well as epidemiological analysis of data, thus streamlining and facilitating the analysis workflow. Data networking facilitates the exchange of information between collaborating centers. The described approach simplifies and improves the genotyping of Campylobacter, allowing cost- and time-efficient routine monitoring.

Genetic relationships of Aeromonas strains inferred from 16S rRNA, gyrB and rpoB gene sequences

Genetic relationships among bacterial strains belonging to the genus Aeromonas were inferred from 16S rRNA, gyrB and rpoB gene sequences. Twenty-eight type or collection strains of the recognized species or subspecies and 33 Aeromonas strains isolated from human and animal specimens as well as from environmental samples were included in the study. As reported previously, the 16S rRNA gene sequence is highly conserved within the genus Aeromonas, having only limited resolution for this very tight group of species. Analysis of a 1.1 kb gyrB sequence confirmed that this gene has high resolving power, with maximal interspecies divergence of 15.2 %. Similar results were obtained by sequencing only 517 bp of the rpoB gene, which showed maximal interspecies divergence of 13 %. The topologies of the gyrB- and rpoB-derived trees were similar. The results confirm the close relationship of species within the genus Aeromonas and show that a phylogenetic approach including several genes is suitable for improving the complicated taxonomy of the genus.

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.

Spectrum and prevalence of mutations involving BrS1- through BrS12-susceptibility genes in a cohort of unrelated patients referred for Brugada syndrome genetic testing: implications for genetic testing

OBJECTIVES The aim of this study was to provide the spectrum and prevalence of mutations in the 12 Brugada syndrome (BrS)-susceptibility genes discovered to date in a single large cohort of unrelated BrS patients. BACKGROUND BrS is a potentially lethal heritable arrhythmia syndrome diagnosed electrocardiographically by coved-type ST-segment elevation in the right precordial leads (V1 to V3; type 1 Brugada electrocardiographic [ECG] pattern) and the presence of a personal/family history of cardiac events. METHODS Using polymerase chain reaction, denaturing high-performance liquid chromatography, and DNA sequencing, comprehensive mutational analysis of BrS1- through BrS12-susceptibility genes was performed in 129 unrelated patients with possible/probable BrS (46 with clinically diagnosed BrS [ECG pattern plus personal/family history of a cardiac event] and 83 with a type 1 BrS ECG pattern only). RESULTS Overall, 27 patients (21%) had a putative pathogenic mutation, absent in 1,400 Caucasian reference alleles, including 21 patients with an SCN5A mutation, 2 with a CACNB2B mutation, and 1 each with a KCNJ8 mutation, a KCND3 mutation, an SCN1Bb mutation, and an HCN4 mutation. The overall mutation yield was 23% in the type 1 BrS ECG pattern-only patients versus 17% in the clinically diagnosed BrS patients and was significantly greater among young men<20 years of age with clinically diagnosed BrS and among patients who had a prolonged PQ interval. CONCLUSIONS We identified putative pathogenic mutations in ∼20% of our BrS cohort, with BrS genes 2 through 12 accounting for <5%. Importantly, the yield was similar between patients with only a type 1 BrS ECG pattern and those with clinically established BrS. The yield approaches 40% for SCN5A-mediated BrS (BrS1) when the PQ interval exceeds 200 ms. Calcium channel-mediated BrS is extremely unlikely in the absence of a short QT interval.

Cardiac channel molecular autopsy: insights from 173 consecutive cases of autopsy-negative sudden unexplained death referred for postmortem genetic testing

OBJECTIVE To perform long QT syndrome and catecholaminergic polymorphic ventricular tachycardia cardiac channel postmortem genetic testing (molecular autopsy) for a large cohort of cases of autopsy-negative sudden unexplained death (SUD). METHODS From September 1, 1998, through October 31, 2010, 173 cases of SUD (106 males; mean ± SD age, 18.4 ± 12.9 years; age range, 1-69 years; 89% white) were referred by medical examiners or coroners for a cardiac channel molecular autopsy. Using polymerase chain reaction, denaturing high-performance liquid chromatography, and DNA sequencing, a comprehensive mutational analysis of the long QT syndrome susceptibility genes (KCNQ1, KCNH2, SCN5A, KCNE1, and KCNE2) and a targeted analysis of the catecholaminergic polymorphic ventricular tachycardia type 1-associated gene (RYR2) were conducted. RESULTS Overall, 45 putative pathogenic mutations absent in 400 to 700 controls were identified in 45 autopsy-negative SUD cases (26.0%). Females had a higher yield (26/67 [38.8%]) than males (19/106 [17.9%]; P<.005). Among SUD cases with exercise-induced death, the yield trended higher among the 1- to 10-year-olds (8/12 [66.7%]) compared with the 11- to 20-year-olds (4/27 [14.8%]; P=.002). In contrast, for those who died during a period of sleep, the 11- to 20-year-olds had a higher yield (9/25 [36.0%]) than the 1- to 10-year-olds (1/24 [4.2%]; P=.01). CONCLUSION Cardiac channel molecular autopsy should be considered in the evaluation of autopsy-negative SUD. Several interesting genotype-phenotype observations may provide insight into the expected yields of postmortem genetic testing for SUD and assist in selecting cases with the greatest potential for mutation discovery and directing genetic testing efforts.

Tularemia in a common marmoset (Callithrix jacchus) diagnosed by 16S rRNA sequencing

We report a case of tularemia in a common marmoset (Callithrix jacchus) diagnosed by determination of the isolate's 16S ribosomal RNA (rRNA) gene sequence. Pathological examination of the animal revealed a multifocal acute necrotizing hepatitis, interstitial nephritis, splenitis, and lymphangitis of the mandibular, retropharyngeal, and cervical and mesenteric lymph nodes. Moreover, multiple foci of acute necrosis were found in the epithelium of the jejunum and the interstitium of the lung. Bacteriological investigations revealed a septicemia. The isolated infectious agent was uncommon, not routinely diagnosed in our laboratory and therefore difficult to identify by conventional tools in a reasonable time and effort. thus, we decided to perform a genetic analysis based on the 16S rRNA gene sequence. Thereby, an infection with Francisella tularensis, the causative agent of tularemia, was unambiguously diagnosed. This shows the great advantage 16S rRNA gene sequencing has as a general identification approach for unusual or rare isolates.

Maine Coon renal screening: ultrasonographical characterisation and preliminary genetic analysis for common genes in cats with renal cysts.

The objective of this study was to assess the prevalence of renal cysts and other renal abnormalities in purebred Maine Coon cats, and to characterise these through genetic typing. Voluntary pre-breeding screening programmes for polycystic kidney disease (PKD) are offered for this breed throughout Switzerland, Germany and other northern European countries. We performed a retrospective evaluation of Maine Coon screening for renal disease at one institution over an 8-year period. Renal ultrasonography was performed in 187 healthy Maine Coon cats. Renal changes were observed in 27 of these cats. Renal cysts were found in seven cats, and were mostly single and unilateral (6/7, 85.7%), small (mean 3.6 mm) and located at the corticomedullary junction (4/6, 66.7%). Sonographical changes indicating chronic kidney disease (CKD) were observed in 10/187 (5.3%) cats and changes of unknown significance were documented in 11/187 (5.9%) cats. All six cats genetically tested for PKD1 were negative for the mutation, and gene sequencing of these cats did not demonstrate any common genetic sequences. Cystic renal disease occurs with a low prevalence in Maine Coons and is unrelated to the PKD observed in Persians and related breeds. Ultrasonographical findings compatible with CKD are not uncommon in juvenile Maine Coons.

Sequencing of TGF-beta pathway genes in familial cases of intracranial aneurysm.

BACKGROUND AND PURPOSE: Familial aggregation of intracranial aneurysms (IA) strongly suggests a genetic contribution to pathogenesis. However, genetic risk factors have yet to be defined. For families affected by aortic aneurysms, specific gene variants have been identified, many affecting the receptors to transforming growth factor-beta (TGF-beta). In recent work, we found that aortic and intracranial aneurysms may share a common genetic basis in some families. We hypothesized, therefore, that mutations in TGF-beta receptors might also play a role in IA pathogenesis. METHODS: To identify genetic variants in TGF-beta and its receptors, TGFB1, TGFBR1, TGFBR2, ACVR1, TGFBR3, and ENG were directly sequenced in 44 unrelated patients with familial IA. Novel variants were confirmed by restriction digestion analyses, and allele frequencies were analyzed in cases versus individuals without known intracranial disease. Similarly, allele frequencies of a subset of known SNPs in each gene were also analyzed for association with IA. RESULTS: No mutations were found in TGFB1, TGFBR1, TGFBR2, or ACVR1. Novel variants identified in ENG (p.A60E) and TGFBR3 (p.W112R) were not detected in at least 892 reference chromosomes. ENG p.A60E showed significant association with familial IA in case-control studies (P=0.0080). No association with IA could be found for any of the known polymorphisms tested. CONCLUSIONS: Mutations in TGF-beta receptor genes are not a major cause of IA. However, we identified rare variants in ENG and TGFBR3 that may be important for IA pathogenesis in a subset of families.

Genetic and physical analysis of the bacteriophage P22 tail protein

A plasmid based genetic system was developed for the tail protein of the Salmonella typhimurium bacteriophage P22 and used to isolate and characterize tail protein mutants. The tail protein is a trimeric structural protein of the phage and an endorhamnosidase whose activity is essential for infection. The gene for the tail protein has previously been cloned into a plasmid expression vector and sequenced. A plate complementation assay for tail protein produced from the cloned gene was developed and used to isolate 27 tail protein mutants following mutagenesis of the cloned gene. These mutations were mapped into 12 deletion intervals using deletions which were made on plasmids in vitro and crossed onto P22. The base substitutions were determined by DNA sequencing. The majority of mutants had missense or nonsense mutations in the protein coding portion of the gene; however four of the mutants were in the putative transcription terminator. The oligomeric state of tail protein from the 15 missense mutants was investigated using SDS and nondenaturing polyacrylamide gel electrophoresis of cell lysates. Wild-type tail protein retains its trimeric structure in SDS gels at room temperature. Two of the mutant proteins also migrated as trimers in SDS gels, yet one of these had a considerably faster mobility than wild-type trimer. Its migration was the same as wild-type in a nondenaturing gel, so it is thought to be a trimer which is partially denatured by SDS. Four of the mutants produced proteins which migrate at the position of a monomer in an SDS gel but cannot be seen on a nondenaturing gel. These proteins are thought to be either monomers or soluble aggregates which cannot enter the nondenaturing gel. The remainder of mutants produce protein which is degraded. The mutant tail protein which had normal trimeric mobility on SDS and nondenaturing gels was purified. This protein has essentially wild-type ability to attach to phage capsids, but its endorhamnosidase activity is only 4% of wild-type. ^

Occurrence and Genetic Characteristics of Third-Generation Cephalosporin-Resistant Escherichia coli in Swiss Retail Meat

Prevalence and genetic relatedness were determined for third-generation cephalosporin-resistant Escherichia coli (3GC-R-Ec) detected in Swiss beef, veal, pork, and poultry retail meat. Samples from meat-packing plants (MPPs) processing 70% of the slaughtered animals in Switzerland were purchased at different intervals between April and June 2013 and analyzed. Sixty-nine 3GC-R-Ec isolates were obtained and characterized by microarray, PCR/DNA sequencing, Multi Locus Sequence Typing (MLST), and plasmid replicon typing. Plasmids of selected strains were transformed by electroporation into E. coli TOP10 cells and analyzed by plasmid MLST. The prevalence of 3GC-R-Ec was 73.3% in chicken and 2% in beef meat. No 3GC-R-Ec were found in pork and veal. Overall, the blaCTX-M-1 (79.4%), blaCMY-2 (17.6%), blaCMY-4 (1.5%), and blaSHV-12 (1.5%) β-lactamase genes were detected, as well as other genes conferring resistance to chloramphenicol (cmlA1-like), sulfonamides (sul), tetracycline (tet), and trimethoprim (dfrA). The 3GC-R-Ec from chicken meat often harbored virulence genes associated with avian pathogens. Plasmid incompatibility (Inc) groups IncI1, IncFIB, IncFII, and IncB/O were the most frequent. A high rate of clonality (e.g., ST1304, ST38, and ST93) among isolates from the same MPPs suggests that strains persist at the plant and spread to meat at the carcass-processing stage. Additionally, the presence of the blaCTX-M-1 gene on an IncI1 plasmid sequence type 3 (IncI1/pST3) in genetically diverse strains indicates interstrain spread of an epidemic plasmid. The blaCMY-2 and blaCMY-4 genes were located on IncB/O plasmids. This study represents the first comprehensive assessment of 3GC-R-Ec in meat in Switzerland. It demonstrates the need for monitoring contaminants and for the adaptation of the Hazard Analysis and Critical Control Point concept to avoid the spread of multidrug-resistant bacteria through the food chain.

Genome and transcriptome sequencing identifies breeding targets in the orphan crop tef (Eragrostis tef)

Background: Tef (Eragrostis tef), an indigenous cereal critical to food security in the Horn of Africa, is rich in minerals and protein, resistant to many biotic and abiotic stresses and safe for diabetics as well as sufferers of immune reactions to wheat gluten. We present the genome of tef, the first species in the grass subfamily Chloridoideae and the first allotetraploid assembled de novo. We sequenced the tef genome for marker-assisted breeding, to shed light on the molecular mechanisms conferring tef's desirable nutritional and agronomic properties, and to make its genome publicly available as a community resource. Results: The draft genome contains 672 Mbp representing 87% of the genome size estimated from flow cytometry. We also sequenced two transcriptomes, one from a normalized RNA library and another from unnormalized RNASeq data. The normalized RNA library revealed around 38000 transcripts that were then annotated by the SwissProt group. The CoGe comparative genomics platform was used to compare the tef genome to other genomes, notably sorghum. Scaffolds comprising approximately half of the genome size were ordered by syntenic alignment to sorghum producing tef pseudo-chromosomes, which were sorted into A and B genomes as well as compared to the genetic map of tef. The draft genome was used to identify novel SSR markers, investigate target genes for abiotic stress resistance studies, and understand the evolution of the prolamin family of proteins that are responsible for the immune response to gluten. Conclusions: It is highly plausible that breeding targets previously identified in other cereal crops will also be valuable breeding targets in tef. The draft genome and transcriptome will be of great use for identifying these targets for genetic improvement of this orphan crop that is vital for feeding 50 million people in the Horn of Africa.

Whole-Genome Sequencing of a Canine Family Trio Reveals a FAM83G Variant Associated with Hereditary Footpad Hyperkeratosis

Over 250 Mendelian traits and disorders, caused by rare alleles have been mapped in the canine genome. Although each disease is rare in the dog as a species, they are collectively common and have major impact on canine health. With SNP-based genotyping arrays, genome-wide association studies (GWAS) have proven to be a powerful method to map the genomic region of interest when 10-20 cases and 10-20 controls are available. However, to identify the genetic variant in associated regions, fine-mapping and targeted re-sequencing is required. Here we present a new approach using whole-genome sequencing (WGS) of a family trio without prior GWAS. As a proof-of-concept, we chose an autosomal recessive disease known as hereditary footpad hyperkeratosis (HFH) in Kromfohrl änder dogs. To our knowledge, this is the first time this family trio WGS-approach, has successfully been used to identify a genetic variant that perfectly segregates with a canine disorder. The sequencing of three Kromfohrl änder dogs from a family trio (an affected offspring and both its healthy parents) resulted in an average genome coverage of 9.2X per individual. After applying stringent filtering criteria for candidate causative coding variants, 527 single nucleotide variants (SNVs) and 15 indels were found to be homozygous in the affected offspring and heterozygous in the parents. Using the computer software packages ANNOVAR and SIFT to functionally annotate coding sequence differences and to predict their functional effect, resulted in seven candidate variants located in six different genes. Of these, only FAM83G:c155G>C (p.R52P) was found to be concordant in eight additional cases and 16 healthy Kromfohrl änder dogs.

Complex genetic background in a large family with Brugada syndrome.

The Brugada syndrome (BrS) is an inherited arrhythmia characterized by ST-segment elevation in V1-V3 leads and negative T wave on standard ECG. BrS patients are at risk of sudden cardiac death (SCD) due to ventricular tachyarrhythmia. At least 17 genes have been proposed to be linked to BrS, although recent findings suggested a polygenic background. Mutations in SCN5A, the gene coding for the cardiac sodium channel Nav1.5, have been found in 15-30% of index cases. Here, we present the results of clinical, genetic, and expression studies of a large Iranian family with BrS carrying a novel genetic variant (p.P1506S) in SCN5A. By performing whole-cell patch-clamp experiments using HEK293 cells expressing wild-type (WT) or p.P1506S Nav1.5 channels, hyperpolarizing shift of the availability curve, depolarizing shift of the activation curve, and hastening of the fast inactivation process were observed. These mutant-induced alterations lead to a loss of function of Nav1.5 and thus suggest that the p.P1506S variant is pathogenic. In addition, cascade familial screening found a family member with BrS who did not carry the p.P1506S mutation. Additional next generation sequencing analyses revealed the p.R25W mutation in KCNH2 gene in SCN5A-negative BrS patients. These findings illustrate the complex genetic background of BrS found in this family and the possible pathogenic role of a new SCN5A genetic variant.

Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^

IDENTIFYING GENETIC VARIANTS AND CHARACTERIZING THEIR ROLE IN CLUBFOOT

Clubfoot is a common, complex birth defect affecting 4,000 newborns in the United States and 135,000 world-wide each year. The clubfoot deformity is characterized by inward and rigid downward displacement of one or both feet, along with persistent calf muscle hypoplasia. Despite strong evidence for a genetic liability, there is a limited understanding of the genetic and environmental factors contributing to the etiology of clubfoot. The studies described in this dissertation were performed to identify variants and/or genes associated with clubfoot. Genome-wide linkage scan performed on ten multiplex clubfoot families identified seven new chromosomal regions that provide new areas to search for clubfoot genes. Troponin C (TNNC2) the strongest candidate gene, located in 20q12-q13.11, is involved in muscle contraction. Exon sequencing of TNNC2 did not identify any novel coding variants. Interrogation of fifteen muscle contraction genes found strong associations with SNPs located in potential regulatory regions of TPM1 (rs4075583 and rs3805965), TPM2 (rs2025126 and rs2145925) and TNNC2 (rs383112 and rs437122). In previous studies, a strong association was found with rs3801776 located in the basal promoter of HOXA9, a gene also involved in muscle development and patterning. Altogether, this data suggests that SNPs located in potential regulatory regions of genes involved in muscle development and function could alter transcription factor binding leading to changes in gene expression. Functional analysis of 3801776/HOXA9, rs2025126/TPM2 and rs2145925/TPM2 showed altered protein binding, which significantly influenced promoter activity. Although the ancestral allele (G) of rs4075583/TPM1 creates a DNA-protein complex, it did not affect TPM1 promoter activity. However and importantly, in the context of a haplotype, rs4075583/G significantly decreased TPM1 promoter activity. These results suggest dysregulation of multiple skeletal muscle genes, TPM1, TPM2, TNNC2 and HOXA9, working in concert may contribute to clubfoot. However, specific allelic combinations involving these four regulatory SNPs did not confer a significantly higher risk for clubfoot. Other combinations of these variants are being evaluated. Moreover, these variants may interact with yet to be discovered variants in other genes to confer a higher clubfoot risk. Collectively, we show novel evidence for the role of skeletal muscle genes in clubfoot indicating that there are multiple genetic factors contributing to this complex birth defect.