Longitudinal copy number, whole exome and targeted deep sequencing of 'good risk' IGHV-mutated CLL patients with progressive disease.

The biological features of IGHV-M chronic lymphocytic leukemia responsible for disease progression are still poorly understood. We undertook a longitudinal study close to diagnosis, pre-treatment and post relapse in 13 patients presenting with cMBL or Stage A disease and good-risk biomarkers (IGHV-M genes, no del(17p) or del(11q) and low CD38 expression) who nevertheless developed progressive disease, of whom 10 have required therapy. Using cytogenetics, fluorescence in situ hybridisation, genome-wide DNA methylation and copy number analysis together with whole exome, targeted deep- and Sanger sequencing at diagnosis, we identified mutations in established chronic lymphocytic leukemia driver genes in nine patients (69%), non-coding mutations (PAX5 enhancer region) in three patients and genomic complexity in two patients. Branching evolutionary trajectories predominated (n=9/13), revealing intra-tumoural epi- and genetic heterogeneity and sub-clonal competition before therapy. Of the patients subsequently requiring treatment, two had sub-clonal TP53 mutations that would not be detected by standard methodologies, three qualified for the very-low-risk category defined by integrated mutational and cytogenetic analysis and yet had established or putative driver mutations and one patient developed progressive, therapy-refractory disease associated with the emergence of an IGHV-U clone. These data suggest that extended genomic and immunogenetic screening may have clinical utility in patients with apparent good-risk disease.

Profound parental bias associated with chromosome 14 acquired uniparental disomy indicates targeting of an imprinted locus.

Acquired uniparental disomy (aUPD) is a common finding in myeloid malignancies and typically acts to convert a somatically acquired heterozygous mutation to homozygosity. We sought to identify the target of chromosome 14 aUPD (aUPD14), a recurrent abnormality in myeloid neoplasms and population cohorts of elderly individuals. We identified 29 cases with aUPD14q that defined a minimal affected region (MAR) of 11.2 Mb running from 14q32.12 to the telomere. Exome sequencing (n=7) did not identify recurrently mutated genes, but methylation-specific PCR at the imprinted MEG3-DLK1 locus located within the MAR demonstrated loss of maternal chromosome 14 and gain of paternal chromosome 14 (P<0.0001), with the degree of methylation imbalance correlating with the level of aUPD (r=0.76; P=0.0001). The absence of driver gene mutations in the exomes of three individuals with aUPD14q but no known haematological disorder suggests that aUPD14q may be sufficient to drive clonal haemopoiesis. Analysis of cases with both aUPD14q and JAK2 V617F (n=11) indicated that aUPD14q may be an early event in some cases but a late event in others. We conclude that aUPD14q is a recurrent abnormality that targets an imprinted locus and may promote clonal haemopoiesis either as an initiating event or as a secondary change.

Ethnicity and outcome of young breast cancer patients in the United Kingdom: the POSH study.

BACKGROUND: Black ethnic groups have a higher breast cancer mortality than Whites. American studies have identified variations in tumour biology and unequal health-care access as causative factors. We compared tumour pathology, treatment and outcomes in three ethnic groups in young breast cancer patients treated in the United Kingdom. METHODS: Women aged ≤ 40 years at breast cancer diagnosis were recruited to the POSH national cohort study (MREC: 00/06/69). Personal characteristics, tumour pathology and treatment data were collected at diagnosis. Follow-up data were collected annually. Overall survival (OS) and distant relapse-free survival (DRFS) were assessed using Kaplan-Meier curves, and multivariate analyses were performed using Cox regression. RESULTS: Ethnicity data were available for 2915 patients including 2690 (91.0%) Whites, 118 (4.0%) Blacks and 87 (2.9%) Asians. Median tumour diameter at presentation was greater in Blacks than Whites (26.0 mm vs 22.0 mm, P=0.0103), and multifocal tumours were more frequent in both Blacks (43.4%) and Asians (37.0%) than Whites (28.9%). ER/PR/HER2-negative tumours were significantly more frequent in Blacks (26.1%) than Whites (18.6%, P=0.043). Use of chemotherapy was similarly high in all ethnic groups (89% B vs 88.6% W vs 89.7% A). A 5-year DRFS was significantly lower in Blacks than Asians (62.8% B vs 77.0% A, P=0.0473) or Whites (62.8 B% vs 77.0% W, P=0.0053) and a 5-year OS for Black patients, 71.1% (95% CI: 61.0-79.1%), was significantly lower than that of Whites (82.4%, 95% CI: 80.8-83.9%, W vs B: P=0.0160). In multivariate analysis, Black ethnicity had an effect on DRFS in oestrogen receptor (ER)-positive patients that is independent of body mass index, tumour size, grade or nodal status, HR: 1.60 (95% CI: 1.03-2.47, P=0.035). CONCLUSION: Despite equal access to health care, young Black women in the United Kingdom have a significantly poorer outcome than White patients. Black ethnicity is an independent risk factor for reduced DRFS particularly in ER-positive patients.

A novel prognostic model in myeloma based on co-segregating adverse FISH lesions and the ISS: analysis of patients treated in the MRC Myeloma IX trial.

The association of genetic lesions detected by fluorescence in situ hybridization (FISH) with survival was analyzed in 1069 patients with newly presenting myeloma treated in the Medical Research Council Myeloma IX trial, with the aim of identifying patients associated with the worst prognosis. A comprehensive FISH panel was performed, and the lesions associated with short progression-free survival and overall survival (OS) in multivariate analysis were +1q21, del(17p13) and an adverse immunoglobulin heavy chain gene (IGH) translocation group incorporating t(4;14), t(14;16) and t(14;20). These lesions frequently co-segregated, and there was an association between the accumulation of these adverse FISH lesions and a progressive impairment of survival. This observation was used to define a series of risk groups based on number of adverse lesions. Taking this approach, we defined a favorable risk group by the absence of adverse genetic lesions, an intermediate group with one adverse lesion and a high-risk group defined by the co-segregation of >1 adverse lesion. This genetic grouping was independent of the International Staging System (ISS) and so was integrated with the ISS to identify an ultra-high-risk group defined by ISS II or III and >1 adverse lesion. This group constituted 13.8% of patients and was associated with a median OS of 19.4 months.

A map of the human genome in linkage disequilibrium units.

Two genetic maps with additive distances contribute information about recombination patterns, recombinogenic sequences, and discovery of genes affecting a particular phenotype. Recombination is measured in morgans (w) over a single generation in a linkage map but may cover thousands of generations in a linkage disequilibrium (LD) map measured in LD units (LDU). We used a subset of single nucleotide polymorphisms from the HapMap Project to create a genome-wide map in LDU. Recombination accounts for 96.8% of the LDU variance in chromosome arms and 92.4% in their deciles. However, deeper analysis shows that LDU/w, an estimate of the effective bottleneck time (t), is significantly variable among chromosome arms because (i) the linkage map is approximated from the Haldane function, then adjusted toward the Kosambi function that is more accurate but still exaggerates w for all chromosomes, especially shorter ones; (ii) the non-pseudoautosomal region of the X chromosome is subject to hemizygous selection; and (iii) at resolution less than approximately 40,000 markers per w, there are indeterminacies (holes) in the LD map reflecting intervals of very high recombination. Selection and stochastic variation in small regions must have effects, which remain to be investigated by comparisons among populations. These considerations suggest an optimal strategy to eliminate holes quickly, greatly enhance the resolution of sex-specific linkage maps, and maximize the gain in association mapping by using LD maps.

A metric linkage disequilibrium map of a human chromosome.

We used LDMAP (Maniatis et al. 2002) to analyse SNP data spanning chromosome 22 (Dawson et al. 2002), to obtain a whole-chromosome metric LD map. The LD map, with map distances analogous to the centiMorgan scale of linkage maps, identifies regions of high LD as plateaus ('blocks') and characterises steps which define the relationship between these regions. From this map we estimate that block regions comprise between 32% and 55% of the euchromatic portion of chromosome 22 and that increasing marker density within steps may increase block coverage. Steps are regions of low LD which correspond to areas of variable recombination intensity. The intensity of recombination is related to the height of the step and thus intense recombination hot-spots can be distinguished from more randomly distributed historical events. The LD maps are more closely related to the high-resolution linkage map (Kong et al. 2002) than average measures of rho with recombination accounting for between 34% and 52% of the variance in patterns of LD (r=0.58 - 0.71, p=0.0001). Step regions are closely correlated with a range of sequence motifs including GT/CA repeats. The LD map identifies holes in which greater marker density is required and defines the optimal SNP spacing for positional cloning, which suggests that some multiple of around 50,000 SNPs will be required to efficiently screen Caucasian genomes. Further analyses which investigate selection of informative SNPs and the effect of SNP allele frequency and marker density will refine this estimate.

Recombination, interference and sequence: comparison of chromosomes 21 and 22.

The euchromatic regions of chromosomes 21 and 22 are almost completely sequenced and have similar lengths (33.7-34.6 Mb). This similarity effectively controls for the influence of length, making comparisons of recombination and interference interesting. For both chromosomes, there is less male than female recombination, and male recombination is associated with GT/CA repeats. The striking sex difference may result from greater condensation of chromosomes in paternal meiosis, possibly restricting recombination to regions with longer repeat tracts and/or higher repeat densities. Chiasma interference in both sexes for chromosome 22 and in females for chromosome 21 is close to the genome average. Chromosome 21 is significantly different in male meiosis, with near complete interference, suggesting that even when double recombinants occur they are widely spaced. We propose that this difference is related to the different distribution of GT/CA dinucleotides. These repeats are widely distributed on chromosome 22, perhaps offering greater opportunities for double recombinants to occur within smaller regions, whereas they are largely subtelomeric in distribution on chromosome 21.

The first linkage disequilibrium (LD) maps: delineation of hot and cold blocks by diplotype analysis.

Linkage disequilibrium (LD) provides information about positional cloning, linkage, and evolution that cannot be inferred from other evidence, even when a correct sequence and a linkage map based on more than a handful of families become available. We present theory to construct an LD map for which distances are additive and population-specific maps are expected to be approximately proportional. For this purpose, there is only a modest difference in relative efficiency of haplotypes and diplotypes: resolving the latter into 2-locus haplotypes has significant cost or error and increases information by about 50%. LD maps for a cold spot in 19p13.3 and a more typical region in 3q21 are optimized by interval estimates. For a random sample and trustworthy map the value of LD at large distance can be predicted reliably from information over a small distance and does not depend on the evolutionary variance unless the sample size approaches the population size. Values of the association probability that can be distinguished from the value at large distance are determined not by population size but by time since a critical bottleneck. In these examples, omission of markers with significant Hardy-Weinberg disequilibrium does not improve the map, and widely discrepant draft sequences have similar estimates of the genetic parameters. The LD cold spot in 19p13.3 gives an unusually high estimate of time, supporting an argument that this relationship is general. As predicted for a region with ancient haplotypes or uniformly high recombination, there is no clear evidence of LD clustering. On the contrary, the 3q21 region is resolved into alternating blocks of stable and decreasing LD, as expected from crossover clustering. Construction of a genomewide LD map requires data not yet available, which may be complemented but not replaced by a catalog of haplotypes.

A tournament of linkage tests in complex inheritance.

The performance of some weakly parametric linkage tests in common use was compared on 200 replicates of oligogenic inheritance from Genetic Analysis Workshop 10. Each random sample for the quantitative trait was dichotomized at different thresholds and also selected through 2 affected sibs, generating 8 combinations of sample and variable. The variance component program SOLAR performed best with a continuous trait, even in selected samples, when the population mean was used. The sib-pair program SIBPAL2 was best in most other cases when the phenotype product, population mean, and empirical estimates of pair correlations were used. The BETA program that introduced phenotype products was slightly more powerful than maximum likelihood scores under the null hypothesis and approached but did not exceed SIBPAL2 under its optimal conditions. Type I errors generally exceeded expectations from a chi(2) test, but were conservative with respect to bounds on lods. All methods can be improved by use of the population mean, empirical correlations, logistic representation for affection status, and correct lods for samples that favour the null hypothesis. It remains uncertain whether all information can be extracted by weakly parametric methods and whether correction for ascertainment bias demands a strongly parametric model. Performance on a standard set of simulated data is indispensable for recognising optimal methods.

A sequence-based integrated map of chromosome 22.

The near-completion of the sequence for chromosome 22q revolutionizes map integration. We describe a sequence-based integrated map containing 968 loci including 516 known or predicted gene sequences, 317 STSs not included in these sequences, and 135 nonexpressed multinucleotide polymorphisms. The published sequence spans 34.6 Mb, inclusive of gaps estimated to total 1.1 Mb, compared with a top-down estimate of 43 Mb. This discrepancy is discussed, but will not be resolved until more of the genome is analyzed. The radiation hybrid map has 5% error in order and 34% error in location exceeding 1 Mb. The utility of a composite location based on evidence other than sequence is limited to regions not yet sequenced. A genetic map conditional on sequence order was constructed from pairwise lods. Its length of 74.8 cM in males and 80.2 cM in females is slightly less than the previous estimate not constrained by sequence order. Five recombination hot spots are detected, with differences in location between the sexes. Male recombination correlates with repetitive DNA, whereas female recombination does not. It remains to be seen whether this is true for other human chromosomes. An algorithm to improve the fit of cytogenetic bands sequence location reduces the discrepancies in cytogenetic assignment from 61 to 38. This sequence-based integrated map is represented in the genetic location database (LDB2000), which is available at http://cedar.genetics.soton.ac.uk/public_html/LDB2000.html.

LDB2000: sequence-based integrated maps of the human genome.

MOTIVATION: Integrated maps are useful for gene mapping and establishing the relationship between recombination and sequence. In this paper we describe algorithms and their implementation for constructing sequence-based integrated maps of the human chromosomes, which are presented in LDB2000, a web based resource. Gene mapping efforts are now focussing on linkage disequilibrium mapping and extension of the integrated map to represent the extent of linkage disequilibrium in different genomic regions would further increase the utility of these maps. RESULTS: Sequence-based integrated maps have been completed for chromosomes 21 and 22. These maps provide locations for genes and polymorphic markers in sequence and on genetic linkage, radiation hybrid and cytogenetic scales. Single nucleotide polymorphisms associated with genes in the maps are also included and their sequence locations indicated. Related locus information, such as aliases and expression information, can be searched on the WWW site.

Meta-analysis and retrospective collaboration: two methods to map oligogenes for atopy and asthma.

Combination of evidence over samples, each of which is too small to be conclusive, is the central problem in complex inheritance. There are three approaches: meta-analysis, prospective collaboration, and retrospective collaboration. Our experience with the first and last, which are the most cost-effective, is discussed.