Genomic characterisation of hormone receptor-positive breast cancer arising in very young women.

BACKGROUND: Very young premenopausal women diagnosed with hormone receptor-positive, human epidermal growth factor receptor 2-negative (HR+HER2-) early breast cancer (EBC) have higher rates of recurrence and death for reasons that remain largely unexplained. PATIENTS AND METHODS: Genomic sequencing was applied to HR+HER2- tumours from patients enrolled in the Suppression of Ovarian Function Trial (SOFT) to determine genomic drivers that are enriched in young premenopausal women. Genomic alterations were characterised using next-generation sequencing from a subset of 1276 patients (deep targeted sequencing, n = 1258; whole-exome sequencing in a young-age, case-control subsample, n = 82). We defined copy number (CN) subgroups and assessed for features suggestive of homologous recombination deficiency (HRD). Genomic alteration frequencies were compared between young premenopausal women (<40 years) and older premenopausal women (≥40 years), and assessed for associations with distant recurrence-free interval (DRFI) and overall survival (OS). RESULTS: Younger women (<40 years, n = 359) compared with older women (≥40 years, n = 917) had significantly higher frequencies of mutations in GATA3 (19% versus 16%) and CN amplifications (CNAs) (47% versus 26%), but significantly lower frequencies of mutations in PIK3CA (32% versus 47%), CDH1 (3% versus 9%), and MAP3K1 (7% versus 12%). Additionally, they had significantly higher frequencies of features suggestive of HRD (27% versus 21%) and a higher proportion of PIK3CA mutations with concurrent CNAs (23% versus 11%). Genomic features suggestive of HRD, PIK3CA mutations with CNAs, and CNAs were associated with significantly worse DRFI and OS compared with those without these features. These poor prognostic features were enriched in younger patients: present in 72% of patients aged <35 years, 54% aged 35-39 years, and 40% aged ≥40 years. Poor prognostic features [n = 584 (46%)] versus none [n = 692 (54%)] had an 8-year DRFI of 84% versus 94% and OS of 88% versus 96%. Younger women (<40 years) had the poorest outcomes: 8-year DRFI 74% versus 85% and OS 80% versus 93%, respectively. CONCLUSION: These results provide insights into genomic alterations that are enriched in young women with HR+HER2- EBC, provide rationale for genomic subgrouping, and highlight priority molecular targets for future clinical trials.

Determining epithelial contribution to in vivo mesenchymal tumour expression signature using species-specific microarray profiling analysis of xenografts.

Gene expression profiling using microarrays and xenograft transplants of human cancer cell lines are both popular tools to investigate human cancer. However, the undefined degree of cross hybridization between the mouse and human genomes hinders the use of microarrays to characterize gene expression of both the host and the cancer cell within the xenograft. Since an increasingly recognized aspect of cancer is the host response (or cancer-stroma interaction), we describe here a bioinformatic manipulation of the Affymetrix profiling that allows interrogation of the gene expression of both the mouse host and the human tumour. Evidence of microenvironmental regulation of epithelial mesenchymal transition of the tumour component in vivo is resolved against a background of mesenchymal gene expression. This tool could allow deeper insight to the mechanism of action of anti-cancer drugs, as typically novel drug efficacy is being tested in xenograft systems.

Estimation and assessment of raw copy numbers at the single locus level.

MOTIVATION: Although copy-number aberrations are known to contribute to the diversity of the human DNA and cause various diseases, many aberrations and their phenotypes are still to be explored. The recent development of single-nucleotide polymorphism (SNP) arrays provides researchers with tools for calling genotypes and identifying chromosomal aberrations at an order-of-magnitude greater resolution than possible a few years ago. The fundamental problem in array-based copy-number (CN) analysis is to obtain CN estimates at a single-locus resolution with high accuracy and precision such that downstream segmentation methods are more likely to succeed. RESULTS: We propose a preprocessing method for estimating raw CNs from Affymetrix SNP arrays. Its core utilizes a multichip probe-level model analogous to that for high-density oligonucleotide expression arrays. We extend this model by adding an adjustment for sequence-specific allelic imbalances such as cross-hybridization between allele A and allele B probes. We focus on total CN estimates, which allows us to further constrain the probe-level model to increase the signal-to-noise ratio of CN estimates. Further improvement is obtained by controlling for PCR effects. Each part of the model is fitted robustly. The performance is assessed by quantifying how well raw CNs alone differentiate between one and two copies on Chromosome X (ChrX) at a single-locus resolution (27kb) up to a 200kb resolution. The evaluation is done with publicly available HapMap data. AVAILABILITY: The proposed method is available as part of an open-source R package named aroma.affymetrix. Because it is a bounded-memory algorithm, any number of arrays can be analyzed.

FIRMA: a method for detection of alternative splicing from exon array data.

MOTIVATION: Analyses of EST data show that alternative splicing is much more widespread than once thought. The advent of exon and tiling microarrays means that researchers now have the capacity to experimentally measure alternative splicing on a genome wide level. New methods are needed to analyze the data from these arrays. RESULTS: We present a method, finding isoforms using robust multichip analysis (FIRMA), for detecting differential alternative splicing in exon array data. FIRMA has been developed for Affymetrix exon arrays, but could in principle be extended to other exon arrays, tiling arrays or splice junction arrays. We have evaluated the method using simulated data, and have also applied it to two datasets: a panel of 11 human tissues and a set of 10 pairs of matched normal and tumor colon tissue. FIRMA is able to detect exons in several genes confirmed by reverse transcriptase PCR. AVAILABILITY: R code implementing our methods is contributed to the package aroma.affymetrix.

Replication of KIAA0350, IL2RA, RPL5 and CD58 as multiple sclerosis susceptibility genes in Australians.

A recent genome-wide association study (GWAS) conducted by the International Multiple Sclerosis Genetics Consortium (IMSGC) identified a number of putative MS susceptibility genes. Here we have performed a replication study in 1134 Australian MS cases and 1265 controls for 17 risk-associated single nucleotide polymorphisms (SNPs) reported by the IMSGC. Of 16 SNPs that passed quality control filters, four, each corresponding to a different non-human leukocyte antigen (HLA) gene, were associated with disease susceptibility: KIAA0350 (rs6498169) P=0.001, IL2RA (rs2104286) P=0.033, RPL5 (rs6604026) P=0.041 and CD58 (rs12044852) P=0.042. There was no association (P=0.58) between rs6897932 in the IL7R gene and the risk of MS. No interactions were detected between the replicated IMSGC SNPs and HLA-DRB1*15, gender, disease course, disease progression or age-at-onset. We used a novel Bayesian approach to estimate the extent to which our data increased or decreased evidence for association with the six most-associated IMSGC loci. These analyses indicated that even modest P-values, such as those reported here, can contribute markedly to the posterior probability of 'true' association in replication studies. In conclusion, these data provide support for the involvement of four non-HLA genes in the pathogenesis of MS, and combined with previous data, increase to genome-wide significance (P=3 x 10(-8)) evidence of an association between KIAA0350 and risk of disease.

A genetic screen for behavioral mutations that perturb dopaminergic homeostasis in mice.

Disruption of dopaminergic (DA) systems is thought to play a central role in the addictive process and in the pathophysiology of schizophrenia. Although inheritance plays an important role in the predisposition to these disorders, the genetic basis of this is not well understood. To provide additional insight, we have performed a modifier screen in mice designed to identify mutations that perturb DA homeostasis. With a genetic background sensitized by a mutation in the dopamine transporter (DAT), we used random chemical mutagenesis and screened for mutant mice with locomotor abnormalities. Four mutant lines were identified with quantitatively elevated levels of locomotor activity. Mapping of mutations in these lines identified two loci that alter activity only when dopamine levels are elevated by a DAT mutation and thus would only have been uncovered by this type of approach. One of these quantitative trait loci behaves as an enhancer of DA neurotransmission, whereas the other may act as a suppressor. In addition, we also identified three loci which are not dependent on the sensitized background but which also contribute to the overall locomotor phenotype.

A score test for the linkage analysis of qualitative and quantitative traits based on identity by descent data from sib-pairs.

We propose a general likelihood-based approach to the linkage analysis of qualitative and quantitative traits using identity by descent (IBD) data from sib-pairs. We consider the likelihood of IBD data conditional on phenotypes and test the null hypothesis of no linkage between a marker locus and a gene influencing the trait using a score test in the recombination fraction theta between the two loci. This method unifies the linkage analysis of qualitative and quantitative traits into a single inferential framework, yielding a simple and intuitive test statistic. Conditioning on phenotypes avoids unrealistic random sampling assumptions and allows sib-pairs from differing ascertainment mechanisms to be incorporated into a single likelihood analysis. In particular, it allows the selection of sib-pairs based on their trait values and the analysis of only those pairs having the most informative phenotypes. The score test is based on the full likelihood, i.e. the likelihood based on all phenotype data rather than just differences of sib-pair phenotypes. Considering only phenotype differences, as in Haseman and Elston (1972) and Kruglyak and Lander (1995), may result in important losses in power. The linkage score test is derived under general genetic models for the trait, which may include multiple unlinked genes. Population genetic assumptions, such as random mating or linkage equilibrium at the trait loci, are not required. This score test is thus particularly promising for the analysis of complex human traits. The score statistic readily extends to accommodate incomplete IBD data at the test locus, by using the hidden Markov model implemented in the programs MAPMAKER/SIBS and GENEHUNTER (Kruglyak and Lander, 1995; Kruglyak et al., 1996). Preliminary simulation studies indicate that the linkage score test generally matches or outperforms the Haseman-Elston test, the largest gains in power being for selected samples of sib-pairs with extreme phenotypes.

Statistical analysis of ordered tetrads.

Ordered tetrad data yield information on chromatid interference, chiasma interference, and centromere locations. In this article, we show that the assumption of no chromatid interference imposes certain constraints on multilocus ordered tetrad probabilities. Assuming no chromatid interference, these constraints can be used to order markers under general chiasma processes. We also derive multilocus tetrad probabilities under a class of chiasma interference models, the chi-square models. Finally, we compare centromere map functions under the chi-square models with map functions proposed in the literature. Results in this article can be applied to order genetic markers and map centromeres using multilocus ordered tetrad data.

Statistical analysis of half-tetrads.

Half-tetrads, where two meiotic products from a single meiosis are recovered together, arise in different forms in a variety of organisms. Closely related to ordered tetrads, half-tetrads yield information on chromatid interference, chiasma interference, and centromere positions. In this article, for different half-tetrad types and different marker configurations, we derive the relations between multilocus half-tetrad probabilities and multilocus ordered tetrad probabilities. These relations are used to obtain equality and inequality constraints among multilocus half-tetrad probabilities that are imposed by the assumption of no chromatid interference. We illustrate how to apply these results to study chiasma interference and to map centromeres using multilocus half-tetrad data.

Modeling interference in genetic recombination.

In analyzing genetic linkage data it is common to assume that the locations of crossovers along a chromosome follow a Poisson process, whereas it has long been known that this assumption does not fit the data. In many organisms it appears that the presence of a crossover inhibits the formation of another nearby, a phenomenon known as "interference." We discuss several point process models for recombination that incorporate position interference but assume no chromatid interference. Using stochastic simulation, we are able to fit the models to a multilocus Drosophila dataset by the method of maximum likelihood. We find that some biologically inspired point process models incorporating one or two additional parameters provide a dramatically better fit to the data than the usual "no-interference" Poisson model.

On genetic map functions.

Various genetic map functions have been proposed to infer the unobservable genetic distance between two loci from the observable recombination fraction between them. Some map functions were found to fit data better than others. When there are more than three markers, multilocus recombination probabilities cannot be uniquely determined by the defining property of map functions, and different methods have been proposed to permit the use of map functions to analyze multilocus data. If for a given map function, there is a probability model for recombination that can give rise to it, then joint recombination probabilities can be deduced from this model. This provides another way to use map functions in multilocus analysis. In this paper we show that stationary renewal processes give rise to most of the map functions in the literature. Furthermore, we show that the interevent distributions of these renewal processes can all be approximated quite well by gamma distributions.

Statistical analysis of crossover interference using the chi-square model.

The chi-square model (also known as the gamma model with integer shape parameter) for the occurrence of crossovers along a chromosome was first proposed in the 1940's as a description of interference that was mathematically tractable but without biological basis. Recently, the chi-square model has been reintroduced into the literature from a biological perspective. It arises as a result of certain hypothesized constraints on the resolution of randomly distributed crossover intermediates. In this paper under the assumption of no chromatid interference, the probability for any single spore or tetrad joint recombination pattern is derived under the chi-square model. The method of maximum likelihood is then used to estimate the chi-square parameter m and genetic distances among marker loci. We discuss how to interpret the goodness-of-fit statistics appropriately when there are some recombination classes that have only a small number of observations. Finally, comparisons are made between the chi-square model and some other tractable models in the literature.

Statistical analysis of chromatid interference.

The nonrandom occurrence of crossovers along a single strand during meiosis can be caused by either chromatid interference, crossover interference or both. Although crossover interference has been consistently observed in almost all organisms since the time of the first linkage studies, chromatid interference has not been as thoroughly discussed in the literature, and the evidence provided for it is inconsistent. In this paper with virtually no restrictions on the nature of crossover interference, we describe the constraints that follow from the assumption of no chromatid interference for single spore data. These constraints are necessary consequences of the assumption of no chromatid interference, but their satisfaction is not sufficient to guarantee no chromatid interference. Models can be constructed in which chromatid interference clearly exists but is not detectable with single spore data. We then extend our analysis to cover tetrad data, which permits more powerful tests of no chromatid interference. We note that the traditional test of no chromatid interference based on tetrad data does not make full use of the information provided by the data, and we offer a statistical procedure for testing the no chromatid interference constraints that does make full use of the data. The procedure is then applied to data from several organisms. Although no strong evidence of chromatid interference is found, we do observe an excess of two-strand double recombinations, i.e., negative chromatid interference.

A comparison of match-only algorithms for the analysis of Plasmodium falciparum oligonucleotide arrays.

This study is motivated by two data sets which employ a custom Plasmodium falciparum version of the Affymetrix GeneChip, containing only perfect match (PM) oligonucleotides. A PM-only chip cannot be analysed using the standard Affymetrix-supplied software. We compared the performance of three match-only algorithms on these data: the Match Only Integral Distribution (MOID) algorithm, Robust Multichip Analysis (RMA), and the Model Based Expression Index (MBEI). We validated the differential expression of several genes using quantitative reverse transcriptase-PCR. We also performed a comparison using two publicly available 'benchmarking' data sets: the Latin Square spike-in data set generated by Affymetrix, and the Gene Logic dilution series. Since we know what the true fold changes are in these special data sets, they are helpful for assessment of expression algorithms.

Factors associated with human immunodeficiency virus seroconversion in homosexual men in three San Francisco cohort studies, 1984-1989.

A total of 83 HIV seroconversions occurred between 1984 and 1989 in three San Francisco cohorts of homosexual and bisexual men. A nested case-control analysis was performed to assess the risk of seroconversion associated with sexual practices. Strong associations were found with total number of intercourse partners and receptive anal intercourse. Weaker, but significant associations were found with receptive oral intercourse. Individuals reporting condom use with some partners were more likely to become infected than those reporting no condom use. Some of this difference was due to increased numbers of sexual partners among those reporting some condom use, but the association remained significant in multivariate analyses. This observation indicates there is some unmeasured risk associated with those using condoms, such as more HIV seropositive partners.

Phat--a gene finding program for Plasmodium falciparum.

We describe and assess the performance of the gene finding program pretty handy annotation tool (Phat) on sequence from the malaria parasite Plasmodium falciparum. Phat is based on a generalized hidden Markov model (GHMM) similar to the models used in GENSCAN, Genie and HMMgene. In a test set of 44 confirmed gene structures Phat achieves nucleotide-level sensitivity and specificity of greater than 95%, performing as well as the other P. falciparum gene finding programs Hexamer and GlimmerM. Phat is particularly useful for P. falciparum and other eukaryotes for which there are few gene finding programs available as it is distributed with code for retraining it on new organisms. Moreover, the full source code is freely available under the GNU General Public License, allowing for users to further develop and customize it.

An algorithm for haplotype analysis.

This paper proposes an algorithm for haplotype analysis based on a Monte Carlo method. Haplotype configurations are generated according to the distribution of joint haplotypes of individuals in a pedigree given their phenotype data, via a Markov chain Monte Carlo algorithm. The haplotype configuration which maximizes this conditional probability distribution can thus be estimated. In addition, the set of haplotype configurations with relatively high probabilities can also be estimated as possible alternatives to the most probable one. This flexibility enables geneticists to choose the haplotype configurations which are most reasonable to them, allowing them to include their knowledge of the data under analysis.

DNA sequencing with transposons.

The use of transposons offers the possibility of a directed approach to DNA sequencing, where a target DNA up to about 6kb in length can be sequenced quickly and with minimal redundancy. Transposons are mobile DNA elements which can be inserted in a reasonably random fashion into the target DNA. An important part of this process is the location of the transposon insertions (known as mapping) and the selection of a sensible subset of transposons to use as priming sites for sequencing reactions. This paper presents a probabilistic method of scoring selected subsets of transposons and a graph-theoretic algorithm for selection of a subset of maximal score.

Over- and underrepresentation of short DNA words in herpesvirus genomes.

The relative abundance and rarity of DNA words have been recognized in previous biological studies to have implications for the regulation, repair, and evolutionary mechanisms of a genome. In this paper, we review several different measures of abundance and rarity of DNA words, including z-scores, representation ratios, and cross-ratios, that have appeared in the recent literature, and examine the concordance among them using the human cytomegalovirus genome sequence. We then rank all words of length k = 2, ..., 5 of seven herpesvirus genomes according to their abundance, as measured by one of the z-scores based upon a stationary Markov model of order k-2. Using a simple metric on the ranks of 2-words of the seven herpesvirus sequences, we construct an evolutionary tree. Several 3-words are observed to be consistently over- or underrepresented in all seven herpesviruses. Furthermore, clusters of some of the most over- and underrepresented 4- and 5-words in the genomes are identified with functional sites such as the origins of replication and regulatory signals of individual viruses.

Alveolar lining layer is thin and continuous: low-temperature scanning electron microscopy of rat lung.

The low-temperature electron microscope, which preserves aqueous structures as solid water at liquid nitrogen temperature, was used to image the alveolar lining layer, including surfactant and its aqueous subphase, of air-filled lungs frozen in anesthetized rats at 15-cmH2O transpulmonary pressure. Lining layer thickness was measured on cross fractures of walls of the outermost subpleural alveoli that could be solidified with metal mirror cryofixation at rates sufficient to limit ice crystal growth to 10 nm and prevent appreciable water movement. The thickness of the liquid layer averaged 0.14 micron over relatively flat portions of the alveolar walls, 0.89 micron at the alveolar wall junctions, and 0.09 micron over the protruding features (9 rats, 20 walls, 16 junctions, and 146 areas), for an area-weighted average thickness of 0.2 micron. The alveolar lining layer appears continuous, submerging epithelial cell microvilli and intercellular junctional ridges; varies from a few nanometers to several micrometers in thickness, and serves to smooth the alveolar air-liquid interface in lungs inflated to zone 1 or 2 conditions.