Natural variation in autumn expression is the major adaptive determinant distinguishing Arabidopsis FLC haplotypes

In Arabidopsis thaliana, winter is registered during vernalization through the temperature-dependent repression and epigenetic silencing of floral repressor FLOWERING LOCUS C (FLC). Natural Arabidopsis accessions show considerable variation in vernalization. However, which aspect of the FLC repression mechanism is most important for adaptation to different environments is unclear. By analysing FLC dynamics in natural variants and mutants throughout winter in three field sites, we find that autumnal FLC expression, rather than epigenetic silencing, is the major variable conferred by the distinct Arabidopsis FLChaplotypes. This variation influences flowering responses of Arabidopsis accessions resulting in an interplay between promotion and delay of flowering in different climates to balance survival and, through a post-vernalization effect, reproductive output. These data reveal how expression variation through non-coding cis variation at FLC has enabled Arabidopsis accessions to adapt to different climatic conditions and year-on-year fluctuations.

Evolution of P2A and P5A ATPases: ancient gene duplications and the red algal connection to green plants revisited

In a search for slowly evolving nuclear genes that may cast light on the deep evolution of plants, we carried out phylogenetic analyses of two well-characterized subfamilies of P-type pumps (P2A and P5A ATPases) from representative branches of the eukaryotic tree of life. Both P-type ATPase genes were duplicated very early in eukaryotic evolution and before the divergence of the present eukaryotic supergroups. Synapomorphies identified in the sequences provide evidence that green plants and red algae are more distantly related than are green plants and eukaryotic supergroups in which secondary or tertiary plastids are common, such as several groups belonging to the clade that includes Stramenopiles, Alveolata, Rhizaria, Cryptophyta and Haptophyta (SAR). We propose that red algae branched off soon after the first photosynthesizing eukaryote had acquired a primary plastid, while in another lineage that led to SAR, the primary plastid was lost but, in some cases, regained as a secondary or tertiary plastid.

High Diversity of Microcystin Chemotypes within a Summer Bloom of the Cyanobacterium Microcystis botrys

The fresh-water cyanobacterium Microcystis is known to form blooms world-wide, and is often responsible for the production of microcystins found in lake water. Microcystins are non-ribosomal peptides with toxic effects, e.g. on vertebrates, but their function remains largely unresolved. Moreover, not all strains produce microcystins, and many different microcystin variants have been described. Here we explored the diversity of microcystin variants within Microcystis botrys, a common bloom-former in Sweden. We isolated a total of 130 strains through the duration of a bloom in eutrophic Lake Vomb, and analyzed their microcystin profiles with tandem mass spectrometry (LC-MS/MS). We found that microcystin producing (28.5%) and non-producing (71.5%) M. botrys strains, co-existed throughout the bloom. However, microcystin producing strains were more prevalent towards the end of the sampling period. Overall, 26 unique M. botrys chemotypes were identified, and while some chemotypes re-occurred, others were found only once. The M. botrys chemotypes showed considerable variation both in terms of number of microcystin variants, as well as in what combinations the variants occurred. To our knowledge, this is the first report on microcystin chemotype variation and dynamics in M. botrys. In addition, our study verifies the co-existence of microcystin and non-microcystin producing strains, and we propose that environmental conditions may be implicated in determining their composition.

RNAi as a Tool to Study Virulence in the Pathogenic Yeast Candida glabrata

The yeast Candida glabrata is a major opportunistic pathogen causing mucosal and systemic infections in humans. Systemic infections caused by this yeast have high mortality rates and are difficult to treat due to this yeast’s intrinsic and frequently adapting antifungal resistance. To understand and treat C. glabrata infections, it is essential to investigate the molecular basis of C. glabrata virulence and resistance. We established an RNA interference (RNAi) system in C. glabrata by expressing the Dicer and Argonaute genes from Saccharomyces castellii (a budding yeast with natural RNAi). Our experiments with reporter genes and putative virulence genes showed that the introduction of RNAi resulted in 30 and 70% gene-knockdown for the construct-types antisense and hairpin, respectively. The resulting C. glabrata RNAi strain was used for the screening of a gene library for new virulence-related genes. Phenotypic profiling with a high-resolution quantification of growth identified genes involved in the maintenance of cell integrity, antifungal drugs, and ROS resistance. The genes identified by this approach are promising targets for the treatment of C. glabrata infections.

Neofunctionalization of Mitochondrial Proteins and Incorporation into Signaling Networks in Plants

Because of their symbiotic origin, many mitochondrial proteins are well conserved across eukaryotic kingdoms. It is however less obvious how specific lineages have obtained novel nuclear-encoded mitochondrial proteins. Here, we report a case of mitochondrial neofunctionalization in plants. Phylogenetic analysis of genes containing the Domain of Unknown Function 295 (DUF295) revealed that the domain likely originated in Angiosperms. The C-terminal DUF295 domain is usually accompanied by an N-terminal F-box domain, involved in ubiquitin ligation via binding with ASK1/SKP1-type proteins. Due to gene duplication, the gene family has expanded rapidly, with 94 DUF295-related genes in Arabidopsis thaliana alone. Two DUF295 family subgroups have uniquely evolved and quickly expanded within Brassicaceae. One of these subgroups has completely lost the F-box, but instead obtained strongly predicted mitochondrial targeting peptides. We show that several representatives of this DUF295 Organellar group are effectively targeted to plant mitochondria and chloroplasts. Furthermore, many DUF295 Organellar genes are induced by mitochondrial dysfunction, whereas F-Box DUF295 genes are not. In agreement, several Brassicaceae-specific DUF295 Organellar genes were incorporated in the evolutionary much older ANAC017-dependent mitochondrial retrograde signaling pathway. Finally, a representative set of DUF295 T-DNA insertion mutants was created. No obvious aberrant phenotypes during normal growth and mitochondrial dysfunction were observed, most likely due to the large extent of gene duplication and redundancy. Overall, this study provides insight into how novel mitochondrial proteins can be created via “intercompartmental” gene duplication events. Moreover, our analysis shows that these newly evolved genes can then be specifically integrated into relevant, pre-existing coexpression networks.

Temperature Sensing Is Distributed throughout the Regulatory Network that Controls FLC Epigenetic Silencing in Vernalization

Many organisms need to respond to complex, noisy environmental signals for developmental decision making. Here, we dissect how Arabidopsis plants integrate widely fluctuating field temperatures over month-long timescales to progressively upregulate VERNALIZATION INSENSITIVE3 (VIN3) and silence FLOWERING LOCUS C (FLC), aligning flowering with spring. We develop a mathematical model for vernalization that operates on multiple timescales-long term (month), short term (day), and current (hour)-and is constrained by experimental data. Our analysis demonstrates that temperature sensing is not localized to specific nodes within the FLC network. Instead, temperature sensing is broadly distributed, with each thermosensory process responding to specific features of the plants' history of exposure to warm and cold. The model accurately predicts FLC silencing in new field data, allowing us to forecast FLC expression in changing climates. We suggest that distributed thermosensing may be a general property of thermoresponsive regulatory networks in complex natural environments.

Absence of warmth permits epigenetic memory of winter in Arabidopsis

Plants integrate widely fluctuating temperatures to monitor seasonal progression. Here, we investigate the temperature signals in field conditions that result in vernalisation, the mechanism by which flowering is aligned with spring. We find that multiple, distinct aspects of the temperature profile contribute to vernalisation. In autumn, transient cold temperatures promote transcriptional shutdown of Arabidopsis FLOWERING LOCUS C (FLC), independently of factors conferring epigenetic memory. As winter continues, expression of VERNALIZATION INSENSITIVE3 (VIN3), a factor needed for epigenetic silencing, is upregulated by at least two independent thermosensory processes. One integrates long-term cold temperatures, while the other requires the absence of daily temperatures above 15 °C. The lack of spikes of high temperature, not just prolonged cold, is thus the major driver for vernalisation. Monitoring of peak daily temperature is an effective mechanism to judge seasonal progression, but is likely to have deleterious consequences for vernalisation as the climate becomes more variable.

Large deletions of the PROS1 gene in a large fraction of mutation-negative patients with protein S deficiency

Protein S deficiency is an autosomal dominant disorder that results from mutations in the PROS1 gene. Conventional mutation detection techniques fail to detect a pathogenic PROS1 mutation in approximately 50% of cases. The present study investigates whether large deletions of PROS1 are found in families where mutations in the PROS1 gene have not been found despite sequencing. For this purpose, a dense set of SNP and microsatellite markers were used in segregation analysis to identify deletions. Large deletions were identified by this technique in three out of eight investigated families (38%). The deletions encompassed at least 35 kb, 437 kb and 449 kb respectively. The deletions were confirmed by quantitative PCR. Haplotype analysis showed that the three large deletions and the five other disease haplotypes were all different. All of the eight disease haplotypes co-segregated with protein S deficiency, but each of the five non-deletion haplotypes were present also in normal individuals. In conclusion: Large deletions of PROS1 are relatively common in protein S deficiency patients and screening for large deletions in PROS1 mutation-negative individuals are therefore warranted.

Genome Sequencing Reveals the Origin of the Allotetraploid Arabidopsis suecica

Polyploidy is an example of instantaneous speciation when it involves the formation of a new cytotype that is incompatible with the parental species. Because new polyploid individuals are likely to be rare, establishment of a new species is unlikely unless polyploids are able to reproduce through self-fertilization (selfing), or asexually. Conversely, selfing (or asexuality) makes it possible for polyploid species to originate from a single individual-a bona fide speciation event. The extent to which this happens is not known. Here, we consider the origin of Arabidopsis suecica, a selfing allopolyploid between Arabidopsis thaliana and Arabidopsis arenosa, which has hitherto been considered to be an example of a unique origin. Based on whole-genome re-sequencing of 15 natural A. suecica accessions, we identify ubiquitous shared polymorphism with the parental species, and hence conclusively reject a unique origin in favor of multiple founding individuals. We further estimate that the species originated after the last glacial maximum in Eastern Europe or central Eurasia (rather than Sweden, as the name might suggest). Finally, annotation of the self-incompatibility loci in A. suecica revealed that both loci carry non-functional alleles. The locus inherited from the selfing A. thaliana is fixed for an ancestral non-functional allele, whereas the locus inherited from the outcrossing A. arenosa is fixed for a novel loss-of-function allele. Furthermore, the allele inherited from A. thaliana is predicted to transcriptionally silence the allele inherited from A. arenosa, suggesting that loss of self-incompatibility may have been instantaneous.

AS3MT-mediated tolerance to arsenic evolved by multiple independent horizontal gene transfers from bacteria to eukaryotes

Organisms have evolved the ability to tolerate toxic substances in their environments, often by producing metabolic enzymes that efficiently detoxify the toxicant. Inorganic arsenic is one of the most toxic and carcinogenic substances in the environment, but many organisms, including humans, metabolise inorganic arsenic to less toxic metabolites. This multistep process produces mono-, di-, and trimethylated arsenic metabolites, which the organism excretes. In humans, arsenite methyltransferase (AS3MT) appears to be the main metabolic enzyme that methylates arsenic. In this study, we examined the evolutionary origin of AS3MT and assessed the ability of different genotypes to produce methylated arsenic metabolites. Phylogenetic analysis suggests that multiple, independent horizontal gene transfers between different bacteria, and from bacteria to eukaryotes, increased tolerance to environmental arsenic during evolution. These findings are supported by the observation that genetic variation in AS3MT correlates with the capacity to methylate arsenic. Adaptation to arsenic thus serves as a model for how organisms evolve to survive under toxic conditions.

Genetic variation of the Toll-like receptors in a Swedish allergic rhinitis case population

BACKGROUND

Variation in the 10 toll-like receptor (TLR) genes has been significantly associated with allergic rhinitis (AR) in several candidate gene studies and three large genome-wide association studies. These have all investigated common variants, but no investigations for rare variants (MAF ≤ 1%) have been made in AR. The present study aims to describe the genetic variation of the promoter and coding sequences of the 10 TLR genes in 288 AR patients.

METHODS

Sanger sequencing and Ion Torrent next-generation sequencing was used to identify polymorphisms in a Swedish AR population and these were subsequently compared and evaluated using 1000Genomes and Exome Aggregation Consortium (ExAC) data.

RESULTS

The overall level of genetic variation was clearly different among the 10 TLR genes. The TLR10-TLR1-TLR6 locus was the most variable, while the TLR7-TLR8 locus was consistently showing a much lower level of variation. The AR patients had a total of 37 promoter polymorphisms with 14 rare (MAF ≤ 1%) and 14 AR-specific polymorphisms. These numbers were highly similar when comparing the AR and the European part of the 1000Genomes populations, with the exception of TLR10 where a significant (P = 0.00009) accumulation of polymorphisms were identified. The coding sequences had a total of 119 polymorphisms, 68 were rare and 43 were not present in the European part of the 1000Genomes population. Comparing the numbers of rare and AR-specific SNPs in the patients with the European part of the 1000Genomes population it was seen that the numbers were quite similar both for individual genes and for the sum of all 10 genes. However, TLR1, TLR5, TLR7 and TLR9 showed a significant excess of rare variants in the AR population when compared to the non-Finnish European part of ExAC. In particular the TLR1 S324* nonsense mutation was clearly overrepresented in the AR population.

CONCLUSIONS

Most TLR genes showed a similar level of variation between AR patients and public databases, but a significant excess of rare variants in AR patients were detected in TLR1, TLR5, TLR7, TLR9 and TLR10. This further emphasizes the frequently reproduced TLR10-TLR1-TLR6 locus as being involved in the pathogenesis of allergic rhinitis.

Sequencing of the genus Arabidopsis identifies a complex history of nonbifurcating speciation and abundant trans-specific polymorphism

The notion of species as reproductively isolated units related through a bifurcating tree implies that gene trees should generally agree with the species tree and that sister taxa should not share polymorphisms unless they diverged recently and should be equally closely related to outgroups. It is now possible to evaluate this model systematically. We sequenced multiple individuals from 27 described taxa representing the entire Arabidopsis genus. Cluster analysis identified seven groups, corresponding to described species that capture the structure of the genus. However, at the level of gene trees, only the separation of Arabidopsis thaliana from the remaining species was universally supported, and, overall, the amount of shared polymorphism demonstrated that reproductive isolation was considerably more recent than the estimated divergence times. We uncovered multiple cases of past gene flow that contradict a bifurcating species tree. Finally, we showed that the pattern of divergence differs between gene ontologies, suggesting a role for selection.

Chronic Rhinosinusitis Patients Show Accumulation of Genetic Variants in PARS2

Genetic studies of chronic rhinosinusitis (CRS) have identified a total of 53 CRS-associated SNPs that were subsequently evaluated for their reproducibility in a recent study. The rs2873551 SNP in linkage disequilibrium with PARS2 showed the strongest association signal. The present study aims to comprehensively screen for rare variants in PARS2 and evaluate for accumulation of such variants in CRS-patients. Sanger sequencing and long-range PCR were used to screen for rare variants in the putative promoter region and coding sequence of 310 CRS-patients and a total of 21 variants were detected. The mutation spectrum was then compared with data from European populations of the 1000Genomes project (EUR) and the Exome Aggregation Consortium (ExAC). The CRS population showed a significant surplus of low-frequency variants compared with ExAC data. Haplotype analysis of the region showed a significant excess of rare haplotypes in the CRS population compared to the EUR population. Two missense mutations were also genotyped in the 310 CRS patients and 372 CRS-negative controls, but no associations with the disease were found. This is the first re-sequencing study in CRS research and also the first study to show an association of rare variants with the disease.

Characterization of genetic variation in TLR8 in relation to allergic rhinitis

BACKGROUND

A previous investigation of all 10 TLR genes for associations with allergic rhinitis (AR) detected a number of significant SNPs in the TLR8 locus. The associations indicated that an accumulation of rare variants could explain the signal. This study therefore searches for rare variants in the TLR8 region and also investigates the reproducibility of previous SNP associations.

METHODS

The TLR8 gene was resequenced in 288 AR patients from Malmö and the data were compared with publically available data. Seven previously AR-associated SNPs from TLR8 were analyzed for AR associations in 422 AR patients and 859 controls from the BAMSE cohort. The associations detected in present and previous studies were compared.

RESULTS

Sequencing detected 13 polymorphisms (three promotor and 10 coding) among 288 AR patients. Four of the coding polymorphisms were rare (MAF < 1%) and three of those were novel. Two coding polymorphisms were benign missense mutations and the rest were synonymous. Comparison with 1000Genomes and Exome Aggregation Consortium data revealed no accumulation of rare variants in the AR cases. The AR association tests made using the BAMSE cohort yielded five P-values <0.05. Tests of IgE levels yielded four significant SNP associations to birch pollen. Comparing results between different populations revealed opposing risk alleles, different gender effects, and response to different allergens in the different populations.

CONCLUSIONS

Rare variants in TLR8 are not associated with AR. Comparison of present and previous association studies reveals contradictory results for common variants. Thus, no associations exist between genetic variation in TLR8 and AR.

Replication of genomewide associations with allergic sensitization and allergic rhinitis

BACKGROUND

Three genomewide metastudies have recently reported associations with self-reported allergic rhinitis and allergic sensitization. The three studies together identified a set of 37 loci but showed low concordance. This study investigates the reproducibility of the detected single nucleotide polymorphism (SNP) associations in an extensively characterized longitudinal cohort, BAMSE.

METHODS

Phenotypic evaluation of allergic rhinitis (AR) and allergic sensitization was performed on 2153 children from BAMSE at 8 and 16 years of age. Allele frequencies of 39 SNPs were investigated for association with the exact allergic phenotypes of the metastudies. Odds ratios and false discovery rates were calculated, and the impact of asthma was evaluated. The cases were also evaluated for age at onset effects (≤ or >8 years of age).

RESULTS

Association tests of the 39 SNPs identified 12 SNPs with P-values < 0.05 and Q-values < 0.10. Two of the four loci (TLR6-TLR1 and HLA-DQA1-HLA-DQB1) identified in all three original studies were also identified in this study. Three SNPs located in the TLR6-TLR1 locus had the lowest P-values and Q-values < 0.1 when using a well-defined AR phenotype. Two loci showed significant age at onset effects, but the effect of asthma on the associations was very limited.

CONCLUSION

The TLR6-TLR1 locus is likely to have a central role in the development of allergic disease. The association between genetic variation in the SSTR1-MIPOL1 and TSLP-SLC25A46 loci and age at onset is the first report of age at onset effects in allergic rhinitis.

Analysis of the yeast short-term Crabtree effect and its origin

The short-term Crabtree effect is defined as the immediate occurrence of aerobic alcoholic fermentation in response to provision of a pulse of excess sugar to sugar-limited yeast cultures. Here we have characterized ten yeast species with a clearly defined phylogenetic relationship. Yeast species were cultivated under glucose-limited conditions, and we studied their general carbon metabolism in response to a glucose pulse. We generated an extensive collection of data on glucose and oxygen consumption, and ethanol and carbon dioxide generation. We conclude that the Pichia, Debaryomyces, Eremothecium and Kluyveromyces marxianus yeasts do not exhibit any significant ethanol formation, while Kluyveromyces lactis behaves as an intermediate yeast, and Lachancea, Torulaspora, Vanderwaltozyma and Saccharomyces yeasts exhibit rapid ethanol accumulation. Based on the present data and our previous data relating to the presence of the long-term Crabtree effect in over 40 yeast species, we speculate that the origin of the short-term effect may coincide with the origin of the long-term Crabtree effect in the Saccharomycetales lineage, occurring ~ 150 million years ago.

Genetic variation in KLK2 and KLK3 is associated with concentrations of hK2 and PSA in serum and seminal plasma in young men

BACKGROUND

Genetic variants in KLK2 and KLK3 have been associated with increased serum concentrations of their encoded proteins, human kallikrein-related peptidase 2 (hK2) and prostate-specific antigen (PSA), and with prostate cancer in older men. Low PSA concentrations in seminal plasma (SP) have been associated with low sperm motility. To evaluate whether KLK2 and KLK3 genetic variants affect physiological prostatic secretion, we studied the association of SNPs with hK2 and PSA concentrations in SP and serum of young, healthy men.

METHODS

Leukocyte DNA was extracted from 303 male military conscripts (median age 18.1 years). Nine SNPs across KLK2-KLK3 were genotyped. We measured PSA and hK2 in SP and serum using immunofluorometric assays. The association of genotype frequencies with hK2 and PSA concentrations was tested with the Kruskal-Wallis test.

RESULTS

Four KLK2 SNPs (rs198972, rs198977, rs198978, and rs80050017) were strongly associated with hK2 concentrations in SP and serum, with individuals homozygous for the major alleles having 3- to 7-fold higher concentrations than the intermediate concentrations found in other homozygotes and heterozygotes (all P < 0.001). Three of these SNPs were significantly associated with percentage of free PSA (%fPSA) in serum (all P < 0.007). Three KLK3 SNPs showed associations with PSA in SP, and the rs1058205 SNP was associated with total PSA in serum (P = 0.001) and %fPSA (P = 0.015).

CONCLUSIONS

Associations observed in young, healthy men between the SP and serum concentrations of hK2 and PSA and several genetic variants in KLK2 and KLK3 could be useful to refine models of PSA cutoff values in prostate cancer testing.

Origin of Swedish hemophilia B mutations

BACKGROUND

More than 1100 mutations that cause hemophilia B (HB) have been identified. At the same time, specific F9 mutations are present at high frequencies in certain populations, which raise questions about the origin of HB mutations.

OBJECTIVES

To describe the mutation spectrum of all HB families in Sweden and investigate if mutations appearing in several families are due to independent recurrent mutations (RMs) or to a common mutation event (i.e. are identical by descent (IBD)).

PATIENTS/METHODS

The registered Swedish HB population consists of patients from 86 families. Mutations were identified by resequencing and identical haplotypes were defined using 74 markers and a control population of 285 individuals. The ages of IBD mutations were estimated using ESTIAGE.

RESULTS

Out of 77 presumably unrelated patients with substitution mutations, 47 patients (61%) had mutations in common with other patients. Haplotyping of the 47 patients showed that 24 patients had IBD mutations (51%) with estimated ages of between two and 23 generations. A majority of these patients had mild disease. Eight of the 15 mutations observed in more than one family were C>T transitions in CpG sites and all eight were RMs.

CONCLUSIONS

The association of IBD mutations with a mild phenotype is similar to what has been previously observed in hemophilia A. Noteworthy features of the mutations that are common to more than one family are the equal proportions of patients with RM and IBD mutations and the correlation between the occurrence of RMs and C>T transitions at CpG sites.

Copy number variants in the kallikrein gene cluster

The kallikrein gene family (KLK1-KLK15) is the largest contiguous group of protease genes within the human genome and is associated with both risk and outcome of cancer and other diseases. We searched for copy number variants in all KLK genes using quantitative PCR analysis and analysis of inheritance patterns of single nucleotide polymorphisms. Two deletions were identified: one 2235-bp deletion in KLK9 present in 1.2% of alleles, and one 3394-bp deletion in KLK15 present in 4.0% of alleles. Each deletion eliminated one complete exon and created out-of-frame coding that eliminated the catalytic triad of the resulting truncated gene product, which therefore likely is a non-functional protein. Deletion breakpoints identified by DNA sequencing located the KLK9 deletion breakpoint to a long interspersed element (LINE) repeated sequence, while the deletion in KLK15 is located in a single copy sequence. To search for an association between each deletion and risk of prostate cancer (PC), we analyzed a cohort of 667 biopsied men (266 PC cases and 401 men with no evidence of PC at biopsy) using short deletion-specific PCR assays. There was no association between evidence of PC in this cohort and the presence of either gene deletion. Haplotyping revealed a single origin of each deletion, with most recent common ancestor estimates of 3000-8000 and 6000-14 000 years for the deletions in KLK9 and KLK15, respectively. The presence of the deletions on the same haplotypes in 1000 Genomes data of both European and African populations indicate an early origin of both deletions. The old age in combination with homozygous presence of loss-of-function variants suggests that some kallikrein-related peptidases have non-essential functions.

Yeast "make-accumulate-consume" life strategy evolved as a multi-step process that predates the whole genome duplication

When fruits ripen, microbial communities start a fierce competition for the freely available fruit sugars. Three yeast lineages, including baker’s yeast Saccharomyces cerevisiae, have independently developed the metabolic activity to convert simple sugars into ethanol even under fully aerobic conditions. This fermentation capacity, named Crabtree effect, reduces the cell-biomass production but provides in nature a tool to out-compete other microorganisms. Here, we analyzed over forty Saccharomycetaceae yeasts, covering over 200 million years of the evolutionary history, for their carbon metabolism. The experiments were done under strictly controlled and uniform conditions, which has not been done before. We show that the origin of Crabtree effect in Saccharomycetaceae predates the whole genome duplication and became a settled metabolic trait after the split of the S. cerevisiae and Kluyveromyces lineages, and coincided with the origin of modern fruit bearing plants. Our results suggest that ethanol fermentation evolved progressively, involving several successive molecular events that have gradually remodeled the yeast carbon metabolism. While some of the final evolutionary events, like gene duplications of glucose transporters and glycolytic enzymes, have been deduced, the earliest molecular events initiating Crabtree effect are still to be determined.

Investigating highly replicated asthma genes as candidate genes for allergic rhinitis

BACKGROUND

Asthma genetics has been extensively studied and many genes have been associated with the development or severity of this disease. In contrast, the genetic basis of allergic rhinitis (AR) has not been evaluated as extensively. It is well known that asthma is closely related with AR since a large proportion of individuals with asthma also present symptoms of AR, and patients with AR have a 5-6 fold increased risk of developing asthma. Thus, the relevance of asthma candidate genes as predisposing factors for AR is worth investigating. The present study was designed to investigate if SNPs in highly replicated asthma genes are associated with the occurrence of AR.

METHODS

A total of 192 SNPs from 21 asthma candidate genes reported to be associated with asthma in 6 or more unrelated studies were genotyped in a Swedish population with 246 AR patients and 431 controls. Genotypes for 429 SNPs from the same set of genes were also extracted from a Singapore Chinese genome-wide dataset which consisted of 456 AR cases and 486 controls. All SNPs were subsequently analyzed for association with AR and their influence on allergic sensitization to common allergens.

RESULTS

A limited number of potential associations were observed and the overall pattern of P-values corresponds well to the expectations in the absence of an effect. However, in the tests of allele effects in the Chinese population the number of significant P-values exceeds the expectations. The strongest signals were found for SNPs in NPSR1 and CTLA4. In these genes, a total of nine SNPs showed P-values <0.001 with corresponding Q-values <0.05. In the NPSR1 gene some P-values were lower than the Bonferroni correction level. Reanalysis after elimination of all patients with asthmatic symptoms excluded asthma as a confounding factor in our results. Weaker indications were found for IL13 and GSTP1 with respect to sensitization to birch pollen in the Swedish population.

CONCLUSIONS

Genetic variation in the majority of the highly replicated asthma genes were not associated to AR in our populations which suggest that asthma and AR could have less in common than previously anticipated. However, NPSR1 and CTLA4 can be genetic links between AR and asthma and associations of polymorphisms in NPSR1 with AR have not been reported previously.

Poor reproducibility of allergic rhinitis SNP associations

Replication of reported associations is crucial to the investigation of complex disease. More than 100 SNPs have previously been reported as associated with allergic rhinitis (AR), but few of these have been replicated successfully. To investigate the general reproducibility of reported AR-associations in candidate gene studies, one Swedish (352 AR-cases, 709 controls) and one Singapore Chinese population (948 AR-cases, 580 controls) were analyzed using 49 AR-associated SNPs. The overall pattern of P-values indicated that very few of the investigated SNPs were associated with AR. Given published odds ratios (ORs) most SNPs showed high power to detect an association, but no correlations were found between the ORs of the two study populations or with published ORs. None of the association signals were in common to the two genome-wide association studies published in AR, indicating that the associations represent false positives or have much lower effect-sizes than reported.

Origin of Swedish hemophilia A mutations

BACKGROUND

Hemophilia A (HA) has a high level of variation within the disease class, with more than 1000 mutations being listed in the HAMSTeRS database. At the same time a number of F8 mutations are present in specific populations at high frequencies.

OBJECTIVES

The simultaneous presence of large numbers of rare mutations and a small number of high-frequency mutations raises questions about the origins of HA mutations. The present study was aimed at describing the origins of HA mutations in the complete Swedish population. The primary issue was to determine what proportion of identical mutations are identical by descent (IBD) and what proportion are attributable to recurrent mutation events. The age of IBD mutations was also determined.

PATIENTS/METHODS

In Sweden, the care of HA is centralized, and the Swedish HA population consists of ~ 750 patients from > 300 families (35% severe, 15% moderate, and 50% mild). Identical haplotypes were defined by single-nucleotide polymorphism and microsatellite haplotyping, and the ages of the mutations were estimated with estiage.

RESULTS

Among 212 presumably unrelated patients with substitution mutations, 97 (46%) had mutations in common with other patients. Haplotyping of the 97 patients showed that 47 had IBD mutations (22%) with estimated ages of between two and 35 generations. The frequency of mild disease increased with an increasing number of patients sharing the mutations.

CONCLUSIONS

A majority of the IBD mutations are mild and have age estimates of a few hundred years, but some could date back to the Middle Ages.

A high density RFLP linkage map of sugar beet

A high density sugar beet RFLP map with an average distance of 1.5 cM between markers has been constructed. The map covers 621 cM and includes 413 markers distributed over the nine linkage groups of sugar beet. The map is based on two F2 populations representing two different pairs of parents. The two sets of data were integrated into a single map using 90 markers that were common to both data sets. The quality of the map was assessed in several ways. The common markers were used to investigate how often the loci had been mapped in the same order in the two F2 populations. For closely situated markers (<1.5 cM) the order specified in the map is uncertain, but for markers separated by more than 2 cM the locus order is highly reliable. The error rate of the overall process was estimated at 0.3% by independently repeating the analysis of 41 markers. The map is comparatively short, with a map length corresponding to approximately 1.4 crossovers per bivalent. Another feature of the map is a high degree of clustering of markers along the linkage groups. With the possible exception of linkage group 2, each linkage group shows one major cluster, which in most cases is situated in the centre of the linkage group. Our interpretation is that sugar beet, in comparison with most other species, has an extreme localization of recombination. Key words : sugar beet, linkage, RFLP, clustering.

Estimating genetic variation in sugar beets and wild beets using pools of individuals

The study describes the genetic structure in sugar beets and in wild beets (Beta vulgaris) using 30 RFLP markers. Samples consisting of pooled plant material of 100 individuals from each line and population were used to analyse 120 sugar beet breeding lines and 91 wild beet populations. Greater variation was found among the wild populations than among the breeding lines. Although the two major groups of breeding lines, monogerm and multigerm, had approximately equal amounts of genetic variation, in the monogerm group more of this variation was partitioned among the lines than within the lines. Furthermore, despite most of the variation being shared by the two groups, the two groups were found to be separated along the first two components in a principal component analysis. Computer simulations were carried out to evaluate the usefulness of the pooled-sample strategy employed in the investigation. These simulations showed the use of pooled samples to be a better alternative than that of analysing a few plants individually.

Toll-like receptor gene polymorphisms are associated with allergic rhinitis: a case control study

Background

The Toll-like receptor proteins are important in host defense and initiation of the innate and adaptive immune responses. A number of studies have identified associations between genetic variation in the Toll-like receptor genes and allergic disorders such as asthma and allergic rhinitis. The present study aim to search for genetic variation associated with allergic rhinitis in the Toll-like receptor genes.

Methods

A first association analysis genotyped 73 SNPs in 182 cases and 378 controls from a Swedish population. Based on these results an additional 24 SNPs were analyzed in one Swedish population with 352 cases and 709 controls and one Chinese population with 948 cases and 580 controls.

Results

The first association analysis identified 4 allergic rhinitis-associated SNPs in the TLR7-TLR8 gene region. Subsequent analysis of 24 SNPs from this region identified 7 and 5 significant SNPs from the Swedish and Chinese populations, respectively. The corresponding risk-associated haplotypes are significant after Bonferroni correction and are the most common haplotypes in both populations. The associations are primarily detected in females in the Swedish population, whereas it is seen in males in the Chinese population. Further independent support for the involvement of this region in allergic rhinitis was obtained from quantitative skin prick test data generated in both populations.

Conclusions

Haplotypes in the TLR7-TLR8 gene region were associated with allergic rhinitis in one Swedish and one Chinese population. Since this region has earlier been associated with asthma and allergic rhinitis in a Danish linkage study this speaks strongly in favour of this region being truly involved in the development of this disease.

CLC- a novel susceptibility gene for allergic rhinitis?

BACKGROUND

Studies of the nasal lavage fluid proteome have previously identified proteins differently expressed in patients with symptomatic allergic rhinitis, e.g. S100A7, prolactin-inducible protein (PIP), wingless-type MMTV integration site family, member 2B (WNT2B), Charcot-Leyden crystal protein (CLC) and palate lung nasal epithelial clone (PLUNC). The aim of the present study was to investigate if genetic variation associated with allergic rhinitis can be found in these genes.

METHODS

Peripheral blood was collected from 251 patients with birch and/or grass pollen-induced allergic rhinitis and 386 nonatopic healthy controls. A total of 39 single nucleotide polymorphisms (SNPs) distributed over the genes PIP, WNT2B, CLC and PLUNC were selected from dbSNP, genotyped and investigated for associations with allergic rhinitis. Twelve additional SNPs were subsequently analysed for CLC.

RESULTS

All 22 investigated SNPs in CLC were polymorphic. Ten SNPs yielded significant differences between cases and controls with respect to genotype frequencies. Homozygotes for the minor allele were more common in allergic individuals compared to healthy controls. The minor alleles of these SNPs were all located on the same haplotype. Furthermore, homozygotes for the minor allele of two of the promoter SNPs had higher average scores for birch in skin prick test. In contrast, for seven SNPs within the gene, heterozygotes and homozygotes for the major allele had higher average scores for grass. None of the other three genes showed association.

CONCLUSION

Genetic variation in CLC was found to be associated with allergic rhinitis. The pattern of variation is compatible with a recessive inheritance model and the previously observed altered protein levels detected in patients with allergic rhinitis.

A haplotype in the inducible T-cell tyrosine kinase is a risk factor for seasonal allergic rhinitis

BACKGROUND

Identification of disease-associated single nucleotide polymorphisms (SNPs) in seasonal allergic rhinitis (SAR) may be facilitated by focusing on genes in a disease-associated pathway.

OBJECTIVE

To search for SNPs in genes that belong to the T-cell receptor (TCR) pathway and that change in expression in allergen-challenged CD4+ cells from patients with SAR.

METHODS

CD4+ cells from patients with SAR were analysed with gene expression microarrays. Allele, genotype and haplotype frequencies were compared in 251 patients and 386 healthy controls.

RESULTS

Gene expression microarray analysis of allergen-challenged CD4+ cells from patients with SAR showed that 25 of 38 TCR pathway genes were differentially expressed. A total of 62 SNPs were analysed in eight of the 25 genes; ICOS, IL4, IL5, IL13, CSF2, CTLA4, the inducible T-cell tyrosine kinase (ITK) and CD3D. Significant chi-squared values were identified for several markers in the ITK kinase gene region. A total of five SNPs were nominally significant at the 5% level. Haplotype analysis of the five significant SNPs showed increased frequency of a haplotype that covered most of the coding part of ITK. The functional relevance of ITK was supported by analysis of an independent material, which showed increased expression of ITK in allergen-challenged CD4+ cells from patients, but not from controls.

CONCLUSION

Analysis of SNPs in TCR pathway genes revealed that a haplotype that covers a major part of the coding sequence of ITK is a risk factor for SAR.

Comprehensive evaluation of genetic variation in S100A7 suggests an association with the occurrence of allergic rhinitis

Background

S100A7 is a calcium-binding protein with chemotactic and antimicrobial properties. S100A7 protein levels are decreased in nasal lavage fluid from individuals with ongoing allergic rhinitis, suggesting a role for S100A7 in allergic airway inflammation. The aims of this study were to describe genetic variation in S100A7 and search for associations between this variation and allergic rhinitis.

Methods

Peripheral blood was collected from 184 atopic patients with a history of pollen-induced allergic rhinitis and 378 non-atopic individuals, all of Swedish origin. DNA was extracted and the S100A7 gene was resequenced in a subset of 47 randomly selected atopic individuals. Nine polymorphisms were genotyped in 184 atopic and 378 non-atopic individuals and subsequently investigated for associations with allergic rhinitis as well as skin prick test results. Haplotypes were estimated and compared in the two groups.

Results

Thirteen polymorphisms were identified in S100A7, of which 7 were previously undescribed. rs3014837 (G/C), which gives rise to an Asp → Glu amino acid shift, had significantly increased minor allele frequency in atopic individuals. The major haplotype, containing the major allele at all sites, was more common in non-atopic individuals, while the haplotype containing the minor allele at rs3014837 was equally more common among the atopic individuals. Additionally, heterozygotes at this site had significantly higher scores in skin prick tests for 9 out of 11 tested allergens, compared to homozygotes.

Conclusion

This is the first study describing genetic variation, associated with allergy, in S100A7. The results indicate that rs3014837 is linked to allergic rhinitis in our Swedish population and render S100A7 a strong candidate for further investigations regarding its role in allergic inflammation.

The evolutionary history of the common chloroplast genome of Arabidopsis thaliana and A. suecica

The evolutionary history of the common chloroplast (cp) genome of the allotetraploid Arabidopsis suecica and its maternal parent A. thaliana was investigated by sequencing 50 fragments of cpDNA, resulting in 98 polymorphic sites. The variation in the A. suecica sample was small, in contrast to that of the A. thaliana sample. The time to the most recent common ancestor (T(MRCA)) of the A. suecica cp genome alone was estimated to be about one 37th of the T(MRCA) of both the A. thaliana and A. suecica cp genomes. This corresponds to A. suecica having a MRCA between 10 000 and 50 000 years ago, suggesting that the entire species originated during, or before, this period of time, although the estimates are sensitive to assumptions made about population size and mutation rate. The data was also consistent with the hypothesis of A. suecica being of single origin. Isolation-by-distance and population structure in A. thaliana depended upon the geographical scale analysed; isolation-by-distance was found to be weak on the global scale but locally pronounced. Within the genealogical cp tree of A. thaliana, there were indications that the root of the A. suecica species is located among accessions of A. thaliana that come primarily from central Europe. Selective neutrality of the cp genome could not be rejected, despite the fact that it contains several completely linked protein-coding genes.

Evolution of chloroplast mononucleotide microsatellites in Arabidopsis thaliana

The level of variation and the mutation rate were investigated in an empirical study of 244 chloroplast microsatellites in 15 accessions of Arabidopsis thaliana. In contrast to SNP variation, microsatellite variation in the chloroplast was found to be common, although less common than microsatellite variation in the nucleus. No microsatellite variation was found in coding regions of the chloroplast. To evaluate different models of microsatellite evolution as possible explanations for the observed pattern of variation, the length distribution of microsatellites in the published DNA sequence of the A. thaliana chloroplast was subsequently used. By combining information from these two analyses we found that the mode of evolution of the chloroplast mononucleotide microsatellites was best described by a linear relation between repeat length and mutation rate, when the repeat lengths exceeded about 7 bp. This model can readily predict the variation observed in non-coding chloroplast DNA. It was found that the number of uninterrupted repeat units had a large impact on the level of chloroplast microsatellite variation. No other factors investigated--such as the position of a locus within the chromosome, or imperfect repeats--appeared to affect the variability of chloroplast microsatellites. By fitting the slippage models to the Genbank sequence of chromosome 1, we show that the difference between microsatellite variation in the nucleus and the chloroplast is largely due to differences in slippage rate.

Mode of reproduction in Arabidopsis suecica

The breeding system of Arabidopsis suecica was investigated through genetic analysis of microsatellite segregation patterns in five controlled crosses as well as in 16 single-mother families collected in the wild. Analysis of single and two-locus segregations in the F2 generation following a cross clearly shows that A. suecica is reproduces sexually. The single-mother families show a high level of homozygosity corroborating earlier results indicating a high level of inbreeding. The high level of individual homozygosity is due both to a high level of selfing and to the underlying population structure.

A unique recent origin of the allotetraploid species Arabidopsis suecica: Evidence from nuclear DNA markers

A coalescent-based method was used to investigate the origins of the allotetraploid Arabidopsis suecica, using 52 nuclear microsatellite loci typed in eight individuals of A. suecica and 14 individuals of its maternal parent Arabidopsis thaliana, and four short fragments of genomic DNA sequenced in a sample of four individuals of A. suecica and in both its parental species A. thaliana and Arabidopsis arenosa. All loci were variable in A. thaliana but only 24 of the 52 microsatellite loci and none of the four sequence fragments were variable in A. suecica. We explore a number of possible evolutionary scenarios for A. suecica and conclude that it is likely that A. suecica has a recent, unique origin between 12,000 and 300,000 years ago. The time estimates depend strongly on what is assumed about population growth and rates of mutation. When combined with what is known about the history of glaciations, our results suggest that A. suecica originated south of its present distribution in Sweden and Finland and then migrated north, perhaps in the wake of the retreating ice.

Statistical behavior of complex cancer karyotypes

Epithelial tumors commonly show complex and variable karyotypes that obscure the identification of general patterns of the karyotypic evolution. To overcome some of these problems, we previously systematically analyzed the accumulated cytogenetic data from individual tumor types by using various statistical means. In the present study, we compare previous results obtained for nine tumor types and perform several meta-analyses of data obtained from a number of epithelial tumors, including head and neck, kidney, bladder, breast, colorectal, ovarian, and lung cancer, as well as from malignant melanoma and Wilms tumor, with the specific aim of discovering common patterns of karyotypic evolution. We show that these tumors frequently develop through a hypo- or a hyperdiploid pathway and progress by an increasing number of alternative imbalances through at least two karyotypic phases, Phases I and II, and possibly through a third, Phase III. During Phase I, the karyotypes exhibited a power law distribution of both the number of changes per tumor and the frequency distribution at which bands were involved in breaks. At the transition from Phase I to Phase II/III, the observed power law distributions were lost, indicating a transition from an ordered and highly structured process to a disordered and chaotic pattern. The change in karyotypic orderliness at the transition from Phase I to Phase II/III was also shown by a drastic difference in karyotypic entropy.

Identification of cytogenetic subgroups and karyotypic pathways of clonal evolution in follicular lymphomas

Follicular lymphoma (FL) is characterized by the activation of BCL2 through t(14;18)(q32;q21). Additional acquired mutations are necessary to generate a fully malignant clonal proliferation. Many of these secondary genetic alterations are visible in the clonal karyotype; however, the sequence by which they arise and their influence on clinical behavior have not been determined. The ability to address these issues has been hampered by the lack of computational methods to manipulate complex chromosomal data in a sufficiently large cohort of cases. In the present investigation, we analyzed secondary karyotypic alterations in 336 cases of FL with t(14;18) to identify the most common regions of recurrent chromosomal gain or loss. This revealed 29 recurrent changes present in more than 5% of the tumors. Each tumor karyotype was then assessed for the presence or absence of each of these 29 specific changes. By statistical means, we show that the chromosomal changes arise in an apparent temporal order, with distinct early and late changes. We identify, by principal-components analysis, four possible cytogenetic pathways that characterize the early stages of clonal evolution, which converge to a common route at later stages. We show that FLs with t(14;18) may be classified into cytogenetic subgroups determined by the presence or absence of 6q-, +7, or der(18)t(14;18). Correlation with clinical outcomes in a subset of cases with clinical data revealed del(17p) and +12 to be correlated with an adverse clinical outcome. The clinical implications of these pathways of clonal evolution need to be examined on a prospective basis in a large cohort of FLs.

Dissecting karyotypic patterns in malignant melanomas: Temporal clustering of losses and gains in melanoma karyotypic evolution

Malignant melanomas can be divided into two major subtypes, involving either the skin or eye melanomas. Both tumor forms exhibit highly complex karyotypes with nonrandom recurrent chromosomal imbalances. Loss of chromosome 3, the short arm of chromosome 1, and gain of 8q have been suggested to be associated with eye melanomas, whereas gain of 6p and loss of 6q have been more often seen in skin melanomas. Imbalances implicated in tumor progression include among others, -10 and +7. In spite of the abundance of cytogenetic information, with more than 300 published karyotypes, very little is known about the mode of karyotypic evolution or of the presence of possible cytogenetic pathways. In our investigation, we have used 362 melanoma karyotypes, including both the skin and eye subtypes, to identify the most frequently occurring imbalances. Tumor cases were then classified with respect to the presence or absence of these imbalances and statistically analyzed in order to assess the order of appearance of chromosomal imbalances, the presence of karyotypic pathways, as well as possible cytogenetic subtypes. We show that the melanomas develop through one mode of karyotypic evolution, common to both low and high complexity karyotypes, and we establish the temporal order by which the different imbalances occur. By applying several statistical methods, we show that at least two cytogenetic pathways of clonal evolution exist in malignant melanomas, one initiated with -3 and one with +6p, and that these pathways operate in both skin and eye melanomas.

Chloroplast DNA indicates a single origin of the allotetraploid Arabidopsis suecica

DNA sequencing was performed on up to 12 chloroplast DNA regions [giving a total of 4288 base pairs (bp) in length] from the allopolyploid Arabidopsis suecica (48 accessions) and its two parental species, A. thaliana (25 accessions) and A. arenosa (seven accessions). Arabidopsis suecica was identical to A. thaliana at all 93 sites where A. thaliana and A. arenosa differed, thus showing that A. thaliana is the maternal parent of A. suecica. Under the assumption that A. thaliana and A. arenosa separated 5 million years ago, we estimated a substitution rate of 2.9 x 10(-9) per site per year in noncoding single copy sequence. Within A. thaliana we found 12 substitution (single bp) and eight insertion/deletion (indel) polymorphisms, separating the 25 accessions into 15 haplotypes. Eight of the A. thaliana accessions from central Sweden formed one cluster, which was separated from a cluster consisting of central European and extreme southern Swedish accessions. This latter cluster also included the A. suecica accessions, which were all identical except for one 5 bp indel. We interpret this low level of variation as a strong indication that A. suecica effectively has a single origin, which we dated at 20 000 years ago or more.

Ovarian carcinoma develops through multiple modes of chromosomal evolution

Ovarian carcinoma has the highest mortality of all of the gynecologic cancers. The chromosomal changes in this tumor type are highly complex, and the karyotypes typically show severe aneuploidy. Despite the abundance of cytogenetic information, with approximately 400 published karyotypes, very little is known about the mode of karyotypic evolution and the possible presence of cytogenetic pathways related to tumor development. In the present investigation we used 387 ovarian carcinoma karyotypes to identify the most frequent genomic imbalances. Tumor cases were then classified with respect to the presence or absence of these imbalances and statistically analyzed to assess the order of appearance of chromosomal imbalances, as well as possible karyotypic pathways and cytogenetic subtypes. We establish the temporal order by which the different imbalances occur and show that at least two cytogenetic pathways exist, one characterized by +7, +8q, and +12, and one by 6q- and 1q-. We show that ovarian carcinomas develop through at least three phases of karyotypic evolution. At the early stages, Phase I, the karyotypic evolution seems to proceed though step-wise acquisition of changes. The transition to Phase II showed signs of an increased chromosomal instability, most probably caused by extensive telomere crisis and the onset of breakage-fusion-bridge cycles. This process was linked to the presence of imbalances characteristic for the 6q-/1q- pathway. The transition to Phase III involved triploidization and was also linked to the presence of the 6q-/1q- pathway.

Dissecting karyotypic patterns in colorectal tumors: two distinct but overlapping pathways in the adenoma-carcinoma transition

More than 500 colorectal tumors with clonal chromosomal abnormalities have been reported. Although the pattern of aberrations is nonrandom, no specific primary or secondary karyotypic abnormality has been identified. Also, the chronological order in which the aberrations appear during disease progression is not well known. One reason why our understanding of the cytogenetic evolution is unclear is the high degree of karyotypic complexity seen in these tumors. To overcome some of these difficulties we have previously used several statistical methods that allow identification and interpretation of karyotypic pathways as well as establishment of a temporal order of appearance of the imbalances. These methods were applied on 531 colorectal tumor karyotypes. By using a resampling strategy, 1p-, +7, 7q-, and +12p were identified as early events. Two major and two minor cytogenetic pathways were identified by means of principal component analysis. The two major pathways were initiated with 1p- and +7, and the minor pathways were initiated with +12p and 7q-. The +7/+12p tumors were found to be hyperdiploid, whereas those with 1p-/7q- were pseudodiploid. We also show that the adenoma-carcinoma transition in the 1p- pathway is strongly linked to karyoytypic evolution, whereas the +7 pathway is not, and that the cytogenetic pathways are separated at both early and late stages.

Genetic variation in Arabidopsis suecica and its parental species A. arenosa and A. thaliana

Random amplified polymorphic DNA (RAPD) markers were used to estimate the level of genetic variation in Swedish accessions of the allopolyploid Arabidopsis suecica and its parental species A. thaliana and A. arenosa. The results showed clear differences among the three species with respect to the level of variation. A. arenosa was highly variable, A. thaliana showed a moderate level of variation whereas A. suecica was much less variable than the two other species. An extended analysis covering 19 Swedish populations of A. suecica corroborated the low level of variation in this species, yet 16 unique phenotypes were observed. No isolation by distance was observed. When the genetic variation was partitioned among and within populations of A. suecica, the results showed that the majority of the variation (81%) occurred among populations. This result is interpreted as a strong indication that A. suecica is autogamous in nature.

Coping with complexity. multivariate analysis of tumor karyotypes

Human cancers are characterized by chromosomal aberrations, and an increasing number of specific balanced rearrangements have been found among malignant hematologic disorders. Most solid tumors, however, exhibit a much more complex cytogenetic pattern. Although these chromosome changes show a nonrandom distribution, tumor-specific aberrations are uncommon, and the solid tumors often contain a large number of abnormalities and also display extensive cytogenetic variability. The high level of karyotypic complexity has made a systematic characterization of the chromosomal patterns difficult. In order to better understand the biological relevance of highly abnormal karyotypes in tumor cell populations, novel statistical strategies are needed. We have developed and adapted several methods that may be useful for the evaluation of general patterns of karyotypic complexity, including distribution analysis of cytogenetic imbalances, temporal analysis for time of occurrence of aberrations, and principal component analysis for reconstructing karyotypic pathways. By applying these methods on the chromosomal changes presently known, distinct subgroups have been identified among breast, kidney, bladder, colon, and brain tumors.

Multivariate analysis of chromosomal imbalances in breast cancer delineates cytogenetic pathways and reveals complex relationships among imbalances

More than 550 breast adenocarcinomas with clonal chromosomal abnormalities have been reported. Although the aberration pattern is clearly nonrandom, no specific primary or secondary karyotypic abnormality has been identified, and furthermore the chronological order in which the aberrations appear during disease progression is not well known. The high degree of karyotypic complexity in epithelial tumors such as breast cancer is one reason why our understanding of the sequential order of cytogenetic evolution is unclear. To overcome some of these difficulties, we have used several statistical methods that allow identification and interpretation of karyotypic pathways. These methods were applied on 538 breast cancer karyotypes. The distribution of the number of imbalances/tumor showed a monomodal appearance, indicating that one single mode of karyotypic evolution is operating in this tumor type. We show that there exists a temporal order with respect to the appearance of chromosomal imbalances. The imbalances +1pq, 1q-, 3p-, and +7 appear earlier than expected from random events, and two cytogenetic pathways, one initiated by +1q and followed by 11q- and -22, the other initiated by either 3p- or 1q- and followed by 1p-, 3q-, and 6q-, can be discerned. We also show that +7 and +8q behave independently of the other imbalances and cannot, by simple means, be incorporated in the identified pathway scheme. Although the cytogenetic pathways are well separated at earlier stages, they later converge and include a common set of late imbalances.