Identification of Two Novel Mutations in PKHD1 Gene from Two Families with Polycystic Kidney Disease by Whole Exome Sequencing

Background Polycystic kidney disease (PKD) is an autosomal recessive disorder resulting from mutations in the PKHD1 gene on chromosome 6 (6p12), a large gene spanning 470 kb of genomic DNA.ObjectiveThe aim of the present study was to report newly identified mutations in the PKHD1 gene in two Iranian families with PKD.Materials and Methods Genetic alterations of a 3-month-old boy and a 27-year-old girl with PKD were evaluated using whole-exome sequencing. The PCR direct sequencing was performed to analyse the co-segregation of the variants with the disease in the family. Finally, the molecular function of the identified novel mutations was evaluated by in silico study. ResultsIn the 3 month-old boy, a novel homozygous frameshift mutation was detected in the PKHD1 gene, which can cause PKD. Moreover, we identified three novel heterozygous missense mutations in ATIC, VPS13B, and TP53RK genes. In the 27-year-old woman, with two recurrent abortions history and two infant mortalities at early weeks due to metabolic and/or renal disease, we detected a novel missense mutation on PKHD1 gene and a novel mutation in ETFDH gene.Conclusion In general, we have identified two novel mutations in the PKHD1 gene. These molecular findings can help accurately correlate genotype and phenotype in families with such disease in order to reduce patient births through preoperative genetic diagnosis or better management of disorders.     

In Silico Prediction of Epitopes in Virulence Proteins of Mycobacterium ulcerans for Vaccine Designing

Background Mycobacterium ulcerans is the fundamental agent of the third most common Mycobacterial disease known as Buruli Ulcer (BU). It is an infection of the skin and soft tissue affecting the human population worldwide. Presently, the vaccine is not available against BU.ObjectiveThis study aimed to investigate the vaccine potential of virulence proteins of M. ulcerans computationally.MethodsChromosome encoded virulence proteins of Mycobacterium ulcerans strain Agy99 were selected, which were available at the VFDB database. These proteins were analyzed for their subcellular localization, antigenicity, and human non-homology analysis. Ten virulence factors were finally chosen and analyzed for further study. Three-dimensional structures for selected proteins were predicted using Phyre2. B cell and T cell epitope analysis was done using methods available at Immune Epitope Database and Analysis Resource. Antigenicity, allergenicity, and toxicity analysis were also done to predict epitopes. Molecular docking analysis was done for T cell epitopes, those showing overlap with B cell epitopes.ResultsSelected virulence proteins were predicted with B cell and T cell epitopes. Some of the selected proteins were found to be already reported as antigenic in other mycobacteria. Some of the predicted epitopes also had similarities with experimentally identified epitopes of M. ulcerans and M. tuberculosis which further supported our predictions.ConclusionIn-silico approach used for the vaccine candidate identification predicted some virulence proteins that could be proved important in future vaccination strategies against this chronic disease. Predicted epitopes require further experimental validation for their potential use as peptide vaccines.     

Identification of Potential Pleiotropic Genes for Immune and Skeletal Diseases Using Multivariate MetaCCA Analysis

BackgroundImmune and skeletal systems physiologically and pathologically interact with each other. Immune and skeletal diseases may share potential pleiotropic genetics factors, but the shared specific genes are largely unknown.ObjectiveThis study aimed to investigate the overlapping genetic factors between multiple diseases (including rheumatoid arthritis (RA), psoriasis, osteoporosis, osteoarthritis, sarcopenia, and fracture).MethodsThe canonical correlation analysis (metaCCA) approach was used to identify the shared genes for six diseases by integrating genome-wide association study (GWAS)-derived summary statistics. The versatile Gene-based Association Study (VEGAS2) method was further applied to refine and validate the putative pleiotropic genes identified by metaCCA.ResultsAbout 157 (p<8.19E-6), 319 (p<3.90E-6), and 77 (p<9.72E-6) potential pleiotropic genes were identified shared by two immune diseases, four skeletal diseases, and all of the six diseases, respectively. The top three significant putative pleiotropic genes shared by both immune and skeletal diseases, including HLA-B, TSBP1, and TSBP1-AS1 (p<E-300), were located in the major histocompatibility complex (MHC) region. Nineteen of 77 putative pleiotropic genes identified by metaCCA analysis were associated with at least one disease in the VEGAS2 analysis. Specifically, the majority (18) of these 19 putative validated pleiotropic genes were associated with RA.ConclusionThe metaCCA method identified some pleiotropic genes shared by the immune and skeletal diseases. These findings help to improve our understanding of the shared genetic mechanisms and signaling pathways underlying immune and skeletal diseases.     

The genome sequence of Caenorhabditis briggsae: a platform for comparative genomics

The soil nematodes Caenorhabditis briggsae and Caenorhabditis elegans diverged from a common ancestor roughly 100 million years ago and yet are almost indistinguishable by eye. They have the same chromosome number and genome sizes, and they occupy the same ecological niche. To explore the basis for this striking conservation of structure and function, we have sequenced the C. briggsae genome to a high-quality draft stage and compared it to the finished C. elegans sequence. We predict approximately 19,500 protein-coding genes in the C. briggsae genome, roughly the same as in C. elegans. Of these, 12,200 have clear C. elegans orthologs, a further 6,500 have one or more clearly detectable C. elegans homologs, and approximately 800 C. briggsae genes have no detectable matches in C. elegans. Almost all of the noncoding RNAs (ncRNAs) known are shared between the two species. The two genomes exhibit extensive colinearity, and the rate of divergence appears to be higher in the chromosomal arms than in the centers. Operons, a distinctive feature of C. elegans, are highly conserved in C. briggsae, with the arrangement of genes being preserved in 96% of cases. The difference in size between the C. briggsae (estimated at approximately 104 Mbp) and C. elegans (100.3 Mbp) genomes is almost entirely due to repetitive sequence, which accounts for 22.4% of the C. briggsae genome in contrast to 16.5% of the C. elegans genome. Few, if any, repeat families are shared, suggesting that most were acquired after the two species diverged or are undergoing rapid evolution. Coclustering the C. elegans and C. briggsae proteins reveals 2,169 protein families of two or more members. Most of these are shared between the two species, but some appear to be expanding or contracting, and there seem to be as many as several hundred novel C. briggsae gene families. The C. briggsae draft sequence will greatly improve the annotation of the C. elegans genome. Based on similarity to C. briggsae, we found strong evidence for 1,300 new C. elegans genes. In addition, comparisons of the two genomes will help to understand the evolutionary forces that mold nematode genomes.     

Major structural differences and novel potential virulence mechanisms from the genomes of multiple campylobacter species

Sequencing and comparative genome analysis of four strains of Campylobacter including C. lari RM2100, C. upsaliensis RM3195, and C. coli RM2228 has revealed major structural differences that are associated with the insertion of phage- and plasmid-like genomic islands, as well as major variations in the lipooligosaccharide complex. Poly G tracts are longer, are greater in number, and show greater variability in C. upsaliensis than in the other species. Many genes involved in host colonization, including racR/S, cadF, cdt, ciaB, and flagellin genes, are conserved across the species, but variations that appear to be species specific are evident for a lipooligosaccharide locus, a capsular (extracellular) polysaccharide locus, and a novel Campylobacter putative licABCD virulence locus. The strains also vary in their metabolic profiles, as well as their resistance profiles to a range of antibiotics. It is evident that the newly identified hypothetical and conserved hypothetical proteins, as well as uncharacterized two-component regulatory systems and membrane proteins, may hold additional significant information on the major differences in virulence among the species, as well as the specificity of the strains for particular hosts.     

Genome sequence of Rickettsia bellii illuminates the role of amoebae in gene exchanges between intracellular pathogens

The recently sequenced Rickettsia felis genome revealed an unexpected plasmid carrying several genes usually associated with DNA transfer, suggesting that ancestral rickettsiae might have been endowed with a conjugation apparatus. Here we present the genome sequence of Rickettsia bellii, the earliest diverging species of known rickettsiae. The 1,552,076 base pair–long chromosome does not exhibit the colinearity observed between other rickettsia genomes, and encodes a complete set of putative conjugal DNA transfer genes most similar to homologues found in Protochlamydia amoebophila UWE25, an obligate symbiont of amoebae. The genome exhibits many other genes highly similar to homologues in intracellular bacteria of amoebae. We sought and observed sex pili-like cell surface appendages for R. bellii. We also found that R. bellii very efficiently multiplies in the nucleus of eukaryotic cells and survives in the phagocytic amoeba, Acanthamoeba polyphaga. These results suggest that amoeba-like ancestral protozoa could have served as a genetic “melting pot” where the ancestors of rickettsiae and other bacteria promiscuously exchanged genes, eventually leading to their adaptation to the intracellular lifestyle within eukaryotic cells.     

Importing statistical measures into Artemis enhances gene identification in the <Emphasis Type="Italic">Leishmania</Emphasis> genome project

Background

Seattle Biomedical Research Institute (SBRI) as part of the Leishmania Genome Network (LGN) is sequencing chromosomes of the trypanosomatid protozoan species Leishmania major. At SBRI, chromosomal sequence is annotated using a combination of trained and untrained non-consensus gene-prediction algorithms with ARTEMIS, an annotation platform with rich and user-friendly interfaces.

Results

Here we describe a methodology used to import results from three different protein-coding gene-prediction algorithms (GLIMMER, TESTCODE and GENESCAN) into the ARTEMIS sequence viewer and annotation tool. Comparison of these methods, along with the CODON USAGE algorithm built into ARTEMIS, shows the importance of combining methods to more accurately annotate the L. major genomic sequence.

Conclusion

An improvised and powerful tool for gene prediction has been developed by importing data from widely-used algorithms into an existing annotation platform. This approach is especially fruitful in the Leishmania genome project where there is large proportion of novel genes requiring manual annotation.

     

Adaptive gene expression divergence inferred from population genomics

Detailed studies of individual genes have shown that gene expression divergence often results from adaptive evolution of regulatory sequence. Genome-wide analyses, however, have yet to unite patterns of gene expression with polymorphism and divergence to infer population genetic mechanisms underlying expression evolution. Here, we combined genomic expression data—analyzed in a phylogenetic context—with whole genome light-shotgun sequence data from six Drosophila simulans lines and reference sequences from D. melanogaster and D. yakuba. These data allowed us to use molecular population genetics to test for neutral versus adaptive gene expression divergence on a genomic scale. We identified recent and recurrent adaptive evolution along the D. simulans lineage by contrasting sequence polymorphism within D. simulans to divergence from D. melanogaster and D. yakuba. Genes that evolved higher levels of expression in D. simulans have experienced adaptive evolution of the associated 3′ flanking and amino acid sequence. Concomitantly, these genes are also decelerating in their rates of protein evolution, which is in agreement with the finding that highly expressed genes evolve slowly. Interestingly, adaptive evolution in 5′ cis-regulatory regions did not correspond strongly with expression evolution. Our results provide a genomic view of the intimate link between selection acting on a phenotype and associated genic evolution.     

BEACON: automated tool for Bacterial GEnome Annotation ComparisON

Background

Genome annotation is one way of summarizing the existing knowledge about genomic characteristics of an organism. There has been an increased interest during the last several decades in computer-based structural and functional genome annotation. Many methods for this purpose have been developed for eukaryotes and prokaryotes. Our study focuses on comparison of functional annotations of prokaryotic genomes. To the best of our knowledge there is no fully automated system for detailed comparison of functional genome annotations generated by different annotation methods (AMs).

Results

The presence of many AMs and development of new ones introduce needs to: a/ compare different annotations for a single genome, and b/ generate annotation by combining individual ones. To address these issues we developed an Automated Tool for Bacterial GEnome Annotation ComparisON (BEACON) that benefits both AM developers and annotation analysers. BEACON provides detailed comparison of gene function annotations of prokaryotic genomes obtained by different AMs and generates extended annotations through combination of individual ones. For the illustration of BEACON’s utility, we provide a comparison analysis of multiple different annotations generated for four genomes and show on these examples that the extended annotation can increase the number of genes annotated by putative functions up to 27 %, while the number of genes without any function assignment is reduced.

Conclusions

We developed BEACON, a fast tool for an automated and a systematic comparison of different annotations of single genomes. The extended annotation assigns putative functions to many genes with unknown functions. BEACON is available under GNU General Public License version 3.0 and is accessible at: http://www.cbrc.kaust.edu.sa/BEACON/.

Electronic supplementary material

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

The complete genome sequence and comparative genome analysis of the high pathogenicity Yersinia enterocolitica strain 8081

The human enteropathogen, Yersinia enterocolitica, is a significant link in the range of Yersinia pathologies extending from mild gastroenteritis to bubonic plague. Comparison at the genomic level is a key step in our understanding of the genetic basis for this pathogenicity spectrum. Here we report the genome of Y. enterocolitica strain 8081 (serotype 0:8; biotype 1B) and extensive microarray data relating to the genetic diversity of the Y. enterocolitica species. Our analysis reveals that the genome of Y. enterocolitica strain 8081 is a patchwork of horizontally acquired genetic loci, including a plasticity zone of 199 kb containing an extraordinarily high density of virulence genes. Microarray analysis has provided insights into species-specific Y. enterocolitica gene functions and the intraspecies differences between the high, low, and nonpathogenic Y. enterocolitica biotypes. Through comparative genome sequence analysis we provide new information on the evolution of the Yersinia. We identify numerous loci that represent ancestral clusters of genes potentially important in enteric survival and pathogenesis, which have been lost or are in the process of being lost, in the other sequenced Yersinia lineages. Our analysis also highlights large metabolic operons in Y. enterocolitica that are absent in the related enteropathogen, Yersinia pseudotuberculosis, indicating major differences in niche and nutrients used within the mammalian gut. These include clusters directing, the production of hydrogenases, tetrathionate respiration, cobalamin synthesis, and propanediol utilisation. Along with ancestral gene clusters, the genome of Y. enterocolitica has revealed species-specific and enteropathogen-specific loci. This has provided important insights into the pathology of this bacterium and, more broadly, into the evolution of the genus. Moreover, wider investigations looking at the patterns of gene loss and gain in the Yersinia have highlighted common themes in the genome evolution of other human enteropathogens.     

Draft Genome Sequence of Xylaria grammica EL000614, a Strain Producing Grammicin,a Potent Nematicidal Compound

An endolichenic fungus, Xylaria grammica strain {"type":"entrez-nucleotide","attrs":{"text":"EL000614","term_id":"124922317"}}EL000614, showed strong nematicidal effects against plant pathogenic nematode, Meloidogyne incognita by producing grammicin. We report genome assembly of X. grammica {"type":"entrez-nucleotide","attrs":{"text":"EL000614","term_id":"124922317"}}EL000614 comprised of 25 scaffolds with a total length of 54.73 Mb, N50 of 4.60 Mb, and 99.8% of BUSCO completeness. GC contents of this genome were 44.02%. Gene families associated with biosynthesis of secondary metabolites or regulatory proteins were identified out of 13,730 gene models predicted.     

Africanization of a feral honey bee (Apis mellifera) population in South Texas: does a decade make a difference?

The arrival to the United States of the Africanized honey bee, a hybrid between European subspecies and the African subspecies Apis mellifera scutellata, is a remarkable model for the study of biological invasions. This immigration has created an opportunity to study the dynamics of secondary contact of honey bee subspecies from African and European lineages in a feral population in South Texas. An 11‐year survey of this population (1991–2001) showed that mitochondrial haplotype frequencies changed drastically over time from a resident population of eastern and western European maternal ancestry, to a population dominated by the African haplotype. A subsequent study of the nuclear genome showed that the Africanization process included bidirectional gene flow between European and Africanized honey bees, giving rise to a new panmictic mixture of A. m. scutellata‐ and European‐derived genes. In this study, we examined gene flow patterns in the same population 23 years after the first hybridization event occurred. We found 28 active colonies inhabiting 92 tree cavities surveyed in a 5.14 km² area, resulting in a colony density of 5.4 colonies/km². Of these 28 colonies, 25 were of A. m. scutellata maternal ancestry, and three were of western European maternal ancestry. No colonies of eastern European maternal ancestry were detected, although they were present in the earlier samples. Nuclear DNA revealed little change in the introgression of A. m. scutellata‐derived genes into the population compared to previous surveys. Our results suggest this feral population remains an admixed swarm with continued low levels of European ancestry and a greater presence of African‐derived mitochondrial genetic composition.