Association of CD247 (CD3ζ) gene polymorphisms with T1D and AITD in the population of northern Sweden.

BACKGROUND: T1D and AITD are autoimmune disorders commonly occurring in the same family and even in the same individual. The genetic contribution to these disorders is complex making uncovering of susceptibility genes very challenging. The general aim of this study was to identify loci and genes contributing to T1D/AITD susceptibility. Our strategy was to perform linkage and association studies in the relatively genetically homogenous population of northern Sweden. We performed a GWLS to find genomic regions linked to T1D/AITD in families from northern Sweden and we performed an association study in the families to test for association between T1D/AITD and variants in previously published candidate genes as well as a novel candidate gene, CD247. METHODS: DNA prepared from 459 individuals was used to perform a linkage and an association study. The ABI PRISM Linkage Mapping Set v2.5MD10 was employed for an initial 10-cM GWLS, and additional markers were added for fine mapping. Merlin was used for linkage calculations. For the association analysis, a GoldenGate Custom Panel from Illumina containing 79 SNPs of interest was used and FBAT was used for association calculations. RESULTS: Our study revealed linkage to two previously identified chromosomal regions, 4q25 and 6p22, as well as to a novel chromosomal region, 1q23. The association study replicated association to PTPN22, HLA-DRB1, INS, IFIH1, CTLA4 and C12orf30. Evidence in favor of association was also found for SNPs in the novel susceptibility gene CD247. CONCLUSIONS: Several risk loci for T1D/AITD identified in published association studies were replicated in a family material, of modest size, from northern Sweden. This provides evidence that these loci confer disease susceptibility in this population and emphasizes that small to intermediate sized family studies in this population can be used in a cost-effective manner for the search of genes involved in complex diseases. The linkage study revealed a chromosomal region in which a novel T1D/AITD susceptibility gene, CD247, is located. The association study showed association between T1D/AITD and several variants in this gene. These results suggests that common susceptibility genes act in concert with variants of CD247 to generate genetic risk for T1D/AITD in this population.

The phylogeographic pattern of Francisella tularensis in Sweden indicates a Scandinavian origin of Eurosiberian tularaemia.

Previous studies of the causative agent of tularaemia, Francisella tularensis have identified phylogeographic patterns suggestive of environmental maintenance reservoirs. To investigate the phylogeography of tularaemia in Sweden, we selected 163 clinical isolates obtained during 1995-2009 in 10 counties and sequenced one isolate's genome to identify new genetic markers. An improved typing scheme based on two indels and nine SNPs was developed using hydrolysis or TaqMan MGB probe assays. The results showed that much of the known global genetic diversity of F. tularensis subsp. holarctica is present in Sweden. Thirteen of the 163 isolates belonged to a new genetic group that is basal to all other known members of the major genetic clade B.I, which is spread across the Eurosiberian region. One hundred and twenty-five of the 163 Swedish isolates belonged to B.I, but individual clades' frequencies differed from county to county (P < 0.001). Subsequent analyses revealed a correlation between genotype variation over time and recurrent outbreaks at specific places, supporting the 'maintenance reservoir' environmental maintenance hypothesis. Most importantly, the findings reveal the presence of diverse source populations of F. tularensis subsp. holarctica in Sweden and suggest a historical spread of the disease from Scandinavia to other parts of Eurosiberia.

Genome characterisation of the genus Francisella reveals insight into similar evolutionary paths in pathogens of mammals and fish.

BACKGROUND: Prior to this study, relatively few strains of Francisella had been genome-sequenced. Previously published Francisella genome sequences were largely restricted to the zoonotic agent F. tularensis. Only limited data were available for other members of the Francisella genus, including F. philomiragia, an opportunistic pathogen of humans, F. noatunensis, a serious pathogen of farmed fish, and other less well described endosymbiotic species. RESULTS: We determined the phylogenetic relationships of all known Francisella species, including some for which the phylogenetic positions were previously uncertain. The genus Francisella could be divided into two main genetic clades: one included F. tularensis, F. novicida, F. hispaniensis and Wolbachia persica, and another included F. philomiragia and F. noatunensis.Some Francisella species were found to have significant recombination frequencies. However, the fish pathogen F. noatunensis subsp. noatunensis was an exception due to it exhibiting a highly clonal population structure similar to the human pathogen F. tularensis. CONCLUSIONS: The genus Francisella can be divided into two main genetic clades occupying both terrestrial and marine habitats. However, our analyses suggest that the ancestral Francisella species originated in a marine habitat. The observed genome to genome variation in gene content and IS elements of different species supports the view that similar evolutionary paths of host adaptation developed independently in F. tularensis (infecting mammals) and F. noatunensis subsp. noatunensis (infecting fish).

Increased knowledge of Francisella genus diversity highlights the benefits of optimised DNA-based assays.

BACKGROUND: Recent advances in sequencing technologies offer promising tools for generating large numbers of genomes, larger typing databases and improved mapping of environmental bacterial diversity. However, DNA-based methods for the detection of Francisella were developed with limited knowledge about genetic diversity. This, together with the high sequence identity between several Francisella species, means there is a high risk of false identification and detection of the highly virulent pathogen Francisella tularensis. Moreover, phylogenetic reconstructions using single or limited numbers of marker sequences often result in incorrect tree topologies and inferred evolutionary distances. The recent growth in publicly accessible whole-genome sequences now allows evaluation of published genetic markers to determine optimal combinations of markers that minimise both time and laboratory costs. RESULTS: In the present study, we evaluated 38 previously published DNA markers and the corresponding PCR primers against 42 genomes representing the currently known diversity of the genus Francisella. The results highlight that PCR assays for Francisella tularensis are often complicated by low specificity, resulting in a high probability of false positives. A method to select a set of one to seven markers for obtaining optimal phylogenetic resolution or diagnostic accuracy is presented. CONCLUSIONS: Current multiple-locus sequence-typing systems and detection assays of Francisella, could be improved by redesigning some of the primers and reselecting typing markers. The use of only a few optimally selected sequence-typing markers allows construction of phylogenetic topologies with almost the same accuracy as topologies based on whole-genome sequences.

Extended-release quetiapine fumarate (quetiapine XR) as adjunctive therapy in major depressive disorder (MDD) in patients with an inadequate response to ongoing antidepressant treatment: a multicentre, randomized, double-blind, placebo-controlled study.

This study evaluated once-daily extended-release quetiapine fumarate (quetiapine XR) as adjunctive therapy in patients with major depressive disorder (MDD) with inadequate response to ongoing antidepressant treatment. In this 8-wk (6-wk active treatment/2-wk post-treatment drug-discontinuation/follow-up), multicentre, double-blind, placebo-controlled, Phase III study, 446 patients were randomized to quetiapine XR 150 mg/d, 300 mg/d, or placebo adjunct to ongoing antidepressant treatment. The primary endpoint was the change from randomization to week 6 in Montgomery-Asberg Depression Rating Scale (MADRS) total score. At week 6, MADRS total scores significantly improved with quetiapine XR 300 mg/d vs. placebo (-14.7 vs. -11.7, p<0.01). Quetiapine XR 300 mg/d showed significant improvements vs. placebo for: MADRS total score from week 1 onwards; MADRS response [(> or = 50% total score reduction) 58.9% vs. 46.2%, p<0.05] and remission [(total score < or = 8) 42.5% vs. 24.5%, p<0.01] rates; Hamilton Depression Rating Scale (HAMD) (-13.53 vs. -10.80, p<0.01) and Clinical Global Impression-Severity of illness (CGI-S) change (-1.52 vs. -1.23, p<0.05) at week 6. For quetiapine XR 150 mg/d, improvements were not significantly different vs. placebo, except for MADRS (weeks 1 and 2) and HAMD (week 6) total scores. Withdrawal rates due to adverse events (AEs) were: quetiapine XR 150 mg/d 11.5%, 300 mg/d 19.5%, and placebo 0.7%. The most common AEs (>10%) with quetiapine XR were dry mouth, somnolence, sedation, dizziness, constipation, nausea, insomnia, headache, and fatigue. In this study, quetiapine XR 300 mg/d as adjunctive therapy in patients with MDD with an inadequate response to ongoing antidepressant treatment was effective at week 6. However, the difference from placebo for quetiapine XR 150 mg/d at week 6 was not statistically significant. Both doses studied (150 and 300 mg/d) were effective at week 1 and generally well tolerated.

Biological amplification of low frequency mutations unravels laboratory culture history of the bio-threat agent Francisella tularensis.

Challenges of investigating a suspected bio attack include establishing if microorganisms have been cultured to produce attack material and to identify their source. Addressing both issues, we have investigated genetic variations that emerge during laboratory culturing of the bacterial pathogen Francisella tularensis. Key aims were to identify genetic variations that are characteristic of laboratory culturing and explore the possibility of using biological amplification to identify genetic variation present at exceedingly low frequencies in a source sample. We used parallel serial passage experiments and high-throughput sequencing of F. tularensis to explore the genetic variation. We found that during early laboratory culture passages of F. tularensis, gene duplications emerged in the pathogen genome followed by single-nucleotide polymorphisms in genes for bacterial capsule synthesis. Based on a biological enrichment scheme and the use of high-throughput sequencing, we identified genetic variation that likely pre-existed in a source sample. The results support that capsule synthesis gene mutations are common during laboratory culture, and that a biological amplification strategy is useful for linking a F. tularensis sample to a specific laboratory variant among many highly similar variants.

A likelihood ratio-based approach for improved source attribution in microbiological forensic investigations.

A common objective in microbial forensic investigations is to identify the origin of a recovered pathogenic bacterium by DNA sequencing. However, there is currently no consensus about how degrees of belief in such origin hypotheses should be quantified, interpreted, and communicated to wider audiences. To fill this gap, we have developed a concept based on calculating probabilistic evidential values for microbial forensic hypotheses. The likelihood-ratio method underpinning this concept is widely used in other forensic fields, such as human DNA matching, where results are readily interpretable and have been successfully communicated in juridical hearings. The concept was applied to two case scenarios of interest in microbial forensics: (1) identifying source cultures among series of very similar cultures generated by parallel serial passage of the Tier 1 pathogen Francisella tularensis, and (2) finding the production facilities of strains isolated in a real disease outbreak caused by the human pathogen Listeria monocytogenes. Evidence values for the studied hypotheses were computed based on signatures derived from whole genome sequencing data, including deep-sequenced low-frequency variants and structural variants such as duplications and deletions acquired during serial passages. In the F. tularensis case study, we were able to correctly assign fictive evidence samples to the correct culture batches of origin on the basis of structural variant data. By setting up relevant hypotheses and using data on cultivated batch sources to define the reference populations under each hypothesis, evidential values could be calculated. The results show that extremely similar strains can be separated on the basis of amplified mutational patterns identified by high-throughput sequencing. In the L. monocytogenes scenario, analyses of whole genome sequence data conclusively assigned the clinical samples to specific sources of origin, and conclusions were formulated to facilitate communication of the findings. Taken together, these findings demonstrate the potential of using bacterial whole genome sequencing data, including data on both low frequency SNP signatures and structural variants, to calculate evidence values that facilitate interpretation and communication of the results. The concept could be applied in diverse scenarios, including both epidemiological and forensic source tracking of bacterial infectious disease outbreaks.

Galleria mellonella Reveals Niche Differences Between Highly Pathogenic and Closely Related Strains of Francisella spp.

Francisella tularensis, a highly virulent bacteria that causes the zoonotic disease tularemia, is considered a potential agent of biological warfare and bioterrorism. Although the host range for several species within the Francisella is known, little is known about the natural reservoirs of various Francisella species. The lack of knowledge regarding the environmental fates of these pathogens greatly reduces the possibilities for microbial risk assessments. The greater wax moth (Galleria mellonella) is an insect of the order Lepidoptera that has been used as an alternative model to study microbial infection during recent years. The aim of this study was to evaluate G. mellonella as a model system for studies of human pathogenic and closely related opportunistic and non-pathogenic strains within the Francisella genus. The employed G. mellonella larvae model demonstrated differences in lethality between human pathogenic and human non-pathogenic or opportunistic Francisella species. The F. novicida, F. hispaniensis and F. philomiragia strains were significantly more virulent in the G. mellonella model than the strains of human pathogens F. t. holarctica and F. t. tularensis. Our data show that G. mellonella is a possible in vivo model of insect immunity for studies of both opportunistic and virulent lineages of Francisella spp., that produces inverse results regarding lethality in G. mellonella and incapacitating disease in humans. The results provide insight into the potential host specificity of F. tularensis and closely related members of the same genus, thus increasing our present understanding of Francisella spp. ecology.

Accounting for Bacterial Overlap Between Raw Water Communities and Contaminating Sources Improves the Accuracy of Signature-Based Microbial Source Tracking.

Microbial source tracking (MST) analysis is essential to identifying and mitigating the fecal pollution of water resources. The signature-based MST method uses a library of sequences to identify contaminants based on operational taxonomic units (OTUs) that are unique to a certain source. However, no clear guidelines for how to incorporate OTU overlap or natural variation in the raw water bacterial community into MST analyses exist. We investigated how the inclusion of bacterial overlap between sources in the library affects source prediction accuracy. To achieve this, large-scale sampling - including feces from seven species, raw sewage, and raw water samples from water treatment plants - was followed by 16S rRNA amplicon sequencing. The MST library was defined using three settings: (i) no raw water communities represented; (ii) raw water communities selected through clustering analysis; and (iii) local water communities collected across consecutive years. The results suggest that incorporating either the local background or representative bacterial composition improves MST analyses, as the results were positively correlated to measured levels of fecal indicator bacteria and the accuracy at which OTUs were assigned to the correct contamination source increased fourfold. Using the proportion of OTUs with high source origin probability, underpinning a contaminating signal, is a solid foundation in a framework for further deciphering and comparing contaminating signals derived in signature-based MST approaches. In conclusion, incorporating background bacterial composition of water in MST can improve mitigation efforts for minimizing the spread of pathogenic and antibiotic resistant bacteria into essential freshwater resources.

Genome-wide linkage scan of common stroke in families from northern Sweden.

BACKGROUND AND PURPOSE: Taking advantage of low genetic variations in northern Sweden, we performed a genome-wide linkage scan to investigate the susceptibility loci for common forms of stroke. METHODS: Fifty-six families, containing multiple cases of stroke and whose data had been previously used to replicate linkage to the phosphodiesterase 4D locus on chromosome 5q, were genotyped in a genome-wide scan. Fine mapping was performed, and subsequently 53 additional families from the same region were genotyped over the candidate regions. RESULTS: Linkage calculations were performed by using 3 different disease models, from a very broad (all stroke cases defined by World Health Organization MONICA criteria) to a narrower (ischemic stroke only) stroke phenotype. With all models, nonparametric multipoint linkage analysis yielded allele-sharing log of the odds (LOD) scores >1.2 at 9 locations: 1p34, 5q13, 7q35, 9q22, 9q34, 13q32, 14q32, 18p11, and 20q13. The highest allele-sharing LOD scores were obtained on chromosomes 5q (previously reported), 1p (LOD=2.09), and 18p (LOD=2.14). Fine mapping resulted in increased allele-sharing LOD scores for chromosome 5q (previously reported) and 9q22 (LOD=1.56), but all others decreased. Combining these initial results with a subsequent analysis of 53 additional families, we obtained the highest allele-sharing LOD scores on chromosomes 5q, 13q, and 18p, although none reached the initial genome-wide allele-sharing LOD scores. CONCLUSIONS: Genetic analysis of stroke in families from northern Sweden did not identify any new major stroke loci. This indicates that multiple minor susceptibility loci in addition to the previously known locus on chromosome 5 could contribute to the disease.

TCF7L2 polymorphisms are associated with type 2 diabetes in northern Sweden.

A recent study found association of one microsatellite and five single nucleotide polymorphisms (SNPs) in intron 3 of the TCF7L2 gene with type 2 diabetes (T2D) in the Icelandic, Danish and American populations. The aim of the present study was to investigate if those SNPs were associated to T2D in two (family- and population-based) cohorts from northern Sweden. We genotyped four of the associated SNPs in a case-control cohort consisting of 872 T2D cases and 857 controls matched with respect to age, sex and geographical origin and in a sample of 59 extended families (148 affected and 83 unaffected individuals). Here, we report replication of association between T2D and three SNPs in the case-control (rs7901695, P=0.003; rs7901346, P=0.00002; and rs12255372, P=0.000004) and two SNPs in the family-based (rs7901695, P=0.01 and rs7901346, P=0.04) samples from northern Sweden. This replication strengthens the evidence for involvement of TCF7L2 in T2D.

Linkage but not association of calpain-10 to type 2 diabetes replicated in northern Sweden.

We present data from a genome-wide scan identifying genetic factors conferring susceptibility to type 2 diabetes. The linkage analysis was based on 59 families from northern Sweden, consisting of a total of 129 cases of type 2 diabetes and 19 individuals with impaired glucose tolerance. Model-free linkage analysis revealed a maximum multipoint logarithm of odds score of 3.19 for D2S2987 at 267.7 cM (P=0.00058), suggesting that a gene conferring susceptibility to type 2 diabetes in the northern Swedish population resides in the 2q37 region. These data replicate, in a European population, previously identified linkage of marker loci in this region to type 2 diabetes in Mexican Americans. In contrast, no evidence in support of association to the previously identified single nucleotide polymorphisms in the calpain-10 gene was observed in a case-control cohort derived from the same population.

CT60 genotype does not affect CTLA-4 isoform expression despite association to T1D and AITD in northern Sweden.

BACKGROUND: Polymorphisms in and around the CTLA-4 gene have previously been associated to T1D and AITD in several populations. One such single nucleotide polymorphism (SNP), CT60, has been reported to affect the expression level ratio of the soluble (sCTLA-4) to full length CTLA-4 (flCTLA-4) isoforms. The aims of our study were to replicate the association previously published by Ueda et al. of polymorphisms in the CTLA-4 region to T1D and AITD and to determine whether the CT60 polymorphism affects the expression level ratio of sCTLA-4/flCTLA-4 in our population. METHODS: Three SNPs were genotyped in 253 cases (104 AITD cases and 149 T1D cases) and 865 ethnically matched controls. Blood from 23 healthy individuals was used to quantify mRNA expression of CTLA-4 isoforms in CD4+ cells using real-time PCR. Serum from 102 cases and 59 healthy individuals was used to determine the level of sCTLA-4 protein. RESULTS: Here we show association of the MH30, CT60 and JO31 polymorphisms to T1D and AITD in northern Sweden. We also observed a higher frequency of the CT60 disease susceptible allele in our controls compared to the British, Italian and Dutch populations, which might contribute to the high frequency of T1D in Sweden. In contrast to previously published findings, however, we were unable to find differences in the sCTLA-4/flCTLA-4 expression ratio based on the CT60 genotype in 23 healthy volunteers, also from northern Sweden. Analysis of sCTLA-4 protein levels in serum showed no correlation between sCTLA-4 protein levels and disease status or CT60 genotype. CONCLUSION: Association was found between T1D/AITD and all three polymorphisms investigated. However, in contrast to previous investigations, sCTLA-4 RNA and protein expression levels did not differ based on CT60 genotype. Our results do not rule out the CT60 SNP as an important polymorphism in the development of T1D or AITD, but suggest that further investigations are necessary to elucidate the effect of the CTLA-4 region on the development of T1D and AITD.

Linkage of ischemic stroke to the PDE4D region on 5q in a Swedish population.

BACKGROUND AND PURPOSE: Recent Icelandic studies have demonstrated linkage for common forms of stroke to chromosome 5q12 and association between phosphodiesterase4D (PDE4D) and ischemic stroke. Using a candidate region approach, we wanted to test the validity of these findings in a different population from northern Sweden. METHODS: A total of 56 families with 117 affected individuals were included in the linkage study. Genotyping was performed with polymorphic microsatellite markers with an average distance of 4.5 cM on chromosome 5. In the association study, 275 cases of first-ever stroke were included together with 550 matched community controls. Polymorphisms were tested individually for association of PDE4D to stroke. RESULTS: Maximum allele-sharing lod score in favor of linkage was observed at marker locus D5S424 (lod score=2.06; P=0.0010). Conditional logistic regression calculations revealed no significant association of ischemic stroke to the defined at-risk allele in PDE4D (odds ratio, 1.1; 95% confidence interval, 0.84 to 1.45). A protective effect may though be implied for 2 of the polymorphisms analyzed in PDE4D. CONCLUSIONS: Using a candidate region approach in a set of stroke families from northern Sweden, we have replicated linkage of stroke susceptibility to the PDE4D gene region on chromosome 5q. Association studies in an independent nested case-control sample from the same geographically located population suggested that different alleles confer susceptibility/protection to stroke in the Icelandic and the northern Swedish populations.

Eight new genomes and synthetic controls increase the accessibility of rapid melt-MAMA SNP typing of Coxiella burnetii.

The case rate of Q fever in Europe has increased dramatically in recent years, mainly because of an epidemic in the Netherlands in 2009. Consequently, there is a need for more extensive genetic characterization of the disease agent Coxiella burnetii in order to better understand the epidemiology and spread of this disease. Genome reference data are essential for this purpose, but only thirteen genome sequences are currently available. Current methods for typing C. burnetii are criticized for having problems in comparing results across laboratories, require the use of genomic control DNA, and/or rely on markers in highly variable regions. We developed in this work a method for single nucleotide polymorphism (SNP) typing of C. burnetii isolates and tissue samples based on new assays targeting ten phylogenetically stable synonymous canonical SNPs (canSNPs). These canSNPs represent previously known phylogenetic branches and were here identified from sequence comparisons of twenty-one C. burnetii genomes, eight of which were sequenced in this work. Importantly, synthetic control templates were developed, to make the method useful to laboratories lacking genomic control DNA. An analysis of twenty-one C. burnetii genomes confirmed that the species exhibits high sequence identity. Most of its SNPs (7,493/7,559 shared by >1 genome) follow a clonal inheritance pattern and are therefore stable phylogenetic typing markers. The assays were validated using twenty-six genetically diverse C. burnetii isolates and three tissue samples from small ruminants infected during the epidemic in the Netherlands. Each sample was assigned to a clade. Synthetic controls (vector and PCR amplified) gave identical results compared to the corresponding genomic controls and are viable alternatives to genomic DNA. The results from the described method indicate that it could be useful for cheap and rapid disease source tracking at non-specialized laboratories, which requires accurate genotyping, assay accessibility and inter-laboratory comparisons.

Telomere length and heredity: Indications of paternal inheritance.

Cellular telomere length is linked to replicative life span. Telomere repeats are lost in peripheral blood cells in vivo by age, and women show less telomere attrition than men. Previous reports have indicated that telomere length and chromosome-specific telomere-length patterns partly are inherited. The mode of heredity has not been clarified, but a link to the X chromosome was recently suggested. We analyzed peripheral mononuclear cells from 49 unrelated families for telomere length using a real-time PCR method. Short-term cultured Epstein-Barr virus-transformed lymphoblasts from the same individuals (n = 130) were analyzed for ability to maintain telomere length and possible gender-linked inheritance. A statistically significant association between telomere lengths comparing father-son (P = 0.011, n = 20) and father-daughter (P = 0.005, n = 22) pairs was found. However, no correlation was observed between mother-daughter (P = 0.463, n = 23) or mother-son (P = 0.577, n = 18). The father-offspring correlation was highly significant (P < 0.0001), in contrast to mother-offspring (P = 0.361). Epstein-Barr virus cultures demonstrated in most cases telomere preservation inversely related to initial mononuclear cell telomere length with short telomeres displaying the most pronounced elongation. Telomere length is inherited, and evidence for a father-to-offspring heritage of this trait was obtained, whereas in vitro telomere length maintenance was found to be dependent on the initial telomere length.

Homozygosity mapping of familial glioma in Northern Sweden.

About 5% of glioma cases are familial. Most glioma families are not ascribed to the well-known glioma predisposing syndromes. One segregation analysis has supported an autosomal recessive gene in glioma families, which could be studied by homozygosity mapping. The ancestors of seven glioma families from the northern region of Sweden were traced through genealogical databases. A common ancestor and inbreeding were traced to give support to an autosomal recessive gene. Homozygosity mapping was performed with a genome-wide scan of 811 markers with linkage calculations. The families were geographically mapped to see if familial glioma was more common in northern compared with southern Sweden. Three of the seven families were remotely related. Homozygosity mapping did not reveal any allele homozygous for all three families. However, there was a geographical clustering of glioma families in the northern region of Sweden. A non-parametric analysis showed an allele-sharing LOD score of 1.05 for marker D1S196 on chromosome 1q23. Genealogical studies linking glioma families might be a tool for linkage in a small set of families. This study did not support an autosomal recessive gene, implicating a low penetrant dominant gene as a possible explanation to the glioma family clustering.

A multivariate approach applied to microarray data for identification of genes with cell cycle-coupled transcription.

We have analyzed microarray data using a modeling approach based on the multivariate statistical method partial least squares (PLS) regression to identify genes with periodic fluctuations in expression levels coupled to the cell cycle in the budding yeast, Saccharomyces cerevisiae. PLS has major advantages for analyzing microarray data since it can model data sets with large numbers of variables and with few observations. A response model was derived describing the expression profile over time expected for periodically transcribed genes, and was used to identify budding yeast transcripts with similar profiles. PLS was then used to interpret the importance of the variables (genes) for the model, yielding a ranking list of how well the genes fitted the generated model. Application of an appropriate cutoff value, calculated from randomized data, allows the identification of genes whose expression appears to be synchronized with cell cycling. Our approach also provides information about the stage in the cell cycle where their transcription peaks. Three synchronized yeast cell microarray data sets were analyzed, both separately and combined. Cell cycle-coupled periodicity was suggested for 455 of the 6,178 transcripts monitored in the combined data set, at a significance level of 0.5%. Among the candidates, 85% of the known periodic transcripts were included. Analysis of the three data sets separately yielded similar ranking lists, showing that the method is robust.