Broadcasters, receivers, functional groups of metabolites and the link to heart failure progression using polygenic factors

In a prospective study with records of heart failure (HF) incidence, we present metabolite profiling data from individuals without HF at baseline. We uncovered the interconnectivity of metabolites using data-driven and causal networks augmented with polygenic factors. Exploring the networks, we identified metabolite broadcasters, receivers, mediators, and subnetworks corresponding to functional classes of metabolites, and provided insights into the link between metabolomic architecture and regulation in health. We incorporated the network structure into the identification of metabolites associated with HF to control the effect of confounding metabolites. We identified metabolites associated with higher or lower risk of HF incidence, the associations that were not confounded by the other metabolites, such as glycine, ureidopropionic and glycocholic acids, and LPC 18:2. We revealed the underlying relationships of the findings. For example, asparagine directly influenced glycine, and both were inversely associated with HF. These two metabolites were influenced by polygenic factors and only essential amino acids which are not synthesized in the human body and come directly from the diet. Metabolites may play a critical role in linking genetic background and lifestyle factors to HF progression. Revealing the underlying connectivity of metabolites associated with HF strengthens the findings and facilitates a mechanistic understanding of HF progression.

Quantifying factors that affect polygenic risk score performance across diverse ancestries and age groups for body mass index

Polygenic risk scores (PRS) have led to enthusiasm for precision medicine. However, it is well documented that PRS do not generalize across groups differing in ancestry or sample characteristics e.g., age. Quantifying performance of PRS across different groups of study participants, using genome-wide association study (GWAS) summary statistics from multiple ancestry groups and sample sizes, and using different linkage disequilibrium (LD) reference panels may clarify which factors are limiting PRS transferability. To evaluate these factors in the PRS generation process, we generated body mass index (BMI) PRS (PRSBMI) in the Electronic Medical Records and Genomics (eMERGE) network (N=75,661). Analyses were conducted in two ancestry groups (European and African) and three age ranges (adult, teenagers, and children). For PRSBMI calculations, we evaluated five LD reference panels and three sets of GWAS summary statistics of varying sample size and ancestry. PRSBMI performance increased for both African and European ancestry individuals using cross-ancestry GWAS summary statistics compared to European-only summary statistics (6.3% and 3.7% relative R2 increase, respectively, pAfrican=0.038, pEuropean=6.26x10-4). The effects of LD reference panels were more pronounced in African ancestry study datasets. PRSBMI performance degraded in children; R2 was less than half of teenagers or adults. The effect of GWAS summary statistics sample size was small when modeled with the other factors. Additionally, the potential of using a PRS generated for one trait to predict risk for comorbid diseases is not well understood especially in the context of cross-ancestry analyses - we explored clinical comorbidities from the electronic health record associated with PRSBMI and identified significant associations with type 2 diabetes and coronary atherosclerosis. In summary, this study quantifies the effects that ancestry, GWAS summary statistic sample size, and LD reference panel have on PRS performance, especially in cross-ancestry and age-specific analyses.

Development and validation of a trans-ancestry polygenic risk score for type 2 diabetes in diverse populations

BACKGROUND

Type 2 diabetes (T2D) is a worldwide scourge caused by both genetic and environmental risk factors that disproportionately afflicts communities of color. Leveraging existing large-scale genome-wide association studies (GWAS), polygenic risk scores (PRS) have shown promise to complement established clinical risk factors and intervention paradigms, and improve early diagnosis and prevention of T2D. However, to date, T2D PRS have been most widely developed and validated in individuals of European descent. Comprehensive assessment of T2D PRS in non-European populations is critical for equitable deployment of PRS to clinical practice that benefits global populations.

METHODS

We integrated T2D GWAS in European, African, and East Asian populations to construct a trans-ancestry T2D PRS using a newly developed Bayesian polygenic modeling method, and assessed the prediction accuracy of the PRS in the multi-ethnic Electronic Medical Records and Genomics (eMERGE) study (11,945 cases; 57,694 controls), four Black cohorts (5137 cases; 9657 controls), and the Taiwan Biobank (4570 cases; 84,996 controls). We additionally evaluated a post hoc ancestry adjustment method that can express the polygenic risk on the same scale across ancestrally diverse individuals and facilitate the clinical implementation of the PRS in prospective cohorts.

RESULTS

The trans-ancestry PRS was significantly associated with T2D status across the ancestral groups examined. The top 2% of the PRS distribution can identify individuals with an approximately 2.5-4.5-fold of increase in T2D risk, which corresponds to the increased risk of T2D for first-degree relatives. The post hoc ancestry adjustment method eliminated major distributional differences in the PRS across ancestries without compromising its predictive performance.

CONCLUSIONS

By integrating T2D GWAS from multiple populations, we developed and validated a trans-ancestry PRS, and demonstrated its potential as a meaningful index of risk among diverse patients in clinical settings. Our efforts represent the first step towards the implementation of the T2D PRS into routine healthcare.

Arrhythmia Variant Associations and Reclassifications in the eMERGE-III Sequencing Study

BACKGROUND

Sequencing Mendelian arrhythmia genes in individuals without an indication for arrhythmia genetic testing can identify carriers of pathogenic or likely pathogenic (P/LP) variants. However, the extent to which these variants are associated with clinically meaningful phenotypes before or after return of variant results is unclear. In addition, the majority of discovered variants are currently classified as variants of uncertain significance, limiting clinical actionability.

METHODS

The eMERGE-III study (Electronic Medical Records and Genomics Phase III) is a multicenter prospective cohort that included 21 846 participants without previous indication for cardiac genetic testing. Participants were sequenced for 109 Mendelian disease genes, including 10 linked to arrhythmia syndromes. Variant carriers were assessed with electronic health record-derived phenotypes and follow-up clinical examination. Selected variants of uncertain significance (n=50) were characterized in vitro with automated electrophysiology experiments in HEK293 cells.

RESULTS

As previously reported, 3.0% of participants had P/LP variants in the 109 genes. Herein, we report 120 participants (0.6%) with P/LP arrhythmia variants. Compared with noncarriers, arrhythmia P/LP carriers had a significantly higher burden of arrhythmia phenotypes in their electronic health records. Fifty-four participants had variant results returned. Nineteen of these 54 participants had inherited arrhythmia syndrome diagnoses (primarily long-QT syndrome), and 12 of these 19 diagnoses were made only after variant results were returned (0.05%). After in vitro functional evaluation of 50 variants of uncertain significance, we reclassified 11 variants: 3 to likely benign and 8 to P/LP.

CONCLUSIONS

Genome sequencing in a large population without indication for arrhythmia genetic testing identified phenotype-positive carriers of variants in congenital arrhythmia syndrome disease genes. As the genomes of large numbers of people are sequenced, the disease risk from rare variants in arrhythmia genes can be assessed by integrating genomic screening, electronic health record phenotypes, and in vitro functional studies.

REGISTRATION

URL: https://www.

CLINICALTRIALS

gov; Unique identifier; NCT03394859.

A unified framework identifies new links between plasma lipids and diseases from electronic medical records across large-scale cohorts

Plasma lipids are known heritable risk factors for cardiovascular disease, but increasing evidence also supports shared genetics with diseases of other organ systems. We devised a comprehensive three-phase framework to identify new lipid-associated genes and study the relationships among lipids, genotypes, gene expression and hundreds of complex human diseases from the Electronic Medical Records and Genomics (347 traits) and the UK Biobank (549 traits). Aside from 67 new lipid-associated genes with strong replication, we found evidence for pleiotropic SNPs/genes between lipids and diseases across the phenome. These include discordant pleiotropy in the HLA region between lipids and multiple sclerosis and putative causal paths between triglycerides and gout, among several others. Our findings give insights into the genetic basis of the relationship between plasma lipids and diseases on a phenome-wide scale and can provide context for future prevention and treatment strategies.

Genome-wide Modeling of Polygenic Risk Score in Colorectal Cancer Risk

Accurate colorectal cancer (CRC) risk prediction models are critical for identifying individuals at low and high risk of developing CRC, as they can then be offered targeted screening and interventions to address their risks of developing disease (if they are in a high-risk group) and avoid unnecessary screening and interventions (if they are in a low-risk group). As it is likely that thousands of genetic variants contribute to CRC risk, it is clinically important to investigate whether these genetic variants can be used jointly for CRC risk prediction. In this paper, we derived and compared different approaches to generating predictive polygenic risk scores (PRS) from genome-wide association studies (GWASs) including 55,105 CRC-affected case subjects and 65,079 control subjects of European ancestry. We built the PRS in three ways, using (1) 140 previously identified and validated CRC loci; (2) SNP selection based on linkage disequilibrium (LD) clumping followed by machine-learning approaches; and (3) LDpred, a Bayesian approach for genome-wide risk prediction. We tested the PRS in an independent cohort of 101,987 individuals with 1,699 CRC-affected case subjects. The discriminatory accuracy, calculated by the age- and sex-adjusted area under the receiver operating characteristics curve (AUC), was highest for the LDpred-derived PRS (AUC = 0.654) including nearly 1.2 M genetic variants (the proportion of causal genetic variants for CRC assumed to be 0.003), whereas the PRS of the 140 known variants identified from GWASs had the lowest AUC (AUC = 0.629). Based on the LDpred-derived PRS, we are able to identify 30% of individuals without a family history as having risk for CRC similar to those with a family history of CRC, whereas the PRS based on known GWAS variants identified only top 10% as having a similar relative risk. About 90% of these individuals have no family history and would have been considered average risk under current screening guidelines, but might benefit from earlier screening. The developed PRS offers a way for risk-stratified CRC screening and other targeted interventions.

Abstract PD3-03: Impact of the breast cancer polygenic risk score on preventive endocrine therapy adherence and endocrine therapy usage on quality of life - The Genetic Risk Estimate (GENRE) trial

Background: Studies demonstrate breast cancer risk reduction of 50-65% with the use of endocrine therapy (ET) and yet drug uptake and adherence in this setting is suboptimal even among high risk women. A Polygenic Risk Score (PRS) comprised of 77 BC genetic susceptibility loci (Single Nucleotide Polymorphisms (SNP)) can provide a personalized risk assessment and potentially influence ET adherence. We assessed ET adherence at 1 year comparing women whose risk estimate increased due to PRS versus women whose risk estimate decreased due to PRS. The effect of ET use on quality of life (hot flashes, night sweats, vaginal dryness, weight gain, joint pain) was evaluated.

Methods: Eligible women required either a 5 year Gail Model risk of ≥3% or 10 year IBIS ( International Breast Intervention Study or Tyrer-Cuzick model) of ≥5%. Women with a history of breast cancer (BC) or hereditary BC syndrome were excluded. High risk women were counseled at baseline using their Gail and IBIS risk scores and ET options were discussed including benefits and risks. Participants completed a self-reported questionnaire at baseline to assess their understanding of breast cancer risk and decision to take preventive ET. Blood samples were obtained and genotyped for 77 SNPs, and the odds ratio from the SNPs were used to modify the IBIS and Gail risk estimates. The BC -PRS risk estimate information was shared with study participants that reflected the IBIS and Gail risk estimates for 5 year, 10 year, & lifetime BC risk with and without the PRS. Follow up questionnaires at year 1 were administered to assess drug adherence and ET usage on quality of life.

Results: 151 women were enrolled at Mayo Clinic Rochester and CancerCare Manitoba from 2016 to 2017. The median age was 56.1 (range 36-76.4), 35.6% were premenopausal, 98.7% were Caucasian and 64.7% had>1 family member with BC. Median 5yr, 10yr, & lifetime IBIS- PRS risk estimates were 3.8% (2.0-11.5), 7.9% (5.0-23.1), and 25.3% (5.5 to 92.2). At year 1 (n=112 women) 46 % of those with an increase in risk when considering the BC-PRS score and 16 % with a decrease in risk were taking ET ( p< 0.001). Types of ET taken: tamoxifen-18, raloxifene- 3, exemestane -14 and missing-1. Women taking ET reported weight gain ( 19.4% vs 6.7%, p=0.04) and more joint pain ( 27.8 % vs 12%, p=0.04) when compared to women not taking ET.

Conclusion: In high risk women, BC-PRS risk estimates in addition to standard BC risk calculators had a significant impact on preventive ET adherence. ET use was associated with weight gain and joint pain.

Citation Format: Sandhya Pruthi, Julian O Kim, Daniel Schaid, Andrew Cooke, Christina Kim, Benjamin Goldenberg, Jason Sinnwell, Debjani Grenier, Fergus Couch, Celine Vachon. Impact of the breast cancer polygenic risk score on preventive endocrine therapy adherence and endocrine therapy usage on quality of life - The Genetic Risk Estimate (GENRE) trial [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr PD3-03.

Atlas-CNV: a validated approach to call single-exon CNVs in the eMERGESeq gene panel

PURPOSE

To provide a validated method to confidently identify exon-containing copy-number variants (CNVs), with a low false discovery rate (FDR), in targeted sequencing data from a clinical laboratory with particular focus on single-exon CNVs.

METHODS

DNA sequence coverage data are normalized within each sample and subsequently exonic CNVs are identified in a batch of samples, when the target log₂ ratio of the sample to the batch median exceeds defined thresholds. The quality of exonic CNV calls is assessed by C-scores (Z-like scores) using thresholds derived from gold standard samples and simulation studies. We integrate an ExonQC threshold to lower FDR and compare performance with alternate software (VisCap).

RESULTS

Thirteen CNVs were used as a truth set to validate Atlas-CNV and compared with VisCap. We demonstrated FDR reduction in validation, simulation, and 10,926 eMERGESeq samples without sensitivity loss. Sixty-four multiexon and 29 single-exon CNVs with high C-scores were assessed by Multiplex Ligation-dependent Probe Amplification (MLPA).

CONCLUSION

Atlas-CNV is validated as a method to identify exonic CNVs in targeted sequencing data generated in the clinical laboratory. The ExonQC and C-score assignment can reduce FDR (identification of targets with high variance) and improve calling accuracy of single-exon CNVs respectively. We propose guidelines and criteria to identify high confidence single-exon CNVs.

Detecting potential pleiotropy across cardiovascular and neurological diseases using univariate, bivariate, and multivariate methods on 43,870 individuals from the eMERGE network

The link between cardiovascular diseases and neurological disorders has been widely observed in the aging population. Disease prevention and treatment rely on understanding the potential genetic nexus of multiple diseases in these categories. In this study, we were interested in detecting pleiotropy, or the phenomenon in which a genetic variant influences more than one phenotype. Marker-phenotype association approaches can be grouped into univariate, bivariate, and multivariate categories based on the number of phenotypes considered at one time. Here we applied one statistical method per category followed by an eQTL colocalization analysis to identify potential pleiotropic variants that contribute to the link between cardiovascular and neurological diseases. We performed our analyses on ~530,000 common SNPs coupled with 65 electronic health record (EHR)-based phenotypes in 43,870 unrelated European adults from the Electronic Medical Records and Genomics (eMERGE) network. There were 31 variants identified by all three methods that showed significant associations across late onset cardiac- and neurologic- diseases. We further investigated functional implications of gene expression on the detected "lead SNPs" via colocalization analysis, providing a deeper understanding of the discovered associations. In summary, we present the framework and landscape for detecting potential pleiotropy using univariate, bivariate, multivariate, and colocalization methods. Further exploration of these potentially pleiotropic genetic variants will work toward understanding disease causing mechanisms across cardiovascular and neurological diseases and may assist in considering disease prevention as well as drug repositioning in future research.

Abstract B40: Rare variant discovery in known cancer genes from whole-exome sequencingof African American hereditary prostate cancer families

African American Hereditary Prostate Cancer Study (AAHPC) was developed as a national collaboration to explore the role of genetics in the causation of hereditary prostate cancer (HPC) in African American (AA) men. AAHPC is in partnership with the International Consortium for Prostate Cancer Genetics (ICPCG), which conducts collaborative studies of HPC genetics in multiplex families. As part of an ICPCG sequencing study of 539 affected individuals from 366 HPC pedigrees, we performed whole exome sequencing in 21 ICPCG AA families, of which there were 14 AAHPC affected men from 11 pedigrees. The combined ICPCG AA cohort consisted of N=26 affected members. Post-variant calling quality control (QC) was implemented using Golden Helix SVS 8 software with filters set for removal of variants with Read Depth < 10, Quality Score < 20, Quality Score: Read Depth Ratio < 0.5, Call Rate < 0.75. Variants were additionally filtered by MAF based on the NHLBI ESP650051-V2 exomes variant frequencies for the AA population using a MAF threshold of 5%. Following QC, 176/299 SNVs and 20/39 INDELs remained for further analysis. In these analyses, we focused on 13 known cancer genes (MSR1, AR, BRCA1, BRCA2, BTNL2, EPHB2, CDH1, RNASEL, ELAC2, HOXB13, CHEK2, TP53 and NBN). Three sequenced families had > 1 affected members sequenced (2 or 3 per family) and the remaining 18 families had one member sequenced. Under the dominant model, our preliminary results show that no rare variants in the 13 candidate genes were found in 3/3 affecteds in two families. Rare SNVs in seven candidate genes (AR, CDH1, ELAC2, HOXB13, RNASEL, BRCA2 and EPHB2) were found in 2/3 affecteds for two families and 2/2 affecteds in one family. Several of the remaining 18 affected men (1 sequenced per pedigree) shared the same rare SNV in these candidate genes. For INDELs, rare variants in three candidate genes were found in pedigrees with ≥ 2 affecteds. Several of the remaining 18 affected men (one sequenced per pedigree) shared the same rare INDEL. Additional QC is underway to validate these variants and bioinformatic analyses are being used to predict effects of the variants in an effort to unravel the complex genetic heterogeneity of HPC in AA.

Citation Format: Cheryl D. Cropp, Shannon K. McDonnell, Sumit Middha, Melissa DeRycke, Danielle M. Karyadi, Daniel Schaid, Stephen N. Thibodeau, William B. Isaacs, Elaine A. Ostrander, Janet Stanford, Kathleen A. Cooney, Joan E. Bailey-Wilson, John D. Carpten. Rare variant discovery in known cancer genes from whole-exome sequencingof African American hereditary prostate cancer families. [abstract]. In: Proceedings of the Eighth AACR Conference on The Science of Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; Nov 13-16, 2015; Atlanta, GA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2016;25(3 Suppl):Abstract nr B40.

Fine-mapping the HOXB region detects common variants tagging a rare coding allele: evidence for synthetic association in prostate cancer

The HOXB13 gene has been implicated in prostate cancer (PrCa) susceptibility. We performed a high resolution fine-mapping analysis to comprehensively evaluate the association between common genetic variation across the HOXB genetic locus at 17q21 and PrCa risk. This involved genotyping 700 SNPs using a custom Illumina iSelect array (iCOGS) followed by imputation of 3195 SNPs in 20,440 PrCa cases and 21,469 controls in The PRACTICAL consortium. We identified a cluster of highly correlated common variants situated within or closely upstream of HOXB13 that were significantly associated with PrCa risk, described by rs117576373 (OR 1.30, P = 2.62×10(-14)). Additional genotyping, conditional regression and haplotype analyses indicated that the newly identified common variants tag a rare, partially correlated coding variant in the HOXB13 gene (G84E, rs138213197), which has been identified recently as a moderate penetrance PrCa susceptibility allele. The potential for GWAS associations detected through common SNPs to be driven by rare causal variants with higher relative risks has long been proposed; however, to our knowledge this is the first experimental evidence for this phenomenon of synthetic association contributing to cancer susceptibility.

HOXB13 is a susceptibility gene for prostate cancer: results from the International Consortium for Prostate Cancer Genetics (ICPCG)

Prostate cancer has a strong familial component but uncovering the molecular basis for inherited susceptibility for this disease has been challenging. Recently, a rare, recurrent mutation (G84E) in HOXB13 was reported to be associated with prostate cancer risk. Confirmation and characterization of this finding is necessary to potentially translate this information to the clinic. To examine this finding in a large international sample of prostate cancer families, we genotyped this mutation and 14 other SNPs in or flanking HOXB13 in 2,443 prostate cancer families recruited by the International Consortium for Prostate Cancer Genetics (ICPCG). At least one mutation carrier was found in 112 prostate cancer families (4.6 %), all of European descent. Within carrier families, the G84E mutation was more common in men with a diagnosis of prostate cancer (194 of 382, 51 %) than those without (42 of 137, 30 %), P = 9.9 × 10(-8) [odds ratio 4.42 (95 % confidence interval 2.56-7.64)]. A family-based association test found G84E to be significantly over-transmitted from parents to affected offspring (P = 6.5 × 10(-6)). Analysis of markers flanking the G84E mutation indicates that it resides in the same haplotype in 95 % of carriers, consistent with a founder effect. Clinical characteristics of cancers in mutation carriers included features of high-risk disease. These findings demonstrate that the HOXB13 G84E mutation is present in ~5 % of prostate cancer families, predominantly of European descent, and confirm its association with prostate cancer risk. While future studies are needed to more fully define the clinical utility of this observation, this allele and others like it could form the basis for early, targeted screening of men at elevated risk for this common, clinically heterogeneous cancer.

Abstract 4727: Breast cancer prevention and selective estrogen response modulators (SERMs): Pharmacogenomics and differential estrogen and SERM regulation of BRCA1 and BRCA2 expression

Introduction: The NSABP P-1 and P-2 breast cancer prevention trials have demonstrated that the SERMs tamoxifen and raloxifene can reduce the risk of breast cancer in high risk women.

Methods: We performed a nested matched case-control GWAS utilizing DNA from participants enrolled in P-1 and P-2. 592 participants developed invasive breast cancer or ductal carcinoma in situ, and were matched with 1171 controls who did not. Genome-wide genotyping with the Illumina 610-Quad chip and CYP2D6 genotyping for pertinent alleles were performed. Functional genomic studies involved the use of cultured cells, with siRNA knockdown, followed by qRT-PCR and Western blot analyses.

Results: Eleven SNPs with p-values < 3E-05 were identified. Twenty one additional SNPs with p < 4E-05 were identified by imputation around the genotyped SNPs on chromosomes (Chrs) 4, 8, 9, 13, 16 and 22, which were then genotyped with the Invader platform. Initial functional genomic studies focused on 6 SNPs on Chr16 (p-values 1.81-9.55E-06), all of which were in ZNF423, a gene encoding a putative zinc-finger protein. The Chr16 minor variant SNP sequences had less risk than the common variant (OR=0.7). Our functional studies showed that incubation of Hs578T cells stably transfected with estrogen receptor (ER) α with 0.1 nM estradiol (E2) induced expression of ZNF423, BRCA1/2 for WT but not variant SNP sequences. There was also a striking difference between E2-induced ZNF423, BRCA1/2 expression for variant and WT SNPs in lymphoblastoid cell lines (LCLs) stably transfected with ERα, with only WT showing induction. ChIP assay showed that ZNF423 could bind to the 5’-flanking region of BRCA1, and reporter gene assays showed that ZNF423 could regulate BRCA1/2 transcription. In LCLs stably transfected with ERα and genotyped for ZNF423 SNPs, estrogen-dependent expression of all 3 genes occurred only in the presence of the WT, but not variant ZNF423 SNP sequences. Blockade of ERα with 4-OH tamoxifen prevented induction of expression in cells with WT sequences, but resulted in gene dose-dependent increases in BRCA1/2 expression in cells with variant SNP sequences-suggested that the minor variant has a “protective” effect during clinical SERM therapy and verified by ChIP assays.

Conclusions: The SNPs in ZNF423 identified during GWAS were associated with differential E2-dependent induction of ZNF423, BRCA1 and BRCA2-implying that ZNF423 is “upstream” for the estrogen-dependent induction of BRCA1/2 expression. However, LCLs with variant SNP sequences showed greatly enhanced expression of ZNF423, BRCA1/2 during 4-OH tamoxifen exposure, while cells with WT sequences did not. These results may reveal a novel mechanism for SERM-dependent regulation of BRCA1/2 expression and may have implications for patient selection for SERM therapy.

Funded by U10-CA37377, U10-CA69974 and U14-GM61388.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 4727. doi:10.1158/1538-7445.AM2011-4727

Glycine and a glycine dehydrogenase (GLDC) SNP as citalopram/escitalopram response biomarkers in depression: pharmacometabolomics-informed pharmacogenomics

Major depressive disorder (MDD) is a common psychiatric disease. Selective serotonin reuptake inhibitors (SSRIs) are an important class of drugs used in the treatment of MDD. However, many patients do not respond adequately to SSRI therapy. We used a pharmacometabolomics-informed pharmacogenomic research strategy to identify citalopram/escitalopram treatment outcome biomarkers. Metabolomic assay of plasma samples from 20 escitalopram remitters and 20 nonremitters showed that glycine was negatively associated with treatment outcome (P = 0.0054). This observation was pursued by genotyping tag single-nucleotide polymorphisms (SNPs) for genes encoding glycine synthesis and degradation enzymes, using 529 DNA samples from SSRI-treated MDD patients. The rs10975641 SNP in the glycine dehydrogenase (GLDC) gene was associated with treatment outcome phenotypes. Genotyping for rs10975641 was carried out in 1,245 MDD patients in the Sequenced Treatment Alternatives to Relieve Depression (STAR*D) study, and its presence was significant (P = 0.02) in DNA taken from these patients. These results highlight a possible role for glycine in SSRI response and illustrate the use of pharmacometabolomics to "inform" pharmacogenomics.

A Genome-Wide Association Study in Patients Experiencing Musculoskeletal Adverse Events on Aromatase Inhibitors as Adjuvant Therapy in Early Breast Cancer Entered on NCIC CTG Trial MA.27. A Pharmacogenetics Research Network-RIKEN Collaboration

Background: The third-generation aromatase inhibitors (AIs) are established adjuvant therapy for postmenopausal women with early stage breast cancer. Musculoskeletal adverse events (MS AE) have become increasingly recognized as an issue for many women receiving AIs resulting in decreased quality of life and, in some instances, discontinuation of the AI. MA.27 is an ongoing randomized clinical trial comparing anastrozole with exemestane that accrued 7576 patients (pts), including 6827 from North America, about 79% of whom provided DNA and consent for genetic testing. A review of withdrawals related to adverse events revealed that the most common cause was MS AE. The hypothesis for our study was that a genome-wide association (GWA) case control study would identify single nucleotide polymorphisms (SNPs) associated with MS AE.Methods: Cases were defined as pts with grade 3 or 4 MS AE (NCI CTCAE v3.0) or who went off treatment for any grade MS AE within the first 2 years and were matched to 2 controls based on treatment arm, length of follow-up, prior chemotherapy, age and prior celecoxib (at least 3 months prior to MS AE). The study was restricted to the 94% of pts self-identified as Caucasian. Genotyping was performed with the Illumina Human610-Quad. Genotype quality control prior to analysis required call rates ≥98%. Eigenstrat analyses were performed to control for population stratification.Results: The GWA study included 293 cases and 586 controls with 6 pts (0.7%) (4 cases, 2 controls) excluded for call rates <98%. 551,358 SNPs were used in the analyses after removing 29,439 for a minor allele frequency <1% and 82 with Hardy Weinberg Equilibrium p-values <10-6. Eigenvectors did not impact the results when used as covariates. Four SNPs with the lowest p-values (3.3-8.2 x 10-6, Armitage test), 3 in high LD (r2>.8) on chromosome (ch) 14 (rs7158782, rs7159713, rs2369049) and 1 on the X ch (rs6637820) were identified. The gene closest (4000-7000 bp) to the 3 SNPs on ch 14 was T-Cell Leukemia 1A (TCL1A), encoding a protein known to augment AKT kinase activity. Expression array data from lymphoblastoid cell lines from 300 subjects of 3 ethnicities revealed that all 3 SNPs were associated with decreased TCL1A expression after adjusting for ethnicity using 2 probe sets (p=0.006 to 0.03). Electromobility shift assay (EMSA) showed that 2 of these SNPs (rs7158782, rs7159713) displayed a shift and the former SNP showed substantially less binding for variant than wild type. The gene closest to the X ch SNP was Immunoglobulin Superfamily 1 (IGSF1), but EMSA did not show a shift for this SNP.Conclusions: MS AEs represent a major impediment to optimal use of AIs in women with breast cancer and this GWA study identified SNPs on ch 14 that provide a focus for further research to identify pts at risk for, and means to ameliorate, this adverse event.(Supported in part by NIH grants U01GM61388, U01GM63173, P50CA116201 and U10CA77202)

Citation Information: Cancer Res 2009;69(24 Suppl):Abstract nr 15.

Multiple novel prostate cancer predisposition loci confirmed by an international study: the PRACTICAL Consortium

A recent genome-wide association study found that genetic variants on chromosomes 3, 6, 7, 10, 11, 19 and X were associated with prostate cancer risk. We evaluated the most significant single-nucleotide polymorphisms (SNP) in these loci using a worldwide consortium of 13 groups (PRACTICAL). Blood DNA from 7,370 prostate cancer cases and 5,742 male controls was analyzed by genotyping assays. Odds ratios (OR) associated with each genotype were estimated using unconditional logistic regression. Six of the seven SNPs showed clear evidence of association with prostate cancer (P = 0.0007-P = 10(-17)). For each of these six SNPs, the estimated per-allele OR was similar to those previously reported and ranged from 1.12 to 1.29. One SNP on 3p12 (rs2660753) showed a weaker association than previously reported [per-allele OR, 1.08 (95% confidence interval, 1.00-1.16; P = 0.06) versus 1.18 (95% confidence interval, 1.06-1.31)]. The combined risks associated with each pair of SNPs were consistent with a multiplicative risk model. Under this model, and in combination with previously reported SNPs on 8q and 17q, these loci explain 16% of the familial risk of the disease, and men in the top 10% of the risk distribution have a 2.1-fold increased risk relative to general population rates. This study provides strong confirmation of these susceptibility loci in multiple populations and shows that they make an important contribution to prostate cancer risk prediction.

Gemcitabine and cytosine arabinoside cytotoxicity: association with lymphoblastoid cell expression

Two cytidine analogues, gemcitabine (dFdC) and 1-beta-d-arabinofuranosylcytosine (AraC), show significant therapeutic effect in a variety of cancers. However, response to these drugs varies widely. Evidence from tumor biopsy samples shows that expression levels for genes involved in the cytidine transport, metabolism, and bioactivation pathway contribute to this variation in response. In the present study, we set out to test the hypothesis that variation in gene expression both within and outside of this "pathway" might influence sensitivity to gemcitabine and AraC. Specifically, Affymetrix U133 Plus 2.0 GeneChip and cytotoxicity assays were performed to obtain basal mRNA expression and IC(50) values for both drugs in 197 ethnically defined Human Variation Panel lymphoblastoid cell lines. Genes with a high degree of association with IC(50) values were involved mainly in cell death, cancer, cell cycle, and nucleic acid metabolism pathways. We validated selected significant genes by performing real-time quantitative reverse transcription-PCR and selected two representative candidates, NT5C3 (within the pathway) and FKBP5 (outside of the pathway), for functional validation. Those studies showed that down-regulation of NT5C3 and FKBP5 altered tumor cell sensitivity to both drugs. Our results suggest that cell-based model system studies, when combined with complementary functional characterization, may help to identify biomarkers for response to chemotherapy with these cytidine analogues.

Breast cancer risk reduction and membrane-bound catechol O-methyltransferase genetic polymorphisms

Catechol O-methyltransferase (COMT)-catalyzed methylation of catecholestrogens has been proposed to play a protective role in estrogen-induced genotoxic carcinogenesis. We have taken a comprehensive approach to test the hypothesis that genetic variation in COMT might influence breast cancer risk. Fifteen COMT single nucleotide polymorphisms (SNPs) selected on the basis of in-depth resequencing of the COMT gene were genotyped in 1,482 DNA samples from a Mayo Clinic breast cancer case control study. Two common SNPs in the distal promoter for membrane-bound (MB) COMT, rs2020917 and rs737865, were associated with breast cancer risk reduction in premenopausal women in the Mayo Clinic study, with allele-specific odds ratios (OR) of 0.70 [95% confidence interval (CI), 0.52-0.95] and 0.68 (95% CI, 0.51-0.92), respectively. These two SNPs were then subjected to functional genomic analysis and were genotyped in an additional 3,683 DNA samples from two independent case control studies (GENICA and GESBC). Functional genomic experiments showed that these SNPs could up-regulate transcription and that they altered DNA-protein binding patterns. Furthermore, substrate kinetic and exon array analyses suggested a role for MB-COMT in catecholestrogen inactivation. The GENICA results were similar to the Mayo case control observations, with ORs of 0.85 (95% CI, 0.72-1.00) and 0.85 (95% CI, 0.72-1.01) for the two SNPs. No significant effect was observed in the GESBC study. These studies showed that two SNPs in the COMT distal promoter were associated with breast cancer risk reduction in two of three case control studies, compatible with the results of functional genomic experiments, suggesting a role for MB-COMT in breast cancer risk.

Human aromatase: gene resequencing and functional genomics

Aromatase [cytochrome P450 19 (CYP19)] is a critical enzyme for estrogen biosynthesis, and aromatase inhibitors are of increasing importance in the treatment of breast cancer. We set out to identify and characterize genetic polymorphisms in the aromatase gene, CYP19, as a step toward pharmacogenomic studies of aromatase inhibitors. Specifically, we "resequenced" all coding exons, all upstream untranslated exons plus their presumed core promoter regions, all exon-intron splice junctions, and a portion of the 3'-untranslated region of CYP19 using 240 DNA samples from four ethnic groups. Eighty-eight polymorphisms were identified, resulting in 44 haplotypes. Functional genomic studies were done with the four nonsynonymous coding single nucleotide polymorphisms (cSNP) that we observed, two of which were novel. Those cSNPs altered the following amino acids: Trp39Arg, Thr201Met, Arg264Cys, and Met364Thr. The Cys264, Thr364, and double variant Arg39Cys264 allozymes showed significant decreases in levels of activity and immunoreactive protein when compared with the wild-type (WT) enzyme after transient expression in COS-1 cells. A slight decrease in protein level was also observed for the Arg39 allozyme, whereas Met201 displayed no significant changes in either activity or protein level when compared with the WT enzyme. There was also a 4-fold increase in apparent K(m) value for Thr364 with androstenedione as substrate. Of the recombinant allozymes, only the double mutant (Arg39Cys264) displayed a significant change from the WT enzyme in inhibitor constant for the aromatase inhibitors exemestane and letrozole. These observations indicate that genetic variation in CYP19 might contribute to variation in the pathophysiology of estrogen-dependent disease.

Behavioral risk exposure and host genetics of susceptibility to HIV-1 infection

BACKGROUND

Some individuals are readily infected with low human immunodeficiency virus type 1 (HIV-1) exposure, whereas others appear less susceptible, suggesting that host genetics plays a role in the viral entry pathway. The matched case-control study design with measured risk exposures provides an avenue for discovering genes involved in susceptibility to infection.

METHODS

We conducted a nested case-control study of African Americans (266 HIV-1 seroconverter cases and 532 seronegative controls from the AIDS Link to Intravenous Experience cohort), to examine the association between 50 single-nucleotide polymorphisms (SNPs) in 9 candidate genes (CCR5, CCR2, RANTES, MIP1A, MCP2, IL10, IFNG, MCSF, and IL2) and susceptibility to HIV-1 infection. To account for differential exposure propensities, risk behavior self-reported during semiannual visits was used to estimate a standardized cumulative risk exposure (SCRE). Individual SNPs were evaluated using conditional logistic-regression models, and the inferred haplotypes were assessed in the haplotype trend regression analyses after adjusting for age and SCRE.

RESULTS

Four SNPs (CCR2-V64I, CCR5-2459, MIP1A+954, and IL2+3896) and specific haplotypes in the IL2 and CCR2/CCR5 regions were significantly associated with HIV-1 infection susceptibility in different genetic models.

CONCLUSIONS

Our results suggest that genetic variants in associated host genes may play an important role in susceptibility to HIV-1 infection.

TBX2 is preferentially amplified in BRCA1- and BRCA2-related breast tumors

The chromosome 17q23 region is frequently amplified in breast tumors. Gain of the region is present in 50% of BRCA1-associated breast tumors and 87% of BRCA2-associated breast tumors. The amplification frequency of the RPS6KB1 and TBX2 oncogenes from this amplicon was compared in 27 BRCA1 and BRCA2 mutant breast tumors, 15 breast tumors from high-risk patients with no BRCA1 or BRCA2 mutations, and 62 matched sporadic breast tumor controls. TBX2 was determined to be preferentially amplified and overexpressed in BRCA1 and BRCA2 mutant tumors, whereas RPS6KB1 was not, suggesting a role for TBX2 amplification in the development of BRCA1- and BRCA2-associated breast tumors.

Genetic variation in the transforming growth factor beta1 gene in multiple sclerosis

Transforming growth factor beta1 (TGFbeta1) is a Th2 cytokine encoded on chromosome 19q13, a region possibly linked to multiple sclerosis (MS). TGFbeta1 exerts favorable effects on experimental allergic encephalomyelitis. We performed a comprehensive search for genetic variants in this gene in 122 population-based sporadic cases of MS. We detected six variants, including three missense variants. We tested for association of the variants with susceptibility and course of MS and for linkage and transmission disequilibrium in a family series consisting of 395 samples in 59 pedigrees. Genetic variation in TGFB1 does not appear to contribute in a major way to susceptibility to MS.

A prospective trial of midwest breast cancer patients: a p53 gene mutation is the most important predictor of adverse outcome

Several retrospective studies have suggested p53 gene mutation as an adverse prognostic indicator in breast cancer patients, based on a selective growth advantage of p53 mutant cancer cells and their presumed resistance to current adjuvant therapy regimens. A cohort of 90 Caucasian midwestern breast cancer patients was analyzed prospectively (60 months of follow-up) with a rigorous mutation detection methodology. The presence of a p53 gene mutation was the single most adverse prognostic indicator for recurrence (p = 0.0032) and death (p = 0.0001), and was associated with poor response to both adjuvant (p = 0.0001) and palliative (p = 0.006) therapy. Analysis of the p53 gene with appropriate mutation detection methodology markedly improves the prediction of early recurrence, treatment failure, and death in breast cancer patients.

BRCA1 and BRCA2 have a limited role in familial prostate cancer

Epidemiological studies have suggested that the breast cancer susceptibility genes, BRCA1 and BRCA2, may be involved in the development of prostate cancer. Several studies have screened prostate cancer populations for the presence of BRCA1 and BRCA2 mutations, with few mutations identified. In this study, 22 high-risk prostate cancer families (at least three cases of prostate cancer) were screened by conformation-sensitive gel electrophoresis (CSGE) for mutations in BRCA1 and BRCA2. To maximize the chance of finding mutations in these two genes, families were also selected for the presence of at least two cases of breast and/or ovarian cancer. We identified one previously reported BRCA2 missense mutation and two previously unreported BRCA2 intron polymorphisms. No BRCA1 or BRCA2 truncating mutations were detected. Thus, BRCA1 and BRCA2 appear to have a limited role in familial prostate cancer, and families with both prostate and breast cancer may result from mutations in other predisposition genes.

A prospective trial of midwest breast cancer patients: a p53 gene mutation is the most important predictor of adverse outcome

Several retrospective studies have suggested p53 gene mutation as an adverse prognostic indicator in breast cancer patients, based on a selective growth advantage of p53 mutant cancer cells and their presumed resistance to current adjuvant therapy regimens. A cohort of 90 Caucasian midwestern breast cancer patients was analyzed prospectively (60 months of follow-up) with a rigorous mutation detection methodology. The presence of a p53 gene mutation was the single most adverse prognostic indicator for recurrence (p = 0.0032) and death (p = 0.0001), and was associated with poor response to both adjuvant (p = 0.0001) and palliative (p = 0.006) therapy. Analysis of the p53 gene with appropriate mutation detection methodology markedly improves the prediction of early recurrence, treatment failure, and death in breast cancer patients.

Applicability of the case-parent design in the etiological research of Type 1 diabetes in Chile and other genetically mixed populations

In case-control studies, spurious associations between Human Leukocyte Antigen (HLA) alleles and Type 1 diabetes could arise as a result of population stratification, if there are ethnic differences between cases and non-related controls. The Chilean population has several unique features which make it ideal for the study of the effect of stratification by ethnicity on genetic epidemiological research. The incidence rates of Type 1 diabetes in Chilean Aboriginal populations are very low compared to Caucasian populations, while the frequency of the alleles in HLA loci also vary across ethnic groups. In order to avoid the confounding effect of ethnicity, one possible remedy would be the use of cases and their parents in place of non-related controls. The case-parent design offers an adequate framework for the study of the association between HLA polymorphisms and Type 1 diabetes in the Chilean population and can also be applicable to other genetically mixed populations especially in the Americas.

The association between HLA class I alleles and measles vaccine-induced antibody response: evidence of a significant association

While the Moraten strain measles vaccine is an excellent, safe, and immunogenic vaccine, vaccine failure occurs, presumably when an individual develops an inadequate immune response. In this study, we examined the association of HLA class I genes and measles vaccine-induced antibody levels. We found that the allele distribution of HLA-B alleles differed between non-responders and hyper-responders (p = 0.002). Several class I alleles were associated with non-response (HLA-B13, -B44, and -C5); whereas several other alleles were associated with hyper-response (HLA-B7 and -B51). In addition, non-responders were more likely to be HLA-B homozygous than normal responders (odds ratio 2.1), and more likely to be homozygous than hyper-responders (odds ratio 3.7, p = 0.031 Mantel-Haenzel for trend). Finally, we found evidence of an allele dose-response phenomenon for HLA-B7. We conclude that there are important associations between class I HLA genes and measles antibody levels following immunization.

Evidence for a prostate cancer susceptibility locus on the X chromosome

Over 200,000 new prostate cancer cases are diagnosed in the United States each year, accounting for more than 35% of all cancer cases affecting men, and resulting in 40,000 deaths annually. Attempts to characterize genes predisposing to prostate cancer have been hampered by a high phenocopy rate, the late age of onset of the disease and, in the absence of distinguishing clinical features, the inability to stratify patients into subgroups relative to suspected genetic locus heterogeneity. We previously performed a genome-wide search for hereditary prostate cancer (HPC) genes, finding evidence of a prostate cancer susceptibility locus on chromosome 1 (termed HPC1; ref. 2). Here we present evidence for the location of a second prostate cancer susceptibility gene, which by heterogeneity estimates accounts for approximately 16% of HPC cases. This HPC locus resides on the X chromosome (Xq27-28), a finding consistent with results of previous population-based studies suggesting an X-linked mode of HPC inheritance. Linkage to Xq27-28 was observed in a combined study population of 360 prostate cancer families collected at four independent sites in North America, Finland and Sweden. A maximum two-point lod score of 4.60 was observed at DXS1113, theta=0.26, in the combined data set. Parametric multipoint and non-parametric analyses provided results consistent with the two-point analysis. Significant evidence for genetic locus heterogeneity was observed, with similar estimates of the proportion of linked families in each separate family collection. Genetic mapping of the locus represents an important initial step in the identification of an X-linked gene implicated in the aetiology of HPC.

Major histocompatibility complex class II alleles and the course and outcome of MS: a population-based study

BACKGROUND

The major histocompatibility complex (MHC) has been consistently associated with susceptibility to MS and the course of several other human autoimmune diseases. A putative association between the course and severity of MS and the MHC remains controversial.

METHODS

DR and DQ genotyping by either restriction fragment length polymorphism or sequence-specific PCR-based typing in 119 patients representing 73.4% of the population with MS evaluated in a cross-sectional disability survey and 100 healthy controls from Olmsted County, Minnesota.

RESULTS

We found a positive association between MS susceptibility and the DR15-DQ6 and DR13-DQ7 haplotypes, and we found a negative association with the DR1-DQ5 haplotype. We found a trend to a positive association of primary progressive MS with DR4-DQ8 and DR1-DQ5 and an association of "bout onset" MS with DR17-DQ2. We did not find an association with disease severity, as defined by EDSS/duration.

CONCLUSION

Lack of consistency between different studies may be due to regional variation in MS and limitations of power but likely indicate a minor effect of MHC class II genes on the course and severity of MS.

Expansion of unusual CD4+ T cells in severe rheumatoid arthritis

OBJECTIVE

The repertoire of T cells in patients with rheumatoid arthritis (RA) is characterized by clonal expansion of selected CD4+ T cells, which are autoreactive and lack the expression of the functionally important CD28 molecule. The goal of this study was to determine the contribution of these unusual lymphocytes to the disease process.

METHODS

RA patients (n = 108) and normal controls (n = 53) were examined for the expression of CD4+ CD28- T cells by 2-color fluorescence-activated cell sorter analysis. Clinical data were ascertained by retrospective chart review.

RESULTS

The frequencies of CD4+ CD28- T cells displayed a bimodal distribution, defining carriers and noncarriers in normal subjects and RA patients. In longitudinal studies, the noncarrier and carrier phenotypes were stable over time. Carriers of CD4+ CD28- T cells accumulated in the RA population (64% versus 45%; P = 0.02). The expansion of CD4+ CD28- T cells correlated with extraarticular involvement, but not with disease duration, antirheumatic treatment, or severity of joint destruction. The patient subsets with nodular disease (P = 0.02) and rheumatoid organ disease (P = 0.04) had the highest proportion of CD4+ CD28- T cell carriers. The size of the CD4+ CD28- compartment correlated with extraarticular progression of RA (P = 0.001 in nodular RA, P = 0.003 in rheumatoid organ disease).

CONCLUSION

The bimodality of distribution of CD4+ CD28- T cell frequencies is compatible with genetic control of the generation of these unusual T cells. In RA patients, CD4+ CD28- T cells are not an epiphenomenon of the disease process, but predispose patients to developing inflammatory lesions in extraarticular tissues.

Mutation detection by highly sensitive methods indicates that p53 gene mutations in breast cancer can have important prognostic value

Human cancer cells with a mutated p53 tumor-suppressor gene have a selective growth advantage and may exhibit resistance to ionizing radiation and certain chemotherapeutic agents. To examine the prognostic value of mutations in the p53 gene, a cohort of 90 Midwestern Caucasian breast cancer patients were analyzed with methodology that detects virtually 100% of all mutations. The presence of a p53 gene mutation was by far the single most predictive indicator for recurrence and death (relative risks of 4.7 and 23.2, respectively). Direct detection of p53 mutations had substantially greater prognostic value than immunohistochemical detection of p53 overexpression. Analysis of p53 gene mutations may permit identification of a subset of breast cancer patients who, despite lack of conventional indicators of poor prognosis, are at high risk of early recurrence and death.

Frequency of familial dilated cardiomyopathy

Based on evaluation of 59 probands and their families, we previously demonstrated that over 20% of patients with idiopathic dilated cardiomyopathy (DCM) have familial disease. We acquired 36 additional probands and found familial disease in 24.2% of the 95 probands in this expanded consecutively ascertained cohort. The family history, as reported by the patient and relatives, is often an unreliable indicator of familial disease because patients may be unaware of the significance of a family history of sudden death, arrhythmia or stroke. We demonstrate that careful family history with review of medical records identifies more familial cases than merely asking the patient if there is a family history of DCM. However, even such a careful family history does not identify all familial cases. Some familial cases are identified only by echocardiographic investigation of asymptomatic relatives. We found no clinical attributes of probands, other than family history, which predicted familial disease.

p53 gene mutations in breast cancers in midwestern US women: null as well as missense-type mutations are associated with poor prognosis

We determined the pattern of mutations in exons 2-11 and adjacent intronic regions in breast cancers from Midwestern US white women. Twenty-one mutations were detected in 53 tumors (39.6%). Comparisons of the pattern of mutations within exons 5-9 showed that the frequency of missense mutations (44%) was lower in breast cancers of US Midwestern women than in most tumor types including breast cancers in other populations. Compared to breast cancers reported in a Scottish population, US women had a high frequency of G:C-->T:A transversions (P = 0.046). These findings suggest that environmental or endogenous factors contribute to p53 mutagenesis in mammary tissue to different extents among different populations. With a median follow-up of 19 months, the presence of a mutation was associated with shorter time to disease recurrence (P = 0.05) and shorter survival (P = 0.003). Putative dominant negative missense-type mutations (missense and in-frame microdeletions; P = 0.001) and null mutations (hemizygous nonsense and frameshift mutations; P = 0.007) were equally ominous. Thus, tumors with missense p53 mutations resulting in over-expression of a dysfunctional but otherwise intact protein have a clinical outcome similar to tumors with null mutations resulting in a truncated or garbled protein.

Novel pattern of p53 gene mutations in an American black cohort with high mortality from breast cancer

The pattern of acquired mutations in the p53 gene can be used to study differences in factors contributing to carcinogenesis. We investigated mutations in exons 5-9 and adjacent intronic regions in 47 breast cancers of black women from Michigan, a population with the highest breast-cancer mortality in the US. The 16 mutations detected differed from those of other populations. In particular, the black women had an excess of A:T-->G:C transitions compared with rural white US midwest women. While the causes of the different pattern of acquired mutation remain to be determined, this molecular epidemiological approach detects the consequences of mutagenic processes in specific populations. Mutation patterns will constrain hypotheses to mechanisms consistent with the observed biochemical alterations.

Microsatellite instability in cancer of the proximal colon

Colorectal tumor DNA was examined for somatic instability at (CA)n repeats on human chromosomes 5q, 15q, 17p, and 18q. Differences between tumor and normal DNA were detected in 25 of the 90 (28 percent) tumors examined. This instability appeared as either a substantial change in repeat length (often heterogeneous in nature) or a minor change (typically two base pairs). Microsatellite instability was significantly correlated with the tumor's location in the proximal colon (P = 0.003), with increased patient survival (P = 0.02), and, inversely, with loss of heterozygosity for chromosomes 5q, 17p, and 18q. These data suggest that some colorectal cancers may arise through a mechanism that does not necessarily involve loss of heterozygosity.

Genetic and physical mapping of a novel region close to the fragile X site on the human X chromosome

We report the isolation and characterization of a novel DNA marker (1A1) in Xqter in the region of the fragile X. Genetic studies in families segregating for the fragile X syndrome suggest that 1A1 lies between the disease mutation and the distal locus, DXS52. Studies in normal and fragile X families show that 1A1 is tightly linked to DXS52 (Zmax = 17.20; theta max = 0.03) and F8 (Zmax = 7.01; theta max = 0.08). Multipoint mapping of families supports the order Xcen-DXS105-FRAXA-1A1-DXS52-(F8, DXS115)-Xqter. Pulsed-field gel electrophoresis (PFGE) studies demonstrate that 1A1 defines a new region of at least 2 Mb of DNA not physically linked to DXS52 or F8, thus extending the physical map of Xq27-qter to over 4 Mb. Complex partial digestion PFGE patterns, probably due to differing degrees of methylation, are observed with 1A1 in unrelated normal and fragile-X-positive individuals, whereas other distal markers give uniform digestion profiles. Physical data suggest that 1A1 lies in a region less CpG rich than other distal markers in Xq27-qter.

An intronic region within the human factor VIII gene is duplicated within Xq28 and is homologous to the polymorphic locus DXS115 (767)

The genomic sequences recognized by the anonymous probe 767 (DXS115) are localized to two sites within Xq28. One site lies within intron 22 of the factor VIII gene (FBC). Physical mapping suggests that the second site lies within 1.2 megabases of the F8C gene. The RFLPs detected by 767 are located within the second site. Genetic data suggest that F8C and DXS115 are tightly linked (theta max = .04; Zmax = 8.30). Recombination events in meioses informative for DXS52 (St14), DXS115, and F8C suggest that DXS115 and F8C lie distal to DXS52.

Androgen stimulation and laryngeal development

Androgen-induced changes in laryngeal growth patterns were studied using a sheep animal model. Forty-eight lambs were divided into eight treatment groups. Lambs in seven of the groups were castrated at birth, while lambs in the eighth group served as an intact (noncastrated) control. Six groups were then treated with varying doses of testosterone and dihydrotestosterone, while the seventh served as a castrated, nontreated control. All animals were killed and gross dissections of the larynges were performed. Thirty-four linear and angular measurements were obtained from each larynx. The mean superior thyroid horn separation showed the most dramatic androgen-induced effect (p = 0.023). Laryngeal anterior-posterior diameter, superior thyroid horn height, posterior thyroid cartilage width, thyroid cartilage angle, and vocal process to arytenoid base distances all demonstrated positive dose-response relationships. Hypoandrogenic levels appeared to have an inhibitory effect upon laryngeal growth when compared to castrated controls.