Hidden Markov Model Search Research Articles

Search for gravitational waves from Scorpius X-1 with a hidden Markov model in O3 LIGO data

Results are presented for a semicoherent search for continuous gravitational waves from the low-mass x-ray binary Scorpius X-1, using a hidden Markov model (HMM) to allow for spin wandering. This search improves on previous HMM-based searches of Laser Interferometer Gravitational-Wave Observatory data by including the orbital period in the search template grid, and by analyzing data from the latest (third) observing run. In the frequency range searched, from 60 to 500 Hz, we find no evidence of gravitational radiation. This is the most sensitive search for Scorpius X-1 using a HMM to date. For the most sensitive subband, starting at 256.06 Hz, we report an upper limit on gravitational wave strain (at 95% confidence) of h095%=6.16×10−26, assuming the orbital inclination angle takes its electromagnetically restricted value ι=44°. The upper limits on gravitational wave strain reported here are on average a factor of ∼3 lower than in the second observing run HMM search. This is the first Scorpius X-1 HMM search with upper limits that reach below the indirect torque-balance limit for certain subbands, assuming ι=44°.Received 26 January 2022Accepted 26 August 2022DOI:https://doi.org/10.1103/PhysRevD.106.062002© 2022 American Physical SocietyPhysics Subject Headings (PhySH)Research AreasGravitational wavesPhysical SystemsBinary starsNeutron stars & pulsarsGravitation, Cosmology & Astrophysics

Genome-wide identification and characterization of the JAZ gene family and its expression patterns under various abiotic stresses in Sorghum bicolor

The jasmonate ZIM domain (JAZ) protein belongs to the TIFY ((TIF[F/Y]XG) domain protein) family, which is composed of several plant-specific proteins that play important roles in plant growth, development, and defense responses. However, the mechanism of the sorghum JAZ family in response to abiotic stress remains unclear. In the present study, a total of 17 JAZ genes were identified in sorghum using a Hidden Markov Model search. In addition, real-time quantification polymerase chain reaction (RT-qPCR) was used to analyze the gene expression patterns under abiotic stress. Based on phylogenetic tree analysis, the sorghum JAZ proteins were mainly divided into nine subfamilies. A promoter analysis revealed that the SbJAZ family contains diverse types of promoter cis-acting elements, indicating that JAZ proteins function in multiple pathways upon stress stimulation in plants. According to RT-qPCR, SbJAZ gene expression is tissue-specific. Additionally, under cold, hot, polyethylene glycol, jasmonic acid, abscisic acid, and gibberellin treatments, the expression patterns of SbJAZ genes were distinctly different, indicating that the expression of SbJAZ genes may be coordinated with different stresses. Furthermore, the overexpression of SbJAZ1 in Escherichia coli was found to promote the growth of recombinant cells under abiotic stresses, such as PEG 6000, NaCl, and 40°C treatments. Altogether, our findings help us to better understand the potential molecular mechanisms of the SbJAZ family in sorghum in response to abiotic stresses.

Unveiling the structural features of CysE: a novel target for therapeutic interventions against persistent mycobacteria

World Health Organization (WHO) reports that one-third of the world’s population is infected with a persistent form of Mycobacterium tuberculosis (M.tb), the causative bacterium responsible for causing the dreaded tuberculosis disease. Targeting mycobacterial persisters is important for achieving WHO’s End TB target. The de-novo cysteine biosynthetic pathway is a novel target for addressing M.tb persistence. The two-step pathway comprises of serine acetyltransferase/CysE and O-acetyl-serine-sulfhydrylase/OASS/CysK. The present study is an attempt to understand the structural features of mycobacterial CysE by investigating the divergence amongst orthologous through phylogenetic analysis. Mapping of mycobacterial CysE sequences on the whole orthologous (COG1045) tree segregated the species into four clusters and several isoforms leading to their descendants identification. Interestingly the analysis revealed that the extended C-terminal α-helix believed unique to M.tb is also present in other organisms such as: Campylobacter ureolyticus, Bacillus cereus, Geminocystis herdmanii and Paenibacillus borealis. Further, the Hidden Markov model search against the whole Uniprot database suggests a plausible role of C-terminal α-helix of CysE in strengthening the substrate and/or co-factor binding. In addition, phylogenetic analysis of CysE sequences from the Mycobacteriaceae family facilitates grouping them under ten well-formed and six monophyletic clades, each based on characteristic features with respect to domain architecture, oligomeric assembly, C-terminal tetra-peptide tail, regulatory and feedback mechanism etc. Employing molecular phylogeny in conjunction with structural analysis has provided detailed insights for mycobacterial CysEs as drug target.

HMS-S-S: A tool for the identification of Sulphur metabolism-related genes and analysis of operon structures in genome and metagenome assemblies.

Sulphur compounds are used in a variety of biological processes including respiration and photosynthesis. Sulphide and sulphur compounds of intermediary oxidation state can serve as electron donors for lithotrophic growth while sulphate, thiosulphate and sulphur are used as electron acceptors in anaerobic respiration. The biochemistry underlying the manifold transformations of inorganic sulphur compounds occurring in sulphur metabolizing prokaryotes is astonishingly complex and knowledge about it has immensely increased over the last years. The advent of next-generation sequencing approaches as well as the significant increase of data availability in public databases has driven focus of environmental microbiology to probing the metabolic capacity of microbial communities by analysis of this sequence data. To facilitate these analyses, we created HMS-S-S, a comprehensive equivalogous hidden Markov model (HMM)-supported tool. Protein sequences related to sulphur compound oxidation, reduction, transport and intracellular transfer are efficiently detected and related enzymes involved in dissimilatory sulphur oxidation as opposed to sulphur compound reduction can be confidently distinguished. HMM search results are coupled to corresponding genes, which allows analysis of co-occurrence, synteny and genomic neighbourhood. The HMMs were validated on an annotated test data set and by cross-validation. We also proved its performance by exploring meta-assembled genomes isolated from samples from environments with active sulphur cycling, including members of the cable bacteria, novel Acidobacteria and assemblies from a sulphur-rich glacier, and were able to replicate and extend previous reports.

The STI1-domain is a flexible alpha-helical fold with a hydrophobic groove.

STI1-domains are present in a variety of co-chaperone proteins and are required for the transfer of hydrophobic clients in various cellular processes. The domains were first identified in the yeast Sti1 protein where they were referred to as DP1 and DP2. Based on hidden Markov model searches, this domain had previously been found in other proteins including the mammalian co-chaperone SGTA, the DNA damage response protein Rad23, and the chloroplast import protein Tic40. Here, we refine the domain definition and carry out structure-based sequence alignment of STI1-domains showing conservation of five amphipathic helices. Upon examinations of these identified domains, we identify a preceding helix 0 and unifying sequence properties, determine new molecular models, and recognize that STI1-domains nearly always occur in pairs. The similarity at the sequence, structure, and molecular levels likely supports a unified functional role.

Genome-Wide Mining of MYB Transcription Factors in the Anthocyanin Biosynthesis Pathway of Gossypium Hirsutum.

The MYB family, one of the largest transcription factor (TF) families, plays an important role in plant growth, development, and stress response. Although genome-wide analysis of the MYB family has been performed in many species based on sequence similarity, predicting the potential functions of the MYB genes and classifying the regulators into specific metabolic pathways remains difficult. In this study, using a hidden Markov model search and co-expression regulatory network analysis, we demonstrated a process to screen and identify potential MYB TFs in the anthocyanin biosynthesis pathway of Gossypium hirsutum. As a result, we identified 617 and 784 MYB genes (812 in total) from the previously reported and recently released genomes, respectively. Using 126 structural genes involved in the anthocyanin biosynthesis pathway as targets for several co-expression network analyses, we sorted out 31 R2R3-MYB genes, which are potential regulators in the specific pathway. Phylogenetic and collinearity analyses indicated that 83.9% of the 31 MYB genes originated from whole genome duplication or polyploidization. In addition, we revealed relatively specific regulatory relationships between the MYB TFs and their target structural genes. Approximately, 71% of the MYBs could regulate only a single anthocyanin-related structural gene. Moreover, we found that the A- and D- subgenome homoeologs of MYB TFs in G. hirsutum rarely co-regulate the same target gene. The current study not only demonstrated an easy method to rapidly predict potential TFs in a specific metabolic pathway, but also enhanced our understanding of the evolution, gene characteristics, expression, and regulatory pattern of MYB TFs in G. hirsutum.

Completion of the gut microbial epi-bile acid pathway

ABSTRACT Bile acids are detergent molecules that solubilize dietary lipids and lipid-soluble vitamins. Humans synthesize bile acids with α-orientation hydroxyl groups which can be biotransformed by gut microbiota to toxic, hydrophobic bile acids, such as deoxycholic acid (DCA). Gut microbiota can also convert hydroxyl groups from the α-orientation through an oxo-intermediate to the β-orientation, resulting in more hydrophilic, less toxic bile acids. This interconversion is catalyzed by regio- (C-3 vs. C-7) and stereospecific (α vs. β) hydroxysteroid dehydrogenases (HSDHs). So far, genes encoding the urso- (7α-HSDH & 7β-HSDH) and iso- (3α-HSDH & 3β-HSDH) bile acid pathways have been described. Recently, multiple human gut clostridia were reported to encode 12α-HSDH, which interconverts DCA and 12-oxolithocholic acid (12-oxoLCA). 12β-HSDH completes the epi-bile acid pathway by converting 12-oxoLCA to the 12β-bile acid denoted epiDCA; however, a gene(s) encoding this enzyme has yet to be identified. We confirmed 12β-HSDH activity in cultures of Clostridium paraputrificum ATCC 25780. From six candidate C. paraputrificum ATCC 25780 oxidoreductase genes, we discovered the first gene (DR024_RS09610) encoding bile acid 12β-HSDH. Phylogenetic analysis revealed unforeseen diversity for 12β-HSDH, leading to validation of two additional bile acid 12β-HSDHs through a synthetic biology approach. By comparison to a previous phylogenetic analysis of 12α-HSDH, we identified the first potential C-12 epimerizing strains: Collinsella tanakaei YIT 12063 and Collinsella stercoris DSM 13279. A Hidden Markov Model search against human gut metagenomes located putative 12β-HSDH genes in about 30% of subjects within the cohorts analyzed, indicating this gene is relevant in the human gut microbiome.

Search for gravitational waves from 12 young supernova remnants with a hidden Markov model in Advanced LIGO’s second observing run

Persistent gravitational waves from rapidly rotating neutron stars, such as those found in some young supernova remnants, may fall in the sensitivity band of the advanced Laser Interferometer Gravitational-wave Observatory (aLIGO). Searches for these signals are computationally challenging, as the frequency and frequency derivative are unknown and evolve rapidly due to the youth of the source. A hidden Markov model (HMM), combined with a maximum-likelihood matched filter, tracks rapid frequency evolution semi-coherently in a computationally efficient manner. We present the results of an HMM search targeting 12 young supernova remnants in data from Advanced LIGO's second observing run. Six targets produce candidates that are above the search threshold and survive pre-defined data quality vetoes. However, follow-up analyses of these candidates show that they are all consistent with instrumental noise artefacts.

Bringing to light the molecular evolution of GUX genes in plants.

Hemicellulose and cellulose are essential polysaccharides for plant development and major components of cell wall. They are also an important energy source for the production of ethanol from plant biomass, but their conversion to fermentable sugars is hindered by the complex structure of cell walls. The glucuronic acid substitution of xylan (GUX) enzymes attach glucuronic acid to xylan, a major component of hemicellulose, decreasing the efficiency of enzymes used for ethanol production. Since loss-of-function gux mutants of Arabidopsis thaliana enhance enzyme accessibility and cell wall digestion without adverse phenotypes, GUX genes are potential targets for genetically improving energy crops. However, comprehensive identification of GUX in important species and their evolutionary history are largely lacking. Here, we identified putative GUX proteins using hidden Markov model searches with the GT8 domain and a GUX-specific motif, and inferred the phylogenetic relationship of 18 species with Maximum likelihood and Bayesian approaches. Each species presented a variable number of GUX, and their evolution can be explained by a mixture of divergent, concerted and birth-and-death evolutionary models. This is the first broad insight into the evolution of GUX gene family in plants and will potentially guide genetic and functional studies in species used for biofuel production.

Distribution and functional analysis of the phosphopantetheinyl transferase superfamily in Actinomycetales microorganisms

Phosphopantetheinyl transferases (PPTases) are a superfamily of essential enzymes required for the synthetic processes of many compounds including fatty acid, polyketide, and nonribosomal peptide metabolites. These enzymes activate carrier proteins in specific biosynthetic pathways via the transfer of a phosphopantetheinyl moiety to a serine residue in the conserved motif of carrier proteins. Since many Actinomycetales microorganisms produce a number of polyketide and nonribosomal peptide metabolites, the distribution of PPTase genes was investigated in these microorganisms. PPTases were found in bacterial protein databases using a hidden Markov model search with the PF01648 (4'-phosphopantetheinyl transferase superfamily) model. Actinomycetales microorganisms harbor several genes encoding AcpS-type and Sfp-type PPTases in individual genomes, many of which were associated with the biosynthetic gene cluster for polyketide or nonribosomal peptide metabolites. The properties of these PPTases were evaluated in the heterologous expression system using the biosynthetic gene clusters and genes encoding PPTases found in the present study. Sfp-type PPTases were classified into two subgroups, and although the substrate specificities of the enzymes in one subgroup were wide, the catalytic activities of enzymes in the other subgroup were low. SAV_1784 of Streptomyces avermitilis possessed the most characteristic broad-range activity against several type I polyketide synthases and nonribosomal peptide synthetases.

Transcriptome-wide identification and screening of WRKY factors involved in the regulation of taxol biosynthesis in Taxus chinensis

WRKY, a plant-specific transcription factor family, plays important roles in pathogen defense, abiotic cues, phytohormone signaling, and regulation of plant secondary metabolism. However, little is known about the roles, functions, and mechanisms of WRKY in taxane biosynthesis in Taxus spp. In this study, 61 transcripts were identified from Taxus chinensis transcriptome datasets by using hidden Markov model search. All of these transcripts encoded proteins containing WRKY domains, which were designated as TcWRKY1–61. After phylogenetic analysis of the WRKY domains of TcWRKYs and AtWRKYs, 16, 8, 10, 14, 5, 7, and 1 TcWRKYs were cladded into Group I, IIa–IIe, and III, respectively. Then, six representative TcWRKYs were selected to classify their effects on taxol biosynthesis. After MeJA (methyl jasmonate acid) and SA (salicylic acid) treatments, all of the six TcWRKYs were upregulated by MeJA treatment. TcWRKY44 (IId) and TcWRKY47 (IIa) were upregulated, whereas TcWRKY8 (IIc), TcWRKY20 (III), TcWRKY26 (I), TcWRKY41 (IIe), and TcWRKY52 (IIb) were downregulated by SA treatment. Overexpression experiments showed that the six selected TcWRKYs exerted different effects on taxol biosynthesis. In specific, TcWRKY8 and TcWRKY47 significantly improved the expression levels of taxol-biosynthesis-related genes. Transcriptome-wide identification of WRKY factors in Taxus not only enhances our understanding of plant WRKY factors but also identifies candidate regulators of taxol biosynthesis.

Evolutionary history and functional divergence of the cytochrome P450 gene superfamily between Arabidopsis thaliana and Brassica species uncover effects of whole genome and tandem duplications

BackgroundThe cytochrome P450 monooxygenase (P450) superfamily is involved in the biosynthesis of various primary and secondary metabolites. However, little is known about the effects of whole genome duplication (WGD) and tandem duplication (TD) events on the evolutionary history and functional divergence of P450s in Brassica after splitting from a common ancestor with Arabidopsis thaliana.ResultsUsing Hidden Markov Model search and manual curation, we detected that Brassica species have nearly 1.4-fold as many P450 members as A. thaliana. Most P450s in A. thaliana and Brassica species were located on pseudo-chromosomes. The inferred phylogeny indicated that all P450s were clustered into two different subgroups. Analysis of WGD event revealed that different P450 gene families had appeared after evolutionary events of species. For the TD event analyses, the P450s from TD events in Brassica species can be divided into ancient and recent parts. Our comparison of influence of WGD and TD events on the P450 gene superfamily between A. thaliana and Brassica species indicated that the family-specific evolution in the Brassica lineage can be attributed to both WGD and TD, whereas WGD was recognized as the major mechanism for the recent evolution of the P450 super gene family. Expression analysis of P450s from A. thaliana and Brassica species indicated that WGD-type P450s showed the same expression pattern but completely different expression with TD-type P450s across different tissues in Brassica species. Selection force analysis suggested that P450 orthologous gene pairs between A. thaliana and Brassica species underwent negative selection, but no significant differences were found between P450 orthologous gene pairs in A. thaliana–B. rapa and A. thaliana–B. oleracea lineages, as well as in different subgenomes in B. rapa or B. oleracea compared with A. thaliana.ConclusionsThis study is the first to investigate the effects of WGD and TD on the evolutionary history and functional divergence of P450 gene families in A. thaliana and Brassica species. This study provides a biology model to study the mechanism of gene family formation, particularly in the context of the evolutionary history of angiosperms, and offers novel insights for the study of angiosperm genomes.

GsCHX19.3, a member of cation/H+ exchanger superfamily from wild soybean contributes to high salinity and carbonate alkaline tolerance

Cation/H+ exchangers (CHX) are characterized to be involved in plant growth, development and stress responses. Although soybean genome sequencing has been completed, the CHX family hasn’t yet been systematically analyzed, especially in wild soybean. Here, through Hidden Markov Model search against Glycine soja proteome, 34 GsCHXs were identified and phylogenetically clustered into five groups. Members within each group showed high conservation in motif architecture. Interestingly, according to our previous RNA-seq data, only Group IVa members exhibited highly induced expression under carbonate alkaline stress. Among them, GsCHX19.3 displayed the greatest up-regulation in response to carbonate alkaline stress, which was further confirmed by quantitative real-time PCR analysis. We also observed the ubiquitous expression of GsCHX19.3 in different tissues and its localization on plasma membrane. Moreover, we found that GsCHX19.3 expression in AXT4K, a yeast mutant lacking four ion transporters conferred resistance to low K+ at alkali pH, as well as carbonate stress. Consistently, in Arabidopsis, GsCHX19.3 overexpression increased plant tolerance both to high salt and carbonate alkaline stresses. Furthermore, we also confirmed that GsCHX19.3 transgenic lines showed lower Na+ concentration but higher K+/Na+ values under salt-alkaline stress. Taken together, our findings indicated that GsCHX19.3 contributed to high salinity and carbonate alkaline tolerance.

Genome-wide analysis of UDP-glycosyltransferase super family in Brassica rapa and Brassica oleracea reveals its evolutionary history and functional characterization

BackgroundGlycosyltransferases comprise a highly divergent and polyphyletic multigene family that is involved in widespread modification of plant secondary metabolites in a process called glycosylation. According to conserved domains identified in their amino acid sequences, these glycosyltransferases can be classified into a single UDP-glycosyltransferase (UGT) 1 superfamily.ResultsWe performed genome-wide comparative analysis of UGT genes to trace evolutionary history in algae, bryophytes, pteridophytes, and angiosperms; then, we further investigated the expansion mechanisms and function characterization of UGT gene families in Brassica rapa and Brassica oleracea. Using Hidden Markov Model search, we identified 3, 21, 140, 200, 115, 147, and 147 UGTs in Chlamydomonas reinhardtii, Physcomitrella patens, Selaginella moellendorffii, Oryza sativa, Arabidopsis thaliana, B. rapa, and B. oleracea, respectively. Phylogenetic analysis revealed that UGT80 gene family is an ancient gene family, which is shared by all plants and UGT74 gene family is shared by ferns and angiosperms, but the remaining UGT gene families were shared by angiosperms. In dicot lineage, UGTs among three species were classified into three subgroups containing 3, 6, and 12 UGT gene families. Analysis of chromosomal distribution indicates that 98.6 and 71.4% of UGTs were located on B. rapa and B. oleracea pseudo-molecules, respectively. Expansion mechanism analyses uncovered that whole genome duplication event exerted larger influence than tandem duplication on expansion of UGT gene families in B. rapa, and B. oleracea. Analysis of selection forces of UGT orthologous gene pairs in B. rapa, and B. oleracea compared to A. thaliana suggested that orthologous genes in B. rapa, and B. oleracea have undergone negative selection, but there were no significant differences between A. thaliana –B. rapa and A. thaliana –B. oleracea lineages. Our comparisons of expression profiling illustrated that UGTs in B. rapa performed more discrete expression patterns than these in B. oleracea indicating stronger function divergence. Combing with phylogeny and expression analysis, the UGTs in B. rapa and B. oleracea experienced parallel evolution after they diverged from a common ancestor.ConclusionWe first traced the evolutionary history of UGT gene families in plants and revealed its evolutionary and functional characterization of UGTs in B. rapa, and B. oleracea. This study provides novel insights into the evolutionary history and functional divergence of important traits or phenotype-related gene families in plants.

The mRNA-bound proteome of the human malaria parasite Plasmodium falciparum.

BackgroundGene expression is controlled at multiple levels, including transcription, stability, translation, and degradation. Over the years, it has become apparent that Plasmodium falciparum exerts limited transcriptional control of gene expression, while at least part of Plasmodium’s genome is controlled by post-transcriptional mechanisms. To generate insights into the mechanisms that regulate gene expression at the post-transcriptional level, we undertook complementary computational, comparative genomics, and experimental approaches to identify and characterize mRNA-binding proteins (mRBPs) in P. falciparum.ResultsClose to 1000 RNA-binding proteins are identified by hidden Markov model searches, of which mRBPs encompass a relatively large proportion of the parasite proteome as compared to other eukaryotes. Several abundant mRNA-binding domains are enriched in apicomplexan parasites, while strong depletion of mRNA-binding domains involved in RNA degradation is observed. Next, we experimentally capture 199 proteins that interact with mRNA during the blood stages, 64 of which with high confidence. These captured mRBPs show a significant overlap with the in silico identified candidate RBPs (p < 0.0001). Among the experimentally validated mRBPs are many known translational regulators active in other stages of the parasite’s life cycle, such as DOZI, CITH, PfCELF2, Musashi, and PfAlba1–4. Finally, we also detect several proteins with an RNA-binding domain abundant in Apicomplexans (RAP domain) that is almost exclusively found in apicomplexan parasites.ConclusionsCollectively, our results provide the most complete comparative genomics and experimental analysis of mRBPs in P. falciparum. A better understanding of these regulatory proteins will not only give insight into the intricate parasite life cycle but may also provide targets for novel therapeutic strategies.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-016-1014-0) contains supplementary material, which is available to authorized users.

Core cell cycle regulatory genes in rice and their expression profiles across the growth zone of the leaf.

Rice (Oryza sativa L.) as a model and crop plant with a sequenced genome offers an outstanding experimental system for discovering and functionally analyzing the major cell cycle control elements in a cereal species. In this study, we identified the core cell cycle genes in the rice genome through a hidden Markov model search and multiple alignments supported with the use of short protein sequence probes. In total we present 55 rice putative cell cycle genes with locus identity, chromosomal location, approximate chromosome position and EST accession number. These cell cycle genes include nine cyclin dependent-kinase (CDK) genes, 27 cyclin genes, one CKS gene, two RBR genes, nine E2F/DP/DEL genes, six KRP genes, and one WEE gene. We also provide characteristic protein sequence signatures encoded by CDK and cyclin gene variants. Promoter analysis by the FootPrinter program discovered several motifs in the regulatory region of the core cell cycle genes. As a first step towards functional characterization we performed transcript analysis by RT-PCR to determine gene specific variation in transcript levels along the rice leaves. The meristematic zone of the leaves where cells are actively dividing was identified based on kinematic analysis and flow cytometry. As expected, expression of the majority of cell cycle genes was exclusively associated with the meristematic region. However genes such as different D-type cyclins, DEL1, KRP1/3, and RBR2 were also expressed in leaf segments representing the transition zone in which cells start differentiation.

Phylogenomics reveals surprising sets of essential and dispensable clades of MIKC(c)-group MADS-box genes in flowering plants.

MIKC(C)-group MADS-box genes are involved in the control of many developmental processes in flowering plants. All of these genes are members of one of 17 clades that had already been established in the most recent common ancestor (MRCA) of extant angiosperms. These clades trace back to 11 seed plant-specific superclades that were present in the MRCA of extant seed plants. Due to their important role in plant development and evolution, the origin of the clades of MIKC(C)-group genes has been studied in great detail. In contrast, whether any of these ancestral clades has ever been lost completely in any species has not been investigated so far. Here, we determined the presence of these clades by BLAST, PSI-BLAST, and Hidden Markov Model searches and by phylogenetic methods in the whole genomes of 27 flowering plants. Our data suggest that there are only three superclades of which all members have been lost in at least one of the investigated flowering plant species, and only few additional losses of angiosperm-specific MIKC(C)-group gene clades could be identified. Remarkably, for one seed plant superclade (TM8-like genes) and one angiosperm clade (FLC-like genes), multiple losses were identified, suggesting that the function of these genes is dispensable or that gene loss might have even been adaptive. The clades of MIKC(C)-group genes that have never been wiped out in any of the investigated species comprises, in addition to the expected floral organ identity genes, also TM3-like (SOC1-like), StMADS11-like (SVP-like), AGL17-like and GGM13-like (Bsister) genes, suggesting that these genes are more important for angiosperm development and evolution than has previously been appreciated.

Autoregressive and Iterative Hidden Markov Models for Periodicity Detection and Solenoid Structure Recognition in Protein Sequences

Traditional signal processing methods cannot detect interspersed repeats and generally cannot handle nonstationary signals. In this paper, we propose a new method for periodicity detection in protein sequences to locate interspersed repeats. We first apply the autoregressive model with a sliding window to find possible repeating subsequences within a protein sequence. Then, we utilize an iterative hidden Markov model (HMM) to count the number of subsequences similar to each of the possible repeating subsequences. An iterative HMM search of the potential repeating subsequences can help identify interspersed repeats. Finally, the numbers of repeating subsequences are aggregated together as a feature and used in the classification process. Experiment results show that our method improves the performance of solenoid protein recognition substantially.

Response to Zarsky et al.

Summary Mitochondria evolved from an α-proteobacterial endosymbiont and recent phylogenetic and function-based research has demonstrated that the major pieces of the protein transport machinery were inherited from the symbiont. This includes the SAM machinery for assembly of outer membrane proteins and the TIM machinery for protein transport across, and assembly into, the mitochondrial inner membrane [1–3]. Hidden Markov model (HMM) analysis, which enables a broad, all-encompassing approach for identifying protein homologies, has been very important in detecting members of protein families that are not easily recognized by simple BLAST-based comparisons [1]; HMM searches initially failed to find a Tom40 protein in one group of eukaryotes, the kinetoplastids. These organisms, which include the experimentally-tractable Trypanosoma brucei , have highly developed mitochondria that have evolved from the same ancestor as mitochondria in other eukaryotes. The initial failure to identify a Tom40 homolog in T. brucei was both surprising and exciting.

Generation, annotation, analysis and database integration of 16,500 white spruce EST clusters

BackgroundThe sequencing and analysis of ESTs is for now the only practical approach for large-scale gene discovery and annotation in conifers because their very large genomes are unlikely to be sequenced in the near future. Our objective was to produce extensive collections of ESTs and cDNA clones to support manufacture of cDNA microarrays and gene discovery in white spruce (Picea glauca [Moench] Voss).ResultsWe produced 16 cDNA libraries from different tissues and a variety of treatments, and partially sequenced 50,000 cDNA clones. High quality 3' and 5' reads were assembled into 16,578 consensus sequences, 45% of which represented full length inserts. Consensus sequences derived from 5' and 3' reads of the same cDNA clone were linked to define 14,471 transcripts. A large proportion (84%) of the spruce sequences matched a pine sequence, but only 68% of the spruce transcripts had homologs in Arabidopsis or rice. Nearly all the sequences that matched the Populus trichocarpa genome (the only sequenced tree genome) also matched rice or Arabidopsis genomes. We used several sequence similarity search approaches for assignment of putative functions, including blast searches against general and specialized databases (transcription factors, cell wall related proteins), Gene Ontology term assignation and Hidden Markov Model searches against PFAM protein families and domains. In total, 70% of the spruce transcripts displayed matches to proteins of known or unknown function in the Uniref100 database (blastx e-value < 1e-10). We identified multigenic families that appeared larger in spruce than in the Arabidopsis or rice genomes. Detailed analysis of translationally controlled tumour proteins and S-adenosylmethionine synthetase families confirmed a twofold size difference. Sequences and annotations were organized in a dedicated database, SpruceDB. Several search tools were developed to mine the data either based on their occurrence in the cDNA libraries or on functional annotations.ConclusionThis report illustrates specific approaches for large-scale gene discovery and annotation in an organism that is very distantly related to any of the fully sequenced genomes. The ArboreaSet sequences and cDNA clones represent a valuable resource for investigations ranging from plant comparative genomics to applied conifer genetics.