High resolution assembly and characterization of genomes of Canadian isolates of Salmonella Enteritidis

Background

There is a need to characterize genomes of the foodborne pathogen, Salmonella enterica serovar Enteritidis (SE) and identify genetic information that could be ultimately deployed for differentiating strains of the organism, a need that is yet to be addressed mainly because of the high degree of clonality of the organism. In an effort to achieve the first characterization of the genomes of SE of Canadian origin, we carried out massively parallel sequencing of the nucleotide sequence of 11 SE isolates obtained from poultry production environments (n = 9), a clam and a chicken, assembled finished genomes and investigated diversity of the SE genome.

Results

The median genome size was 4,678,683 bp. A total of 4,833 chromosomal genes defined the pan genome of our field SE isolates consisting of 4,600 genes present in all the genomes, i.e., core genome, and 233 genes absent in at least one genome (accessory genome). Genome diversity was demonstrable by the presence of 1,360 loci showing single nucleotide polymorphism (SNP) in the core genome which was used to portray the genetic distances by means of a phylogenetic tree for the SE isolates. The accessory genome consisted mostly of previously identified SE prophage sequences as well as two, apparently full- sized, novel prophages namely a 28 kb sequence provisionally designated as SE-OLF-10058 (3) prophage and a 43 kb sequence provisionally designated as SE-OLF-10012 prophage.

Conclusions

The number of SNPs identified in the relatively large core genome of SE is a reflection of substantial diversity that could be exploited for strain differentiation as shown by the development of an informative phylogenetic tree. Prophage sequences can also be exploited for SE strain differentiation and lineage tracking. This work has laid the ground work for further studies to develop a readily adoptable laboratory test for the subtyping of SE.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-713) contains supplementary material, which is available to authorized users.     

Genetic association of single nucleotide polymorphisms in dystrobrevin binding protein 1 gene with schizophrenia in a Malaysian population

Dystrobrevin binding protein 1 (DTNBP1) gene is pivotal in regulating the glutamatergic system. Genetic variants of the DTNBP1 affect cognition and thus may be particularly relevant to schizophrenia. We therefore evaluated the association of six single nucleotide polymorphisms (SNPs) with schizophrenia in a Malaysian population (171 cases; 171 controls). Associations between these six SNPs and schizophrenia were tested in two stages. Association signals with p < 0.05 and minor allele frequency > 0.05 in stage 1 were followed by genotyping the SNPs in a replication phase (stage 2). Genotyping was performed with sequenced specific primer (PCR-SSP) and restriction fragment length polymorphism (PCR-RFLP). In our sample, we found significant associations between rs2619522 (allele p = 0.002, OR = 1.902, 95%CI = 1.266 – 2.859; genotype p = 0.002) and rs2619528 (allele p = 0.008, OR = 1.606, 95%CI = 1.130 – 2.281; genotype p = 6.18 × 10⁻⁵) and schizophrenia. Given that these two SNPs may be associated with the pathophysiology of schizophrenia, further studies on the other DTNBP1 variants are warranted.     

Comparative genomics of koala,cattle and sheep strains of Chlamydia pecorum

Background

Chlamydia pecorum is an important pathogen of domesticated livestock including sheep, cattle and pigs. This pathogen is also a key factor in the decline of the koala in Australia. We sequenced the genomes of three koala C. pecorum strains, isolated from the urogenital tracts and conjunctiva of diseased koalas. The genome of the C. pecorum VR629 (IPA) strain, isolated from a sheep with polyarthritis, was also sequenced.

Results

Comparisons of the draft C. pecorum genomes against the complete genomes of livestock C. pecorum isolates revealed that these strains have a conserved gene content and order, sharing a nucleotide sequence similarity > 98%. Single nucleotide polymorphisms (SNPs) appear to be key factors in understanding the adaptive process. Two regions of the chromosome were found to be accumulating a large number of SNPs within the koala strains. These regions include the Chlamydia plasticity zone, which contains two cytotoxin genes (toxA and toxB), and a 77 kbp region that codes for putative type III effector proteins. In one koala strain (MC/MarsBar), the toxB gene was truncated by a premature stop codon but is full-length in IPTaLE and DBDeUG. Another five pseudogenes were also identified, two unique to the urogenital strains C. pecorum MC/MarsBar and C. pecorum DBDeUG, respectively, while three were unique to the koala C. pecorum conjunctival isolate IPTaLE. An examination of the distribution of these pseudogenes in C. pecorum strains from a variety of koala populations, alongside a number of sheep and cattle C. pecorum positive samples from Australian livestock, confirmed the presence of four predicted pseudogenes in koala C. pecorum clinical samples. Consistent with our genomics analyses, none of these pseudogenes were observed in the livestock C. pecorum samples examined. Interestingly, three SNPs resulting in pseudogenes identified in the IPTaLE isolate were not found in any other C. pecorum strain analysed, raising questions over the origin of these point mutations.

Conclusions

The genomic data revealed that variation between C. pecorum strains were mainly due to the accumulation of SNPs, some of which cause gene inactivation. The identification of these genetic differences will provide the basis for further studies to understand the biology and evolution of this important animal pathogen.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-667) contains supplementary material, which is available to authorized users.     

Molecular characterization of porcine MMP19 and MMP23B genes and its association with immune traits

MMP19 and MMP23B belong to the Matrix metalloproteases (MMPs) family, which are zinc-binding endopeptidases that are capable of degrading various components of the extracellular matrix. They are thought to play important roles in embryonic development, reproduction and tissue remodeling, as well as in cell proliferation, differentiation, migration, angiogenesis, apoptosis and host defense. However, they are poorly understood in pigs. Here, we obtained the full length coding region sequence and genomic sequence of the porcine MMP19 and MMP23B genes and analyzed their genomic structures. The deduced amino acid sequence shares similar precursor protein domains with human and mouse MMP19 and MMP23B protein, respectively. Using IMpRH panel, MMP19 was mapped to SSC5p12-q11 (closely linked to microsatellite DK) and MMP23B was mapped to SSC8q11-q12 (linked to microsatellite Sw2521). Quantitative real-time PCR showed that MMP19 was abundantly expressed in the liver, while MMP23B was strongly expressed in the ovarian and heart. Furthermore, both genes were all expressed increasingly in prenatal skeletal muscle during development. Three SNPs were detected by sequencing and PCR-RFLP methods, and association analysis indicated that C203T at exon 5 of MMP19 has a significant association with the blood parameters WBC (G/L) and IgG2 (mg/mL) (P<0.05), SNP C131T at exon 3 of MMP23B is significantly associated with the blood parameters HGB (g/L) and MCH (P<0.05), and A150G in exon 4 has no significant association with the economic traits in pigs.     

Identification of calreticulin as a marker for phagocytosis of apoptotic cells in Drosophila

Apoptotic cell phagocytosis is initiated through the specific interaction between markers for phagocytosis present at the surface of targets and their receptors of phagocytes. Although many molecules have been proposed to be phagocytosis markers and receptors in mammals, information as to the identity of those molecules is limited for invertebrate animals. Calreticulin, a molecular chaperone that functions in the lumen of the endoplasmic reticulum, was recently reported to be the second general marker, the membrane phospholipid phosphatidylserine being the first, for mammalian apoptotic cells to be recognized by phagocytes. We here asked whether or not calreticulin serves as a marker for phagocytosis in Drosophila. Phagocytosis of apoptotic S2 cells by Drosophila hemocyte-derived l(2)mbn cells, which we previously showed to occur independent of phosphatidylserine, was inhibited by the addition of anti-calreticulin antibody. This inhibition was observed when the target cells, but not phagocytes, were pre-incubated with the antibody. In addition, RNA interference-mediated reduction of calreticulin expression in apoptotic S2 cells, but not in l(2)mbn cells, reduced the level of phagocytosis. An immunocytochemical analysis revealed that calreticulin is widely distributed at the surface of viable S2 cells. After the induction of apoptosis, cell surface calreticulin seemed to form aggregates, with no change in its amount. Furthermore, in embryos of a mutant Drosophila strain that expresses calreticulin at a reduced level, the level of phagocytosis of apoptotic cells was about a half of that observed in embryos of a wild-type strain. These results collectively indicate that calreticulin is the first molecule to be identified as a marker for phagocytosis of apoptotic cells by Drosophila phagocytes.     

Genome-wide association study for flowering time,maturity dates and plant height in early maturing soybean (Glycine max) germplasm

Background

Soybean (Glycine max) is a photoperiod-sensitive and self-pollinated species. Days to flowering (DTF) and maturity (DTM), duration of flowering-to-maturity (DFTM) and plant height (PH) are crucial for soybean adaptability and yield. To dissect the genetic architecture of these agronomically important traits, a population consisting of 309 early maturity soybean germplasm accessions was genotyped with the Illumina Infinium SoySNP50K BeadChip and phenotyped in multiple environments. A genome-wide association study (GWAS) was conducted using a mixed linear model that involves both relative kinship and population structure.

Results

The linkage disequilibrium (LD) decayed slowly in soybean, and a substantial difference in LD pattern was observed between euchromatic and heterochromatic regions. A total of 27, 6, 18 and 27 loci for DTF, DTM, DFTM and PH were detected via GWAS, respectively. The Dt1 gene was identified in the locus strongly associated with both DTM and PH. Ten candidate genes homologous to Arabidopsis flowering genes were identified near the peak single nucleotide polymorphisms (SNPs) associated with DTF. Four of them encode MADS-domain containing proteins. Additionally, a pectin lyase-like gene was also identified in a major-effect locus for PH where LD decayed rapidly.

Conclusions

This study identified multiple new loci and refined chromosomal regions of known loci associated with DTF, DTM, DFTM and/or PH in soybean. It demonstrates that GWAS is powerful in dissecting complex traits and identifying candidate genes although LD decayed slowly in soybean. The loci and trait-associated SNPs identified in this study can be used for soybean genetic improvement, especially the major-effect loci associated with PH could be used to improve soybean yield potential. The candidate genes may serve as promising targets for studies of molecular mechanisms underlying the related traits in soybean.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1441-4) contains supplementary material, which is available to authorized users.     

Cloning HSP70 and HSP90 genes of kaluga (<Emphasis Type="Italic">Huso dauricus</Emphasis>) and the effects of temperature and salinity stress on their gene expression

The genes encoding HSP70 and HSP90 proteins were isolated from kaluga by homologous cloning and rapid amplification of complementary DNA (cDNA) ends (RACE). HSP70 (GenBank accession no. {"type":"entrez-nucleotide","attrs":{"text":"KP050541","term_id":"756558159","term_text":"KP050541"}}KP050541) and HSP90 (GenBank accession no. {"type":"entrez-nucleotide","attrs":{"text":"KP050542","term_id":"756558161","term_text":"KP050542"}}KP050542) cDNAs were composed of 2275 and 2718 bp and encoded polypeptides of 650 and 725 amino acids, respectively. Basic Local Alignment Search Tool (BLAST) analysis showed that HSP70 and HSP90 of kaluga shared high identities with those of Acipenser ruthenus, Acipenser schrenckii, and Acipenser baerii (98–99 %). Fluorescent real-time RT-PCR under unstressed conditions revealed that HSP70 and HSP90 were expressed in 11 different tissues of kaluga. Messenger RNA (mRNA) expressions of both HSP70 and HSP90 were highest in the intestine and lowest in the muscle. In addition, the patterns of mRNA expression of HSP70 and HSP90 were similar, although the level of expression was more in HSP90 than in HSP70 (P < 0.05).We also analyzed patterns of HSP70 and HSP90 expression in the muscle, gill, and liver of kaluga under different combinations of temperature and salinity stress, including temperatures of 4,10, 25, and 28 °C at 0 ppt salinity, and salinities of 10, 20, 30, and 40 ppt at 16 °C, where 16 °C at 0 ppt (parts per thousand) served as the control. We found that levels of mRNA expression of both HSP70 and HSP90 were highest at 4 °C in the muscle, gill, and liver and changed little with salinity stress. These results increase understanding of the mechanisms of stress response of cold freshwater fish.     

Expression of calreticulin P-domain results in impairment of secretory pathway in Leishmania donovani and reduced parasite survival in macrophages

The secretory proteins of Leishmania are thought to be involved in the parasite survival inside the insect vector or mammalian host. It is clear from studies in higher eukaryotes that proper folding in the endoplasmic reticulum and targeting out of the endoplasmic reticulum is critical for the function of secretory proteins. The endoplasmic reticulum chaperones such as calreticulin play an important role in the quality control of secretory proteins. However, very little is known about the secretory pathway of trypanosomatid parasites such as Leishmania. In the present study, we show that overexpression of the P-domain of Leishmania donovani calreticulin in transfected L. donovani resulted in a significant reduction in the secretion of the parasite secretory acid phosphatases. This effect is associated with an intracellular accumulation of active enzyme in these transfected parasites. In addition, parasites expressing the P-domain calreticulin showed a significant decrease in survival inside human macrophages. This study suggests that altering the function of an endoplasmic reticulum chaperone such as calreticulin in Leishmania may affect the targeting of proteins that are associated with the virulence of the parasite during their trafficking through the parasite secretory pathway.     

Variant rs2237892 of KCNQ1 Is Potentially Associated with Hypertension and Macrovascular Complications in Type 2 Diabetes Mellitus in A Chinese Han Population

KCNQ1 has been identified as a susceptibility gene of type 2 diabetes mellitus(T2DM) in Asian populations through genome-wide association studies. However, studies on the association between gene polymorphism of KCNQ1 and T2 DM complications remain unclear. To further analyze the association between different alleles at the single nucleotide polymorphism(SNP) rs2237892 within KCNQ1 and TD2 M and its complications, we conducted a case-control study in a Chinese Han population. The C allele of rs2237892 variant contributed to susceptibility to T2DM(odds ratio [OR], 1.45; 95% confidence interval [CI], 1.20–1.75). Genotypes CT(OR, 1.97; 95% CI,1.24–3.15) and CC(OR, 2.49; 95% CI, 1.57–3.95) were associated with an increased risk of T2 DM. Multivariate regression analysis was performed with adjustment of age, gender, and body mass index. We found that systolic blood pressure(P = 0.015), prevalence of hypertension(P = 0.037), and risk of macrovascular disease(OR, 2.10; CI, 1.00–4.45) were significantly higher in subjects with the CC genotype than in the combined population with genotype either CT or     

Cloning and characterization of Giardia intestinalis cyclophilin

The cyclophilins (Cyps) are family members of proteins that exhibit peptidylprolyl cis-trans isomerase (PPIase, EC 5.2.1.8) activity and bind the immunosuppressive agent cyclosprin A (CsA) in varying degrees. During the process of random sequencing of a cDNA library made from Giardia intestinalis WB strain, the cyclophilin gene (gicyp 1) was isolated. An open reading frame of gicyp 1 gene was 576 nucleotides, which corresponded to a translation product of 176 amino acids (Gicyp 1). The identity with other Cyps was about 58-71%. The 13 residues that constituted the CsA binding site of human cyclophilin were also detected in the amino acid sequence of Gicyp 1, including tryptophan residue essential for the drug binding. The single copy of the gicyp 1 gene was detected in the G. intestinalis chromosome by southern hybridization analysis. Recombinant Gicyp 1 protein clearly accelerated the rate of cis-->trans isomerization of the peptide substrate and the catalysis was completely inhibited by the addition of 0.5 microM CsA.     

Low salinity stress experienced by larvae does not affect post-metamorphic growth or survival in three calyptraeid gastropods

Marine larvae that experience some sub-lethal stresses can show effects from those stresses after metamorphosis, even when they seem to recover from those stresses before metamorphosis. In this study we investigated the short and long-term effects of exposing the larvae of three calyptraeid gastropods (Crepidula fornicata, Crepidula onyx, and Crepipatella fecunda) to temporary reductions in salinity. Larvae of all three species showed slower larval growth rates, longer time to metamorphic competence, and substantial mortality after being stressed in seawater at salinities of 10, 15, and 20 for less than 48 h. Larval tolerance to low salinities varied widely within and among species, but longer stresses at lower salinities were generally more harmful to larvae. However, larvae in nearly all experiments that were able to metamorphose survived and grew normally as juveniles; there were no documented “latent effects.” For all three species, starving larvae in full-strength seawater was not as harmful as exposing larvae to low salinity stress, indicating that detrimental effects on larvae were caused by the salinity stress per se, rather than by an indirect effect of salinity stress on feeding. C. fornicata that were stressed with low salinity as juveniles were more tolerant of the stress than larvae: all stressed juveniles lived and showed reduced growth rates for no more than 3 days. Our data suggest that even though reduced salinity is clearly stressful to the larvae of these 3 gastropod species, metamorphosis seems to generally provide individuals with a fresh start.     

A marine killer yeast against the pathogenic yeast strain in crab (Portunus trituberculatus) and an optimization of the toxin production

A pathogenic yeast strain WCY which could cause milky disease in Portunus trituberculatus was identified to be Metschnikowia bicuspidate according to the results of routine yeast identification and 18S rDNA and ITS sequences. After screening of more than 300 yeast strains from different sources in marine environments, it was found that strain YF07b had the highest ability to produce killer toxin against the pathogenic yeast. Strain YF07b was identified to be Pichia anomala according to the results of routine yeast identification and 18S rDNA and ITS sequences. The optimal conditions for killer toxin production by strain YF07b were the production medium with 2.0% NaCl, pH 4.5, cultivation temperature of 20 degrees C and the optimal conditions for action of the crude killer toxin against the pathogenic yeast were the assay medium with 6.0% NaCl, pH 4.5 and temperature 15 degrees C.     

Molecular characterization of calreticulin from Anopheles stephensi midgut cells and functional assay of the recombinant calreticulin with Plasmodium berghei ookinetes

Transmission blocking vaccines (TBVs) that target the antigens on the midgut epithelium of Anopheles mosquitoes are among the promising tools for the elimination of the malaria parasite. Characterization and analysis of effective antigens is the first step to design TBVs. Calreticulin (CRT), a lectin-like protein, from Anopheles albimanus midgut, has shown antigenic features, suggesting a promising and novel TBV target. CRT is a highly conserved protein with similar features in vertebrates and invertebrates including anopheline. We cloned the full-length crt gene from malaria vector, Anopheles stephensi (AsCrt) and explored the interaction of recombinant AsCrt protein, expressed in a prokaryotic system (pGEX-6p-1), with surface proteins of Plasmodium berghei ookinetes by immunofluorescence assay. The cellular localization of AsCrt was determined using the baculovirus expression system. Sequence analysis of the whole cDNA of AsCrt revealed that AsCrt contains an ORF of 1221 bp. The amino acid sequence of AsCrt protein obtained in this study showed 64% homology with similar protein in human. The AsCrt shares the most common features of CRTs from other species. This gene encodes a 406 amino-acid protein with a molecular mass of 46 kDa, which contains a predicted 16 amino-acid signal peptides, conserved cysteine residues, a proline-rich region, and highly acidic C-terminal domain with endoplasmic reticulum retrieval sequence HDEL. The production of GST-AsCrt recombinant protein was confirmed by Western blot analysis using an antibody against the GST protein. The FITC-labeled GST-AsCrt exhibited a significant interaction with P. berghei ookinete surface proteins. Purified recombinant GST-AsCrt, labeled with FITC, displayed specific binding to the surface of P. berghei ookinetes in comparison with control. Moreover, the expression of AsCrt in baculovirus expression system indicated that AsCrt was localized on the surface of Sf9 cells. Our results suggest that AsCrt could be utilized as a potential target for future studies in TBV area for malaria control.