Hidden Markov model speed heuristic and iterative HMM search procedure

Background  

Profile hidden Markov models (profile-HMMs) are sensitive tools for remote protein homology detection, but the main scoring algorithms, Viterbi or Forward, require considerable time to search large sequence databases.     

Rapid mapping and identification of mutations in Caenorhabditis elegans by restriction site-associated DNA mapping and genomic interval pull-down sequencing

Forward genetic screens provide a powerful approach for inferring gene function on the basis of the phenotypes associated with mutated genes. However, determining the causal mutation by traditional mapping and candidate gene sequencing is often the rate-limiting step, especially when analyzing many mutants. We report two genomic approaches for more rapidly determining the identity of the affected genes in Caenorhabditis elegans mutants. First, we report our use of restriction site-associated DNA (RAD) polymorphism markers for rapidly mapping mutations after chemical mutagenesis and mutant isolation. Second, we describe our use of genomic interval pull-down sequencing (GIPS) to selectively capture and sequence megabase-sized portions of a mutant genome. Together, these two methods provide a rapid and cost-effective approach for positional cloning of C. elegans mutant loci, and are also applicable to other genetic model systems.     

Src-inducible association of CrkL with procaspase-8 promotes cell migration

Procaspase-8, the zymogen form of the apoptosis-initiator caspase-8, undergoes phosphorylation following integrin-mediated cell attachment to an extracellular matrix substrate. Concordant with cell attachment to fibronectin, a population of procaspase-8 becomes associated with a peripheral insoluble compartment that includes focal complexes and lamellar microfilaments. Phosphorylation of procaspase-8 both impairs its maturation to the proapoptotic form and can promote cell migration. Here we show that the cytoskeletal adaptor protein CrkL promotes caspase-8 recruitment to the peripheral spreading edge of cells, and that the catalytic domain of caspase-8 directly interacts with the SH2 domain of CrkL. We show that the interaction is abolished by shRNA-mediated silencing of Src, in Src-deficient MEFs, and by pharmacologic inhibitors of the kinase. The results provide insight into how tyrosine kinases may act to coordinate the suppression caspase-8 mediated apoptosis, while promoting cell invasion.     

Developmental and dominance-associated differences in mushroom body structure in the paper wasp Mischocyttarus mastigophorus

Primitively eusocial paper wasps exhibit considerable plasticity in their division of labor. Dominance interactions among nest mates play a strong role in determining the task performance patterns of adult females. We asked whether dominance status and task performance differences were associated with the development of subregions of the mushroom bodies (MB) of female Mischocyttarus mastigophorus queens and workers. We found that the MB calycal neuropils were better developed (relative to the Kenyon cell body layer) in the dominant females that spent more time on the nest. Increased MB calyx development was more strongly associated with social dominance than with high rates of foraging. The MB of queens resembled those of dominant workers. The results suggest that social interactions are particularly relevant to M. mastigophorus females' cognition. By examining the MB of newly emerged females, we also found evidence for significant age-related changes in MB structure.     

The role of hemiparasitic plants: influencing tallgrass prairie quality,diversity, and structure

Wood betony, Orobanchaceae (Pedicularis canadensis) and bastard toadflax, Santalaceae (Comandra umbellata) are two root‐hemiparasitic plant species found in tallgrass prairie communities. Natural resource managers are interested in utilizing these species as “pseudograzers” in grasslands to reduce competitively dominant grasses and thereby increase ecological diversity and quality in prairie restorations and urban plantings. We performed an observational field study at 5 tallgrass prairie sites to investigate the association of hemiparasite abundance with metrics of phylogenetic and ecological diversity, as well as floristic quality. Although no reduction in C₄ grasses was detected, there was a significant association between hemiparasite abundance and increased floristic quality at all 5 sites. Hemiparasite abundance and species richness were positively correlated at one restoration site. In a greenhouse mesocosm experiment, we investigated response to parasitism by P. canadensis in 6 species representing different plant functional groups of the tallgrass prairie. The annual legume partridge pea, Fabaceae (Chamaecrista fasciculata) had the greatest significant dry biomass reduction among 6 host species, but the C₄ grass big bluestem, Poaceae (Andropogon gerardii) had significantly greater aboveground biomass when grown with the hemiparasite. Overall, host species biomass as a total community was significantly reduced in mesocosms, consistent with other investigations that demonstrate influence on community structure by hemiparasitic plant species. Although hemiparasites were not acting as pseudograzers, they have the potential to influence community structure in grassland restorations and remnants.     

Monitoring early fusion dynamics of human immunodeficiency virus type 1 at single-molecule resolution

The fusion of human immunodeficiency virus type 1 (HIV-1) to host cells is a dynamic process governed by the interaction between glycoproteins on the viral envelope and the major receptor, CD4, and coreceptor on the surface of the cell. How these receptors organize at the virion-cell interface to promote a fusion-competent site is not well understood. Using single-molecule force spectroscopy, we map the tensile strengths, lifetimes, and energy barriers of individual intermolecular bonds between CCR5-tropic HIV-1 gp120 and its receptors CD4 and CCR5 or CXCR4 as a function of the interaction time with the cell. According to the Bell model, at short times of contact between cell and virion, the gp120-CD4 bond is able to withstand forces up to 35 pN and has an initial lifetime of 0.27 s and an intermolecular length of interaction of 0.34 nm. The initial bond also has an energy barrier of 6.7 k(B)T (where k(B) is Boltzmann's constant and T is absolute temperature). However, within 0.3 s, individual gp120-CD4 bonds undergo rapid destabilization accompanied by a shortened lifetime and a lowered tensile strength. This destabilization is significantly enhanced by the coreceptor CCR5, not by CXCR4 or fusion inhibitors, which suggests that it is directly related to a conformational change in the gp120-CD4 bond. These measurements highlight the instability and low tensile strength of gp120-receptor bonds, uncover a synergistic role for CCR5 in the progression of the gp120-CD4 bond, and suggest that the cell-virus adhesion complex is functionally arranged about a long-lived gp120-coreceptor bond.     

The potential and realized spread of wildfires across Canada

Given that they can burn for weeks or months, wildfires in temperate and boreal forests may become immense (eg., 10⁰ – 10⁴ km²). However, during the period within which a large fire is ‘active’, not all days experience weather that is conducive to fire spread; indeed most of the spread occurs on a small proportion (e.g., 1 – 15 days) of not necessarily consecutive days during the active period. This study examines and compares the Canada‐wide patterns in fire‐conducive weather (‘potential’ spread) and the spread that occurs on the ground (‘realized’ spread). Results show substantial variability in distributions of potential and realized spread days across Canada. Both potential and realized spread are higher in western than in eastern Canada; however, whereas potential spread generally decreases from south to north, there is no such pattern with realized spread. The realized‐to‐potential fire‐spread ratio is considerably higher in northern Canada than in the south, indicating that proportionally more fire‐conducive days translate into fire progression. An exploration of environmental correlates to spread show that there may be a few factors compensating for the lower potential spread in northern Canada: a greater proportion of coniferous (i.e., more flammable) vegetation, lesser human impacts (i.e., less fragmented landscapes), sufficient fire ignitions, and intense droughts. Because a linear relationship exists between the frequency distributions of potential spread days and realized spread days in a fire zone, it is possible to obtain one from the other using a simple conversion factor. Our methodology thus provides a means to estimate realized fire spread from weather‐based data in regions where fire databases are poor, which may improve our ability to predict future fire activity.     

New mutations in cloned Escherichia coli umuDC genes: novel phenotypes of strains carrying a umuC125 plasmid

The umuDC locus of Escherichia coli is required for most mutagenesis by UV and many chemicals. Mutations in E. coli umuDC genes cloned on pBR322-derived plasmids wer e isolated by two methods. First, spontaneously-arising mutant umuDC plasmids that failed to confe cold-sensitive growth on a lexA51(Def) strain were isolated by selection. Second, mutant umuDC plasmids that affected apparent mutant yield after UV-irradiation in a strain carrying umuD⁺C⁺ in the chromosome were isolated by screening hydroxylamine-mutagenized umuD⁺C⁺ plasmids. pBR322-derived umuD⁺C⁺ plasmids inhibited the induction of the SOS response of lexA⁺ strains as measured by expression of din::Mu dl(lac) Ap) fusionsbut most mutant plasmids did not. Mutant plasmids defective in complementation of chromosomal umuD44, umuC36, or both were found among those selected for failure to confer cold-sensitivity, whereas those identified by the screening procedure yielded mostly mutant plasmids with more complex phenotypes. We studied in greater detail a plasmid pLM109, carrying the umuC125 mutation. This plasmid increased the sensitivity of lexA⁺ strainsto killing by UV-irradiation but was able to complement the deficiencies of umuC mutants in UV mutagenesis. pLM109 failed to confer cold-sensitive growth on lexA(Def) strains but inhibited SOS induction in lexA⁺ strains. The effect of pLM109 on the UV sensitivity of lexA(Def)strains was similar to that of the parental umuD⁺C⁺ plasmid. The mutation responsible for the phenotypes of pLM109 was localized to a 615-bp fragment. DNA sequencing revealed that the umuC125 mutation was a G:C → A:T transition that changed codon 39 of umuC from GCC → GTC thus changing Ala39 to Val39. The implications of the umuC125 mutation for umuDC-dependent effects on UV-mutagenesis and cell survival after UV damage are discussed.     

The role of apolipoprotein A-I helix 10 in apolipoprotein-mediated cholesterol efflux via the ATP-binding cassette transporter ABCA1

Recent studies of Tangier disease have shown that the ATP-binding cassette transporter A1 (ABCA1)/apolipoprotein A-I (apoA-I) interaction is critical for high density lipoprotein particle formation, apoA-I integrity, and proper reverse cholesterol transport. However, the specifics of this interaction are unknown. It has been suggested that amphipathic helices of apoA-I bind to a lipid domain created by the ABCA1 transporter. Alternatively, apoA-I may bind directly to ABCA1 itself. To better understand this interaction, we created several truncation mutants of apoA-I and then followed up with more specific point mutants and helix translocation mutants to identify and characterize the locations of apoA-I required for ABCA1-mediated cholesterol efflux. We found that deletion of residues 221-243 (helix 10) abolished ABCA1-mediated cholesterol efflux from cultured RAW mouse macrophages treated with 8-bromo-cAMP. Point mutations in helix 10 that affected the helical charge distribution reduced ABCA1-mediated cholesterol efflux versus the wild type. We noted a strong positive correlation between cholesterol efflux and the lipid binding characteristics of apoA-I when mutations were made in helix 10. However, there was no such correlation for helix translocations in other areas of the protein as long as helix 10 remained intact at the C terminus. From these observations, we propose an alternative model for apolipoprotein-mediated efflux.