Bioinformatics analysis of the Pedobacter sp. NL19 genome

The last decades of the 20th century defined the genetic engineering advent, climaxing in the development of techniques, such as PCR and Sanger sequencing. This, permitted the appearance of new techniques to sequencing whole genomes, identified as next-generation sequencing. One of the many applications of these techniques is the in silico search for new secondary metabolites, synthesized by microorganisms exhibiting antimicrobial properties. The peptide antibiotics compounds can be classified in two classes, according to their biosynthesis, in ribosomal or nonribosomal peptides. Lanthipeptides are the most studied ribosomal peptides and are characterized by the presence of lanthionine and methylanthionine that result from posttranslational modifications. Lanthipeptides are divided in four classes, depending on their biosynthetic machinery. In class I, a LanB enzyme dehydrate serine and threonine residues in the C-terminus precursor peptide. Then, these residues undergo a cyclization step performed by a LanC enzyme, forming the lanthionine rings. The cleavage and the transport of the peptide is achieved by the LanP and LanT enzymes, respectively. Although, in class II only one enzyme, LanM, is responsible for the dehydration and cyclization steps and also only one enzyme performs the cleavage and transport, LanT. Pedobacter sp. NL19 is a Gram-negative bacterium, isolated from sludge of an abandon uranium mine, in Viseu (Portugal). Antibacterial activity in vitro was detected against several Gram-positive and Gram-negative bacteria. Sequencing and in silico analysis of NL19 genome revealed the presence of 21 biosynthetic clusters for secondary metabolites, including nonribosomal and ribosomal peptides biosynthetic clusters. Four lanthipeptides clusters were predicted, comprising the precursor peptides, the modifying enzymes (LanB and LanC), and also a bifunctional LanT. This result revealed the hybrid nature of the clusters, comprising characteristics from two distinct classes, which are poorly described in literature. The phylogenetic analysis of their enzymes showed that they clustered within the bacteroidetes clade. Furthermore, hybrid gene clusters were also found in other species of this phylum, revealing that it is a common characteristic in this group. Finally, the analysis of NL19 colonies by MALDI-TOF MS allowed the identification of a 3180 Da mass that corresponds to the predicted mass of a lanthipeptide encoded in one of the clusters. However, this result is not fully conclusive and further experiments are needed to understand the full potential of the compounds encoded in this type of clusters. In conclusion, it was determined that NL19 strain has the potential to produce diverse secondary metabolites, including lanthipeptides that were not functionally characterized so far.

Genetic variation of temperature-regulated curd induction in cauliflower: elucidation of floral transition by genome-wide association mapping and gene expression analysis

Cauliflower (Brassica oleracea var. botrytis) is a vernalization-responsive crop. High ambient temperatures delay harvest time. The elucidation of the genetic regulation of floral transition is highly interesting for a precise harvest scheduling and to ensure stable market supply. This study aims at genetic dissection of temperature-dependent curd induction in cauliflower by genome-wide association studies and gene expression analysis. To assess temperature dependent curd induction, two greenhouse trials under distinct temperature regimes were conducted on a diversity panel consisting of 111 cauliflower commercial parent lines, genotyped with 14,385 SNPs. Broad phenotypic variation and high heritability (0.93) were observed for temperature-related curd induction within the cauliflower population. GWA mapping identified a total of 18 QTL localized on chromosomes O1, O2, O3, O4, O6, O8, and O9 for curding time under two distinct temperature regimes. Among those, several QTL are localized within regions of promising candidate flowering genes. Inferring population structure and genetic relatedness among the diversity set assigned three main genetic clusters. Linkage disequilibrium (LD) patterns estimated global LD extent of r(2) = 0.06 and a maximum physical distance of 400 kb for genetic linkage. Transcriptional profiling of flowering genes FLOWERING LOCUS C (BoFLC) and VERNALIZATION 2 (BoVRN2) was performed, showing increased expression levels of BoVRN2 in genotypes with faster curding. However, functional relevance of BoVRN2 and BoFLC2 could not consistently be supported, which probably suggests to act facultative and/or might evidence for BoVRN2/BoFLC-independent mechanisms in temperature regulated floral transition in cauliflower. Genetic insights in temperature-regulated curd induction can underpin genetically informed phenology models and benefit molecular breeding strategies toward the development of thermo-tolerant cultivars.

Genome-scale analysis of the non-cultivable Treponema pallidum reveals extensive within-patient genetic variation

Insights into the genomic adaptive traits of Treponema pallidum, the causative bacterium of syphilis, have long been hampered due to the absence of in vitro culture models and the constraints associated with its propagation in rabbits. Here, we have bypassed the culture bottleneck by means of a targeted strategy never applied to uncultivable bacterial human pathogens to directly capture whole-genome T. pallidum data in the context of human infection. This strategy has unveiled a scenario of discreet T. pallidum interstrain single-nucleotide-polymorphism-based microevolution, contrasting with a rampant within-patient genetic heterogeneity mainly targeting multiple phase-variable loci and a major antigen-coding gene (tprK). TprK demonstrated remarkable variability and redundancy, intra- and interpatient, suggesting ongoing parallel adaptive diversification during human infection. Some bacterial functions (for example, flagella- and chemotaxis-associated) were systematically targeted by both inter- and intrastrain single nucleotide polymorphisms, as well as by ongoing within-patient phase variation events. Finally, patient-derived genomes possess mutations targeting a penicillin-binding protein coding gene (mrcA) that had never been reported, unveiling it as a candidate target to investigate the impact on the susceptibility to penicillin. Our findings decode the major genetic mechanisms by which T. pallidum promotes immune evasion and survival, and demonstrate the exceptional power of characterizing evolving pathogen subpopulations during human infection.

Functional XPB/RAD25 redundancy in Arabidopsis genome: characterization of AtXPB2 and expression analysis

The xeroderma pigmentosum complementation group B (XPB) protein is involved in both DNA repair and transcription in human cells. It is a component of the transcription factor IIH (TFIIH) and is responsible for DNA helicase activity during nucleotide (nt) excision repair (NER). Its high evolutionary conservation has allowed identification of homologous proteins in different organisms, including plants. In contrast to other organisms, Arabidopsis thaliana harbors a duplication of the XPB orthologue (AtXPB1 and AtXPB2), and the proteins encoded by the duplicated genes are very similar (95% amino acid identity). Complementation assays in yeast rad25 mutant strains suggest the involvement of AtXPB2 in DNA repair, as already shown for AtXPB1, indicating that these proteins may be functionally redundant in the removal of DNA lesions in A. thaliana. Although both genes are expressed in a constitutive manner during the plant life cycle, Northern blot analyses suggest that light modulates the expression level of both XPB copies, and transcript levels increase during early stages of development. Considering the high similarity between AtXPB1 and AtXPB2 and that both of predicted proteins may act in DNA repair, it is possible that this duplication may confer more flexibility and resistance to DNA damaging agents in thale cress. (C) 2004 Elsevier B.V. All rights reserved.

Genome-wide copy number variation (CNV) detection in Nelore cattle reveals highly frequent variants in genome regions harboring QTLs affecting production traits.

Background: Copy number variations (CNVs) have been shown to account for substantial portions of observed genomic variation and have been associated with qualitative and quantitative traits and the onset of disease in a number of species. Information from high-resolution studies to detect, characterize and estimate population-specific variant frequencies will facilitate the incorporation of CNVs in genomic studies to identify genes affecting traits of importance. Results: Genome-wide CNVs were detected in high-density single nucleotide polymorphism (SNP) genotyping data from 1,717 Nelore (Bos indicus) cattle, and in NGS data from eight key ancestral bulls. A total of 68,007 and 12,786 distinct CNVs were observed, respectively. Cross-comparisons of results obtained for the eight resequenced animals revealed that 92 % of the CNVs were observed in both datasets, while 62 % of all detected CNVs were observed to overlap with previously validated cattle copy number variant regions (CNVRs). Observed CNVs were used for obtaining breed-specific CNV frequencies and identification of CNVRs, which were subsequently used for gene annotation. A total of 688 of the detected CNVRs were observed to overlap with 286 non-redundant QTLs associated with important production traits in cattle. All of 34 CNVs previously reported to be associated with milk production traits in Holsteins were also observed in Nelore cattle. Comparisons of estimated frequencies of these CNVs in the two breeds revealed 14, 13, 6 and 14 regions in high (>20 %), low (<20 %) and divergent (NEL > HOL, NEL < HOL) frequencies, respectively. Conclusions: Obtained results significantly enriched the bovine CNV map and enabled the identification of variants that are potentially associated with traits under selection in Nelore cattle, particularly in genome regions harboring QTLs affecting production traits.

Genome-Wiide estimates of coancestry, in breeding and effective population size in the Spanish Holstein population

Estimates of effective population size in the Holstein cattle breed have usually been low despite the large number of animals that constitute this breed. Effective population size is inversely related to the rates at which coancestry and inbreeding increase and these rates have been high as a consequence of intense and accurate selection. Traditionally, coancestry and inbreeding coefficients have been calculated from pedigree data. However, the development of genome-wide single nucleotide polymorphisms has increased the interest of calculating these coefficients from molecular data in order to improve their accuracy. In this study, genomic estimates of coancestry, inbreeding and effective population size were obtained in the Spanish Holstein population and then compared with pedigree-based estimates. A total of 11,135 animals genotyped with the Illumina BovineSNP50 BeadChip were available for the study. After applying filtering criteria, the final genomic dataset included 36,693 autosomal SNPs and 10,569 animals. Pedigree data from those genotyped animals included 31,203 animals. These individuals represented only the last five generations in order to homogenise the amount of pedigree information across animals. Genomic estimates of coancestry and inbreeding were obtained from identity by descent segments (coancestry) or runs of homozygosity (inbreeding). The results indicate that the percentage of variance of pedigree-based coancestry estimates explained by genomic coancestry estimates was higher than that for inbreeding. Estimates of effective population size obtained from genome-wide and pedigree information were consistent and ranged from about 66 to 79. These low values emphasize the need of controlling the rate of increase of coancestry and inbreeding in Holstein selection programmes.

Enterococcus faecalis prophage dynamics and contributions to pathogenic traits

Polylysogeny is frequently considered to be the result of an adaptive evolutionary process in which prophages confer fitness and/or virulence factors, thus making them important for evolution of both bacterial populations and infectious diseases. The Enterococcus faecalis V583 isolate belongs to the high-risk clonal complex 2 that is particularly well adapted to the hospital environment. Its genome carries 7 prophage-like elements (V583-pp1 to -pp7), one of which is ubiquitous in the species. In this study, we investigated the activity of the V583 prophages and their contribution to E. faecalis biological traits. We systematically analyzed the ability of each prophage to excise from the bacterial chromosome, to replicate and to package its DNA. We also created a set of E. faecalis isogenic strains that lack from one to all six non-ubiquitous prophages by mimicking natural excision. Our work reveals that prophages of E. faecalis V583 excise from the bacterial chromosome in the presence of a fluoroquinolone, and are able to produce active phage progeny. Intricate interactions between V583 prophages were also unveiled: i) pp7, coined EfCIV583 for E. faecalis chromosomal island of V583, hijacks capsids from helper phage 1, leading to the formation of distinct virions, and ii) pp1, pp3 and pp5 inhibit excision of pp4 and pp6. The hijacking exerted by EfCIV583 on helper phage 1 capsids is the first example of molecular piracy in Gram positive bacteria other than staphylococci. Furthermore, prophages encoding platelet-binding-like proteins were found to be involved in adhesion to human platelets, considered as a first step towards the development of infective endocarditis. Our findings reveal not only a role of E. faecalis V583 prophages in pathogenicity, but also provide an explanation for the correlation between antibiotic usage and E. faecalis success as a nosocomial pathogen, as fluoriquinolone may provoke release of prophages and promote gene dissemination among isolates.

The exploration of potato-associated bacteria in the Central Andean Highlands and their application in integrated crop management systems

Potato is the most important food crop after wheat and rice. A changing climate, coupled with a heightened consumer awareness of how food is produced and legislative changes governing the usage of agrochemicals, means that alternative more integrated and sustainable approaches are needed for crop management practices. Bioprospecting in the Central Andean Highlands resulted in the isolation and in vitro screening of 600 bacterial isolates. The best performing isolates, under in vitro conditions, were field trialled in their home countries. Six of the isolates, Pseudomonas sp. R41805 (Bolivia), Pseudomonas palleroniana R43631 (Peru), Bacillus sp. R47065, R47131, Paenibacillus sp. B3a R49541, and Bacillus simplex M3-4 R49538 (Ecuador), showed significant increase in the yield of potato. Using – omic technologies (i.e. volatilomic, transcriptomic, proteomic and metabolomic), the influence of microbial isolates on plant defence responses was determined. Volatile organic compounds of bacterial isolates were identified using GC/MS. RT-qPCR analysis revealed the significant expression of Ethylene Response Factor 3 (ERF3) and the results of this study suggest that the dual inoculation of potato with Pseudomonas sp. R41805 and Rhizophagus irregularis MUCL 41833 may play a part in the activation of plant defence system via ERF3. The proteomic analysis by 2-DE study has shown that priming by Pseudomonas sp. R41805 can induce the expression of proteins related to photosynthesis and protein folding in in vitro potato plantlets. The metabolomics study has shown that the total glycoalkaloid (TGA) content of greenhouse-grown potato tubers following inoculation with Pseudomonas sp. R41805 did not exceed the acceptable safety limit (200 mg kg-1 FW). As a result of this study, a number of bacteria have been identified with commercial potential that may offer sustainable alternatives in both Andean and European agricultural settings.

Progress in understanding of the molecular basis underlying functional diversification of cyclic di-nucleotide turnover proteins

Cyclic di-GMP was the first cyclic di-nucleotide second messenger described, presaging the discovery of additional cyclic di-nucleotide messengers in bacteria and eukaryotes. The GGDEF diguanylate cyclase (DGC) and EAL and HD-GYP phosphodiesterase (PDE) domains conduct the turnover of cyclic di-GMP. These three unrelated domains belong to superfamilies that exhibit significant variations in function, to include both enzymatically active and inactive members with a subset involved in synthesis and degradation of other cyclic di-nucleotides. Here we summarize current knowledge of sequence and structural varitions that underpin the functional diversification of cyclic di-GMP turnover proteins. Moreover, we highlight that superfamily diversification is not restricted to cyclic di-GMP signaling domains, as particular DHH/DHHA1 domain and HD domain proteins have been shown to act as cyclic di-AMP phosphodiesterases. We conclude with a consideration of the current limitations that such diversity of action places on bioinformatic prediction of the roles of GGDEF, EAL and HD-GYP domain proteins.

Heat shock response in bacteria with large genomes: lessons from rhizobia

Rhizobia are important soil bacteria due to their ability to establish nitrogen-fixing symbioses with legume plants. In this dual lifestyle, as free-living bacteria or as plant symbiont, rhizobia are often exposed to different environmental stresses. The present chapter overviews the current knowledge on the heat shock response of rhizobia, highlighting how these large genome bacteria respond to heat from a transcriptional point of view. Response to heat shock in rhizobia involves genome wide changes in the transcriptome that may affect more than 30% of the genome and involve all replicons. In addition to the expected upregulation of genes already known to be involved in stress response (dnaK, groEL, ibpA, clpB), the reports on the heat shock response in rhizobia also showed particular aspects of stress response in these resourceful bacteria. The transcriptional response to heat in rhizobia includes the overexpression of a large number of genes involved in transcription and carbohydrate transport and metabolism. Additional studies are needed in order to better understand the transcriptional regulation of stress response in bacteria with large genomes.

Multicentric genome-wide association study for primary spontaneous pneumothorax

Despite elevated incidence and recurrence rates for Primary Spontaneous Pneumothorax (PSP), little is known about its etiology, and the genetics of idiopathic PSP remains unexplored. To identify genetic variants contributing to sporadic PSP risk, we conducted the first PSP genome-wide association study. Two replicate pools of 92 Portuguese PSP cases and of 129 age- and sex-matched controls were allelotyped in triplicate on the Affymetrix Human SNP Array 6.0 arrays. Markers passing quality control were ranked by relative allele score difference between cases and controls (|RASdiff|), by a novel cluster method and by a combined Z-test. 101 single nucleotide polymorphisms (SNPs) were selected using these three approaches for technical validation by individual genotyping in the discovery dataset. 87 out of 94 successfully tested SNPs were nominally associated in the discovery dataset. Replication of the 87 technically validated SNPs was then carried out in an independent replication dataset of 100 Portuguese cases and 425 controls. The intergenic rs4733649 SNP in chromosome 8 (between LINC00824 and LINC00977) was associated with PSP in the discovery (P = 4.07E-03, ORC[95% CI] = 1.88[1.22–2.89]), replication (P = 1.50E-02, ORC[95% CI] = 1.50[1.08–2.09]) and combined datasets (P = 8.61E-05, ORC[95% CI] = 1.65[1.29–2.13]). This study identified for the first time one genetic risk factor for sporadic PSP, but future studies are warranted to further confirm this finding in other populations and uncover its functional role in PSP pathogenesis.

he genome and genetics of a high oxidative stress tolerant Serratia sp. LCN16 isolated from the plant parasitic nematode Bursaphelenchus xylophilus

Background: Pine wilt disease (PWD) is a worldwide threat to pine forests, and is caused by the pine wood nematode (PWN) Bursaphelenchus xylophilus. Bacteria are known to be associated with PWN and may have an important role in PWD. Serratia sp. LCN16 is a PWN-associated bacterium, highly resistant to oxidative stress in vitro, and which beneficially contributes to the PWN survival under these conditions. Oxidative stress is generated as a part of the basal defense mechanism used by plants to combat pathogenic invasion. Here, we studied the biology of Serratia sp. LCN16 through genome analyses, and further investigated, using reverse genetics, the role of two genes directly involved in the neutralization of H2O2, namely the H2O2 transcriptional factor oxyR; and the H2O2-targeting enzyme, catalase katA. Results: Serratia sp. LCN16 is phylogenetically most closely related to the phytosphere group of Serratia, which includes S. proteamaculans, S. grimessi and S. liquefaciens. Likewise, Serratia sp. LCN16 shares many features with endophytes (plant-associated bacteria), such as genes coding for plant polymer degrading enzymes, iron uptake/ transport, siderophore and phytohormone synthesis, aromatic compound degradation and detoxification enzymes. OxyR and KatA are directly involved in the high tolerance to H2O2 of Serratia sp. LCN16. Under oxidative stress, Serratia sp. LCN16 expresses katA independently of OxyR in contrast with katG which is under positive regulation of OxyR. Serratia sp. LCN16 mutants for oxyR (oxyR::int(614)) and katA (katA::int(808)) were sensitive to H2O2 in relation with wild-type, and both failed to protect the PWN from H2O2-stress exposure. Moreover, both mutants showed different phenotypes in terms of biofilm production and swimming/swarming behaviors. Conclusions: This study provides new insights into the biology of PWN-associated bacteria Serratia sp. LCN16 and its extreme resistance to oxidative stress conditions, encouraging further research on the potential role of this bacterium in interaction with PWN in planta environment.

Evidences for an opportunistic and endophytic lifestyle of the Bursaphelenchus xylophilus -associated bacteria Serratia marcescens PWN146 isolated from wilting Pinus pinaster

Pine wilt disease (PWD) results from the interaction of three elements: the pathogenic nematode, Bursaphelenchus xylophilus; the insect-vector, Monochamus sp.; and the host tree, mostly Pinus species. Bacteria isolated from B. xylophilus may be a fourth element in this complex disease. However, the precise role of bacteria in this interaction is unclear as both plant-beneficial and as plant-pathogenic bacteria may be associated with PWD. Using whole genome sequencing and phenotypic characterization, we were able to investigate in more detail the genetic repertoire of Serratia marcescens PWN146, a bacterium associated with B. xylophilus. We show clear evidence that S. marcescens PWN146 is able to withstand and colonize the plant environment, without having any deleterious effects towards a susceptible host (Pinus thunbergii), B. xylophilus nor to the nematode model C. elegans. This bacterium is able to tolerate growth in presence of xenobiotic/organic compounds, and use phenylacetic acid as carbon source. Furthermore, we present a detailed list of S. marcescens PWN146 potentials to interfere with plant metabolism via hormonal pathways and/or nutritional acquisition, and to be competitive against other bacteria and/or fungi in terms of resource acquisition or production of antimicrobial compounds. Further investigation is required to understand the role of bacteria in PWD. We have now reinforced the theory that B. xylophilus-associated bacteria may have a plant origin.