Tag: HMM

The little skate genome and the evolutionary emergence of wing-like fins

Nakamura, T. et al. Molecular mechanisms underlying the exceptional adaptations of batoid fins. Proc. Natl Acad. Sci. USA 112, 15940–15945 (2015). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  Turner, N. et al. The evolutionary origins and diversity of the neuromuscular system of paired appendages in batoids. Proc. Biol. Sci….

Continue Reading The little skate genome and the evolutionary emergence of wing-like fins

how to make a header name in a haplotyping script of gatk?

how to make a header name in a haplotyping script of gatk? 1 Hi, I want to ask how we can make the header name as per our choice in a haplotyping script of gatk because by default the header name of the output.vcf file is mentioned as sample1? here…

Continue Reading how to make a header name in a haplotyping script of gatk?

An ancient metalloenzyme evolves through metal preference modulation

Jayaraman, V., Toledo-Patino, S., Noda-Garcia, L. & Laurino, P. Mechanisms of protein evolution. Protein Sci. 31, e4362 (2022). Article  CAS  PubMed  Google Scholar  Huang, R. et al. Enzyme functional evolution through improved catalysis of ancestrally nonpreferred substrates. Proc. Natl Acad. Sci. USA 109, 2966–2971 (2012). Article  CAS  PubMed  PubMed Central …

Continue Reading An ancient metalloenzyme evolves through metal preference modulation

CRISPR-resolved virus-host interactions in a municipal landfill include non-specific viruses, hyper-targeted viral populations, and interviral conflicts

Suttle, C. A. Environmental microbiology: Viral diversity on the global stage. Nat. Microbiol. 1, 16205 (2016). Article  CAS  PubMed  Google Scholar  Schulz, F. et al. Giant virus diversity and host interactions through global metagenomics. Nature 578, 432–436 (2020). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  Dutilh, B. E., Reyes,…

Continue Reading CRISPR-resolved virus-host interactions in a municipal landfill include non-specific viruses, hyper-targeted viral populations, and interviral conflicts

Can’t call subsampled bam file with GATK Haplotypecaller with –disable-tool-default-read-filters

I want to simulate variant calling of an ultra-low-coverage >0.005x bam file. I subsampled reads from the (HG02024) sample of the 1KG phase 3 dataset. My code in R to do so is the following (bam and reference are just path extensions, file is the inital bam file): cov_rate <-…

Continue Reading Can’t call subsampled bam file with GATK Haplotypecaller with –disable-tool-default-read-filters

Unexpected genetic and microbial diversity for arsenic cycling in deep sea cold seep sediments

Joye, S. B. The geology and biogeochemistry of hydrocarbon seeps. Annu. Rev. Earth Planet. Sci. 48, 205–231 (2020). Article  CAS  Google Scholar  Feng, D. et al. Cold seep systems in the South China Sea: An overview. J. Asian Earth SCI 168, 3–16 (2018). Article  Google Scholar  Dubilier, N., Bergin, C….

Continue Reading Unexpected genetic and microbial diversity for arsenic cycling in deep sea cold seep sediments

Dealing with “Too many samples were discarded” – StrainPhlAn

I’m attempting to do a strainphlan analysis of Blautia wexlerae, but nothing I’m doing is working. I think the problem might be the reference genomes that I’m using, but I’m using the 2 primary genomes from NCBI. Is there another place to get correctly formatted references? Error Here’s the primary…

Continue Reading Dealing with “Too many samples were discarded” – StrainPhlAn

Hyperactive nanobacteria with host-dependent traits pervade Omnitrophota

Pruesse, E. et al. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 35, 7188–7196 (2007). Article  CAS  PubMed  PubMed Central  Google Scholar  Glöckner, J. et al. Phylogenetic diversity and metagenomics of candidate division OP3. Environ. Microbiol. 12, 1218–1229…

Continue Reading Hyperactive nanobacteria with host-dependent traits pervade Omnitrophota

How to extract phased haplotypes from GATK HaplotypeCaller

I would like to extract the physically phased haplotypes from a VCF file generated by GATK’s HaplotypeCaller on Illumina data of some isolates from different yeast (S. cerevisiae) strains. According to this FAQ: In the format field of a PGT (Pre-Implantation Genetic Testing) VCF, you may find a description similar…

Continue Reading How to extract phased haplotypes from GATK HaplotypeCaller

how to iteratively search sequences against sequence database using profile HMM

how to iteratively search sequences against sequence database using profile HMM 1 Hi all, I wonder whether there is a tool for iterative sequence searching using a profile HMM, the task as like the domain construction in Pfam database. jackhmmer can iteratively search sequence against a sequence database, but it…

Continue Reading how to iteratively search sequences against sequence database using profile HMM

Transcription of MERVL retrotransposons is required for preimplantation embryo development

MERVL exhibits distinct localization in mouse embryos To understand the dynamics of MERVL expression, we first analyzed publicly available single-cell RNA-sequencing (scRNA-seq) datasets from each blastomere at eight representative stages of preimplantation development18 (Fig. 1a). To define regions of nonredundant MERVLs in mouse genome, we used RepeatMasker to annotate the…

Continue Reading Transcription of MERVL retrotransposons is required for preimplantation embryo development

Importance of mobile genetic element immunity in numerically abundant Trichodesmium clades

Moore CM, Mills MM, Arrigo KR, Berman-Frank I, Bopp L, Boyd PW, et al. Processes and patterns of oceanic nutrient limitation. Nat Geosci. 2013;6:701–10. Article  CAS  Google Scholar  Sohm JA, Webb EA, Capone DG. Emerging patterns of marine nitrogen fixation. Nat Rev Microbiol. 2011;9:499–508. Article  CAS  PubMed  Google Scholar  Zehr…

Continue Reading Importance of mobile genetic element immunity in numerically abundant Trichodesmium clades

A novel and diverse group of Candidatus Patescibacteria from bathypelagic Lake Baikal revealed through long-read metagenomics | Environmental Microbiome

Rinke C, Schwientek P, Sczyrba A, Ivanova NN, Anderson IJ, Cheng J-F, et al. Insights into the phylogeny and coding potential of microbial dark matter nature. Nat Publ Group. 2013;499:431–7. CAS  Google Scholar  Brown CT, Hug LA, Thomas BC, Sharon I, Castelle CJ, Singh A, et al. Unusual biology across…

Continue Reading A novel and diverse group of Candidatus Patescibacteria from bathypelagic Lake Baikal revealed through long-read metagenomics | Environmental Microbiome

Metagenomic and machine learning-aided identification of biomarkers driving distinctive Cd accumulation features in the root-associated microbiome of two rice cultivars

Zhang J, Liu Y, Zhang N, Hu B, Jin T, Xu H, et al. NRT1.1B is associated with root microbiota composition and nitrogen use in field-grown rice. Nat Biotechnol. 2019;37:676–84. Article  CAS  PubMed  Google Scholar  Philippot L, Raaijmakers JM, Lemanceau P, van der Putten WH. Going back to the roots:…

Continue Reading Metagenomic and machine learning-aided identification of biomarkers driving distinctive Cd accumulation features in the root-associated microbiome of two rice cultivars

GATK showing error of reference file index

GATK showing error of reference file index 0 Hi, i am using the GATK version gatk/4.2.2.0 for HaplotypeCaller. I have been facing the reference.fa indexing issue. I tried to index the file using the following command samtools faidx ~/path/PitayaGenomic.fa #three difference formats of reference files -rw-r–r– 1 tariqr 1.3G Feb…

Continue Reading GATK showing error of reference file index

Antioxidant enzymes that target hydrogen peroxide are conserved across the animal kingdom, from sponges to mammals

Our metazoan-wide survey has provided the most comprehensive analysis to date of gene number and phylogenetic distribution of three key antioxidant gene families across the animal kingdom. Genes encoding all three families were observed in 18 metazoan species; the exception is the ctenophore Mnemiopsis leidyi that has PRX and GPX,…

Continue Reading Antioxidant enzymes that target hydrogen peroxide are conserved across the animal kingdom, from sponges to mammals

Genetic mapping of microbial and host traits reveals production of immunomodulatory lipids by Akkermansia muciniphila in the murine gut

Animal studies Animal care and study protocols were approved by the AAALAC-accredited Institutional Animal Care and Use Committee of the College of Agricultural Life Sciences at the University of Wisconsin-Madison (UW-Madison). All experiments with mice were performed under protocols approved by the UW-Madison Animal Care and Use Committee (Protocol number…

Continue Reading Genetic mapping of microbial and host traits reveals production of immunomodulatory lipids by Akkermansia muciniphila in the murine gut

How To Install clustalw on Ubuntu 20.04

In this tutorial we learn how to install clustalw on Ubuntu 20.04. clustalw is global multiple nucleotide or peptide sequence alignment 633246bd8fd1b951f15985f7cbfb1909 Introduction In this tutorial we learn how to install clustalw on Ubuntu 20.04. What is clustalw clustalw is: This program performs an alignment of multiple nucleotide or amino…

Continue Reading How To Install clustalw on Ubuntu 20.04

AlphaFold2 reveals commonalities and novelties in protein structure space for 21 model organisms

Proportion of AlphaFold2 models that can be brought into CATH Identification of domains in AF2 protein models We applied an in-house Hidden-Markov Model-based protocol CATH-Resolve-Hits (CRH)26 to assign domain regions in all sequences from the 21 model organisms modelled by AlphaFold2 (AF2) (see “Methods” for description and Fig. 1). CRH identifies…

Continue Reading AlphaFold2 reveals commonalities and novelties in protein structure space for 21 model organisms

Issue with hmmcalibrate during tutorial.

Issue with hmmcalibrate during tutorial. 1 Hi everyone. I’m trying to do the tutorial for hmmer but I seem to be having an issue for hmmcalibrate. I tried to use: hmmcalibrate globin.hmm But it says: Command ‘hmmcalibrate’ not found, did you mean: command ‘hmm2calibrate’ from deb hmmer2 (2.3.2+dfsg-6) Try: sudo…

Continue Reading Issue with hmmcalibrate during tutorial.

Genomics discovery of giant fungal viruses from subsurface oceanic crustal fluids

Orcutt B, D’Angelo T, Jungbluth SP, Huber JA, Sylvan JB. Microbial life in oceanic crust. OSF Preprints, 2020; doi.org/10.31219/osf.io/2wxe6. Koonin EV. On the origin of cells and viruses: primordial virus world scenario. Ann NY Acad Sci. 2009;1178:47–64. Nigro OD, Jungbluth SP, Lin HT, Hsieh CC, Miranda JA, Schvarcz CR, et…

Continue Reading Genomics discovery of giant fungal viruses from subsurface oceanic crustal fluids

Phylogenetic and AlphaFold predicted structure analyses provide insights for A1 aspartic protease family classification in Arabidopsis

Introduction Proteases regulate various biological processes including protein synthesis and maturation, activity modification, degradation and turnover. Depending on their catalytic mechanisms, these proteases are primarily classified into cysteine, metallo-, serine, threonine and aspartic protease family (Beers et al., 2004). The latter protease family is known as acid protease family because they…

Continue Reading Phylogenetic and AlphaFold predicted structure analyses provide insights for A1 aspartic protease family classification in Arabidopsis

Faking hh-suite workflow / alignment output

Faking hh-suite workflow / alignment output 1 Given the following: query.fasta -> single entry reference.fasta -> multiple entries I now want to ‘fake’ (or just be able to get) an output that looks like a proper *.hhr alignment file, i.e. as if I aligned the query.fasta against the profiles of…

Continue Reading Faking hh-suite workflow / alignment output

A paralog of Pcc1 is the fifth core subunit of the KEOPS tRNA-modifying complex in Archaea

A paralog of Pcc1 is widely distributed among archaea The recent discovery of a fifth subunit of KEOPS in metazoa raised the possibility that a corresponding subunit could exist in archaea. In search for such a subunit, we noticed in the model archaeon Thermococcus kodakarensis KOD1 the existence of two…

Continue Reading A paralog of Pcc1 is the fifth core subunit of the KEOPS tRNA-modifying complex in Archaea

Issue with VCF format while using Pharmcat

Hello everybody, I am using pharmcat tool’s prerprocessor feature to preprocessmy vcf file using the command > python3 pharmcat_vcf_preprocessor.py -vcf sample.vcf But I think there is some issue with my vcf file as this command outputs an error > Reading samples from sample.vcf … Saving output to . > >…

Continue Reading Issue with VCF format while using Pharmcat

Isolation and infection cycle of a polinton-like virus virophage in an abundant marine alga

Koonin, E. V. & Dolja, V. V. Virus world as an evolutionary network of viruses and capsidless selfish elements. Microbiol. Mol. Biol. Rev. 78, 278–303 (2014). Article  CAS  Google Scholar  Pritham, E. J., Putliwala, T. & Feschotte, C. Mavericks, a novel class of giant transposable elements widespread in eukaryotes and…

Continue Reading Isolation and infection cycle of a polinton-like virus virophage in an abundant marine alga

Annelid functional genomics reveal the origins of bilaterian life cycles

Hall, B. K. & Wake, M. H. in The Origin and Evolution of Larval Forms (eds Hall, B. K. & Wake, M. H.) 1–19 (Academic Press, 1999). Nielsen, C. Animal phylogeny in the light of the trochaea theory. Biol. J. Linn. Soc. 25, 243–299 (2008). Article  Google Scholar  Garstang, W….

Continue Reading Annelid functional genomics reveal the origins of bilaterian life cycles

TCR sequence analysis

TCR sequence analysis 2 Hi, I’m doing single-cell TCR sequencing using Mark Davis’ protocol. Here is the link to Davis paper: www.nature.com/nbt/journal/v32/n7/abs/nbt.2938.html The paper used vdjfasta for TCR analysis. But the vdjfasta was set up for antibody sequence analysis. When I use it for TCR sequences, I don’t get exactly…

Continue Reading TCR sequence analysis

CusProSe: a customizable protein annotation software with an application to the prediction of fungal secondary metabolism genes

Development of the CusProSe software CustomProteinSearch (CusProSe) is a generic genome mining software, consisting of two distinct but complementary customizable programs: IterHMMBuild and ProSeCDA. IterHMMBuild is an HMM profile building tool based on an iterative learning process. ProSeCDA is a protein search and annotation tool based on user-defined domain architectures….

Continue Reading CusProSe: a customizable protein annotation software with an application to the prediction of fungal secondary metabolism genes

Finding common genes by taxa on Genbank? : bioinformatics

Hmm, I am assuming that since you are referring to canidae then you want “family level” information. I know you can do this bioinformaticly by getting a large list of different species that are also from different genus (one taxa from each known species /genus). Add in a few taxa…

Continue Reading Finding common genes by taxa on Genbank? : bioinformatics

Scatter Gather principle by chromosome on Gatk

Scatter Gather principle by chromosome on Gatk 0 Hi all, On a quest to optimize gatk pipeline, I met scatter gather principle, so I did following, pids= for chr in chr1 chr2 chr3 chr4 chr5 chr6 chr7 chr8 chr9 chr10 chr11 chr12 chr13 chr14 chr15 chr16 chr17 chr18 chr19 chr20…

Continue Reading Scatter Gather principle by chromosome on Gatk

Unexpected absence of ribosomal protein genes from metagenome-assembled genomes

Hug LA, Baker BJ, Anantharaman K, Brown CT, Probst AJ, Castelle CJ, et al. A new view of the tree of life. Nat Microbiol. 2016;1:16048. Article  PubMed  Google Scholar  Castelle CJ, Wrighton KC, Thomas BC, Hug LA, Brown CT, Wilkins MJ, et al. Genomic expansion of domain archaea highlights roles…

Continue Reading Unexpected absence of ribosomal protein genes from metagenome-assembled genomes

Ancient origin and constrained evolution of the division and cell wall gene cluster in Bacteria

Miyakawa, T., Matsuzawa, H., Matsuhashi, M. & Sugino, Y. Cell wall peptidoglycan mutants of Escherichia coli K-12: existence of two clusters of genes, mra and mrb, for cell wall peptidoglycan biosynthesis. J. Bacteriol. 112, 950–958 (1972). Article  CAS  PubMed  PubMed Central  Google Scholar  Ayala, J. A., Garrido, T., De Pedro,…

Continue Reading Ancient origin and constrained evolution of the division and cell wall gene cluster in Bacteria

Visualization of Pool-HMM results

Visualization of Pool-HMM results 0 Hello everyone, I have run the Pool-HMM tool on my pool-seq data and got final results in this format Chr Start End Pool-HMM 1 312212 413291 11.60071199 1 595314 676207 11.40558651 1 3136358 3180593 11.65338609 1 3788166 3813896 4.561544648 1 5202297 5218014 1.620942559 1 5448808…

Continue Reading Visualization of Pool-HMM results

The genus Serratia revisited by genomics

Merlino, C. P. Bartolomeo Bizio’s letter to the most eminent priest, Angelo Bellani, concerning the phenomenon of the red-colored polenta [translated from the Italian]. J. Bacteriol. 9, 527–543 (1924). Grimont, P. A. D. & Dulong de Rosnay, H. L. C. Numerical study of 60 strains of Serratia. J. Gen. Microbiol….

Continue Reading The genus Serratia revisited by genomics

The evolutionary origin of host association in the Rickettsiales

Salje, J. Cells within cells: Rickettsiales and the obligate intracellular bacterial lifestyle. Nat. Rev. Microbiol. 19, 375–390 (2021). CAS  PubMed  Article  Google Scholar  Wang, S. & Luo, H. Dating Alphaproteobacteria evolution with eukaryotic fossils. Nat. Commun. 12, 3324 (2021). CAS  PubMed  PubMed Central  Article  Google Scholar  Strassert, J. F. H.,…

Continue Reading The evolutionary origin of host association in the Rickettsiales

Allelic expression imbalance of PIK3CA mutations is frequent in breast cancer and prognostically significant

Subjects Normal breast and tumor samples were obtained with the written informed consent from donors and appropriate approval from local ethical committees, with the detailed information described in the respective original publications: normal tissue9, METABRIC14, TCGA35. Differential allelic expression analysis DNA and total RNA from 64 samples of normal breast…

Continue Reading Allelic expression imbalance of PIK3CA mutations is frequent in breast cancer and prognostically significant

Top gru open source projects

Pytorch Seq2seq Tutorials on implementing a few sequence-to-sequence (seq2seq) models with PyTorch and TorchText. Rnn ctc Recurrent Neural Network and Long Short Term Memory (LSTM) with Connectionist Temporal Classification implemented in Theano. Includes a Toy training example. Haste Haste: a fast, simple, and open RNN library Eeg Dl A Deep…

Continue Reading Top gru open source projects

Pangenome-based genome inference allows efficient and accurate genotyping across a wide spectrum of variant classes

Sequencing data We used publicly available sequencing data from the GIAB consortium45, 1000 Genomes Project high-coverage data46 and Human Genome Structural Variation Consortium (HGSVC)4. All datasets include only samples consented for public dissemination of the full genomes. Statistics and reproducibility For generating the assemblies, we used all 14 samples for…

Continue Reading Pangenome-based genome inference allows efficient and accurate genotyping across a wide spectrum of variant classes

HSP object suitable for describing WABA alignments

Bio::Search::HSP::WABAHSP(3) HSP object suitable for describing WABA alignments SYNOPSIS # use this object as you would a GenericHSP # a few other methods have been added including state DESCRIPTION This object implements a few of the extra methods such as hmmstate_string which returns the HMM state representation for the WABA…

Continue Reading HSP object suitable for describing WABA alignments

What is the best reply to hmm?

Hmm 0 views I like this I dislike this Related questions What do u mean by hug? Why is hidden Markov used in speech recognition? What is profile in bioinformatics? What is HMM in bioinformatics? What is hidden Markov in speech recognition? What is hmm called? What Mmmm means? What…

Continue Reading What is the best reply to hmm?

NCBI HMM accession NF040677

NCBI HMM accession Source identifier Product name Label Gene symbol Family type EC number(s) GO term(s) HMM length Sequence cutoff Domain cutoff Number of RefSeq protein hits counting … HMM profile HMM seed Named by this model (…) Other hits (…) No RefSeq protein is named by this evidence. Failed…

Continue Reading NCBI HMM accession NF040677

GATK HaplotypeCaller with interval list

I am trying to use the -L option of GATK HaplotypeCaller to call SNPs and short InDels with in an interval list. My interval list file (top8snp.interval_list) content is as follows: 12 33029845 33030845 + rs24767598 13 40586682 40587682 + rs24748362 18 24373857 24374857 + rs8856159 21 50381146 50382146 +…

Continue Reading GATK HaplotypeCaller with interval list

Petabase-scale sequence alignment catalyses viral discovery

Serratus alignment architecture Serratus (v0.3.0) (github.com/ababaian/serratus) is an open-source cloud-infrastructure designed for ultra-high-throughput sequence alignment against a query sequence or pangenome (Extended Data Fig. 1). Serratus compute costs are dependent on search parameters (expanded discussion available: github.com/ababaian/serratus/wiki/pangenome_design). The nucleotide vertebrate viral pangenome search (bowtie2, database size: 79.8 MB) reached processing rates…

Continue Reading Petabase-scale sequence alignment catalyses viral discovery

Time-course RNASeq of Camponotus floridanus forager and nurse ant brains indicate links between plasticity in the biological clock and behavioral division of labor | BMC Genomics

1. Sharma VK. Adaptive significance of circadian clocks. Chronobiol Int. 2003;20(6):901–19. PubMed  Google Scholar  2. Paranjpe DA, Sharma VK. Evolution of temporal order in living organisms. J Circadian Rhythms. 2005;3(1):7. PubMed  PubMed Central  Google Scholar  3. Yerushalmi S, Green RM. Evidence for the adaptive significance of circadian rhythms. Ecol Lett….

Continue Reading Time-course RNASeq of Camponotus floridanus forager and nurse ant brains indicate links between plasticity in the biological clock and behavioral division of labor | BMC Genomics

Integrating cultivation and metagenomics for a multi-kingdom view of skin microbiome diversity and functions

1. Oh, J. et al. Biogeography and individuality shape function in the human skin metagenome. Nature 514, 59–64 (2014). 2. Byrd, A. L., Belkaid, Y. & Segre, J. A. The human skin microbiome. Nat. Rev. Microbiol. 16, 143–155 (2018). CAS  PubMed  Google Scholar  3. Oh, J. et al. Temporal stability…

Continue Reading Integrating cultivation and metagenomics for a multi-kingdom view of skin microbiome diversity and functions

Towards the biogeography of prokaryotic genes

1. Sunagawa, S. et al. Structure and function of the global ocean microbiome. Science 348, 1261359 (2015). PubMed  Google Scholar  2. Zou, Y. et al. 1,520 reference genomes from cultivated human gut bacteria enable functional microbiome analyses. Nat. Biotechnol. 37, 179–185 (2019). CAS  PubMed  PubMed Central  Google Scholar  3. Mohammad,…

Continue Reading Towards the biogeography of prokaryotic genes

alphafold2: HHblits failed – githubmemory

I’ve tried using the standard alphafold2 setup via docker (converted to a singularity container) via the setup described at github.com/kalininalab/alphafold_non_docker, and both result in the following error: […] E1210 12:01:01.009660 22603932526400 hhblits.py:141] – 11:49:18.512 INFO: Iteration 1 E1210 12:01:01.009703 22603932526400 hhblits.py:141] – 11:49:19.070 INFO: Prefiltering database E1210 12:01:01.009746 22603932526400 hhblits.py:141]…

Continue Reading alphafold2: HHblits failed – githubmemory

pfam_scan.pl can’t find the pfamdb

pfam_scan.pl can’t find the pfamdb 1 I am trying to run pfam_scan.pl script which keep generating this error below though both Pfam-A.hmm and Pfam-A.hmm.dat files are in /pfamdb. Can someone please help me identify the errors and resolve this? perl /media/owner/b45f8e7a-003c-4573-8841-bcb5f76f281f/sn/rgaugury/PfamScan/pfam_scan.pl -fasta Hannuus_494_r1.2.protein.fa -dir /media/owner/b45f8e7a-003c-4573-8841-bcb5f76f281f/sn/rgaugury/database/pfamdb FATAL: can’t find “Pfam-A.hmm” and/or…

Continue Reading pfam_scan.pl can’t find the pfamdb

NC 002528 | Virtual Laboratory Wiki

LOCUS NC_002528 640681 bp DNA circular BCT 07-JUN-2007 DEFINITION Buchnera aphidicola str. APS (Acyrthosiphon pisum), complete genome. ACCESSION NC_002528 VERSION NC_002528.1 GI:15616630 KEYWORDS . SOURCE Buchnera aphidicola str. APS (Acyrthosiphon pisum) ORGANISM Buchnera aphidicola str. APS (Acyrthosiphon pisum) Bacteria; Proteobacteria; Gammaproteobacteria; Enterobacteriales; Enterobacteriaceae; Buchnera. REFERENCE 1 (sites) AUTHORS Shigenobu,S., Watanabe,H.,…

Continue Reading NC 002528 | Virtual Laboratory Wiki

Mg2+-dependent conformational rearrangements of CRISPR-Cas12a R-loop complex are mandatory for complete double-stranded DNA cleavage

Significance CRISPR-Cas12a has emerged as attractive molecular scissors alternative to Cas9 owing to its unique features including fewer off-target effects, an alternative protospacer-adjacent motif sequence, pre-CRISPR RNA processing activity, and indiscriminate single-stranded DNase activity. However, despite these advantages, Cas12a has not been well utilized as recently reported base and prime…

Continue Reading Mg2+-dependent conformational rearrangements of CRISPR-Cas12a R-loop complex are mandatory for complete double-stranded DNA cleavage

Snakemake: MissingInputException

Snakemake: MissingInputException 0 Hello, I am trying to create a simple Snakemake workflow and I am having some issues. My file looks like this: ——————– ARCHIVE_FILE = ‘output.tar.gz’ **a single output file** OUTPUT_FILE = ‘output/{species}.out’ **a single input file** INPUT_FILE = ‘proteins/{species}.fasta’ **Build the list of input files.** INP =…

Continue Reading Snakemake: MissingInputException

Aspergillus fumigatus pan-genome analysis identifies genetic variants associated with human infection

1. Latgé, J. P. and Chamilos, G. Aspergillus fumigatus and Aspergillosis in 2019. Clin. Microbiol. Rev. doi.org/10.1128/CMR.00140-18 (2019). 2. Invasive Aspergillosis. LIFE www.life-worldwide.org/fungal-diseases/invasive-aspergillosis (2020). 3. Harrison, N. et al. Incidence and characteristics of invasive fungal diseases in allogeneic hematopoietic stem cell transplant recipients: a retrospective cohort study. BMC Infect. Dis….

Continue Reading Aspergillus fumigatus pan-genome analysis identifies genetic variants associated with human infection

minimac4: autopkgtest regression: *** stack smashing detected ***: terminated

Source: minimac4 Version: 1.0.2-3 X-Debbugs-CC: debian…@lists.debian.org Severity: serious User: debian…@lists.debian.org Usertags: regression Dear maintainer(s), With a recent upload of minimac4 the autopkgtest of minimac4 fails in testing when that autopkgtest is run with the binary packages of minimac4 from unstable. It passes when run with only packages from testing. In…

Continue Reading minimac4: autopkgtest regression: *** stack smashing detected ***: terminated

simonsvaerd/pytorch-struct – Giters

A library of tested, GPU implementations of core structured prediction algorithms for deep learning applications. HMM / LinearChain-CRF HSMM / SemiMarkov-CRF Dependency Tree-CRF PCFG Binary Tree-CRF … Designed to be used as efficient batched layers in other PyTorch code. Tutorial paper describing methodology. Getting Started !pip install -qU git+https://github.com/harvardnlp/pytorch-struct #…

Continue Reading simonsvaerd/pytorch-struct – Giters

Index of /~psgendb/local/biopython-1.64.old/Bio

Name Last modified Size Description Parent Directory   –   Affy/ 2014-05-29 05:25 –   Align/ 2014-06-11 10:27 –   AlignIO/ 2014-06-11 10:27 –   Alphabet/ 2014-06-11 10:27 –   Application/ 2014-05-29 05:25 –   Blast/ 2014-05-29 05:25 –   CAPS/ 2014-05-29 05:25 –   Cluster/ 2014-05-29 05:25 –  …

Continue Reading Index of /~psgendb/local/biopython-1.64.old/Bio

(PDF) Predicting MoRFs in protein sequences using HMM profiles | Shiu Kumar

(PDF) Predicting MoRFs in protein sequences using HMM profiles | Shiu Kumar – Academia.edu Academia.edu no longer supports Internet Explorer. To browse Academia.edu and the wider internet faster and more securely, please take a few seconds to upgrade your browser. Academia.edu uses cookies to personalize content, tailor ads and improve the…

Continue Reading (PDF) Predicting MoRFs in protein sequences using HMM profiles | Shiu Kumar

How to use output from GeneMark-ES to identify function?

How to use output from GeneMark-ES to identify function? 2 Hello everyone, I would like to ask anyone who know about How to use result from GeneMark-ES program to identify function? For now I already have a result like below of this post. And It’s include Nucleotide sequences output but…

Continue Reading How to use output from GeneMark-ES to identify function?

Single cell genomics reveals plastid-lacking Picozoa are close relatives of red algae

1. Strassert, J. F. H., Irisarri, I., Williams, T. A. & Burki, F. A molecular timescale for eukaryote evolution with implications for the origin of red algal-derived plastids. Nat. Commun. 12, 1879 (2021). ADS  CAS  PubMed  PubMed Central  Google Scholar  2. Burki, F., Roger, A. J., Brown, M. W. &…

Continue Reading Single cell genomics reveals plastid-lacking Picozoa are close relatives of red algae

python – Error while parsing gene bank file using Biopython

This question was migrated from Unix & Linux Stack Exchange because it can be answered on Bioinformatics Stack Exchange. Migrated 8 hours ago. I am trying to extract the protein sequence of specific genes from gene bank like format file obtained from antismash part of which looks like…

Continue Reading python – Error while parsing gene bank file using Biopython

How to extract protein sequence of selected genes from gene bank like format

Hello everyone, I used this program called antismash to predict the secondary metabolite clusters in my genome of interest. The output I got is gene bank like format whish is shown below: LOCUS scaffold_10 47160 bp DNA linear UNK 01-JAN-1980 DEFINITION scaffold_10. ACCESSION scaffold_10 VERSION scaffold_10 KEYWORDS . SOURCE ORGANISM…

Continue Reading How to extract protein sequence of selected genes from gene bank like format

Annotation tools (prokaryotes): prokka vs eggnog

Annotation tools (prokaryotes): prokka vs eggnog 1 Are there any obvious advantages/disadvantages to using one over the other? Both use HMM (I think), both are hierarchical (starting with EggNOG 4.5). Are they simply competitors? Perhaps there is a trade-off between database size (Prokka is smaller) and the quality of curation?…

Continue Reading Annotation tools (prokaryotes): prokka vs eggnog

Venatorbacter cucullus gen. nov sp. nov a novel bacterial predator

1. Pérez, J., Moraleda-Muñoz, A., Marcos-Torres, F. J. & Muñoz-Dorado, J. Bacterial predation: 75 years and counting!. Environ. Microbiol. 18, 766–779 (2016). PubMed  Article  PubMed Central  Google Scholar  2. Linares-Otoya, L. et al. Diversity and antimicrobial potential of predatory bacteria from the Peruvian coastline. Mar. Drugs. 15, E308. doi.org/10.3390/md15100308 (2017)….

Continue Reading Venatorbacter cucullus gen. nov sp. nov a novel bacterial predator

Update protein annotation on databases

Update protein annotation on databases 1 Hello all, I’m a bit new to the computational biology scene, and I would like information regarding how protein annotations are curated and how individuals can submit protein annotations. I ask because my team and I recently published a paper that shows the function…

Continue Reading Update protein annotation on databases

Bioconductor – Bioconductor 3.14 Released

Home Bioconductor 3.14 Released October 27, 2021 Bioconductors: We are pleased to announce Bioconductor 3.14, consisting of 2083 software packages, 408 experiment data packages, 904 annotation packages, 29 workflows and 8 books. There are 89 new software packages, 13 new data experiment packages, 10 new annotation packages, 1 new workflow,…

Continue Reading Bioconductor – Bioconductor 3.14 Released

How to train annotations tools

How to train annotations tools 1 I would like to train Augustus, SNAP and GlimmerHMM. I found protein sequences in GenBank and in orthodb.org. Furthermore, I found HMM files on busco-data.ezlab.org. $ wget -c busco-data.ezlab.org/v5/data/lineages/viridiplantae_odb10.2020-09-10.tar.gz $ wget -c v100.orthodb.org/download/odb10_plants_fasta.tar.gz Are there any instructions on how to train those annotations tools?…

Continue Reading How to train annotations tools

Bioconductor – ADaCGH2

    This package is for version 2.12 of Bioconductor; for the stable, up-to-date release version, see ADaCGH2. Analysis of data from aCGH experiments using parallel computing and ff objects Bioconductor version: 2.12 Analysis and plotting of array CGH data. Allows usage of Circular Binary Segementation, wavelet-based smoothing (both as…

Continue Reading Bioconductor – ADaCGH2

Sample not found in BAM header

GATK mutech2: Sample not found in BAM header 0 I get the following error when running mutech2, any idea what the reason is: A USER ERROR has occurred: Bad input: Sample N-PANCNGS-006 is not in BAM header: [] gatk Mutect2 –native-pair-hmm-threads 30 -R ~/genomes/BWA/Homo_sapiens.GRCh38.dna.primary_assembly.fa -I T-PANCNGS-006.bam -I N-PANCNGS-006.bam -normal N-PANCNGS-006…

Continue Reading Sample not found in BAM header

Runs of homozygosity in Plink

❯ plink1.9 –homozyg –help PLINK v1.90b6.22 64-bit (3 Nov 2020) www.cog-genomics.org/plink/1.9/ (C) 2005-2020 Shaun Purcell, Christopher Chang GNU General Public License v3 –help present, ignoring other flags. –homozyg [{group | group-verbose}] [‘consensus-match’] [‘extend’] [‘subtract-1-from-lengths’] –homozyg-snp <min var count> –homozyg-kb <min length> –homozyg-density <max inverse density (kb/var)> –homozyg-gap <max internal gap…

Continue Reading Runs of homozygosity in Plink

Prokka Annotation or NCBI Annotation

Prokka Annotation or NCBI Annotation 1 Dear All, I have two sets of annotation files for 15 Bacterial genomes. One set from NCBI annotations (from RefSeq) and the other from Prokka (I have run it in my local machine). Which one is advisable to use between the two for all…

Continue Reading Prokka Annotation or NCBI Annotation

How to merge HMM files for HMMSearch?

I’m trying to add several mult-HMM files together but I haven’t had much luck. First I tried to just cat them, then I realized that it won’t work if there is an overlap of HMMs within the files. I then made a script that only concatenates unique ones but once…

Continue Reading How to merge HMM files for HMMSearch?

RNAmmer running

RNAmmer running 1 Dear, I am Kishor from Shanghai. Recently I have been trying to use RNAmmer. But yet to successfully run it. I have made two changes in the rnammer according to instructions, like as follows: my $INSTALL_PATH = “/mnt/genome3/Lab_Users/Kishor/DISK_2/softwares/RNAmmer” **for linux HMMSEARCHBINARY=”/mnt/genome3/LabUsers/Kishor/DISK2/softwares/hmmer2/hmmer−2.3.2/src/” $PERL = “/usr/bin/perl” I also changed…

Continue Reading RNAmmer running

IMPUTE2 – problem with MCMC

IMPUTE2 – problem with MCMC 1 Dear, I’ve run an imputation using the fallowing command Impute -use_prephased_g -known_haps_g output_phased.haps -h 1000GP_Phase3.hap.gz -l 1000GP_Phase3.legend.gz -m genetic_map.txt -int 1 1000000 -iter 30 -o Part1_impute I do not understand why the MCMC iteration are not 30, but 1. (see output) —————- Run parameters…

Continue Reading IMPUTE2 – problem with MCMC

blast protein alignment

28 set blast protein alignment Posted at 20:44h in Sem categoria by BLAST applied the standard genetic code for Query, translating GTG into valine (V). The BLAST is a set of algorithms that attempt to find a short fragment of a query sequence that aligns perfectly with a fragment of…

Continue Reading blast protein alignment

Dice Bioinformatics Engineer Vacancy 2021

Dice Bioinformatics Engineer Vacancy 2021 – Candidates Apply Online Job Title: Bioinformatics Engineer Vertical: Pharmaceutical Address City: San Diego, CA Country: US Position starting out 100% remote due to COVID-19. Anticipated duties. Work closely with world-class scientists, bioinformaticians, and developers to advance the discovery of new biologic medicine. Develop and…

Continue Reading Dice Bioinformatics Engineer Vacancy 2021

bad file format in HMM file

bad file format in HMM file 1 Hi people, I am trying to run hmmpress hmmpress ./data/dbCAN-HMMdb-V8.txt But i have the following error: Working… bad file format in HMM file data/dbCAN-HMMdb-V8.txt Please, help me, Thanks hmmpress hmmer • 29 views A file ending in .txt may not be a HMM…

Continue Reading bad file format in HMM file

BLAST versus HMM search

BLAST versus HMM search 0 Can someone please provide (1) definition of BLAST and HMM searches and (2) describe/explain what the difference is between BLAST and HMM search? PS. This is NOT a homework help. I am trying to understand the concepts behind BLAST and HMM. Thanks in advance! python…

Continue Reading BLAST versus HMM search

Transitional genomes and nutritional role reversals identified for dual symbionts of adelgids (Aphidoidea: Adelgidae)

1. Szathmáry E, Smith JM. The major evolutionary transitions. Nature 1995;374:227–32. PubMed  Google Scholar  2. West SA, Fisher RM, Gardner A, Kiers ET. Major evolutionary transitions in individuality. Proc Natl Acad Sci USA. 2015;112:10112–9. CAS  PubMed  PubMed Central  Google Scholar  3. Moran NA. The coevolution of bacterial endosymbionts and phloem-feeding…

Continue Reading Transitional genomes and nutritional role reversals identified for dual symbionts of adelgids (Aphidoidea: Adelgidae)

Comparative genomic analysis of Methanimicrococcus blatticola provides insights into host adaptation in archaea and the evolution of methanogenesis

1. Hackstein JH, Stumm CK. Methane production in terrestrial arthropods. Proc Natl Acad Sci USA. 1994;91:5441–5. CAS  PubMed  PubMed Central  Article  Google Scholar  2. Hackstein JHP, van Alen TA. Fecal methanogens and vertebrate evolution. Evolution. 1996;50:559–72. PubMed  Article  PubMed Central  Google Scholar  3. Borrel G, McCann A, Deane J, Neto…

Continue Reading Comparative genomic analysis of Methanimicrococcus blatticola provides insights into host adaptation in archaea and the evolution of methanogenesis

How are the HMM cutoff scores in TIGRFAMs determined?

How are the HMM cutoff scores in TIGRFAMs determined? 0 Hello everyone, As indicated in the title, I am wondering how are the HMM cutoff scores in TIGRFAMs determined? Where did I find the score? (If you search a TIGRFAM accession in the NCBI protein family model database, e.g. TIGR02064.1,…

Continue Reading How are the HMM cutoff scores in TIGRFAMs determined?

Can HMMs detect partial proteins sequences?

Can HMMs detect partial proteins sequences? 0 Hello everyone, I am trying to determine the abundance of a gene in a set of metagenomic short reads using the hmm of this gene, then use this method to compare the abundance of the gene across several samples after normalizing for data…

Continue Reading Can HMMs detect partial proteins sequences?

get only one representative fasta sequence per family

Pfam – get only one representative fasta sequence per family 2 Hey can u help me with getting only one representative fasta sequence per family? Is there way to simply do that? cheers X pfam fasta protein • 186 views It’s not trivial. You could use the sequences from the…

Continue Reading get only one representative fasta sequence per family

MAKER genome annotation error with SNAP ab initio prediction

I am trying to do a second round of maker genome annotation with ab initio prediction by snap. The error I am getting is as follows: error: unknown command “genome.hmm”, see ‘snap help’. ERROR: Snap failed –> rank=NA, hostname=bioinformatics ERROR: Failed while preparing ab-inits ERROR: Chunk failed at level:0, tier_type:2…

Continue Reading MAKER genome annotation error with SNAP ab initio prediction

Comment: alphafold online availability and use case

Not my area of expertise particularly but; 1. I don’t think you can use a structure prediction tool to really ‘validate’ HMMER predictions. I’m pretty sure most structure predictors are relying on HMMER or similar HMM based approaches (Martin told me AlphaFold leans on HHBlits API calls for example). I…

Continue Reading Comment: alphafold online availability and use case

So many variants detected.

So many variants detected. 0 Dear All, I have done variant calling in Germline data that has single sample of each individual and two genes. I did following steps, but after checking results I found too many variants. After Haplotypecaller (the step 6) I found 140900 known variants, and the…

Continue Reading So many variants detected.

Help speeding up HMMER’s HMMSearch algorithm for large fasta file with GNU Parallel

I’ve seen that HMMER can be sped up with GNU Parallel: Speed of hmmsearch I have around 100,000 sequences and a HMMER database of around 300 HMM profiles. I’m running everything at once but I’m wondering if it’ll be faster to split up the sequences and/or split up the jobs….

Continue Reading Help speeding up HMMER’s HMMSearch algorithm for large fasta file with GNU Parallel