The impact of genotyping error on haplotype reconstruction and frequency estimation

The choice of genotyping families vs unrelated individuals is a critical factor in any large-scale linkage disequilibrium (LD) study. The use of unrelated individuals for such studies is promising, but in contrast to family designs, unrelated samples do not facilitate detection of genotyping errors, which have been shown to be of great importance for LD and linkage studies and may be even more important in genotyping collaborations across laboratories. Here we employ some of the most commonly-used analysis methods to examine the relative accuracy of haplotype estimation using families vs unrelateds in the presence of genotyping error. The results suggest that even slight amounts of genotyping error can significantly decrease haplotype frequency and reconstruction accuracy, that the ability to detect such errors in large families is essential when the number/complexity of haplotypes is high (low LD/common alleles). In contrast, in situations of low haplotype complexity (high LD and/or many rare alleles) unrelated individuals offer such a high degree of accuracy that there is little reason for less efficient family designs. Moreover, parent-child trios, which comprise the most popular family design and the most efficient in terms of the number of founder chromosomes per genotype but which contain little information for error detection, offer little or no gain over unrelated samples in nearly all cases, and thus do not seem a useful sampling compromise between unrelated individuals and large families. The implications of these results are discussed in the context of large-scale LD mapping projects such as the proposed genome-wide haplotype map.

Equine respiratory viruses in foals in New Zealand

AIMS: To identify the respiratory viruses that are present among foals in New Zealand and to establish the age at which foals first become infected with these viruses. METHODS: Foals were recruited to the study in October/ November 1995 at the age of 1 month (Group A) or in March/ April 1996 at the age of 4-6 months (Groups B and C). Nasal swabs and blood samples were collected at monthly intervals. Nasal swabs and peripheral blood leucocytes (PBL) harvested from heparinised blood samples were used for virus isolation; serum harvested from whole-blood samples was used for serological testing for the presence of antibodies against equine herpesvirus (EHV)-1 or -4, equine rhinitis-A virus (ERAV), equine rhinitis-B virus (ERBV), equine adenovirus 1 (EAdV-1), equine arteritis virus (EAV), reovirus 3 and parainfluenza virus type 3 (PIV3). Twelve foals were sampled until December 1996; the remaining 19 foals were lost from the study at various times prior to this date. RESULTS: The only viruses isolated were EHV 2 and EHV 5. EHV 2 was isolated from 155/157 PBL samples collected during the period of study and from 40/172 nasal swabs collected from 18 foals. All isolations from nasal swabs, except one, were made over a period of 2-4 months from January to April (Group A), March to April (Group B) or May, to July (Group C). EHV 5 was isolated from either PBL, nasal swabs, or both, from 15 foals on 32 occasions. All foals were positive for antibodies to EHV 1 or EHV 4, as tested by serum neutralisation (SN), on at least one sampling occasion and all but one were positive for EHV 1 antibodies measured by enzyme-linked immunosorbent assay (ELISA) on at least one sampling occasion. Recent EHV 1 infection was evident at least once during the period of study in 18/23 (78%) foals for which at least two samples were collected. SN antibodies to ERBV were evident in 19/23 (83%) foals on at least one sampling occasion and 15/23 foals showed evidence of seroconversion to ERBV Antibodies to ERAV were only detected in serum samples collected from foals in Group A and probably represented maternally-derived antibodies. Haemagglutination inhibition (HI) antibody titres greater than or equal to 1:10 to EAdV-1 were evident in 21/23 (91%) foals on at least one sampling occasion and 16/23 foals showed serological evidence of recent EAdV-1 infection. None of the 67 serum samples tested were positive for antibodies to EAV, reovirus 3 or PIV3. There was no clear association between infection with any of the viruses isolated or tested for and the presence of overt clinical signs of respiratory disease. CONCLUSIONS: There was serological and/or virological evidence that EHV-1, EHV-2, EHV-5, EAdV-1 and ERBV infections were present among foals in New Zealand. EHV-2 infection was first detected in foals as young as 3 months of age. The isolation of EHV-2 from nasal swabs preceded serological evidence of infection with other respiratory viruses, suggesting that EHV-2 may predispose foals to other viral infections.

Viruses associated with outbreaks of equine respiratory disease in New Zealand

AIM: To identify viruses associated with respiratory disease in young horses in New Zealand. METHODS: Nasal swabs and blood samples were collected from 45 foals or horses from five separate outbreaks of respiratory disease that occurred in New Zealand in 1996, and from 37 yearlings at the time of the annual yearling sales in January that same year. Virus isolation from nasal swabs and peripheral blood leukocytes (PBL) was undertaken and serum samples were tested for antibodies against equine herpesviruses (EHV-1, EHV-2, EHV-4 and EHV-5), equine rhinitis-A virus (ERAV), equine rhinitis-B virus (ERBV), equine adenovirus 1 (EAdV-1), equine arteritis virus (EAV), reovirus 3 and parainfluenza virus type 3 (PIV3). RESULTS: Viruses were isolated from 24/94 (26%) nasal swab samples and from 77/80 (96%) PBL samples collected from both healthy horses and horses showing clinical signs of respiratory disease. All isolates were identified as EHV-2, EHV-4, EHV-5 or untyped EHV Of the horses and foals tested, 59/82 (72%) were positive for EHV-1 and/or EHV-4 serum neutralising (SN) antibody on at least one sampling occasion, 52/82 (63%) for EHV-1-specific antibody tested by enzyme-linked immunosorbent assay (ELISA), 10/80 (13%) for ERAV SN antibody, 60/80 (75%) for ERBV SN antibody, and 42/80 (53%) for haemagglutination inhibition (HI) antibody to EAdV-1. None of the 64 serum samples tested were positive for antibodies to EAV, reovirus 3 or PIV3. Evidence of infection with all viruses tested was detected in both healthy horses and in horses showing clinical signs of respiratory disease. Recent EHV 2 infection was associated with the development of signs of respiratory disease among yearlings [relative risk (RR) = 2.67, 95% CI = 1.59-4.47, p = 0.0171]. CONCLUSIONS: Of the equine respiratory viruses detected in horses in New Zealand during this study, EHV 2 was most likely to be associated with respiratory disease. However, factors other than viral infection are probably important in the development of clinical signs of disease.

Predictor-corrector methods of Runge-Kutta type for stochastic differential equations

In this paper we construct predictor-corrector (PC) methods based on the trivial predictor and stochastic implicit Runge-Kutta (RK) correctors for solving stochastic differential equations. Using the colored rooted tree theory and stochastic B-series, the order condition theorem is derived for constructing stochastic RK methods based on PC implementations. We also present detailed order conditions of the PC methods using stochastic implicit RK correctors with strong global order 1.0 and 1.5. A two-stage implicit RK method with strong global order 1.0 and a four-stage implicit RK method with strong global order 1.5 used as the correctors are constructed in this paper. The mean-square stability properties and numerical results of the PC methods based on these two implicit RK correctors are reported.

An atomistic simulation of solid state sintering using Monte Carlo methods

This paper presents results on the simulation of the solid state sintering of copper wires using Monte Carlo techniques based on elements of lattice theory and cellular automata. The initial structure is superimposed onto a triangular, two-dimensional lattice, where each lattice site corresponds to either an atom or vacancy. The number of vacancies varies with the simulation temperature, while a cluster of vacancies is a pore. To simulate sintering, lattice sites are picked at random and reoriented in terms of an atomistic model governing mass transport. The probability that an atom has sufficient energy to jump to a vacant lattice site is related to the jump frequency, and hence the diffusion coefficient, while the probability that an atomic jump will be accepted is related to the change in energy of the system as a result of the jump, as determined by the change in the number of nearest neighbours. The jump frequency is also used to relate model time, measured in Monte Carlo Steps, to the actual sintering time. The model incorporates bulk, grain boundary and surface diffusion terms and includes vacancy annihilation on the grain boundaries. The predictions of the model were found to be consistent with experimental data, both in terms of the microstructural evolution and in terms of the sintering time. (C) 2002 Elsevier Science B.V. All rights reserved.

Estimation of total body water from bioelectrical impedance analysis in patients with pancreatic cancer - agreement between three methods of prediction

Introduction Bioelectrical impedance analysis (BIA) is a useful field measure to estimate total body water (TBW). No prediction formulae have been developed or validated against a reference method in patients with pancreatic cancer. The aim of this study was to assess the agreement between three prediction equations for the estimation of TBW in cachectic patients with pancreatic cancer. Methods Resistance was measured at frequencies of 50 and 200 kHz in 18 outpatients (10 males and eight females, age 70.2 +/- 11.8 years) with pancreatic cancer from two tertiary Australian hospitals. Three published prediction formulae were used to calculate TBW - TBWs developed in surgical patients, TBWca-uw and TBWca-nw developed in underweight and normal weight patients with end-stage cancer. Results There was no significant difference in the TBW estimated by the three prediction equations - TBWs 32.9 +/- 8.3 L, TBWca-nw 36.3 +/- 7.4 L, TBWca-uw 34.6 +/- 7.6 L. At a population level, there is agreement between prediction of TBW in patients with pancreatic cancer estimated from the three equations. The best combination of low bias and narrow limits of agreement was observed when TBW was estimated from the equation developed in the underweight cancer patients relative to the normal weight cancer patients. When no established BIA prediction equation exists, practitioners should utilize an equation developed in a population with similar critical characteristics such as diagnosis, weight loss, body mass index and/or age. Conclusions Further research is required to determine the accuracy of the BIA prediction technique against a reference method in patients with pancreatic cancer.

Control of Aedes vectors of dengue in three provinces of Vietnam by use of Mesocyclops (Copepoda) and community-based methods validated by entomologic, clinical, and serological surveillance

We describe remarkable success in controlling dengue vectors, Aedes aegypti (L.) and Aedes albopictus (Skuse), in 6 communes with 11,675 households and 49,647 people in the northern provinces of Haiphong, Hung Yen, and Nam Dinh in Vietnam. The communes were selected for high-frequency use of large outdoor concrete tanks and wells. These were found to be the source of 49.6-98.4% of Ae. aegypti larvae, which were amenable to treatment with local Mesocyclops, mainly M. woutersi Van der Velde, M. aspericornis (Daday) and M. thermocyclopoides Harada. Knowledge, attitude, and practice surveys were performed to determine whether the communities viewed dengue and dengue hemorrhagic fever as a serious health threat; to determine their knowledge of the etiology, attitudes, and practices regarding control methods including Mesocyclops; and to determine their receptivity to various information methods. On the basis of the knowledge, attitude, and practice data, the community-based dengue control program comprised a system of local leaders, health volunteer teachers, and schoolchildren, supported by health professionals. Recycling of discards for economic gain was enhanced, where appropriate, and this, plus 37 clean-up campaigns, removed small containers unsuitable for Mesocyclops treatment. A previously successful eradication at Phan Boi village (Hung Yen province) was extended to 7 other villages forming Di Su commune (1,750 households) in the current study. Complete control was also achieved in Nghia Hiep (Hung Yen province) and in Xuan Phong (Nam Dinh province); control efficacy was greater than or equal to 99.7% in the other 3 communes (Lac Vien in Haiphong, Nghia Dong, and Xuan Kien in Nam Dinh). Although tanks and wells were the key container types of Ae. aegypti productivity, discarded materials were the source of 51% of the standing crop of Ae. albopictus. Aedes albopictus larvae were eliminated from the 3 Nam Dinh communes, and 86-98% control was achieved in the other 3 communes. Variable dengue attack rates made the clinical and serological comparison of control and untreated communes problematic, but these data indicate that clinical surveillance by itself is inadequate to monitor dengue transmission.

Fitting genetic models to twin data with binary and ordered categorical responses: A comparison of structural equation modelling and Bayesian hierarchical models

We compare Bayesian methodology utilizing free-ware BUGS (Bayesian Inference Using Gibbs Sampling) with the traditional structural equation modelling approach based on another free-ware package, Mx. Dichotomous and ordinal (three category) twin data were simulated according to different additive genetic and common environment models for phenotypic variation. Practical issues are discussed in using Gibbs sampling as implemented by BUGS to fit subject-specific Bayesian generalized linear models, where the components of variation may be estimated directly. The simulation study (based on 2000 twin pairs) indicated that there is a consistent advantage in using the Bayesian method to detect a correct model under certain specifications of additive genetics and common environmental effects. For binary data, both methods had difficulty in detecting the correct model when the additive genetic effect was low (between 10 and 20%) or of moderate range (between 20 and 40%). Furthermore, neither method could adequately detect a correct model that included a modest common environmental effect (20%) even when the additive genetic effect was large (50%). Power was significantly improved with ordinal data for most scenarios, except for the case of low heritability under a true ACE model. We illustrate and compare both methods using data from 1239 twin pairs over the age of 50 years, who were registered with the Australian National Health and Medical Research Council Twin Registry (ATR) and presented symptoms associated with osteoarthritis occurring in joints of the hand.

Acquisition and loss of Porphyromonas gingivalis, Actinobacillus actinomycetemcomitans and Prevotella intermedia over a 5-year period: effect of a triclosan/copolymer dentifrice

Objectives: The present study describes the natural history of Porphyromonas gingivalis, Actinobacillus actinomycetemcomitans and Prevotella intermedia over a 5-year period and the effect of a triclosan/copolymer dentifrice on these organisms in a normal adult population. Material and Methods: Subgingival plaque samples were collected from 504 adult volunteers. Probing pocket depths (PPD) and relative attachment levels were measured using an automated probe. Participants were matched for disease status (CPI), plaque index, age and gender, and allocated to receive either a triclosan/copolymer or placebo dentifrice. Re-examination and subgingival plaque sampling was repeated after 1, 2, 3, 4 and 5 years. P. gingivalis, A. actinomycetemcomitans and P. intermedia were detected and quantitated using an enzyme linked immunosorbent assay. Logistic regression and generalised linear modelling were used to analyse the data. Results: This 5-year longitudinal study showed considerable volatility in acquisition and loss (below the level of detection) of all three organisms in this population. Relatively few subjects had these organisms on multiple occasions. While P. gingivalis was related to loss of attachment and to PPD greater than or equal to3.5 mm, there was no relationship between A. actinomycetemcomitans or P. intermedia and disease progression over the 5 years of the study. Smokers with P. gingivalis had more PPD greater than or equal to3.5 mm than smokers without this organism. There was no significant effect of the triclosan dentifrice on P. gingivalis or A. actinomycetemcomitans . Subjects using triclosan were more likely to have P. intermedia than those not using the dentifrice; however this did not translate into these subjects having higher levels of P. intermedia and its presence was uniform showing no signs of increasing over the course of the study. Conclusion: The present 5-year longitudinal study has shown the transient nature of colonisation with P. gingivalis , A. actinomycetemcomitans and P. intermedia in a normal adult population. The use of a triclosan-containing dentifrice did not lead to an overgrowth of these organisms. The clinical effect of the dentifrice would appear to be independent of its antimicrobial properties.

Impulse-response function of splanchnic circulation with model-independent constraints: theory and experimental validation

Modeling physiological processes using tracer kinetic methods requires knowledge of the time course of the tracer concentration in blood supplying the organ. For liver studies, however, inaccessibility of the portal vein makes direct measurement of the hepatic dual-input function impossible in humans. We want to develop a method to predict the portal venous time-activity curve from measurements of an arterial time-activity curve. An impulse-response function based on a continuous distribution of washout constants is developed and validated for the gut. Experiments with simultaneous blood sampling in aorta and portal vein were made in 13 anesthetized pigs following inhalation of intravascular [O-15] CO or injections of diffusible 3-O[ C-11] methylglucose (MG). The parameters of the impulse-response function have a physiological interpretation in terms of the distribution of washout constants and are mathematically equivalent to the mean transit time ( T) and standard deviation of transit times. The results include estimates of mean transit times from the aorta to the portal vein in pigs: (T) over bar = 0.35 +/- 0.05 min for CO and 1.7 +/- 0.1 min for MG. The prediction of the portal venous time-activity curve benefits from constraining the regression fits by parameters estimated independently. This is strong evidence for the physiological relevance of the impulse-response function, which includes asymptotically, and thereby justifies kinetically, a useful and simple power law. Similarity between our parameter estimates in pigs and parameter estimates in normal humans suggests that the proposed model can be adapted for use in humans.

Distribution of Chlamydia pneumoniae DNA in atherosclerotic carotid arteries: significance for sampling procedures

Despite extensive efforts to confirm a direct association between Chlamydia pneumoniae and atherosclerosis, different laboratories continue to report a large variability in detection rates. In this study, we analyzed multiple sections from atherosclerotic carotid arteries from 10 endartectomy patients to determine the location of C. pneumoniae DNA and the number of sections of the plaque required for analysis to obtain a 95% confidence of detecting the bacterium. A sensitive nested PCR assay detected C. pneumoniae DNA in all patients at one or more locations within the plaque. On average, 42% (ranging from 5 to 91%) of the sections from any single patient had C. pneumoniae DNA present. A patchy distribution of C. pneumoniae in the atherosclerotic lesions was observed, with no area of the carotid having significantly more C. pneumoniae DNA present. If a single random 30-mum-thick section was tested, there was only a 35.6 to 41.6% (95% confidence interval) chance of detecting C. pneumoniae DNA in a patient with carotid artery disease. A minimum of 15 sections would therefore be required to obtain a 95% chance of detecting all true positives. The low concentration and patchy distribution of C. pneumoniae DNA in atherosclerotic plaque appear to be among the reasons for inconsistency between laboratories in the results reported.