Rebecca A Betensky

Rebecca Betensky
Rebecca Betensky
Scroll

Chair of the Department of Biostatistics

Professor of Biostatistics

Professional overview

Prior to NYU, Dr. Betensky was Professor of Biostatistics at the Harvard T.H. Chan School of Public Health. She was director of the Harvard Catalyst (Clinical and Translational Science Award) Biostatistics Program; director of the Data and Statistics Core for the Massachusetts Alzheimer’s Disease Research Center; and director of the Biostatistics Neurology Core at Massachusetts General Hospital. Previously, she was the Biostatistics Program Leader for the Dana-Farber/Harvard Cancer Center.

Dr. Betensky’s research focuses on methods for the analysis of censored and truncated outcomes and covariates, which frequently arise from the subsampling of cohort studies. She has a long-time interest in clinical trials, and has written on the evaluation of biomarkers and the use and interpretation of p-values. She has collaborated extensively in studies in neurologic diseases, and serves as statistical editor for Annals of Neurology.

Dr. Betensky was awarded, and directed for 15 years, an NIH T32 training program in neurostatistics and neuroepidemiology for pre- and post-doctoral students in biostatistics and epidemiology and for clinician-scientists. She previously directed Harvard’s Biostatistics programs to promote and support diversity at all levels in the field of quantitative public health. She was also a member of the BMRD Study Section for review of NIH statistical methodology grants; on committees for the Institute of Medicine; and a co-chair of the technical advisory committee for the scientific registry of transplant recipients.

Dr. Betensky an elected Fellow of the American Statistical Association and of the International Statistical Institute, and is a past recipient of the Spiegelman Award from the American Public Health Association. She currently serves as a member of the Board of Scientific Counselors for Clinical Science and Epidemiology at the National Cancer Institute.

Education

AB, Mathematics, Harvard University, Cambridge, MA
PhD, Statistics, Stanford University, Stanford, CA

Areas of research and study

Biology
Biostatistics
Neuroepidemiology
Neurology
Neurostatistics
Translational science

Publications

Publications

Loss of the mismatch repair protein MSH6 in human glioblastomas is associated with tumor progression during temozolomide treatment

Cahill, D. P., Levine, K. K., Betensky, R. A., Codd, P. J., Romany, C. A., Reavie, L. B., Batchelor, T. T., Futreal, P. A., Stratton, M. R., Curry, W. T., Lafrate, A. J., & Louis, D. N. (n.d.).

Publication year

2007

Journal title

Clinical Cancer Research

Volume

13

Issue

7

Page(s)

2038-2045
Abstract
Abstract
Purpose: Glioblastomas are treated by surgical resection followed by radiotherapy [X-ray therapy (XRT)] and the alkylating chemotherapeutic agent temozolomide. Recently, inactivating mutations in the mismatch repair gene MSH6 were identified in two glioblastomas recurrent post-temozolomide. Because mismatch repair pathway inactivation is a known mediator of alkylator resistance in vitro, these findings suggested that MSH6 inactivation was causally linked to these two recurrences. However, the extent of involvement of MSH6 in glioblastoma is unknown. We sought to determine the overall frequency and clinical relevance of MSH6 alterations in glioblastomas. Experimental Design: The MSH6 gene was sequenced in 54 glioblastomas. MSH6 and O6- methylguanine methyltransferase (MGMT) immunohistochemistry was systematically scored in a panel of 46 clinically well-characterized glioblastomas, and the corresponding patient response to treatment evaluated. Results: MSH6 mutation was not observed in any pretreatment glioblastoma (0 of 40), whereas 3 of 14 recurrent cases had somatic mutations (P = 0.015). MSH6 protein expression was detected in all pretreatment (17 of 17) cases examined but, notably, expression was lost in 7 of 17 (41%) recurrences from matched post-XRT + temozolomide cases (P = 0.016). Loss of MSH6 was not associated with O6-methylguanine methyltransferase status. Measurements of in vivo tumor growth using three-dimensional reconstructed magnetic resonance imaging showed that MSH6-negative glioblastomas had a markedly increased rate of growth while under temozolomide treatment (3.17 versus 0.04 cc/mo for MSH6-positive tumors; P = 0.020). Conclusions: Loss of MSH6 occurs in a subset of post-XRT + temozolomide glioblastoma recurrences and is associated with tumor progression during temozolomide treatment, mirroring the alkylator resistance conferred by MSH6 inactivation in vitro. MSH6 deficiency may therefore contribute to the emergence of recurrent glioblastomas during temozolomide treatment.

Multivariate logistic regression for familial aggregation in age at disease onset

Matthews, A. G., Finkelstein, D. M., & Betensky, R. A. (n.d.).

Publication year

2007

Journal title

Lifetime Data Analysis

Volume

13

Issue

2

Page(s)

191-209
Abstract
Abstract
Familial aggregation studies seek to identify diseases that cluster in families. These studies are often carried out as a first step in the search for hereditary factors affecting the risk of disease. It is necessary to account for age at disease onset to avoid potential misclassification of family members who are disease-free at the time of study participation or who die before developing disease. This is especially true for late-onset diseases, such as prostate cancer or Alzheimer's disease. We propose a discrete time model that accounts for the age at disease onset and allows the familial association to vary with age and to be modified by covariates, such as pedigree relationship. The parameters of the model have interpretations as conditional log-odds and log-odds ratios, which can be viewed as discrete time conditional cross hazard ratios. These interpretations are appealing for cancer risk assessment. Properties of this model are explored in simulation studies, and the method is applied to a large family study of cancer conducted by the National Cancer Institute-sponsored Cancer Genetics Network (CGN).

Predicting short-term disability in multiple sclerosis

Gauthier, S. A., Mandel, M., Guttmann, C. R., Glanz, B. I., Khoury, S. J., Betensky, R. A., & Weiner, H. L. (n.d.).

Publication year

2007

Journal title

Neurology

Volume

68

Issue

24

Page(s)

2059-2065
Abstract
Abstract
OBJECTIVE: To develop covariate specific short-term disability curves to demonstrate the probability of progressing by Expanded Disability Status Scale (EDSS) at semiannual visits. METHODS: Semiannual EDSS scores were prospectively collected in 218 relapsing-remitting (RR) and clinically isolated syndrome (CIS) patients as part of the Comprehensive Longitudinal Investigation of Multiple Sclerosis at the Brigham and Women's Hospital (CLIMB) study. Baseline brain parenchymal fraction (BPF) and T2 lesion volume were available on 205 patients. A partial proportional odds model determined the influence of covariates on the change in EDSS score at subsequent visits. A discrete second order Markov transitional model was fit and generated a probability matrix for each subject; the 6-month probabilities of EDSS change were graphically represented. RESULTS: The univariate analysis demonstrated the lowest baseline BPF quartile (OR 1.99; p = 0.0203) and the highest T2 lesion volume quartile (OR 2.19; p = 0.0130) were associated with progression in EDSS. Covariate specific disability curves demonstrated the effect of BPF and T2 lesion volume on short-term progression. In subjects with a 6-month EDSS of 2, the probability of a sustained progression of an EDSS of 3 within 3 years was 0.277 for a subject with low BPF and a high T2 lesion volume vs 0.055 for a subject with high BPF and a low T2 lesion volume. CONCLUSIONS: Markov transitional models allow for the comparison of covariate specific short-term disability changes among groups of patients with multiple sclerosis.

Testing goodness of fit of a uniform truncation model

Mandel, M., & Betensky, R. A. (n.d.).

Publication year

2007

Journal title

Biometrics

Volume

63

Issue

2

Page(s)

405-412
Abstract
Abstract
Several goodness-of-fit tests of a lifetime distribution have been suggested in the literature; many take into account censoring and/or truncation of event times. In some contexts, a goodness-of-fit test for the truncation distribution is of interest. In particular, better estimates of the lifetime distribution can be obtained when knowledge of the truncation law is exploited. In cross-sectional sampling, for example, there are theoretical justifications for the assumption of a uniform truncation distribution, and several studies have used it to improve the efficiency of their survival estimates. The duality of lifetime and truncation in the absence of censoring enables methods for testing goodness of fit of the lifetime distribution to be used for testing goodness of fit of the truncation distribution. However, under random censoring, this duality does not hold and different tests are required. In this article, we introduce several goodness-of-fit tests for the truncation distribution and investigate their performance in the presence of censored event times using simulation. We demonstrate the use of our tests on two data sets.

Tests of association under misclassification: Application to histological sampling in oncology

Betensky, R. A., Louis, D. N., & Cairncross, J. G. (n.d.).

Publication year

2007

Journal title

Statistics in Medicine

Volume

26

Issue

26

Page(s)

4808-4816
Abstract
Abstract
Subjects in tumour studies are often misclassified with respect to histologic features that are not routinely recorded in diagnostic reports and that display heterogeneity within tumours. Pathologic analysis of the tumours may miss the feature of interest if the pathologist was not alerted to detail the microscopic feature of interest or if it is not present in the selected specimens. In this setting, only the subjects for whom the outcome is not found are potentially misclassified. Analyses of associations between the observed, potentially misclassified, outcome and a second outcome are invalid if the probability of misclassification depends on the second outcome. Three natural tests of association based on the observed data depend on different numbers of nuisance parameters. Most promising is a test based on the ratio of proportions of the observed feature. We illustrate this test using a study of the association of imaging parameters with genetic features in subjects with oligodendroglioma, a common brain tumour. In this study, calcification, a feature related to the imaging parameters, was potentially misclassified as not present.

A computationally simple bivariate survival estimator for efficacy and safety

Scholtens, D., & Betensky, R. A. (n.d.).

Publication year

2006

Journal title

Lifetime Data Analysis

Volume

12

Issue

3

Page(s)

365-387
Abstract
Abstract
Both treatment efficacy and safety are typically the primary endpoints in Phase II, and even in some Phase III, clinical trials. Efficacy is frequently measured by time to response, death, or some other milestone event and thus is a continuous, possibly censored, outcome. Safety, however, is frequently measured on a discrete scale; in Eastern Cooperative Oncology Group clinical trial E2290, it was measured as the number of weekly rounds of chemotherapy that were tolerable to colorectal cancer patients. For the joint analysis of efficacy and safety, we propose a non-parametric, computationally simple estimator for the bivariate survival function when one time-to-event is continuous, one is discrete, and both are subject to right-censoring. The bivariate censoring times may depend on each other, but they are assumed to be independent of both event times. We derive a closed-form covariance estimator for the survivor function which allows for inference to be based on any of several possible statistics of interest. In addition, we derive its covariance with respect to calendar time of analysis, allowing for its use in sequential studies.

A computationally tractable multivariate random effects model for clustered binary data

Coull, B. A., Houseman, E. A., & Betensky, R. A. (n.d.).

Publication year

2006

Journal title

Biometrika

Volume

93

Issue

3

Page(s)

587-599
Abstract
Abstract
We consider a multivariate random effects model for clustered binary data that is useful when interest focuses on the association structure among clustered observations. Based on a vector of gamma random effects and a complementary log-log link function, the model yields a likelihood that has closed form, making a frequentist approach to model-fitting straightforward. This closed form yields several advantages over existing methods, including easy inspection of model identifiability and straightforward adjustment for nonrandom ascertainment of subjects, such as that which occurs in family studies of disease aggregation. We use the proposed model to analyse two different binary datasets concerning disease outcome data from a familial aggregation study of breast and ovarian cancer in women and loss of heterozygosity outcomes from a brain tumour study.

A pseudolikelihood approach for simultaneous analysis of array comparative genomic hybridizations

Engler, D. A., Mohapatra, G., Louis, D. N., & Betensky, R. A. (n.d.).

Publication year

2006

Journal title

Biostatistics

Volume

7

Issue

3

Page(s)

399-421
Abstract
Abstract
DNA sequence copy number has been shown to be associated with cancer development and progression. Array-based comparative genomic hybridization (aCGH) is a recent development that seeks to identify the copy number ratio at large numbers of markers across the genome. Due to experimental and biological variations across chromosomes and hybridizations, current methods are limited to analyses of single chromosomes. We propose a more powerful approach that borrows strength across chromosomes and hybridizations. We assume a Gaussian mixture model, with a hidden Markov dependence structure and with random effects to allow for intertumoral variation, as well as intratumoral clonal variation. For ease of computation, we base estimation on a pseudolikelihood function. The method produces quantitative assessments of the likelihood of genetic alterations at each clone, along with a graphical display for simple visual interpretation. We assess the characteristics of the method through simulation studies and analysis of a brain tumor aCGH data set. We show that the pseudolikelihood approach is superior to existing methods both in detecting small regions of copy number alteration and in accurately classifying regions of change when intratumoral clonal variation is present. Software for this approach is available at http://www.biostat.harvard.edu/∼betensky/papers.html.

Activation of STAT3, MAPK, and AKT in malignant astrocytic gliomas: Correlation with EGFR status, tumor grade, and survival

Mizoguchi, M., Betensky, R. A., Batchelor, T. T., Bernay, D. C., Louis, D. N., & Nutt, C. L. (n.d.).

Publication year

2006

Journal title

Journal of Neuropathology and Experimental Neurology

Volume

65

Issue

12

Page(s)

1181-1188
Abstract
Abstract
Diffuse astrocytic gliomas are the most common human glial tumors with glioblastoma being the most malignant form. Epidermal growth factor receptor (EGFR) gene amplification is one of the most common genetic changes in glioblastoma and can lead to the activation of various downstream signaling molecules, including STAT3, MAPK, and AKT. In this study, we investigated the activation status of these 3 signaling molecules as well as wild-type (EGFRwt) and mutant (EGFRvIII) EGFR in 82 malignant astrocytic gliomas (55 glioblastomas and 27 anaplastic astrocytomas) using immunohistochemistry. The presence of EGFRwt, but not EGFRvIII, immunopositivity correlated significantly with prevalent EGFR gene amplification in glioblastomas. STAT3 and AKT activation correlated significantly with EGFR status, although the correlation for p-STAT3 was attributed exclusively to EGFRvIII. The distribution of these 3 activated molecules varied significantly with tumor grade; although activation of STAT3 was essentially identical between anaplastic astrocytomas and glioblastomas, an increase in the activation of MAPK and AKT appeared to correlate with the progression of anaplastic astrocytoma to glioblastoma. Finally, activated STAT3 and AKT were marginally predictive of improved and worse prognosis, respectively. Taken together, these findings begin to elucidate the interrelationship between these signaling pathways in astrocytic gliomas in vivo.

AKT activation in human glioblastomas enhances proliferation via TSC2 and S6 kinase signaling

Riemenschneider, M. J., Betensky, R. A., Pasedag, S. M., & Louis, D. N. (n.d.).

Publication year

2006

Journal title

Cancer Research

Volume

66

Issue

11

Page(s)

5618-5623
Abstract
Abstract
Aberrant AKT (protein kinase B) signaling is common in many cancers, including glioblastoma. Current models suggest that AKT acts directly, or indirectly via the TSC complex, to activate the mammalian target of rapamycin (mTOR) as the main downstream mediator of AKT signaling. mTOR activation results in subsequent activation of S6K and STAT3, as well as suppression (i.e., phosphorylation) of 4E-BP1, leading to cell cycle progression and inhibition of apoptosis. Most studies of this pathway have used in vitro systems or tumor lysate-based approaches. We aimed to delineate these pathways in a primarily in situ manner using immunohistochemistry in a panel of 29 glioblastomas, emphasizing the histologie distribution of molecular changes. Within individual tumors, increased expression levels of p-TSC2, p-mTOR, P-4E-BP1, p-S6K, p-S6, and p-STAT3 were found in regions defined by elevated AKT activation. However, only TSC2, S6K, and S6 activation levels correlated significantly with AKT activation and clustered together in multidimensional scaling analyses. Ki-67 proliferation indices were significantly elevated in p-AKT-overexpressing regions, whereas expression of the apoptosis marker cleaved caspase 3 was generally low and not significantly different between the regions. These findings provide the first in vivo evidence for a close correlation between AKT and TSC2 phosphorylation levels in glioblastoma. Moreover, they suggest that downstream p-AKT effects are primarily mediated by S6 kinase signaling, thus enhancing proliferation rather than inhibiting apoptosis.

Analysis of co-aggregation of cancer based on registry data

Matthews, A. G., Betensky, R. A., Anton-Culver, H., Bowen, D., Griffin, C., Isaacs, C., Kasten, C., Mineau, G., Nayfield, S., Schildkraut, J., Strong, L., Weber, B., & Finkelstein, D. M. (n.d.).

Publication year

2006

Journal title

Community Genetics

Volume

9

Issue

2

Page(s)

87-92
Abstract
Abstract
Objective: An exploratory analysis of co-aggregation of cancers using registry-based data. Methods: We utilized sibships from over 18,000 families who had been recruited to the NCI-sponsored multi-institutional Cancer Genetics Network. The analysis assesses co-aggregation at the individual and family level and adjusts for ascertainment. Results: We found statistically significant familial co-aggregation of lung cancer with pancreatic (adjusted p < 0.001), prostate (adjusted p < 0.003), and colorectal cancers (adjusted p = 0.004). In addition, we found significant familial co-aggregation of pancreatic and colorectal cancers (adjusted p = 0.018), and co-aggregation of hematopoietic and (non-ovarian) gynecologic cancers (adjusted p = 0.01). Conclusion: This analysis identified familial aggregation of cancers for which a genetic component has yet to be established.

Characterization of amyloid deposition in the APPswe/PS1dE9 mouse model of Alzheimer disease

Garcia-Alloza, M., Robbins, E. M., Zhang-Nunes, S. X., Purcell, S. M., Betensky, R. A., Raju, S., Prada, C., Greenberg, S. M., Bacskai, B. J., & Frosch, M. P. (n.d.).

Publication year

2006

Journal title

Neurobiology of Disease

Volume

24

Issue

3

Page(s)

516-524
Abstract
Abstract
Transgenic mice carrying disease-linked forms of genes associated with Alzheimer disease often demonstrate deposition of the β-amyloid as senile plaques and cerebral amyloid angiopathy. We have characterized the natural history of β-amyloid deposition in APPswe/PS1dE9 mice, a particularly aggressive transgenic mouse model generated with mutant transgenes for APP (APPswe: KM594/5NL) and PS1 (dE9: deletion of exon 9). Ex vivo histochemistry showed Aβ deposition by 4 months with a progressive increase in plaque number up to 12 months and a similar increase of Aβ levels. In vivo multiphoton microscopy at weekly intervals showed increasing β-amyloid deposition as CAA and plaques. Although first appearing at an early age, CAA progressed at a significantly slower rate than in the Tg2576 mice. The consistent and early onset of β-amyloid accumulation in the APPswe/PS1dE9 model confirms its utility for studies of biochemical and pathological mechanisms underlying β-amyloid deposition, as well as exploring new therapeutic treatments.

Effects of unmeasured heterogeneity in the linear transformation model for censored data

Zhang, B., Li, Y., & Betensky, R. A. (n.d.).

Publication year

2006

Journal title

Lifetime Data Analysis

Volume

12

Issue

2

Page(s)

191-203
Abstract
Abstract
We investigate the effect of unobserved heterogeneity in the context of the linear transformation model for censored survival data in the clinical trials setting. The unobserved heterogeneity is represented by a frailty term, with unknown distribution, in the linear transformation model. The bias of the estimate under the assumption of no unobserved heterogeneity when it truly is present is obtained. We also derive the asymptotic relative efficiency of the estimate of treatment effect under the incorrect assumption of no unobserved heterogeneity. Additionally we investigate the loss of power for clinical trials that are designed assuming the model without frailty when, in fact, the model with frailty is true. Numerical studies under a proportional odds model show that the loss of efficiency and the loss of power can be substantial when the heterogeneity, as embodied by a frailty, is ignored.

Expression of oligodendroglial and astrocytic lineage markers in diffuse gliomas: Use of YKL-40, ApoE, ASCL1, and NKX2-2

Rousseau, A., Nutt, C. L., Betensky, R. A., Iafrate, A. J., Han, M., Ligon, K. L., Rowitch, D. H., & Louis, D. N. (n.d.).

Publication year

2006

Journal title

Journal of Neuropathology and Experimental Neurology

Volume

65

Issue

12

Page(s)

1149-1156
Abstract
Abstract
The phenotypic heterogeneity of astrocytic and oligodendroglial tumor cells complicates establishing accurate diagnostic criteria, and lineage-specific markers would facilitate diagnosis of glioma subtypes. Based on data from the literature and from expression microarrays, we selected molecules relevant to gliogenesis and glial lineage specificity and then used immunohistochemistry to assess expression of these molecules in 55 diffuse gliomas, including 8 biphasic oligoastrocytomas, 21 oligodendrogliomas (all with 1p/19qloss), 21 astrocytomas, and 5 glioblastomas. For the astrocytic lineage markers (GFAP, YKL-40, and ApoE), GFAP expression was significantly higher in the astrocytic component of oligoastrocytomas compared with the oligodendroglial part; similar patterns were detected for YKL-40 and ApoE, although the differences were not significant. GFAP, YKL-40, and ApoE reliably distinguished grade II-III oligodendrogliomas from grade II-IV astrocytomas (p < 0.0001, p = 0.002, and p < 0.0001, respectively). Among the oligodendroglial lineage markers (Olig2, Sox10, ASCL1, and NKX2-2), ASCL1 and NKX2-2 displayed significantly different immunostaining between oligodendrogliomas and astrocytomas (p = 0.017 and 0.004, respectively), but none clearly differentiated between the 2 glial populations of oligoastrocytomas. In addition to GFAP, therefore, YKL-40, ApoE, ASCL1, and NKX2-2 represent promising tumor cell markers to distinguish oligodendrogliomas from astrocytomas.

Family History of Diabetes Is a Major Determinant of Endothelial Function

Goldfine, A. B., Beckman, J. A., Betensky, R. A., Devlin, H., Hurley, S., Varo, N., Schonbeck, U., Patti, M. E., & Creager, M. A. (n.d.).

Publication year

2006

Journal title

Journal of the American College of Cardiology

Volume

47

Issue

12

Page(s)

2456-2461
Abstract
Abstract
Objectives: We evaluated whether endothelial dysfunction was present in nondiabetic persons with a family history (FH) of diabetes and assessed its relationship with insulin resistance and atherosclerosis risk factors. Background: Atherosclerosis is frequently present when type 2 diabetes (T2D) is first diagnosed. Endothelial dysfunction contributes to atherogenesis. Methods: Oral glucose tolerance and brachial artery flow-mediated, endothelium-dependent vasodilation (EDV) were assessed in 38 nondiabetic subjects; offspring of two parents with T2D (FH+) or with no first-degree relative with diabetes (FH-). Results: Although fasting glucose was higher in FH+ than FH- (5.3 ± 0.1 mmol/l vs. 4.9 ± 0.1 mmol/l, p < 0.03), glycemic burden assessed as 2-h or area-under-the-curve glucose after glucose load or glycosylated hemoglobin (HbA1c), and measures of insulin sensitivity or inflammation did not differ. Brachial artery flow-mediated EDV was reduced in FH+ (7.1 ± 0.9% vs. 11.7 ± 1.6%, p < 0.02), with no difference in nitroglycerin-induced endothelium-independent vasodilatation. In the combined cohort, only FH+ (r2 = 0.12, p < 0.02) and HbA1c (r2 = 0.14, p < 0.02) correlated with EDV. Insulin resistance, assessed by tertile of homeostasis model assessment of insulin resistance (HOMA-IR), was associated with impaired endothelium-dependent vasodilatation in FH- (p < 0.03, analysis of variance), but not in FH+, as even the most insulin-sensitive FH+ offspring had diminished endothelial function. In multiple regression analysis, including established cardiac risk factors, blood pressure and lipids, HbA1c, and HOMA-IR, FH remained a significant determinant of EDV (p = 0.04). Conclusions: Bioavailability of nitric oxide is lower in persons with a strong FH of T2D. Glycemic burden, even in the nondiabetic range, can contribute to endothelial dysfunction. Abnormalities of endothelial function may contribute to atherosclerosis before development of overt diabetes.

Feature-specific penalized latent class analysis for genomic data

Houseman, E. A., Coull, B. A., & Betensky, R. A. (n.d.).

Publication year

2006

Journal title

Biometrics

Volume

62

Issue

4

Page(s)

1062-1070
Abstract
Abstract
Genomic data are often characterized by a moderate to large number of categorical variables observed for relatively few subjects. Some of the variables may be missing or noninformative. An example of such data is loss of heterozygosity (LOH), a dichotomous variable, observed on a moderate number of genetic markers. We first consider a latent class model where, conditional on unobserved membership in one of k classes, the variables are independent with probabilities determined by a regression model of low dimension q. Using a family of penalties including the ridge and LASSO, we extend this model to address higher-dimensional problems. Finally, we present an orthogonal map that transforms marker space to a space of "features" for which the constrained model has better predictive power. We demonstrate these methods on LOH data collected at 19 markers from 93 brain tumor patients. For this data set, the existing unpenalized latent class methodology does not produce estimates. Additionally, we show that posterior classes obtained from this method are associated with survival for these patients.

Glioma test array for use with formalin-fixed, paraffin-embedded tissue: Array comparative genomic hybridization correlates with loss of heterozygosity and fluorescence in situ hybridization

Mohapatra, G., Betensky, R. A., Miller, E. R., Carey, B., Gaumont, L. D., Engler, D. A., & Louis, D. N. (n.d.).

Publication year

2006

Journal title

Journal of Molecular Diagnostics

Volume

8

Issue

2

Page(s)

268-276
Abstract
Abstract
Array-based comparative genomic hybridization (aCGH) is a powerful, high-throughput tool for whole genome analysis. Until recently, aCGH could only be reproducibly performed on frozen tissue samples and with significant tissue amounts. For brain tumors however, paraffin-embedded tissue blocks from small stereotactic biopsies may be the only tissue routinely available. The development of methods to analyze formalin-fixed, paraffin-embedded (FFPE) material therefore has the potential to impact molecular diagnosis in a significant way. To this end, we constructed a BAC array representing chromosomes 1, 7, 19, and X because 1p/19q deletion and EGFR gene amplification provide clinically relevant information for glioma diagnosis. We also optimized a two-step labeling procedure using an amine-modified nucleotide for generating aCGH probes. Using this approach, we analyzed a series of 28 FFPE oligodendroglial tumors for alterations of chromosomes 1, 7, and 19. We also independently assayed these tumors for 1p/19q deletion by fluorescence in situ hybridization and by loss of heterozygosity analyses. The concordance between aCGH, standard loss of heterozygosity and fluorescence in situ hybridization was nearly 100% for the chromosomes analyzed. These results suggest that aCGH could offer an improved molecular diagnostic approach for gliomas because of its ability to detect clinically relevant molecular alterations in small FFPE specimens.

Hospital volume versus outcome: An unusual example of bivariate association

Betensky, R. A., Christian, C. K., Gustafson, M. L., Daley, J., & Zinner, M. J. (n.d.).

Publication year

2006

Journal title

Biometrics

Volume

62

Issue

2

Page(s)

598-604
Abstract
Abstract
The Leapfrog Group, a consortium of more than 100 large employers, purchasing coalitions, and states that collectively provide health insurance to more than 33 million people, convened in 2000 with the goal of using market forces to improve the quality of healthcare. The resulting Leapfrog initiative suggested selective referral of complex procedures to high-volume hospitals and set volume thresholds for five procedures. This was based on the hypothesis that low-volume hospitals have higher mortality, which can be viewed in simplified statistical terms as the hypothesis that the binomial p is a decreasing function of n. The analysis of the correlation between hospitals' standardized mortality ratios (SMR, i.e., the ratio of observed to expected deaths) and hospitals' procedural volumes is revealing about the volume/mortality hypothesis. This presents an unusual pedagogic example in which the detection of correlation in the presence of nonlinear dependence is of primary interest, and thus the Pearson correlation is ideally suited. The frequently preferred nonparametric measures of bivariate association are inappropriate as they are unable to discriminate between correlation and dependence.

Immunoglobulin gene rearrangement analysis in cerebrospinal fluid of patients with lymphoproliferative processes

Baehring, J. M., Hochberg, F. H., Betensky, R. A., Longtine, J., & Sklar, J. (n.d.).

Publication year

2006

Journal title

Journal of the Neurological Sciences

Volume

247

Issue

2

Page(s)

208-216
Abstract
Abstract
Objective: To determine the sensitivity and specificity of clonal immunoglobulin heavy chain gene rearrangement (IGHR) analysis in the distinction of benign and malignant lymphoproliferative diseases. Methods: A retrospective analysis was conducted of patients in whom a malignant lymphoproliferative process was suspected. Cells of CSF samples were collected by centrifugation, resuspended in 100 μl of the supernatant and boiled. A 10 μl aliquot of this lysate served as template for semi-nested polymerase chain reaction using variable and joining region consensus primers. PCR products were analyzed by polyacrylamide gel electrophoresis. Cytopathological diagnosis and flow cytometry results were recorded. Sensitivity and specificity of IGHR analysis, cytopathology and flow cytometry were calculated. Results: Eleven patients (12 specimens) had involvement of leptomeninges at the time of lumbar puncture. Another 25 cases (27 specimens) had normal CSF findings or were diagnosed with benign lymphoproliferative conditions. Sensitivity of CSF cytopathology, flow cytometry and IGHR analysis were 0.27 [95% confidence interval 0.06, 0.61], 0.1 [0.003, 0.45] and 0.58 [0.28, 0.85]. Specificity was 1 [0.86, 1], 0.95 [0.77, 1.0] and 0.85 [0.66, 0.96]. Interpretation: IGHR analysis appears to be a useful addition to morphological and flow cytometry analysis of cerebrospinal fluid in the evaluation of CNS lymphoproliferative processes.

Kinetics of cerebral amyloid angiopathy progression in a transgenic mouse model of Alzheimer disease

Robbins, E. M., Betensky, R. A., Domnitz, S. B., Purcell, S. M., Garcia-Alloza, M., Greenberg, C., Rebeck, G. W., Hyman, B. T., Greenberg, S. M., Frosch, M. P., & Bacskai, B. J. (n.d.).

Publication year

2006

Journal title

Journal of Neuroscience

Volume

26

Issue

2

Page(s)

365-371
Abstract
Abstract
Cerebral amyloid angiopathy (CAA), the deposition of cerebrovascular β-amyloid (Aβ) in the walls of arterial vessels, has been implicated in hemorrhagic stroke and is present in most cases of Alzheimer disease. Previous studies of the progression of CAA in humans and animal models have been limited to the comparison of pathological tissue from different brains at single time points. Our objective was to visualize in real time the initiation and progression of CAA in Tg2576 mice by multiphoton microscopy through cranial windows. Affected vessels were labeled by methoxy-X04, a fluorescent dye that selectively binds cerebrovascular β-amyloid and plaques. With serial imaging sessions spaced at weekly intervals, we were able to observe the earliest appearance of CAA in leptomeningeal arteries as multifocal deposits of band-like Aβ. Over subsequent imaging sessions, we were able to identify growth of these deposits (propagation), as well as appearance of new bands (additional initiation events). Statistical modeling of the data suggested that as the extent of CAA progressed in this vascular bed, there was increased prevalence of propagation over initiation. During the early phases of CAA development, the overall pathology burden progressed at a rate of 0.35% of total available vessel area per day (95% confidence interval, 0.3-0.4%). The consistent rate of disease progression implies that this model is amenable to investigations of therapeutic interventions.

Methods to classify familial relationships in the presence of laboratory errors, without parental data

Zhang, B., & Betensky, R. A. (n.d.).

Publication year

2006

Journal title

Human Genetics

Volume

119

Issue

6

Page(s)

642-648
Abstract
Abstract
We consider the problem of accurate classification of family relationship in the presence of laboratory error without parental data. We first propose an adjusted version of the test statistic proposed by Ehm and Wagner based on the summation over a large number of genetics markers. We then propose use of the Bayes factor as a classification rule. We prove theoretically that the Bayes factor is the optimal classification rule in that the total classification error is minimized. We show via simulations that both the adjusted Ehm and Wagner method and Bayes factor classification rule reduce misclassification errors, and that the Bayes factor classification rule is robust against under-estimation or over-estimation of laboratory errors. For monozygotic twins versus dizygotic twins, the correct classification rate of the Bayes rule is over 99%. For full-siblings versus half-siblings, the Bayes factor classification rule generally outperforms Ehm and Wagner's method (in Am J Hum Genet 62:181-188, 1998, especially when full-sibling proportion is high.

Outcomes in a series of 103 retroperitoneal sarcomas

Pierie, J. P., Betensky, R. A., Choudry, U., Willett, C. G., Souba, W. W., & Ott, M. J. (n.d.).

Publication year

2006

Journal title

European Journal of Surgical Oncology

Volume

32

Issue

10

Page(s)

1235-1241
Abstract
Abstract
Aims: To report the effect on outcome of selection in patients receiving intra-operative electron beam radiation (IOERT) and external beam radiation therapy (EBRT). Methods: One hundred and three patients treated for primary RS were studied. Median follow-up was 27 months. Clinical presentation, tumor characteristics, and treatment methods were analyzed to determine impact on survival and recurrence and if selection was occurring. Results: Mean age was 55 ± 17 years. Mean tumor size was 15 ± 6 cm and 88 were high-grade. Complete gross tumor resection (CR) occurred in 62 patients and improved survival vs. both debulking (p = 0.0005) and biopsy (p < 0.0001). The 5- and 10-year survival rates were 62% and 52% for those with CR vs. 29% and 20% after incomplete resection. Among the 62 CR patients, there was selection to receive additional EBRT ± IOERT in patients with high-grade tumors (p = 0.005) and/or microscopically positive margins (p = 0.011). In these high-risk patients there was a trend for IOERT to further augment survival vs. EBRT alone and to increase the time to both local and distant recurrences (p = 0.036). Conclusions: Complete gross resection is the primary form of curative treatment for retroperitoneal sarcomas. Selection led to patients with high-risk tumors receiving additional radiation therapy. There appears to be a beneficial effect of IOERT plus EBRT in these high-risk patients after complete tumor resection.

Analysis of clonal immunoglobulin heavy chain rearrangements in ocular lymphoma

Baehring, J. M., Androudi, S., Longtine, J. J., Betensky, R. A., Sklar, J., Foster, C. S., & Hochberg, F. H. (n.d.).

Publication year

2005

Journal title

Cancer

Volume

104

Issue

3

Page(s)

591-597
Abstract
Abstract
BACKGROUND. The morphologic diagnosis of primary and metastatic intraocular lymphoma (IOL) was made difficult by the paucicellular specimens with fragile populations of lymphocytes retrieved through pars plana vitrectomy (PPV). The analysis of immunoglobulin heavy chain (IgH) gene rearrangements (AIGHR) was used as an adjunct to cytopathology and flow cytometry in systemic lymphoma. In IOL, the sensitivity and specificity of AIGHR are unknown. METHODS. The authors reviewed the clinical records of patients who underwent PPV for suspicion of IOL at the Immunology and Uveitis Service of the Massachusetts Eye and Ear Infirmary (Boston, MA) between 2000 and 2002. AIGHR was performed as a routine diagnostic test on cell lysates isolated from < 0.5 mL of vitreous fluid. The authors used seminested polymerase chain reaction (PCR) with consensus primers for the VDJ-region of the IgH gene. PCR products were analyzed by polyacrylamide gel electrophoresis. RESULTS. Thirty patients (37 specimens) with chronic vitritis and 17 patients (23 specimens) with IOL were included. The specificity of vitreous fluid cytopathology, flow cytometry, and AIGHR was 1.0, and the sensitivity values were 0.24, 0.36, and 0.64, respectively. AIGHR was negative in two patients for whom cytopathology or flow cytometry revealed the diagnosis of lymphoma. Clonal IGHR was found in four specimens classified as negative for lymphoma based on cytopathology and flow cytometry. CONCLUSIONS. AIGHR supplemented cytopathology and flow cytometry to increase the diagnostic yield in IOL.

Analysis of familial aggregation in the presence of varying family sizes

Matthews, A. G., Finkelstein, D. M., & Betensky, R. A. (n.d.).

Publication year

2005

Journal title

Journal of the Royal Statistical Society. Series C: Applied Statistics

Volume

54

Issue

5

Page(s)

847-862
Abstract
Abstract
Family studies are frequently undertaken as the first step in the search for genetic and/or environmental determinants of disease. Significant familial aggregation of disease is suggestive of a genetic aetiology for the disease and may lead to more focused genetic analysis. Of course, it may also be due to shared environmental factors. Many methods have been proposed in the literature for the analysis of family studies. One model that is appealing for the simplicity of its computation and the conditional interpretation of its parameters is the quadratic exponential model. However, a limiting factor in its application is that it is not reproducible, meaning that all families must be of the same size. To increase the applicability of this model, we propose a hybrid approach in which analysis is based on the assumption of the quadratic exponential model for a selected family size and combines a missing data approach for smaller families with a marginalization approach for larger families. We apply our approach to a family study of colorectal cancer that was sponsored by the Cancer Genetics Network of the National Institutes of Health. We investigate the properties of our approach in simulation studies. Our approach applies more generally to clustered binary data.

Cell-surface MuSK self-association: A crucial role for the putative signal sequence

Bianchetta, M. J., Betensky, R. A., & Cohen, J. B. (n.d.).

Publication year

2005

Journal title

Biochemistry

Volume

44

Issue

49

Page(s)

16229-16238
Abstract
Abstract
The receptor tyrosine kinase MuSK plays a crucial role-both as a signaling molecule and structurally-in the process of clustering nicotinic acetylcholine receptors at the neuromuscular junction. Immunofluorescence microscopy of transiently transfected fibroblasts has been used to visualize the cell-surface distribution of MuSK, which is found in discrete, punctate clusters. This distribution does not result from targeting of MuSK to identified plasma membrane subdomains, and MuSK's association with itself is specific, as MuSK clusters at the cell surface are segregated from clusters of other cotransfected receptor tyrosine kinases. A mutational analysis, using coexpressed pairs of MuSK mutants and chimeras, demonstrates that the putative signal peptide is both necessary and sufficient for association with MuSK. Removal of the intracellular domain or most of the extracellular domain, or replacement of the transmembrane domain, does not abolish MuSK self-association. The N-terminus of the MuSK protein, however, is sufficient to recruit another receptor tyrosine kinase to MuSK clusters. Quantitation and statistical analysis of the amount of colocalization between coexpressed MuSK mutants and chimeras confirm these results.

Contact

rebecca.betensky@nyu.edu 708 Broadway New York, NY, 10003