Genetic association studies under the population stratification, family pedigree and application to genome-wide association studies

This dissertation has three separate parts: the first part deals with the general pedigree association testing incorporating continuous covariates; the second part deals with the association tests under population stratification using the conditional likelihood tests; the third part deals with the genome-wide association studies based on the real rheumatoid arthritis (RA) disease data sets from Genetic Analysis Workshop 16 (GAW16) problem 1. Many statistical tests are developed to test the linkage and association using either case-control status or phenotype covariates for family data structure, separately. Those univariate analyses might not use all the information coming from the family members in practical studies. On the other hand, the human complex disease do not have a clear inheritance pattern, there might exist the gene interactions or act independently. In part I, the new proposed approach MPDT is focused on how to use both the case control information as well as the phenotype covariates. This approach can be applied to detect multiple marker effects. Based on the two existing popular statistics in family studies for case-control and quantitative traits respectively, the new approach could be used in the simple family structure data set as well as general pedigree structure. The combined statistics are calculated using the two statistics; A permutation procedure is applied for assessing the p-value with adjustment from the Bonferroni for the multiple markers. We use simulation studies to evaluate the type I error rates and the powers of the proposed approach. Our results show that the combined test using both case-control information and phenotype covariates not only has the correct type I error rates but also is more powerful than the other existing methods. For multiple marker interactions, our proposed method is also very powerful. Selective genotyping is an economical strategy in detecting and mapping quantitative trait loci in the genetic dissection of complex disease. When the samples arise from different ethnic groups or an admixture population, all the existing selective genotyping methods may result in spurious association due to different ancestry distributions. The problem can be more serious when the sample size is large, a general requirement to obtain sufficient power to detect modest genetic effects for most complex traits. In part II, I describe a useful strategy in selective genotyping while population stratification is present. Our procedure used a principal component based approach to eliminate any effect of population stratification. The paper evaluates the performance of our procedure using both simulated data from an early study data sets and also the HapMap data sets in a variety of population admixture models generated from empirical data. There are one binary trait and two continuous traits in the rheumatoid arthritis dataset of Problem 1 in the Genetic Analysis Workshop 16 (GAW16): RA status, AntiCCP and IgM. To allow multiple traits, we suggest a set of SNP-level F statistics by the concept of multiple-correlation to measure the genetic association between multiple trait values and SNP-specific genotypic scores and obtain their null distributions. Hereby, we perform 6 genome-wide association analyses using the novel one- and two-stage approaches which are based on single, double and triple traits. Incorporating all these 6 analyses, we successfully validate the SNPs which have been identified to be responsible for rheumatoid arthritis in the literature and detect more disease susceptibility SNPs for follow-up studies in the future. Except for chromosome 13 and 18, each of the others is found to harbour susceptible genetic regions for rheumatoid arthritis or related diseases, i.e., lupus erythematosus. This topic is discussed in part III.

Prevalence of cam and pincer-type deformities on hip MRI in an asymptomatic young Swiss female population: a cross-sectional study

OBJECTIVES Femoroacetabular impingement is proposed to cause early osteoarthritis (OA) in the non-dysplastic hip. We previously reported on the prevalence of femoral deformities in a young asymptomatic male population. The aim of this study was to determine the prevalence of both femoral and acetabular types of impingement in young females. METHODS We conducted a population-based cross-sectional study of asymptomatic young females. All participants completed a set of questionnaires and underwent clinical examination of the hip. A random sample was subsequently invited to obtain magnetic resonance images (MRI) of the hip. All MRIs were read for cam-type deformities, increased acetabular depths, labral lesions, and impingement pits. Prevalence estimates of cam-type deformities and increased acetabular depths were estimated, and relationships between deformities and signs of joint damage were examined using logistic regression models. RESULTS The study included 283 subjects, and 80 asymptomatic females with a mean age of 19.3 years attended MRI. Fifteen showed some evidence of cam-type deformities, but none were scored to be definite. The overall prevalence was therefore 0% [95% confidence interval (95% CI) 0-5%]. The prevalence of increased acetabular depth was 10% (95% CI 5-19). No association was found between increased acetabular depth and decreased internal rotation of the hip. Increased acetabular depth was not associated with signs of labral damage. CONCLUSIONS Definite cam-type deformities in women are rare compared to men, whereas the prevalence of increased acetabular depth is higher, suggesting that femoroacetabular impingement has different gender-related biomechanical mechanisms.

Lateral incisor agenesis, canine impaction and characteristics of supernumerary teeth in a South European male population

Objective: To assess the prevalence of lateral incisor agenesis impacted canines and supernumerary teeth in a young adult male population. Materials and Methods: The panoramic radiographs of 1745 military students (mean age: 18.6 ± 0.52 years) who attended the Center of Aviation Medicine of the Armed Forces of Greece during the period 1997-2011 were initially analyzed for lateral incisor agenesis by two observers. After exclusion of the known orthodontic cases, a subgroup of 1636 examinees (mean age: 18.6 ± 0.44 years) was evaluated for canine impaction and supernumerary teeth. Results: Twenty-eight missing lateral incisors were observed in 22 military students, indicating an incidence of 1.3% in the investigated population. No lateral incisor agenesis was detected in the mandibular arch. A prevalence rate of 0.8% was determined for canine impaction in the sample of young adults. The majority of impacted teeth (86.7%) were diagnosed in the maxillary arch. Thirty-five supernumerary teeth were observed in 24 examinees (prevalence rate: 1.5%). The ratio of supernumerary teeth located in the maxilla versus the mandible was 2.2:1. The most common type of supernumerary tooth was the upper distomolar. Conclusion: The prevalence of lateral incisor agenesis, canine impaction, and supernumerary teeth ranged from 0.8 to 1.5% in the sample of male Greek military students.

Detection of nonrandom association of alleles from the distribution of the number of heterozygous loci in a sample.

The distribution of the number of heterozygous loci in two randomly chosen gametes or in a random diploid zygote provides information regarding the nonrandom association of alleles among different genetic loci. Two alternative statistics may be employed for detection of nonrandom association of genes of different loci when observations are made on these distributions: observed variance of the number of heterozygous loci (s2k) and a goodness-of-fit criterion (X2) to contrast the observed distribution with that expected under the hypothesis of random association of genes. It is shown, by simulation, that s2k is statistically more efficient than X2 to detect a given extent of nonrandom association. Asymptotic normality of s2k is justified, and X2 is shown to follow a chi-square (chi 2) distribution with partial loss of degrees of freedom arising because of estimation of parameters from the marginal gene frequency data. Whenever direct evaluations of linkage disequilibrium values are possible, tests based on maximum likelihood estimators of linkage disequilibria require a smaller sample size (number of zygotes or gametes) to detect a given level of nonrandom association in comparison with that required if such tests are conducted on the basis of s2k. Summarization of multilocus genotype (or haplotype) data, into the different number of heterozygous loci classes, thus, amounts to appreciable loss of information.

Population amalgamation and genetic variation: observations on artificially agglomerated tribal populations of Central and South America.

The interpretation of data on genetic variation with regard to the relative roles of different evolutionary factors that produce and maintain genetic variation depends critically on our assumptions concerning effective population size and the level of migration between neighboring populations. In humans, recent population growth and movements of specific ethnic groups across wide geographic areas mean that any theory based on assumptions of constant population size and absence of substructure is generally untenable. We examine the effects of population subdivision on the pattern of protein genetic variation in a total sample drawn from an artificial agglomerate of 12 tribal populations of Central and South America, analyzing the pooled sample as though it were a single population. Several striking findings emerge. (1) Mean heterozygosity is not sensitive to agglomeration, but the number of different alleles (allele count) is inflated, relative to neutral mutation/drift/equilibrium expectation. (2) The inflation is most serious for rare alleles, especially those which originally occurred as tribally restricted "private" polymorphisms. (3) The degree of inflation is an increasing function of both the number of populations encompassed by the sample and of the genetic divergence among them. (4) Treating an agglomerated population as though it were a panmictic unit of long standing can lead to serious biases in estimates of mutation rates, selection pressures, and effective population sizes. Current DNA studies indicate the presence of numerous genetic variants in human populations. The findings and conclusions of this paper are all fully applicable to the study of genetic variation at the DNA level as well.

Comparison of pain-resilient working individuals to population-based case controls with/without momentary low back pain

Background: Few studies have examined the 20% of individuals who never experience an episode of low back pain (LBP). To date, no investigation has been undertaken that examines a group who claim to have never experienced LBP in their lifetime in comparison to two population-based case–control groups with and without momentary LBP. This study investigates whether LBP-resilient workers between 50 and 65 years had better general health, demonstrated more positive health behaviour and were better able to achieve routine activities compared with both case–control groups. Methods: Forty-two LBP-resilient participants completed the same pain assessment questionnaire as a population-based LBP sample from a nationwide, large-scale cross-sectional survey in Switzerland. The LBP-resilient participants were pairwise compared to the propensity score-matched case controls by exploring differences in demographic and work characteristics, and by calculating odds ratios (ORs) and effect sizes. A discriminant analysis explored group differences, while the multiple logistic regression analysis specified single indicators which accounted for group differences. Results: LBP-resilient participants were healthier than the case controls with momentary LBP and achieved routine activities more easily. Compared to controls without momentary LBP, LBP-resilient participants had a higher vitality, a lower workload, a healthier attitude towards health and behaved more healthily by drinking less alcohol. Conclusions: By demonstrating a difference between LBP-resilient participants and controls without momentary LBP, the question that arises is what additional knowledge can be attained. Three underlying traits seem to be relevant about LBP-resilient participants: personality, favourable work conditions and subjective attitudes/attributions towards health. These rationales have to be considered with respect to LBP prevention.

Reference Values and Factors Associated With Renal Resistive Index in a Family-Based Population Study

Increased renal resistive index (RRI) has been recently associated with target organ damage and cardiovascular or renal outcomes in patients with hypertension and diabetes mellitus. However, reference values in the general population and information on familial aggregation are largely lacking. We determined the distribution of RRI, associated factors, and heritability in a population-based study. Families of European ancestry were randomly selected in 3 Swiss cities. Anthropometric parameters and cardiovascular risk factors were assessed. A renal Doppler ultrasound was performed, and RRI was measured in 3 segmental arteries of both kidneys. We used multilevel linear regression analysis to explore the factors associated with RRI, adjusting for center and family relationships. Sex-specific reference values for RRI were generated according to age. Heritability was estimated by variance components using the ASSOC program (SAGE software). Four hundred women (mean age±SD, 44.9±16.7 years) and 326 men (42.1±16.8 years) with normal renal ultrasound had mean RRI of 0.64±0.05 and 0.62±0.05, respectively (P<0.001). In multivariable analyses, RRI was positively associated with female sex, age, systolic blood pressure, and body mass index. We observed an inverse correlation with diastolic blood pressure and heart rate. Age had a nonlinear association with RRI. We found no independent association of RRI with diabetes mellitus, hypertension treatment, smoking, cholesterol levels, or estimated glomerular filtration rate. The adjusted heritability estimate was 42±8% (P<0.001). In a population-based sample with normal renal ultrasound, RRI normal values depend on sex, age, blood pressure, heart rate, and body mass index. The significant heritability of RRI suggests that genes influence this phenotype.

The Role of Cell Sterilization in Population Based Studies of Radiogenic Second Cancers Following Radiation Therapy

Advances in radiotherapy have generated increased interest in comparative studies of treatment techniques and their effectiveness. In this respect, pediatric patients are of specific interest because of their sensitivity to radiation induced second cancers. However, due to the rarity of childhood cancers and the long latency of second cancers, large sample sizes are unavailable for the epidemiological study of contemporary radiotherapy treatments. Additionally, when specific treatments are considered, such as proton therapy, sample sizes are further reduced due to the rareness of such treatments. We propose a method to improve statistical power in micro clinical trials. Specifically, we use a more biologically relevant quantity, cancer equivalent dose (DCE), to estimate risk instead of mean absorbed dose (DMA). Our objective was to demonstrate that when DCE is used fewer subjects are needed for clinical trials. Thus, we compared the impact of DCE vs. DMA on sample size in a virtual clinical trial that estimated risk for second cancer (SC) in the thyroid following craniospinal irradiation (CSI) of pediatric patients using protons vs. photons. Dose reconstruction, risk models, and statistical analysis were used to evaluate SC risk from therapeutic and stray radiation from CSI for 18 patients. Absorbed dose was calculated in two ways: with (1) traditional DMA and (2) with DCE. DCE and DMA values were used to estimate relative risk of SC incidence (RRCE and RRMA, respectively) after proton vs. photon CSI. Ratios of RR for proton vs. photon CSI (RRRCE and RRRMA) were then used in comparative estimations of sample size to determine the minimal number of patients needed to maintain 80% statistical power when using DCE vs. DMA. For all patients, we found that protons substantially reduced the risk of developing a second thyroid cancer when compared to photon therapy. Mean RRR values were 0.052±0.014 and 0.087±0.021 for RRRMA and RRRCE, respectively. However, we did not find that use of DCE reduced the number of patents needed for acceptable statistical power (i.e, 80%). In fact, when considerations were made for RRR values that met equipoise requirements and the need for descriptive statistics, the minimum number of patients needed for a micro-clinical trial increased from 17 using DMA to 37 using DCE. Subsequent analyses revealed that for our sample, the most influential factor in determining variations in sample size was the experimental standard deviation of estimates for RRR across the patient sample. Additionally, because the relative uncertainty in dose from proton CSI was so much larger (on the order of 2000 times larger) than the other uncertainty terms, it dominated the uncertainty in RRR. Thus, we found that use of corrections for cell sterilization, in the form of DCE, may be an important and underappreciated consideration in the design of clinical trials and radio-epidemiological studies. In addition, the accurate application of cell sterilization to thyroid dose was sensitive to variations in absorbed dose, especially for proton CSI, which may stem from errors in patient positioning, range calculation, and other aspects of treatment planning and delivery.

The impact of cancer on the population of Salvador-Bahia, Brazil

The impact of cancer on the population of Salvador-Bahia, Brazil was studied using mortality data available from the Brazilian National Bureau of Vital Statistics. Average annual site, age, and gender specific and adjusted cancer mortality rates were determined for the years 1977-83 and contrasted with United States cancer mortality rates for the year of 1977. The accuracy of the cancer mortality rates generated by this research was determined by comparing the underlying causes of death as coded on death certificates to pathology reports and to hospital diagnosis of a sample of 966 deaths occurring in Salvador during the year of 1983. To further explore the information available on the death certificate, a population based decedent case control study was used to determine the relationship between type of occupation (proxy for exposure) and mortality by cancer sites known to be occupationally related.^ The rates in Salvador for cancer of the stomach, oral cavity, and biliary passages are, on average, two fold higher than the U.S. rates.^ The death certificate was found to be accurate for 65 percent of the 485 cancer deaths studied. Thirty five histologically confirmed cancer deaths were found in a random sample of 481 deaths from other causes. This means that, approximately 700 more deaths may be lost among the remainder 10,073 death certificates stating a cause other than cancer.^ In addition, despite the known limitations of decedent case-control studies, cancers of the oral cavity OR = 2.44, CI = 1.17-5.09, stomach OR = 2.31, CI = 1.18-4.52, liver OR = 4.06, CI = 1.27-12.99, bladder OR = 6.77, CI = 1.5-30.67, and lymphoma OR = 2.55, CI = 1.04-6.25 had elevated point estimates, for different age strata indicating an association between these cancers and occupations that led to exposure to petroleum and its derivates. ^

Evaluation of androgenicity, abdominal adiposity, blood pressure and types of ovulatory patterns in a sample of women with menstrual irregularities

Objective: To determine the prevalence of and the relationships between the degree and source of hyperandrogenemia, ovulatory patterns and cardiovascular disease risk indicators (blood pressure, indices or amount of obesity and fat distribution) in women with menstrual irregularities seen at endocrinologists' clinic. Design: A cross-sectional study design. Participants: A sample of 159 women with menstrual irregularities, aged 15-44, seen at endocrinologists' clinic. Main Outcome Measures: androgen levels, body mass index (BMI), waist-hip ratio (WHR), systolic and diastolic blood pressure (SBP & DBP), source of androgens, ovulatory activity. Results: The prevalence of hyperandrogenemia was 54.7% in this study sample. As expected, women with acne or hirsutism had an odds ratio 12.5 (95%CI = 5.2-25.5) times and 36 (95%CI = 12.9-99.5) times more likely to have hyperandrogenemia than those without acne or hirsutism. The main findings of this study were the following: Hyperandrogenemic women were more likely to have oligomenorrheic cycles (OR = 3.8, 95%CI = 1.5-9.9), anovulatory cycles (OR = 6.6, 95%CI = 2.8-15.4), general obesity (BMI $\ge$ 27) (OR = 6.8, 95%CI = 2.2-27.2) and central obesity (WHR $\ge$ 127) (OR = 14.5, 95%CI = 6.1-38.7) than euandrogenemic women. Hyperandrogenemic women with non-suppressible androgens had a higher mean BMI (29.3 $\pm$ 8.9) than those with suppressible androgens (27.9 $\pm$ 7.9); the converse was true for abdominal adiposity (WHR). Hyperandrogenemic women had a 2.4 odds ratio (95%CI = 1.0-6.2) for an elevated SBP and a 2.7 odds ratio (95%CI = 0.8-8.8) for elevated DBP. When age differences were accounted for, this relationship was strengthened and further strengthened when sources of androgens were controlled. When the differences in BMI were controlled, the odds ratio for elevated SBP in hyperandrogenemic women increased to 8.8 (95%CI = 1.1-69.9). When the age, the source of androgens, the amount of obesity and the type of obesity were controlled, hyperandrogenemic women had 13.5 (95%CI = 1.1-158.9) odds ratio for elevated SBP. Conclusions: In this study population, the presence of menstrual irregularities are highly predictive for the presence of elevated androgens. Women with elevated androgens have a high risk for obesity, more specifically for central obesity. The androgenemic status is an independent predictor of blood pressure elevation. It is probable that in the general population, the presence of menstrual irregularities are predictive of hyperandrogenemia. There is a great need for a population study of the prevalence of hyperandrogenemia and for longitudinal studies in hyperandrogenemic women (adrenarche to menopause) to investigate the evolution of these relationships. ^

Twenty-three per cent of the Swiss adult population are using complementary medicine

Aim of the study: This study investigated the use among the Swiss adult population and the regional dissemination in Switzerland of various methods of complementary medicine (CM). It focused on CM methods that required visiting a physician or therapist and excluded e.g. over-the-counter drugs. Data and Methods: Data of the Swiss Health Survey 2007 were obtained from the Swiss Federal Statistical Office. This survey is performed every 5 years in a sample and is representative of the Swiss resident population from the age of 15 on. It consists of a telephone interview followed by a written questionnaire (2007: 18'760 and 14'432 respondents, respectively) and includes questions about people's state of health, general living conditions, lifestyle, health insurance and usage of health services. Users and non-users of CM were compared using logistic regression models. Results: 23.0 % of the Swiss adult population (women: 30.5 %, men: 15.2 %) used CM during the 12 months before the survey. Homeopathy (6.4 %), osteopathy (5.4 %) and acupuncture (4.9 %) were the most popular methods. The average number of treatments within 12 months for these three methods was 3.1 ± 3.6, 3.5 ± 3.3 and 6.6 ± 5.8, respectively. For treatments with homeopathy and acupuncture, medical practitioners were more commonly consulted than non-medical practitioners, for treatments with osteopathy no difference was found. By means of logistic regression, CM users and non-users were compared. There were significant differences in the use of CM between genders, age groups, levels of education and areas of living. Women, people aged 25 to 64 years, and people with higher levels of education used CM more commonly than men, people below 25 or above 64 years of age, or those with poorer education. Lake Geneva region and central Switzerland had a higher proportion of CM users than the other regions. Discussion: Almost one fourth of the Swiss adult population had used CM within 12 months before the survey. User profiles were comparable to those in other countries. Despite a generally lower self-perceived health status, elderly people were less likely to use CM. Reference: Klein SD, Frei-Erb M, Wolf U. Usage of complementary medicine across Switzerland. Results of the Swiss Health Survey 2007. Swiss Med Wkly. 2012;142:w13666.

Twenty-three per cent of the Swiss adult population are using complementary medicine

Aim of the study Various forms of complementary medicine (CM) play an important role in the Swiss health care system, they are appreciated by a majority of the population and mostly used complementarily rather than alternatively to conventional medicine. This study investigates, how many people in Switzerland are actually being treated with CM, and what the most popular methods of CM are. Data Data of the Swiss Health Survey 2007 were obtained from the Swiss Federal Statistical Office. This survey is performed every 5 years amongst a sample of the Swiss resident population above 15 years of age. It consists of a telephone interview followed by a written questionnaire (2007: 18'760 and 14'432 respondents, respectively) and includes questions about people's state of health, general living conditions, lifestyle, health insurance and usage of health services. Results 23.0% of the Swiss adult population (women: 30.5%, men: 15.2%) used CM during the 12 months before the survey. Homeopathy (6.4%), osteopathy (5.4%) and acupuncture (4.9%) were the most popular methods. The average number of treatments within 12 months for these three methods was 3.1, 3.5 and 6.6, respectively. For treatments with homeopathy and acupuncture, medical practitioners were more commonly consulted than non-medical practitioners, for treatments with osteopathy no difference was found. By means of logistic regression, CM users and non-users were compared. There were significant differences in the use of CM between genders, age groups, levels of education and areas of living. Women, people aged 25 to 64 years, and people with higher education used CM more commonly than men, people below 25 or above 64 years of age, or those with poorer education. Lake Geneva region and central Switzerland had a higher proportion of CM users than the other regions. Discussion While 2 years ago, 67.0% of the Swiss population approved a referendum in favour of CM, we find that 23.0% are in fact using it. Current political discussions focus on effectiveness, cost effectiveness and suitability of CM to decide which methods should be permanently covered by the basic health insurance.

Associations between mood, anxiety or substance use disorders and inflammatory markers after adjustment for multiple covariates in a population-based study.

Inflammation is one possible mechanism underlying the associations between mental disorders and cardiovascular diseases (CVD). However, studies on mental disorders and inflammation have yielded inconsistent results and the majority did not adjust for potential confounding factors. We examined the associations of several pro-inflammatory cytokines (IL-1β, IL-6 and TNF-α) and high sensitive C-reactive protein (hsCRP) with lifetime and current mood, anxiety and substance use disorders (SUD), while adjusting for multiple covariates. The sample included 3719 subjects, randomly selected from the general population, who underwent thorough somatic and psychiatric evaluations. Psychiatric diagnoses were made with a semi-structured interview. Major depressive disorder was subtyped into "atypical", "melancholic", "combined atypical-melancholic" and "unspecified". Associations between inflammatory markers and psychiatric diagnoses were assessed using multiple linear and logistic regression models. Lifetime bipolar disorders and atypical depression were associated with increased levels of hsCRP, but not after multivariate adjustment. After multivariate adjustment, SUD remained associated with increased hsCRP levels in men (β = 0.13 (95% CI: 0.03,0.23)) but not in women. After multivariate adjustment, lifetime combined and unspecified depression were associated with decreased levels of IL-6 (β = -0.27 (-0.51,-0.02); β = -0.19 (-0.34,-0.05), respectively) and TNF-α (β = -0.16 (-0.30,-0.01); β = -0.10 (-0.19,-0.02), respectively), whereas current combined and unspecified depression were associated with decreased levels of hsCRP (β = -0.20 (-0.39,-0.02); β = -0.12 (-0.24,-0.01), respectively). Our data suggest that the significant associations between increased hsCRP levels and mood disorders are mainly attributable to the effects of comorbid disorders, medication as well as behavioral and physical CVRFs.

Age matters in the prevalence and clinical significance of ultra-high-risk for psychosis symptoms and criteria in the general population: Findings from the BEAR and BEARS-kid studies

Early detection of psychosis is an important topic in psychiatry. Yet, there is limited information on the prevalence and clinical significance of high-risk symptoms in children and adolescents as compared to adults. We examined ultra-high-risk (UHR) symptoms and criteria in a sample of individuals aged 8-40 years from the general population of Canton Bern, Switzerland, enrolled from June 2011 to May 2014. The current presence of attenuated psychotic symptoms (APS) and brief intermittent psychotic symptoms (BLIPS) and the fulfillment of onset/worsening and frequency requirements for these symptoms in UHR criteria were assessed using the Structured Interview for Psychosis Risk Syndromes. Additionally, perceptive and non-perceptive APS were differentiated. Psychosocial functioning and current non-psychotic DSM-IV axis I disorders were also surveyed. Well-trained psychologists performed assessments. Altogether, 9.9% of subjects reported APS and none BLIPS, and 1.3% met all the UHR requirements for APS. APS were related to more current axis I disorders and impaired psychosocial functioning, indicating some clinical significance. A strong age effect was detected around age 16: compared to older individuals, 8-15-year olds reported more perceptive APS, that is, unusual perceptual experiences and attenuated hallucinations. Perceptive APS were generally less related to functional impairment, regardless of age. Conversely, non-perceptive APS were related to low functioning, although this relationship was weaker in those below age 16. Future studies should address the differential effects of perceptive and non-perceptive APS, and their interaction with age, also in terms of conversion to psychosis.

Unresolved Grief in a National Sample of Bereaved Parents: Impaired Mental and Physical Health 4 to 9 Years Later

PURPOSE: To assess unresolved parental grief, the associated long-term impact on mental and physical health, and health service use. PATIENTS AND METHODS: This anonymous, mail-in questionnaire study was performed as a population-based investigation in Sweden between August 2001 and October 2001. Four hundred forty-nine parents who lost a child as a result of cancer 4 to 9 years earlier completed the survey (response rate, 80%). One hundred ninety-one (43%) of the bereaved parents were fathers, and 251 (56%) were mothers. Bereaved parents were asked whether or not, and to what extent, they had worked through their grief. They were also asked about their physical and psychological well-being. For outcomes of interest, we report relative risk (RR) with 95% CIs as well as unadjusted odds ratios and adjusted odds ratios. RESULTS: Parents with unresolved grief reported significantly worsening psychological health (fathers: RR, 3.6; 95% CI, 2.0 to 6.4; mothers: RR, 2.9; 95% CI, 1.9 to 4.4) and physical health (fathers: RR, 2.8; 95% CI, 1.8 to 4.4; mothers: RR, 2.3; 95% CI, 1.6 to 3.3) compared with those who had worked through their grief. Fathers with unresolved grief also displayed a significantly higher risk of sleep difficulties (RR, 6.7; 95% CI, 2.5 to 17.8). Mothers, however, reported increased visits with physicians during the previous 5 years (RR, 1.7; 95% CI, 1.1 to 2.6) as well as a greater likelihood of taking sick leave when they had not worked through their grief (RR, 2.1; 95% CI, 1.2 to 3.5). CONCLUSION: Parents who have not worked through their grief are at increased risk of long-term mental and physical morbidity, increased health service use, and increased sick leave.