Estimating synonymous and nonsynonymous substitution rates

Partitioning the total substitution rate into synnonymous and nonsynonymous components is a key aspect of many analyses in molecular evolution. Numerous methods exist for estimating these rates. However, until recently none of the estimation procedures were based on a sound statistical footing. In this paper, the evolutionary model of Muse and Gaut (1994) is used as the basis for two sets of parameters quantifying silent and replacement substitution rates. The parameters are shown to be equal when the four nucleotides are equally frequent and unequal otherwise. Maximum-likelihood estimation of these parameters is described, and the performance of these estimates is compared to that of existing estimation procedures. It is shown that the estimates of Nei and Gojobori (1986) are not unbiased for either set of parameters, although they provide very good estimates for one set as long as sequence divergence is not too high. However, some disturbing properties are found for the Nei and Gojobori estimates. In particular, it is shown that the expected value of the Nei and Gojobori estimate of silent substitution rate is a function of both the silent and replacement substitution rates. The maximum-likelihood estimates have no such problems.      

A new species of Heth Cobb, 1898 (Nematoda: Rhigonematida: Hethidae) from an Australian millipede (Myriapoda: Diplopoda: Spirostreptida)

Abstract

A new species of nematode, Heth baudini sp. n. from a diplopod (Spirostreptida: Iulomorphidae Verhoeff, 1924) collected in Queensland, Australia, is described and illustrated. The cephalic and cervical cuticular ornamentation of females of H. baudini sp. n. is similar to those of South-East Asian and Australasian Heth species. Heth baudini sp. n. females are particularly close to Heth taynguyeni from Vietnam but can be distinguished by the shape of the lateral lappets, which in H. taynguyeni limit the trapezium-shaped region of smooth cuticle unlike the elliptical region in H. baudini sp. n., and by the presence of lateral spines only half the size. The cuticle of the H. baudini sp. n. is finely annulated along the entire body, whereas H. taynguyeni has broader rings behind the first pair of lateral spines, each consisting of five or six narrower rings separated from each other by deeper furrows. Males of H. baudini sp. n. are characterised by the presence of a bursa-like fold on the tail and can be distinguished from other species of the genus by the presence of somatic papillae embedded into the bursal fold.     

Blood flow index using near-infrared spectroscopy and indocyanine green as a minimally invasive tool to assess respiratory muscle blood flow in humans

Near-infrared spectroscopy (NIRS) in combination with indocyanine green (ICG) dye has recently been used to measure respiratory muscle blood flow (RMBF) in humans. This method is based on the Fick principle and is determined by measuring ICG in the respiratory muscles using transcutaneous NIRS in relation to the [ICG] in arterial blood as measured using photodensitometry. This method is invasive since it requires arterial cannulation, repeated blood withdrawals, and reinfusions. A less invasive alternative is to calculate a relative measure of blood flow known as the blood flow index (BFI), which is based solely on the NIRS ICG curve, thus negating the need for arterial cannulation. Accordingly, the purpose of this study was to determine whether BFI can be used to measure RMBF at rest and during voluntary isocapnic hyperpnea at 25, 40, 55, and 70% of maximal voluntary ventilation in seven healthy humans. BFI was calculated as the change in maximal [ICG] divided by the rise time of the NIRS-derived ICG curve. Intercostal and sternocleidomastoid muscle BFI were correlated with simultaneously measured work of breathing and electromyography (EMG) data from the same muscles. BFI showed strong relationships with the work of breathing and EMG for both respiratory muscles. The coefficients of determination (R(2)) comparing BFI vs. the work of breathing for the intercostal and sternocleidomastoid muscles were 0.887 (P < 0.001) and 0.863 (P < 0.001), respectively, whereas the R(2) for BFI vs. EMG for the intercostal and sternocleidomastoid muscles were 0.879 (P < 0.001) and 0.930 (P < 0.001), respectively. These data suggest that the BFI closely reflects RMBF in conscious humans across a wide range of ventilations and provides a less invasive and less technically demanding alternative to measuring RMBF.     

Risks and probabilities of breast cancer: short-term versus lifetime probabilities.

OBJECTIVE: To calculate age-specific short-term and lifetime probabilities of breast cancer among a cohort of Canadian women. DESIGN: Double decrement life table. SETTING: Alberta. SUBJECTS: Women with first invasive breast cancers registered with the Alberta Cancer Registry between 1985 and 1987. MAIN OUTCOME MEASURES: Lifetime probability of breast cancer from birth and for women at various ages; short-term (up to 10 years) probability of breast cancer for women at various ages. RESULTS: The lifetime probability of breast cancer is 10.17% at birth and peaks at 10.34% at age 25 years, after which it decreases owing to a decline in the number of years over which breast cancer risk will be experienced. However, the probability of manifesting breast cancer in the next year increases steadily from the age of 30 onward, reaching 0.36% at 85 years. The probability of manifesting the disease within the next 10 years peaks at 2.97% at age 70 and decreases thereafter, again owing to declining probabilities of surviving the interval. CONCLUSIONS: Given that the incidence of breast cancer among Albertan women during the study period was similar to the national average, we conclude that currently more than 1 in 10 women in Canada can expect to have breast cancer at some point during their life. However, risk varies considerably over a woman''s lifetime, with most risk concentrated after age 49. On the basis of the shorter-term age-specific risks that we present, the clinician can put breast cancer risk into perspective for younger women and heighten awareness among women aged 50 years or more.     

c-Abl has high intrinsic tyrosine kinase activity that is stimulated by mutation of the Src homology 3 domain and by autophosphorylation at two distinct regulatory tyrosines

Using the specific Abl tyrosine kinase inhibitor STI 571, we purified unphosphorylated murine type IV c-Abl and measured the kinetic parameters of c-Abl tyrosine kinase activity in a solution with a peptide-based assay. Unphosphorylated c-Abl exhibited substantial peptide kinase activity with K(m) of 204 microm and V(max) of 33 pmol min(-1). Contrary to previous observations using immune complex kinase assays, we found that a transforming c-Abl mutant with a Src homology 3 domain point mutation (P131L) had significantly (about 6-fold) higher intrinsic kinase activity than wild-type c-Abl (K(m) = 91 microm, V(max) = 112 pmol min(-1)). Autophosphorylation stimulated the activity of wild-type c-Abl about 18-fold and c-Abl P131L about 3.6-fold, resulting in highly active kinases with similar catalytic rates. The autophosphorylation rate was dependent on Abl protein concentration consistent with an intermolecular reaction. A tyrosine to phenylalanine mutation (Y412F) at the c-Abl residue homologous to the c-Src catalytic domain autophosphorylation site impaired the activation of wild-type c-Abl by 90% but reduced activation of c-Abl P131L by only 45%. Mutation of a tyrosine (Tyr-245) in the linker region between the Src homology 2 and catalytic domains that is conserved among the Abl family inhibited the autophosphorylation-induced activation of wild-type c-Abl by 50%, whereas the c-Abl Y245F/Y412F double mutant was minimally activated by autophosphorylation. These results support a model where c-Abl is inhibited in part through an intramolecular Src homology 3-linker interaction and stimulated to full catalytic activity by sequential phosphorylation at Tyr-412 and Tyr-245.     

Detailed protein sequence alignment based on Spectral Similarity Score (SSS)

Background  

The chemical property and biological function of a protein is a direct consequence of its primary structure. Several algorithms have been developed which determine alignment and similarity of primary protein sequences. However, character based similarity cannot provide insight into the structural aspects of a protein. We present a method based on spectral similarity to compare subsequences of amino acids that behave similarly but are not aligned well by considering amino acids as mere characters. This approach finds a similarity score between sequences based on any given attribute, like hydrophobicity of amino acids, on the basis of spectral information after partial conversion to the frequency domain.     

Proffered papers 14.00–15.00 Tuesday 16 September 2003 8 Ampullary,pancreatic and biliary duct exfoliative brush cytology prepared by the thinprep method: performance and morphology

Background ERCP‐directed brush cytology is used to sample lesions of the pancreatic and biliary ducts and the ampulla of Vater. With conventional preparations, the sensitivity and specificity range from 44% to 63% and 80% to 98%, respectively, and increased N : C ratio, nuclear molding and loss of honeycombing are reliable features of malignancy. The performance and morphology of specimens prepared by ThinPrep, a liquid‐based cytology technique is mostly unknown. Methods The laboratory information system was searched for all cases prepared by ThinPrep. Patient disease classification of benign or malignant was determined by linkage with the provincial cancer registry and was the gold standard against which sensitivity, specificity, positive (PPV) and negative predictive values (NPV) were calculated. True positives and negatives were reviewed to identify predictive cytomorphologic features. Results Between 1996 and 2001, there were 149 ThinPrep specimens; 55 (37%) were reported as positive for malignancy and 94 (63%) as negative. Disease was classified as malignant in 86 (58%) patients and benign in 63 (42%). There were 42 false negative, 11 false positive, 52 true negative, and 44 true positive cytology results. Sensitivity was 51.2% (CI; 40.2 : 62.0), specificity 82.5% (CI; 70.5 : 90.6), and PPV and NPV 80.0% (CI; 66.6 : 89) and 55.3% (CI; 44.7 : 65.5), respectively. Cell groups with crowded, enlarged, irregular nuclei and nuclear features of vesicular chromatin and large, multiple irregular nucleoli correlated with malignant disease, while monolayered sheets of uniform columnar cells, regular nuclei and a finely granular chromatin correlated with benign disease. Conclusions The performance of ThinPrep brushings from this anatomic site equals conventional preparations. Cytomorphologic features of malignancy are more frequent and pronounced with ThinPrep.     

The structure of mouse HP1 suggests a unique mode of single peptide recognition by the shadow chromo domain dimer

The heterochromatin protein 1 (HP1) family of proteins is involved in gene silencing via the formation of heterochromatic structures. They are composed of two related domains: an N-terminal chromo domain and a C-terminal shadow chromo domain. Present results suggest that chromo domains may function as protein interaction motifs, bringing together different proteins in multi-protein complexes and locating them in heterochromatin. We have previously determined the structure of the chromo domain from the mouse HP1beta protein, MOD1. We show here that, in contrast to the chromo domain, the shadow chromo domain is a homodimer. The intact HP1beta protein is also dimeric, where the interaction is mediated by the shadow chromo domain, with the chromo domains moving independently of each other at the end of flexible linkers. Mapping studies, with fragments of the CAF1 and TIF1beta proteins, show that an intact, dimeric, shadow chromo domain structure is required for complex formation.