Growth Dynamics of Ductal Carcinoma in Situ Recapitulate Normal Breast Development

Ductal carcinoma in situ (DCIS) and invasive breast cancer share many morphologic, proteomic, and genomic alterations. Yet in contrast to invasive cancer, many DCIS tumors do not progress and may remain indolent over decades. To better understand the heterogenous nature of this disease, we reconstructed the growth dynamics of 18 DCIS tumors based on the geo-spatial distribution of their somatic mutations. The somatic mutation topographies revealed that DCIS is multiclonal and consists of spatially discontinuous subclonal lesions. Here we show that this pattern of spread is consistent with a new 'Comet' model of DCIS tumorigenesis, whereby multiple subclones arise early and nucleate the buds of the growing tumor. The discontinuous, multiclonal growth of the Comet model is analogous to the branching morphogenesis of normal breast development that governs the rapid expansion of the mammary epithelium during puberty. The branching morphogenesis-like dynamics of the proposed Comet model diverges from the canonical model of clonal evolution, and better explains observed genomic spatial data. Importantly, the Comet model allows for the clinically relevant scenario of extensive DCIS spread, without being subjected to the selective pressures of subclone competition that promote the emergence of increasingly invasive phenotypes. As such, the normal cell movement inferred during DCIS growth provides a new explanation for the limited risk of progression in DCIS and adds biologic rationale for ongoing clinical efforts to reduce DCIS overtreatment.

Archival single-cell genomics reveals persistent subclones during DCIS progression

Ductal carcinoma in situ (DCIS) is a common precursor of invasive breast cancer. Our understanding of its genomic progression to recurrent disease remains poor, partly due to challenges associated with the genomic profiling of formalin-fixed paraffin-embedded (FFPE) materials. Here, we developed Arc-well, a high-throughput single-cell DNA-sequencing method that is compatible with FFPE materials. We validated our method by profiling 40,330 single cells from cell lines, a frozen tissue, and 27 FFPE samples from breast, lung, and prostate tumors stored for 3-31 years. Analysis of 10 patients with matched DCIS and cancers that recurred 2-16 years later show that many primary DCIS had already undergone whole-genome doubling and clonal diversification and that they shared genomic lineages with persistent subclones in the recurrences. Evolutionary analysis suggests that most DCIS cases in our cohort underwent an evolutionary bottleneck, and further identified chromosome aberrations in the persistent subclones that were associated with recurrence.

Molecular classification and biomarkers of clinical outcome in breast ductal carcinoma in situ: Analysis of TBCRC 038 and RAHBT cohorts

Ductal carcinoma in situ (DCIS) is the most common precursor of invasive breast cancer (IBC), with variable propensity for progression. We perform multiscale, integrated molecular profiling of DCIS with clinical outcomes by analyzing 774 DCIS samples from 542 patients with 7.3 years median follow-up from the Translational Breast Cancer Research Consortium 038 study and the Resource of Archival Breast Tissue cohorts. We identify 812 genes associated with ipsilateral recurrence within 5 years from treatment and develop a classifier that predicts DCIS or IBC recurrence in both cohorts. Pathways associated with recurrence include proliferation, immune response, and metabolism. Distinct stromal expression patterns and immune cell compositions are identified. Our multiscale approach employed in situ methods to generate a spatially resolved atlas of breast precancers, where complementary modalities can be directly compared and correlated with conventional pathology findings, disease states, and clinical outcome.

Abstract PR002: Genomic predictor can discriminate between high- and low-risk DCIS

Introduction: Ductal carcinoma in situ (DCIS) is considered a non-obligate precursor of invasive ductal carcinoma. With the aim of preventing a subsequent invasive cancer, all DCIS lesions are currently treated with surgical excision often supplemented with radiotherapy (RT). To prevent DCIS over- or undertreatment, a reliable marker of DCIS invasiveness risk is urgently needed. Methods: We studied two large DCIS cohorts: the Sloane cohort, a prospective breast screening cohort from the UK (median follow-up of 11 years), and a Dutch population-based cohort (NKI, median follow-up of 13 years). FFPE tissue specimens from patients with pure primary DCIS after breast-conserving surgery (BCS) +/- RT that did develop a subsequent ipsilateral event (DCIS or invasive) were considered as cases, whereas patients that did not develop any form of recurrence up to the last follow-up or death were considered as controls. We performed copy number analysis (CNA) and RNAseq analysis on 229 cases (80 DCIS only recurrences) and 344 controls. Results: DCIS was classified into the PAM50 subtypes using RNAseq data which revealed an enrichment of luminal A phenotype in DCIS that did not recur (P = 0.01, Fisher Exact test). No single copy number aberration was more common in cases compared to controls. RNAseq data did not reveal any genes significantly over/under-expressed in cases versus controls after FDR correction. However, by limiting the analysis to samples that had not had RT and excluding pure DCIS recurrences, we could develop a penalized Cox model from RNAseq data. The model was trained on weighted samples (to correct for the biased sampling of the case-control dataset) from the NKI series with double loop cross-validation. The genes were selected using the Elastic net framework of penalization. Using this predicted hazard ratio, the samples were split into high, medium, and low-risk quantiles, with a recurrence risk of 23%, 7% and 2%, respectively at 5 years (p = 10-10, Wald test). The NKI-trained predictor was independently validated in the Sloane No RT no DCIS recurrence cohort (p = 0.02, Wald test). GSEA analysis revealed proliferation hallmarks enriched in the recurrence predictor (FDR = 0.058). The RNAseq predictor was more predictive of recurrence than PAM50, clinical features (Grade, Her2 and ER) and the 12-gene Oncotype DCIS score (p < 0.001, permutation test using the Wald statistic) in both the NKI and Sloane series. Conclusion: Genomic profiling of two independent series of DCIS with outcome data did not reveal any clear associations with recurrence until analysis was limited to a set of samples who had not had radiotherapy and DCIS recurrences were excluded. We then identified an RNAseq-based classifier that could differentiate primary DCIS in low-, medium-, and high-risk groups, and validated it in an independent cohort. This classifier, if validated in other datasets, will allow us to identify women who do not need intensive treatment for their DCIS.

Citation Format: Maria Roman Escorza, Michael Sheinman, Tycho Bismeijer, Ahmed A. Ahmed, Vandna Shah, Jeffrey R. Marks, Lorraine M. King, Anargyros Megalios, Lindy L. Visser, Marlous Hoogstrat, Helen R. Davies, Tapsi Kumar, Deborah Collyar, Hilary Stobart, Sarah Pinder, Nicholas N. Navin, Andrew Futreal, Serena Nik-Zainal, E. Shelley Hwang, Esther H. Lips, Alastair Thompson, Lodewyk F.A. Wessels, Jelle Wesseling, Elinor J. Sawyer. Genomic predictor can discriminate between high- and low-risk DCIS [abstract]. In: Proceedings of the AACR Special Conference on Rethinking DCIS: An Opportunity for Prevention?; 2022 Sep 8-11; Philadelphia, PA. Philadelphia (PA): AACR; Can Prev Res 2022;15(12 Suppl_1): Abstract nr PR002.

Spatial interplay of tissue hypoxia and T-cell regulation in ductal carcinoma in situ

Hypoxia promotes aggressive tumor phenotypes and mediates the recruitment of suppressive T cells in invasive breast carcinomas. We investigated the role of hypoxia in relation to T-cell regulation in ductal carcinoma in situ (DCIS). We designed a deep learning system tailored for the tissue architecture complexity of DCIS, and compared pure DCIS cases with the synchronous DCIS and invasive components within invasive ductal carcinoma cases. Single-cell classification was applied in tandem with a new method for DCIS ductal segmentation in dual-stained CA9 and FOXP3, whole-tumor section digital pathology images. Pure DCIS typically has an intermediate level of colocalization of FOXP3+ and CA9+ cells, but in invasive carcinoma cases, the FOXP3+ (T-regulatory) cells may have relocated from the DCIS and into the invasive parts of the tumor, leading to high levels of colocalization in the invasive parts but low levels in the synchronous DCIS component. This may be due to invasive, hypoxic tumors evolving to recruit T-regulatory cells in order to evade immune predation. Our data support the notion that hypoxia promotes immune tolerance through recruitment of T-regulatory cells, and furthermore indicate a spatial pattern of relocalization of T-regulatory cells from DCIS to hypoxic tumor cells. Spatial colocalization of hypoxic and T-regulatory cells may be a key event and useful marker of DCIS progression.

Abstract 5108: Copy number analysis of pure DCIS and association with recurrence

With the widespread adoption of breast cancer screening the incidence of pure ductal carcinoma in situ (DCIS) has increased. As DCIS is considered a non-obligate precursor of invasive ductal carcinoma most women with pure DCIS are treated with breast conserving surgery (BCS) +/- radiotherapy. However, for many this is likely to be overtreatment as only a minority will develop a subsequent ipsilateral recurrence. Studies also show that only ~60% of these ipsilateral recurrences are invasive disease with the remainder being pure DCIS. To predict which women are most likely to benefit from interventions, there is a need to identify biomarkers that are associated with invasive recurrence. Our aim was to assess whether copy number aberrations (CNAs) could be used to identify DCIS that was likely to recur as invasive disease or remain recurrence-free during long-time follow up.

We performed somatic copy number profiling on 309 pure DCIS samples that had not developed an ipsilateral event (controls), 198 that had developed subsequent ipsilateral invasive disease (INV-cases) and 58 that had developed subsequent ipsilateral pure DCIS (DCIS-cases). The samples were obtained from two large nation-wide cohorts: the Sloane cohort, a prospective breast screening cohort from the UK with a median follow up of 12.5 years and a Dutch population based cohort, with a median follow up of 13 years. CNAs were assessed using the CytoSNP array or low pass whole genome sequencing and analyzed using GISTIC.

Integrative cluster (IntClust) subtyping revealed that only 5 subtypes were well represented in DCIS compared to 10 in invasive disease and the distribution of clusters between INV-cases and controls was similar with the exception of IntClust 4, which was significantly more common in controls (P= 0.025, Fishers exact test). IntClust 4 is characterized to have low levels of genomic instability and a CNA-devoid. INV-cases were globally more aberrant than controls (P = 0.006, Wilcoxon test) as assessed by the chromosomal instability index (CIN) score. GISTIC identified 17 recurrent amplifications, 21 recurrent gains and 22 recurrent losses in the whole cohort. Six of these regions were more common in INV-cases compared to controls: amplifications at 17q24.1 and 8p11.23, losses at 1p36.13 and 11q23.2 ​and gains at 17q21.33​ and 16p (Nominal P < 0.05 and FDR < 0.1, Fishers exact test). Subgroup analysis of ER+, Her2- INV-cases versus controls revealed an additional differential CNA, amplification at 11q13.3 more common in cases.

DCIS-cases had similar CNAs to INV-cases and were more aberrant than controls in terms of CIN score (P < 0.037, Wilcoxon test) but not as aberrant as INV-cases.

In conclusion, we have identified potential CNAs that are associated with invasive recurrence. Further analysis will integrate gene expression with copy number data to identify which genes are being targeted by these CNAs in order to identify pathways important in progression of DCIS.

Citation Format: Ahmed A. Ahmed, Maria Roman-Escorza, Tycho Bismeijer, Michael Sheinman, Vandna Shah, Rana Shami, Jeffrey R. Marks, Lorraine M. King, Anargyros Megalios, Lindy L. Visser, Marlous Hoogstraat, Helen R. Davies, Tapsi Kumar, Deborah Collyar, Hilary Stobart, Sarah Pinder, Nicholas N. Navin, Andrew Futreal, Serena Nik-Zainal, E. Shelley Hwang, Lodewyk F. Wessels, Esther H. Lips, Alastair Thompson, Jelle Wesseling, Elinor J. Sawyer. Copy number analysis of pure DCIS and association with recurrence [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 5108.

Genomic analysis defines clonal relationships of ductal carcinoma in situ and recurrent invasive breast cancer

Ductal carcinoma in situ (DCIS) is the most common form of preinvasive breast cancer and, despite treatment, a small fraction (5-10%) of DCIS patients develop subsequent invasive disease. A fundamental biologic question is whether the invasive disease arises from tumor cells in the initial DCIS or represents new unrelated disease. To address this question, we performed genomic analyses on the initial DCIS lesion and paired invasive recurrent tumors in 95 patients together with single-cell DNA sequencing in a subset of cases. Our data show that in 75% of cases the invasive recurrence was clonally related to the initial DCIS, suggesting that tumor cells were not eliminated during the initial treatment. Surprisingly, however, 18% were clonally unrelated to the DCIS, representing new independent lineages and 7% of cases were ambiguous. This knowledge is essential for accurate risk evaluation of DCIS, treatment de-escalation strategies and the identification of predictive biomarkers.

Anomaly Detection of Calcifications in Mammography Based on 11,000 Negative Cases

In mammography, calcifications are one of the most common signs of breast cancer. Detection of such lesions is an active area of research for computer-aided diagnosis and machine learning algorithms. Due to limited numbers of positive cases, many supervised detection models suffer from overfitting and fail to generalize. We present a one-class, semi-supervised framework using a deep convolutional autoencoder trained with over 50,000 images from 11,000 negative-only cases. Since the model learned from only normal breast parenchymal features, calcifications produced large signals when comparing the residuals between input and reconstruction output images. As a key advancement, a structural dissimilarity index was used to suppress non-structural noises. Our selected model achieved pixel-based AUROC of 0.959 and AUPRC of 0.676 during validation, where calcification masks were defined in a semi-automated process. Although not trained directly on any cancers, detection performance of calcification lesions on 1,883 testing images (645 malignant and 1238 negative) achieved 75% sensitivity at 2.5 false positives per image. Performance plateaued early when trained with only a fraction of the cases, and greater model complexity or a larger dataset did not improve performance. This study demonstrates the potential of this anomaly detection approach to detect mammographic calcifications in a semi-supervised manner with efficient use of a small number of labeled images, and may facilitate new clinical applications such as computer-aided triage and quality improvement.

Abstract GS4-07: The Breast PreCancer Atlas DCIS genomic signatures define biology and correlate with clinical outcomes: An analysis of TBCRC 038 and RAHBT cohorts

Background. DCIS consists of a molecularly heterogeneous group of premalignant lesions, with variable risk of invasive progression. Understanding biomarkers for invasive progression could help individualize treatment recommendations based upon tumor biology. As part of the NCI Human Tumor Atlas Network (HTAN), we conducted comprehensive genomic analyses on two large DCIS case-control cohorts. Methods. We performed smart3-seq and low-pass whole genome sequencing on two independent, retrospective, longitudinally sampled DCIS case-control cohorts. TBCRC 038 was a multicenter cohort diagnosed with DCIS between 1998 and 2016 at one of the Translational Breast Cancer Research sites; the RAHBT (Resource of Archival Human Breast Tissue) cohort included women identified through the St. Louis Breast Tissue Repository, and the Women’s Health Repository diagnosed between 1997 and 2001. We studied the spectrum of molecular changes present and sought genomic predictors of subsequent ipsilateral breast events (iBEs: DCIS recurrence or invasive progression) in both DCIS epithelium and stroma in formalin fixed paraffin embedded tissue. We generated de novo tumor and stroma-centric subtypes for DCIS that represents fundamental transcriptomic organization. Copy number analysis was performed using low-pass DNA sequencing. Non-negative matrix factorization (NMF) was applied to the RNA expression of all coding genes to identify clusters. A negative-binomial regression model was used to identify differentially expressed genes. Results. We analyzed 677 DCIS samples from 481 patients with 7.1 years median follow-up. In TBCRC samples, we identified three clusters via NMF in TBCRC referred to as ER low, quiescent, and ER high. The ER-low cluster had significantly higher levels of ERBB2 and lower levels of ESR1 compared to quiescent and ER-high clusters. Quiescent cluster lesions were less proliferative and less metabolically active than ER high and ER low subtypes. These findings were replicated in the RAHBT cohort. Focusing on the stromal component of DCIS from laser capture microdissection in RAHBT samples, we identified four distinct DCIS-associated stromal clusters. A “normal-like” stromal cluster with ECM organization and PI3K-AKT signaling; a “collagen-rich” stromal cluster; a “desmoplastic” stromal cluster with high fibroblast and total myeloid abundance, mostly associated with macrophages and myeloid dendritic cells (mDC); and an “immune-dense” stromal cluster. Further, we compared differentially expressed genes in patients with or without subsequent iBEs within 5 years of diagnosis. Hypothesizing that the resulting 812 DE genes (DESeq2) represent multiple routes to subsequent iBEs, we leveraged NMF to identify paths to progression. In both TBCRC and RAHBT cohorts, poor outcome groups exhibited increased ER, MYC signaling, and oxidative phosphorylation, supporting that these pathways are important for DCIS recurrence and progression. Conclusion. Comprehensive genomic profiling in two independent DCIS cohorts with longitudinal outcomes shows distinct DCIS stromal expression patterns and immune cell composition. RNA expression profiles reveal underlying tumor biology that is associated with later iBEs in both cohorts. These studies provide new insight into DCIS biology and will guide the design of diagnostic strategies to prevent invasive progression.

Citation Format: Siri H Strand, Belén Rivero-Gutiérrez, Kathleen E Houlahan, Jose A Seoane, Lorraine M King, Tyler Risom, Lunden Simpson, Sujay Vennam, Aziz Khan, Timothy Hardman, Bryan E Harmon, Fergus J Couch, Kristalyn Gallagher, Mark Kilgore, Shi Wei, Angela DeMichele, Tari King, Priscilla F McAuliffe, Julie Nangia, Joanna Lee, Jennifer Tseng, Anna Maria Storniolo, Alastair Thompson, Gaorav Gupta, Robyn Burns, Deborah J Veis, Katherine DeSchryver, Chunfang Zhu, Magdalena Matusiak, Jason Wang, Shirley X Zhu, Jen Tappenden, Daisy Yi Ding, Dadong Zhang, Jingqin Luo, Shu Jiang, Sushama Varma, Cody Straub, Sucheta Srivastava, Christina Curtis, Rob Tibshirani, Robert Michael Angelo, Allison Hall, Kouros Owzar, Kornelia Polyak, Carlo Maley, Jeffrey R Marks, Graham A Colditz, E Shelley Hwang, Robert B West. The Breast PreCancer Atlas DCIS genomic signatures define biology and correlate with clinical outcomes: An analysis of TBCRC 038 and RAHBT cohorts [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr GS4-07.

Abstract P1-22-05: Identifying predictors of invasive recurrence based on molecular profiles of DCIS lesions

Introduction Ductal carcinoma in situ (DCIS) is a non-obligate precursor of invasive breast cancer. Patients with DCIS are routinely treated by breast-conserving surgery often supplemented by radiotherapy, although many will never develop invasive disease. To date, no robust predictors of invasive breast cancer recurrence following DCIS have been identified. In our efforts to find such predictors, we performed gene expression, copy number and mutation analysis on two large DCIS cohorts with long-term follow-up. Methods Two nested case control series were analyzed, where cases are defined as DCIS with a subsequent invasive breast cancer and controls remained disease free during follow up. Cases and controls were matched on age and on follow up duration and were derived from two nation-wide cohort studies. The Sloane cohort is a prospective breast screening cohort from the UK, median follow up is 6 years (range 1-10). The Dutch cohort is population-based and had a median follow up of 13 years (range 2-23). We performed copy number analysis using CytoSNP array or low pass whole genome sequencing (lpWGS) on 310 controls and 196 cases, and RNA-seq on 295 controls and 206 cases. Results First analyses on the copy number data suggest that cases are genetically more aberrant with multiple regions of amplification compared to controls (p < 0.05). RNA-seq was used to classify DCIS into the PAM50 subtypes which did not appear to be predictive of recurrence. Initial RNA-seq analysis did not show consistent gene expression differences between cases and controls in the Sloane or Dutch cohorts, possibly explained by differences in clinical characteristics of the cohorts. A new computational method has been developed accounting for the differences in follow-up times, results will be presented at SABCS. Targeted sequencing revealed that the most common mutations were in PIK3CA and TP53, but there was no association with recurrence. Conclusion Only small molecular differences were identified between DCIS that recurs as invasive breast cancer and DCIS that remains disease-free. Currently, we are seeking to identify reproducible differences by a combined analysis of two population-based cohorts in a time dependent fashion. These will be presented at the SABCS. This work was supported by Cancer Research UK and by KWF Dutch Cancer Society (ref.C38317/A24043)

Citation Format: Tycho Bismeijer, Ahmed A Ahmed, Michael Sheinman, Maria Roman-Escorza, Vandna Shah, Jeffrey R Marks, Lorraine M King, Anargyros Megalios, Lindy L Visser, Marlous Hoogstraat, Helen R Davies, Tapsi Kumar, Deborah Collyar, Hilary Stobart, Nicholas N Navin, Andrew Futreal, Serena Nik-Zainal, Shelley Hwang, Esther H Lips, Alastair Thompson, Lodewyk FA Wessels, Elinor J Sawyer, Jelle Wesseling, Grand Challenge PRECISION Consortium. Identifying predictors of invasive recurrence based on molecular profiles of DCIS lesions [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P1-22-05.

Prediction of Upstaging in Ductal Carcinoma in Situ Based on Mammographic Radiomic Features

Background Improving diagnosis of ductal carcinoma in situ (DCIS) before surgery is important in choosing optimal patient management strategies. However, patients may harbor occult invasive disease not detected until definitive surgery. Purpose To assess the performance and clinical utility of mammographic radiomic features in the prediction of occult invasive cancer among women diagnosed with DCIS on the basis of core biopsy findings. Materials and Methods In this Health Insurance Portability and Accountability Act-compliant retrospective study, digital magnification mammographic images were collected from women who underwent breast core-needle biopsy for calcifications that was performed at a single institution between September 2008 and April 2017 and yielded a diagnosis of DCIS. The database query was directed at asymptomatic women with calcifications without a mass, architectural distortion, asymmetric density, or palpable disease. Logistic regression with regularization was used. Differences across training and internal test set by upstaging rate, age, lesion size, and estrogen and progesterone receptor status were assessed by using the Kruskal-Wallis or χ² test. Results The study consisted of 700 women with DCIS (age range, 40-89 years; mean age, 59 years ± 10 [standard deviation]), including 114 with lesions (16.3%) upstaged to invasive cancer at subsequent surgery. The sample was split randomly into 400 women for the training set and 300 for the testing set (mean ages: training set, 59 years ± 10; test set, 59 years ± 10; P = .85). A total of 109 radiomic and four clinical features were extracted. The best model on the test set by using all radiomic and clinical features helped predict upstaging with an area under the receiver operating characteristic curve of 0.71 (95% CI: 0.62, 0.79). For a fixed high sensitivity (90%), the model yielded a specificity of 22%, a negative predictive value of 92%, and an odds ratio of 2.4 (95% CI: 1.8, 3.2). High specificity (90%) corresponded to a sensitivity of 37%, positive predictive value of 41%, and odds ratio of 5.0 (95% CI: 2.8, 9.0). Conclusion Machine learning models that use radiomic features applied to mammographic calcifications may help predict upstaging of ductal carcinoma in situ, which can refine clinical decision making and treatment planning. © RSNA, 2022.

A new method to accurately identify single nucleotide variants using small FFPE breast samples

Most tissue collections of neoplasms are composed of formalin-fixed and paraffin-embedded (FFPE) excised tumor samples used for routine diagnostics. DNA sequencing is becoming increasingly important in cancer research and clinical management; however it is difficult to accurately sequence DNA from FFPE samples. We developed and validated a new bioinformatic pipeline to use existing variant-calling strategies to robustly identify somatic single nucleotide variants (SNVs) from whole exome sequencing using small amounts of DNA extracted from archival FFPE samples of breast cancers. We optimized this strategy using 28 pairs of technical replicates. After optimization, the mean similarity between replicates increased 5-fold, reaching 88% (range 0-100%), with a mean of 21.4 SNVs (range 1-68) per sample, representing a markedly superior performance to existing tools. We found that the SNV-identification accuracy declined when there was less than 40 ng of DNA available and that insertion-deletion variant calls are less reliable than single base substitutions. As the first application of the new algorithm, we compared samples of ductal carcinoma in situ of the breast to their adjacent invasive ductal carcinoma samples. We observed an increased number of mutations (paired-samples sign test, P < 0.05), and a higher genetic divergence in the invasive samples (paired-samples sign test, P < 0.01). Our method provides a significant improvement in detecting SNVs in FFPE samples over previous approaches.

Abstract PR02: Inferring the evolutionary dynamics of ductal carcinoma in situ through multi-regional sequencing and mathematical modeling

Introduction. The natural history of preinvasive breast cancer, or ductal carcinoma in situ (DCIS) remains poorly understood. Overcoming this gap would allow risk-appropriate treatment for patients diagnosed with DCIS. We used a multiregional sequencing approach in combination with mathematical modeling to characterize the evolutionary dynamics of DCIS initiation and progression. Methods. We analyzed a cohort of 18 patients diagnosed with DCIS, either with (n=9) or without (n=9) synchronous invasive cancer. Based on whole exome sequencing, tumor-specific mutation panels were constructed, each targeting 29-75 mutations (median: 60). From each tumor, and using selective ultraviolet radiation fractionation (SURF), we microdissected small spots (encompassing 1-3 duct cross-sections) from 3-4 spatially separated microscope sections (mean slide separation: 1.25cm, range: 0.34-6.0cm). Spots were spatially registered and genotyped based on targeted sequencing of the tumor-specific mutation panels. For each tumor, we performed unsupervised clonal deconvolution of the spot genotypes (CloneFinder) and constructed phylogenetic subclone trees. To quantify the spatial patterns of subclonal mutations, we introduced a dispersion index (DI), ranging from low (DI=0%) to high (DI=100%). To provide a spatio-temporal context for the heterogeneity patterns we developed a family of stochastic mathematical models of DCIS initiation and progression. Thereby, we embedded the evolutionary dynamics of tumor cell expansion in the branching topology of mammary ductal trees. Results. A total of 485 microdissected spots (median per tumor: 23, range: 10-50) were spatially registered and sequenced (median depth: 9,000x). All tumors were multiclonal, containing a median of 5 subclones (range: 2-14). Surprisingly, the correlation between spatial and genomic distances of spots was low. Individual subclones were diffusely dispersed across tumors. DCIS with synchronous DCIS and invasive cancer (mixed DCIS) had a higher mutation dispersion (DI=84.7%) than those without (pure DCIS, DI=70.5%; p=0.03, Wilcoxon rank-sum test). Mixed DCIS also had a higher fraction of spots containing more than one subclone than pure DCIS (median: 30.4% vs 0%, p=0.03). Among 7 mixed DCIS with invasive spots, 5 showed evidence of multiclonal invasion, that is more than one invading subclones were found in both in situ and invasive regions of the tumor. Mathematical modeling analyses show that the observed spatial patterns of genetic heterogeneity are consistent with a single expansion of mixing subclones across the ductal tree architecture. Conclusions. Our findings provide novel insights into the early growth and invasion dynamics of DCIS lesions. Further, we identified potential evolutionary markers for the delineation between indolent (pure) and aggressive (mixed) DCIS. This constitutes an important step towards identification of patients with low-risk DCIS who could be appropriately managed with less aggressive treatment.

Citation Format: Marc D. Ryser, Inmaculada C. Sorribes, Matthew Greenwald, Ethan Wu, Allison Hall, Diego Mallo, Lorraine M. King, Timothy Hardman, Lunden Simpson, Carlo C. Maley, Jeffrey R. Marks, Darryl Shibata, E. Shelley Hwang. Inferring the evolutionary dynamics of ductal carcinoma in situ through multi-regional sequencing and mathematical modeling [abstract]. In: Proceedings of the AACR Virtual Special Conference on Tumor Heterogeneity: From Single Cells to Clinical Impact; 2020 Sep 17-18. Philadelphia (PA): AACR; Cancer Res 2020;80(21 Suppl):Abstract nr PR02.

Minimal barriers to invasion during human colorectal tumor growth

Intra-tumoral heterogeneity (ITH) could represent clonal evolution where subclones with greater fitness confer more malignant phenotypes and invasion constitutes an evolutionary bottleneck. Alternatively, ITH could represent branching evolution with invasion of multiple subclones. The two models respectively predict a hierarchy of subclones arranged by phenotype, or multiple subclones with shared phenotypes. We delineate these modes of invasion by merging ancestral, topographic, and phenotypic information from 12 human colorectal tumors (11 carcinomas, 1 adenoma) obtained through saturation microdissection of 325 small tumor regions. The majority of subclones (29/46, 60%) share superficial and invasive phenotypes. Of 11 carcinomas, 9 show evidence of multiclonal invasion, and invasive and metastatic subclones arise early along the ancestral trees. Early multiclonal invasion in the majority of these tumors indicates the expansion of co-evolving subclones with similar malignant potential in absence of late bottlenecks and suggests that barriers to invasion are minimal during colorectal cancer growth.

Prediction of Upstaged Ductal Carcinoma In Situ Using Forced Labeling and Domain Adaptation

OBJECTIVE

The goal of this study is to use adjunctive classes to improve a predictive model whose performance is limited by the common problems of small numbers of primary cases, high feature dimensionality, and poor class separability. Specifically, our clinical task is to use mammographic features to predict whether ductal carcinoma in situ (DCIS) identified at needle core biopsy will be later upstaged or shown to contain invasive breast cancer.

METHODS

To improve the prediction of pure DCIS (negative) versus upstaged DCIS (positive) cases, this study considers the adjunctive roles of two related classes: atypical ductal hyperplasia (ADH), a non-cancer type of breast abnormity, and invasive ductal carcinoma (IDC), with 113 computer vision based mammographic features extracted from each case. To improve the baseline Model A's classification of pure vs. upstaged DCIS, we designed three different strategies (Models B, C, D) with different ways of embedding features or inputs.

RESULTS

Based on ROC analysis, the baseline Model A performed with AUC of 0.614 (95% CI, 0.496-0.733). All three new models performed better than the baseline, with domain adaptation (Model D) performing the best with an AUC of 0.697 (95% CI, 0.595-0.797).

CONCLUSION

We improved the prediction performance of DCIS upstaging by embedding two related pathology classes in different training phases.

SIGNIFICANCE

The three new strategies of embedding related class data all outperformed the baseline model, thus demonstrating not only feature similarities among these different classes, but also the potential for improving classification by using other related classes.

Prediction of Occult Invasive Disease in Ductal Carcinoma in Situ Using Deep Learning Features

PURPOSE

The aim of this study was to determine whether deep features extracted from digital mammograms using a pretrained deep convolutional neural network are prognostic of occult invasive disease for patients with ductal carcinoma in situ (DCIS) on core needle biopsy.

METHODS

In this retrospective study, digital mammographic magnification views were collected for 99 subjects with DCIS at biopsy, 25 of which were subsequently upstaged to invasive cancer. A deep convolutional neural network model that was pretrained on nonmedical images (eg, animals, plants, instruments) was used as the feature extractor. Through a statistical pooling strategy, deep features were extracted at different levels of convolutional layers from the lesion areas, without sacrificing the original resolution or distorting the underlying topology. A multivariate classifier was then trained to predict which tumors contain occult invasive disease. This was compared with the performance of traditional "handcrafted" computer vision (CV) features previously developed specifically to assess mammographic calcifications. The generalization performance was assessed using Monte Carlo cross-validation and receiver operating characteristic curve analysis.

RESULTS

Deep features were able to distinguish DCIS with occult invasion from pure DCIS, with an area under the receiver operating characteristic curve of 0.70 (95% confidence interval, 0.68-0.73). This performance was comparable with the handcrafted CV features (area under the curve = 0.68; 95% confidence interval, 0.66-0.71) that were designed with prior domain knowledge.

CONCLUSIONS

Despite being pretrained on only nonmedical images, the deep features extracted from digital mammograms demonstrated comparable performance with handcrafted CV features for the challenging task of predicting DCIS upstaging.

Abstract GS5-04: Prediction of occult invasive disease in ductal carcinoma in situ using deep learning features

Background: Deep learning, especially deep convolutional neural network (CNN), has emerged as a promising approach for many image recognition or classification tasks, demonstrating human or even superhuman performance. Used as feature extractor, some pre-trained CNN models can match or surpass the performance of domain-specific, “handcrafted” features. In this study, we aim to determine whether deep features extracted from digital mammograms using a pre-trained deep CNN are prognostic of occult invasive disease for patients with ductal carcinoma in situ (DCIS) on core needle biopsy.

Materials and Methods: In this retrospective study, we collected digital mammography magnification views for 99 subjects with DCIS at biopsy, 25 of which were subsequently upstaged to invasive cancer. We utilized a deep CNN model that was pre-trained on non-medical images (e.g., animals, plants, instruments) as the feature extractor. Through a statistical pooling strategy, we extracted deep features at different levels of convolutional layers from the lesion areas, without sacrificing the original resolution or distorting the underlying topology. A multivariate classifier was then trained to predict which tumors contain occult invasive disease. This was compared to the performance of traditional “handcrafted” computer vision (CV) features previously developed specifically to assess mammographic calcifications. The generalization performance was assessed using Monte Carlo cross validation and receiver operating characteristic (ROC) curve analysis.

Results: Deep features were able to distinguish DCIS with occult invasion from pure DCIS, with an area under the ROC curve (AUC-ROC) equal to 0.70 (95% CI: 0.68-0.73). This performance was comparable to the "handcrafted" CV features (AUC-ROC = 0.68, 95% CI: 0.66-0.71) that were designed with prior domain knowledge.

Conclusion: In spite of being pre-trained on only non-medical images, the deep features extracted from digital mammograms demonstrated comparable performance to "handcrafted" CV features for the challenging task of predicting DCIS upstaging.

Acknowledgments: This work was supported in part by NIH/NCI R01-CQA185138 and DOD Breast Cancer Research Program W81XWH-14-1-0473.

Citation Format: Lo JY, Grimm LJ, Mazurowski MA, Baker JA, Marks JR, King LM, Maley CC, Hwang E-SS, Shi B. Prediction of occult invasive disease in ductal carcinoma in situ using deep learning features [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr GS5-04.

Abstract P4-15-10: Second breast events after DCIS: Where, what and when?

Background: Detection and diagnosis of ductal carcinoma in situ has substantially increased since the widespread use of mammographic screening with the incidence of DCIS increasing between 1975 from 5.8 per 100,000 to 33.8 per 100,000 women in 2010. While surgical resection of DCIS remains the main treatment, outcomes in terms of in breast re-occurrence or contralateral disease are unpredictable. Recently, three large clinical trials of active surveillance versus standard of care have been initiated. We sought to compare the clinical, imaging, histopathological features and treatment of DCIS for contemporary patients who subsequently developed a second in situ or invasive breast lesion among women treated within a single institution.

Methods: A contemporary prospective cohort study comprising 2,509 women treated for DCIS between 2000 and 2014 was examined for diagnostic imaging, surgery, pathology, adjuvant treatment and second breast events. Patients with primary bilateral DCIS or a bilateral second event were excluded as were patients with antecedent invasive breast cancer. Patients were grouped into two populations based on ipsilateral and contralateral second breast event.

Results: Amongst 2,509 women treated for DCIS between 2000 and 2014, of whom 146 had a second breast event (5.8%). 77 (52.7% of events, 3.07% overall) developed an ipsilateral second breast event (37 DCIS, 40 invasive breast cancer) and 69 (47.3%, 2.75%) patients developed a contralateral second breast lesion (34 DCIS, 35 invasive breast cancer). Patients who developed a contralateral second event were older than those who developed an ipsilateral second event (58 years vs. 52 years, p=0.003). Patients receiving segmental mastectomies were more likely to develop an ipsilateral second event than a contralateral event (85.7% vs. 60.9%, p=0.001). Patients who developed a second contralateral breast event were significantly older at diagnosis (63 years vs. 58 years, p=0.003). There was no significant difference in the number of (69 vs. 77) or time between contralateral and ipsilateral second breast events (5.2 years and 5.1 years, p=0.8) nor in the number of patients with invasive or in situ second events (75 vs. 71). Patients with ipsilateral second breasts events were more likely to undergo bilateral mastectomy than those with contralateral second events (15 vs. 2, p=0.003).

Conclusions: Following standard of care treatment for DCIS, 1 in 20 women develop further DCIS or invasive disease at a median of 5 years with very similar proportions of ipsilateral or contralateral and DCIS or invasive disease. Current trials of active surveillance should consider a median 5 years follow up as a critical time point for reporting results.

Key Words: Ductal carcinoma in situ, invasive breast cancer, population-based cohort, surgery, radiotherapy, endocrine therapy.

Citation Format: King LM, Yi M, Krishnamurthy S, Hunt KK, Thompson AM. Second breast events after DCIS: Where, what and when? [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr P4-15-10.

Can Occult Invasive Disease in Ductal Carcinoma In Situ Be Predicted Using Computer-extracted Mammographic Features?

RATIONALE AND OBJECTIVES

This study aimed to determine whether mammographic features assessed by radiologists and using computer algorithms are prognostic of occult invasive disease for patients showing ductal carcinoma in situ (DCIS) only in core biopsy.

MATERIALS AND METHODS

In this retrospective study, we analyzed data from 99 subjects with DCIS (74 pure DCIS, 25 DCIS with occult invasion). We developed a computer-vision algorithm capable of extracting 113 features from magnification views in mammograms and combining these features to predict whether a DCIS case will be upstaged to invasive cancer at the time of definitive surgery. In comparison, we also built predictive models based on physician-interpreted features, which included histologic features extracted from biopsy reports and Breast Imaging Reporting and Data System-related mammographic features assessed by two radiologists. The generalization performance was assessed using leave-one-out cross validation with the receiver operating characteristic curve analysis.

RESULTS

Using the computer-extracted mammographic features, the multivariate classifier was able to distinguish DCIS with occult invasion from pure DCIS, with an area under the curve for receiver operating characteristic equal to 0.70 (95% confidence interval: 0.59-0.81). The physician-interpreted features including histologic features and Breast Imaging Reporting and Data System-related mammographic features assessed by two radiologists showed mixed results, and only one radiologist's subjective assessment was predictive, with an area under the curve for receiver operating characteristic equal to 0.68 (95% confidence interval: 0.57-0.81).

CONCLUSIONS

Predicting upstaging for DCIS based upon mammograms is challenging, and there exists significant interobserver variability among radiologists. However, the proposed computer-extracted mammographic features are promising for the prediction of occult invasion in DCIS.

Characterization and clinical validation of MCM2 and TOP2A monoclonal antibodies in the BD ProEx™ C assay: An immunoassay which detects aberrant S-phase induction in cervical tissue

BACKGROUND

The Papanicolaou (Pap) screen has been successful in reducing cervical cancer; but exhibits low sensitivity when detecting cervical dysplasia. Use of molecular biomarkers in Pap tests may improve diagnostic accuracy.

DESIGN

Monoclonal antibodies to Minichromosome Maintenance Protein 2 (MCM2) and DNA Topoisomerase II α (TOP2A) were selected for use in IHC based on their ability to differentiate normal from diseased cervical tissues in tissue microarrays. Enhanced Green Fluorescent Protein Western blot analysis was used to help identify binding epitopes specific to MCM2 and TOP2A antibody clones. Antibody affinity was determined by solution phase affinity measurement and immunohistochemistry was performed using high affinity MCM2 or TOP2A antibodies on serial histological sections.

RESULTS

Antibody clones to MCM2 and TOP2A clones were selected based on their ability to detect over expression in abnormal cervical epithelia. In IHC, MCM2-27C5.6 and MCM2-26H6.19 demonstrated superior staining in abnormal cervical tissue over the MCM2-CRCT2.1 antibody. A combination of MCM2 and TOP2A antibodies showed greater staining when compared to staining with any of the antibodies alone on serial histological sections. Distinct linear epitopes were elucidated for each of the MCM2 and TOP2A clones. Affinity values (Kd) for MCM2 or TOP2A antibodies had a similar range. In a research study, the MCM2 and TOP2A (BD ProEx™ C) antibody cocktail showed increased epithelia staining with increasing dysplasia. The use of BD ProEx™ C in combination with H&E staining enhanced immunohistochemical discrimination of dysplastic and non-dysplastic FFPE cervical tissue specimens.

CONCLUSIONS

BD ProEx™ C containing MCM2 and TOP2A antibodies showed strong specific nuclear staining that correlated with increased dysplasia and lesion severity. Enhanced performance of the antibodies was linked to their unique topography recognition. BD ProEx™ C incorporates antibodies that enhance detection of CIN2+ cervical disease.

Abstract P6-05-03: Genomic diversity of ductal carcinoma in situ (DCIS) as a driver of invasion and metastasis

Background: Recent evidence indicates that breast cancers may have high levels of heterogeneity. Based on principles of tumor evolution, we are investigating whether DCIS with higher levels of intra-tumoral genetic heterogeneity are more likely to progress to invasive and/or metastatic disease.

Methods: Cases (DCIS with co-existing invasive or metastatic cancer) and controls (pure DCIS) are identified from Duke Pathology archives. From cases and controls, we are analyzing two areas of DCIS separated by &gt;1cm with germ line DNA from the same subject to measure allele frequencies of somatic mutations and copy number variation (CNV). Small amounts of FFPE derived DNA are made into NGS libraries for full exome sequencing and hybridization to a 4.2 million element SNP array. Comparison of allele frequencies and CNV is made between specimens from the same cancer to assess heterogeneity.

Results: We identified a series of pure DCIS (controls for this study) and generated high coverage sequencing data from 20ng of FFPE DNA from 12 samples (4 subjects, germline + 2 DCIS containing areas from each subject) as proof of principle. We compared the occurrence of deletions, insertions and SNP's after a 20X coverage filtration. From these data, we identified between 50 to greater than 600 sequence changes in these DCIS compared to normal DNA. Present/absent calls and allele frequencies indicate both significant and variable degrees of heterogeneity even in these pure DCIS samples. Additional data is now being gathered and analyzed based on this established technical approach.

Conclusion: We have demonstrated the feasibility of acquiring high quality NGS data from archival DCIS specimens allowing us to test the hypothesis that genetic heterogeneity is a driving force in breast cancer progression. The degree and nature of this heterogeneity will be assessed in a panel of pure DCIS and DCIS co-existing with invasive and/or metastatic cancer. We are now generating and analyzing these data for Symposium presentation.

Citation Format: King LM, Marks JR, Hall AH, Temko D, Graham TA, Mardis ER, Maley CC, Hwang E. Genomic diversity of ductal carcinoma in situ (DCIS) as a driver of invasion and metastasis. [abstract]. In: Proceedings of the Thirty-Eighth Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2015 Dec 8-12; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2016;76(4 Suppl):Abstract nr P6-05-03.

Hands-free Head-movement Gesture Recognition using Artificial Neural Networks and the Magnified Gradient Function

This paper presents a hands-free head-movement gesture classification system using a Neural Network employing the Magnified Gradient Function (MGF) algorithm. The MGF increases the rate of convergence by magnifying the first order derivative of the activation function, whilst guaranteeing convergence. The MGF is tested on able-bodied and disabled users to measure its accuracy and performance. It is shown that for able-bodied users, a classification improvement from 98.25% to 99.85% is made, and 92.08% to 97.50% for disabled users.

Isolation of RNA from cell lines and cervical cytology specimens stored in BD SurePath preservative fluid and downstream detection of housekeeping gene and HPV E6 expression using real time RT-PCR

This study was performed to demonstrate that RNA isolated from cell lines and cervical cytology specimens stored in SurePath preservative fluid would be functional in real-time RT-PCR assays. RNA was isolated from cervical cell lines or cytology samples stored in SurePath preservative at room temperature for 2-5 weeks using five commercially available RNA purification kits, three of which contain proteinases. The quality of the RNA was assessed by real time RT-PCR amplification of GAPDH, GUSB, U1A, HPV 16 and 18 E6 mRNAs. RNA was isolated successfully from cells that were stored in SurePath preservative fluid with only the three protocols that contained proteinases. GAPDH was amplified in 98-100% of the samples, GUSB in 90-98%, and the least abundant transcript, U1A, was amplified in 81-96% of the samples. HPV 16 and 18 E6 transcripts were detected in 56% of high grade, 39% of low grade and 2% of normal samples, with a concordance between DNA genotype and E6 mRNA expression of 97%. We demonstrated that RNA can be extracted from cervical cell lines and cytology specimens stored in BD SurePath preservative fluid with three different procedures that all contain proteinases. This RNA is suitable for real-time RT-PCR applications.

Early driver fatigue detection from electroencephalography signals using artificial neural networks

This paper describes a driver fatigue detection system using an artificial neural network (ANN). Using electroencephalogram (EEG) data sampled from 20 professional truck drivers and 35 non professional drivers, the time domain data are processed into alpha, beta, delta and theta bands and then presented to the neural network to detect the onset of driver fatigue. The neural network uses a training optimization technique called the magnified gradient function (MGF). This technique reduces the time required for training by modifying the standard back propagation (SBP) algorithm. The MGF is shown to classify professional driver fatigue with 81.49% accuracy (80.53% sensitivity, 82.44% specificity) and non-professional driver fatigue with 83.06% accuracy (84.04% sensitivity and 82.08% specificity).

Real-time head movement system and embedded Linux implementation for the control of power wheelchairs

Mobility has become very important for our quality of life. A loss of mobility due to an injury is usually accompanied by a loss of self-confidence. For many individuals, independent mobility is an important aspect of self-esteem. Head movement is a natural form of pointing and can be used to directly replace the joystick whilst still allowing for similar control. Through the use of embedded LINUX and artificial intelligence, a hands-free head movement wheelchair controller has been designed and implemented successfully. This system provides for severely disabled users an effective power wheelchair control method with improved posture, ease of use and attractiveness.

Single nucleotide polymorphisms in the CYP2J2 and CYP2C8 genes and the risk of hypertension

CYP2J2 and CYP2C8 metabolize arachidonic acid (AA) to cis-epoxyeicosatrienoic acids (EETs), which play a central role in regulating renal tubular fluid-electrolyte transport and vascular tone. We hypothesized that functionally relevant polymorphisms in the CYP2J2 or CYP2C8 genes influence hypertension risk. We examined associations between CYP2J2*7 (G-50 T promoter) and CYP2C8*3 (Arg139Lys and Lys399Arg, which are in 100% linkage disequilibrium) polymorphisms and hypertension in a biethnic population from Tennessee. CYP2J2*7 variant allele frequency was significantly higher in African-Americans versus Caucasians (14.1% versus 7.7%, P=0.01), irrespective of hypertension status. When analysed separately by race, the genotype distribution of the CYP2J2*7 variant allele was not significantly different among African-Americans with/without hypertension, but was significantly different among Caucasians with/without hypertension (P=0.03). Indeed, the odds ratio of having hypertension attributable to carrying the CYP2J2*7 variant allele adjusted for age, gender, body mass index and family history was 0.39 (95% confidence interval 0.17-0.89) among Caucasians, suggesting a protective effect. Additional subgroup analyses revealed a significantly lower CYP2J2*7 variant allele frequency in hypertensive versus normotensive Caucasian males (5.6% versus 12.5%, P=0.02) and in hypertensive versus normotensive Caucasians without a family history of hypertension (1.5% versus 11.0%, P=0.03). With respect to the CYP2C8*3 variant, genotype distribution and allele frequencies were similar between normotensive and hypertensive subjects. This study provides evidence for an association between CYP2J2*7 genotype and hypertension in Caucasian males and Caucasians without a family history of hypertension, but suggests no association between CYP2C8*3 genotype and hypertension. Confirmation of these findings in additional populations is warranted.

Fine mapping of the X-linked split-hand/split-foot malformation (SHFM2) locus to a 5.1-Mb region on Xq26.3 and analysis of candidate genes

Split-hand/split-foot malformation (SHFM) is a genetically heterogeneous disorder, with five known loci, that causes a lack of median digital rays, syndactyly, and aplasia or hypoplasia of the phalanges, metacarpals, and metatarsals. In the only known SHFM2 family, affected males and homozygous females exhibit monodactyly or bidactyly of the hands and lobster-claw feet. This family (1) was revisited to include additional subjects and genealogical data. All 39 affected males and three females fully expressed the SHFM, while 13 carrier females examined exhibited partial expression of SHFM. We narrowed the previously linked 22-Mb genetic interval on Xq24-q26 (2), by analyzing additional family members and typing additional markers. The results define a 5.1-Mb region with a new centromeric boundary at DXS1114 and a telomeric boundary at DXS1192. We did not identify mutations in the exons and exon/intron boundaries of 19 candidate genes. These data suggest that the mutation may lie in a regulatory region of one of these candidate genes or in another gene within the SHFM2 region with unclear role in limb development.

A locus for spondylocarpotarsal synostosis syndrome at chromosome 3p14

Spondylocarpotarsal synostosis syndrome is a rare autosomal recessive disorder characterised by vertebral fusions, frequently manifesting as an unsegmented vertebral bar, as well as fusions of the carpal and tarsal bones. In a study of three consanguineous families and one non-consanguineous family, linkage analysis was used to establish the chromosomal location of the disease gene. Linkage analysis localised the disease gene to chromosome 3p14. A maximum lod score of 6.49 (q = 0) was obtained for the marker at locus D3S3532 on chromosome 3p. Recombination mapping narrowed the linked region to the 5.7 cM genetic interval between the markers at loci D3S3724 and D3S1300. A common region of homozygosity was found between the markers at loci D3S3724 and D3S1300, defining a physical interval of approximately 4 million base pairs likely to contain the disease gene. Identification of the gene responsible for this disorder will provide insight into the genes that play a role in the formation of the vertebral column and joints.

Field testing the influence of sperm competition based on sperm mobility in breeder turkey toms

1. Commercial reproduction of turkeys relies on pooling of semen from multiple males for inseminations. Understanding how sperm characteristics influence paternity under commercial breeding conditions is important to improving production efficiency. 2. The objective of this study was to evaluate progeny production of individual toms following commercial practices of pooling semen to determine if sperm mobility influences progeny production in field conditions. 3. A total of 104 toms were evaluated for sperm mobility. A subset of 10 toms were housed together and semen was collected, pooled and used to inseminate hens (n = 28). Hens were inseminated at 30 weeks of age and weekly thereafter. 4. Ejaculates from each tom were evaluated on two separate days for sperm mobility. Semen from each tom was diluted and layered upon 6% (wt/vol) Accudenz solution. The sperm suspension was incubated at 41 degrees C for 5 min and absorbance was measured with a spectrophotometer. 5. Toms were ranked by absorbance and categorised as high or low if mobility score was +/- 1 SD from the flock mean (average). 6. For parentage determination, DNA was extracted from tom, hen and poult blood. Poult parentage (n = 276) was determined at one day of age or at 14 weeks by analysis of marker genotypes that were generated by polymerase chain reaction (PCR) amplification of genomic DNA with selected microsatellite markers. 7. Sperm mobility differed across males with absorbance values ranging from 0.147 to 0.366. 8. Findings demonstrate differences in poult production among individual toms when semen from multiple males was pooled and inseminated. Toms classified as high, average and low produced 55, 41 and 4% of the offspring, respectively. 9. It appears that sperm mobility is a trait that influences sperm competition among toms under field conditions where sperm numbers inseminated from individual toms are not controlled or constant and that toms with low sperm mobility produce few offspring.

Cloning of CYP2J2 gene and identification of functional polymorphisms

CYP2J2 is abundant in cardiovascular tissue and active in the metabolism of arachidonic acid to eicosanoids that possess potent anti-inflammatory, vasodilatory, and fibrinolytic properties. We cloned and sequenced the entire CYP2J2 gene (approximately 40.3 kb), which contains nine exons and eight introns. We then sequenced the CYP2J2 exons and intron-exon boundaries in 72 healthy persons representing African, Asian, and European/white populations as part of the National Institutes of Health/National Institute of Environmental Health Sciences Environmental Genome Single Nucleotide Polymorphism Program. A variety of polymorphisms were found, four of which resulted in coding changes (Arg158Cys, Ile192Asn, Asp342Asn, and Asn404Tyr). A fifth variant (Thr143Ala) was identified by screening a human heart cDNA library. All five variant cDNAs of CYP2J2 were generated by site-directed mutagenesis and expressed in Sf9 insect cells by using a baculovirus system. The recombinant wild-type and variant CYP2J2 proteins immunoreacted with peptide-based antibodies to CYP2J2 and displayed typical cytochrome P450 (P450) CO-difference spectra; however, the Asn404Tyr and Ile192Asn variants also had prominent spectral peaks at 420 nm. The ability of these variants to metabolize arachidonic acid and linoleic acid was compared with that of wild-type CYP2J2. Three variants (Asn404Tyr, Arg158Cys, and Thr143Ala) showed significantly reduced metabolism of both arachidonic acid and linoleic acid. The Ile192Asn variant showed significantly reduced activity toward arachidonic acid only. The Asp342Asn variant showed similar metabolism to wild-type CYP2J2 for both endogenous substrates. Based on these data, we conclude that allelic variants of the human CYP2J2 gene exist and that some of these variants result in a P450 protein that has reduced catalytic function. Insofar as CYP2J2 products have effects in the cardiovascular system, we speculate that these variants may be functionally relevant.

Segregation of spermatozoa within sperm storage tubules of fowl and turkey hens

In avian species, spermatozoa reside in the oviduct for prolonged periods in specialized structures known as sperm storage tubules, but little is known about the relative distribution of spermatozoa in these tubules after successive inseminations by different males. The staining efficacies of various fluorescent dyes for fowl and turkey spermatozoa were evaluated to investigate one proposed mechanism of sperm competition. Hens were then inseminated at different intervals with stained and unstained spermatozoa to observe the spatial distribution of spermatozoa within the storage tubules. Several novel fluorescent lipophilic tracers that successfully stain mammalian spermatozoa either did not stain fowl or turkey spermatozoa, or greatly impaired sperm motility. In contrast, Hoechst 33342 readily stained sperm nuclei (fowl: 25 nmol l-1; turkey: 77 nmol l-1) within 4 h without inhibiting sperm motility, or affecting fertility or the hatching ability of the eggs. Hens were tandemly inseminated with equal numbers of stained or unstained spermatozoa at 24 h intervals and were killed 24 h after the final insemination to study sperm entry and storage within the tubules. Oviductal mucosa containing sperm storage tubules was removed, and individual tubules were classified as containing stained spermatozoa, unstained spermatozoa, a mixture of stained and unstained spermatozoa, or as not containing spermatozoa. Results from the present study indicate that spermatozoa from two different inseminations generally segregate into different storage tubules in both fowl and turkey hens. Storage tubules containing mixed populations of spermatozoa were found in only 4% of fowl and 12% of turkey storage tubules examined. Thus, the mechanism of last-male precedence does not appear to be due to the stratification of spermatozoa within the tubules.

Segregation of spermatozoa within sperm storage tubules of fowl and turkey hens

In avian species, spermatozoa reside in the oviduct for prolonged periods in specialized structures known as sperm storage tubules, but little is known about the relative distribution of spermatozoa in these tubules after successive inseminations by different males. The staining efficacies of various fluorescent dyes for fowl and turkey spermatozoa were evaluated to investigate one proposed mechanism of sperm competition. Hens were then inseminated at different intervals with stained and unstained spermatozoa to observe the spatial distribution of spermatozoa within the storage tubules. Several novel fluorescent lipophilic tracers that successfully stain mammalian spermatozoa either did not stain fowl or turkey spermatozoa, or greatly impaired sperm motility. In contrast, Hoechst 33342 readily stained sperm nuclei (fowl: 25 nmol l-1; turkey: 77 nmol l-1) within 4 h without inhibiting sperm motility, or affecting fertility or the hatching ability of the eggs. Hens were tandemly inseminated with equal numbers of stained or unstained spermatozoa at 24 h intervals and were killed 24 h after the final insemination to study sperm entry and storage within the tubules. Oviductal mucosa containing sperm storage tubules was removed, and individual tubules were classified as containing stained spermatozoa, unstained spermatozoa, a mixture of stained and unstained spermatozoa, or as not containing spermatozoa. Results from the present study indicate that spermatozoa from two different inseminations generally segregate into different storage tubules in both fowl and turkey hens. Storage tubules containing mixed populations of spermatozoa were found in only 4% of fowl and 12% of turkey storage tubules examined. Thus, the mechanism of last-male precedence does not appear to be due to the stratification of spermatozoa within the tubules.

Identification of human FEM1A, the ortholog of a C. elegans sex-differentiation gene

We report the isolation, genomic structure, chromosomal location, and expression pattern of the FEM1A gene, the human ortholog of the Caenorhabditis elegans fem-1 and mouse Fem1a genes. The coding sequence is 1851 bp and encodes a 617 amino acid protein. The human FEM1A protein has 65% identity with the mouse Fem1a protein and 34% identity with the C. elegans fem-1 protein, indicating conservation of this protein. The N-terminal region of the encoded protein contains six ankyrin repeat elements, a motif found in signaling and transcriptional regulatory molecules such as Notch and glp1. The gene was highly expressed in human kidney and cardiac tissue, and was expressed at lower levels in multiple tissues, including cartilage. FEM1A was localized to chromosome 5q23.1, a region of conserved synteny with a portion of mouse chromosome 17 that contains Fem1a. In C. elegans, fem-1 is involved in a pathway necessary for sex determination. The identification of a human homolog of this conserved gene suggests a potential role for this sex-determining molecule in humans.

Functional characterization of cyclooxygenase-2 polymorphisms

Cyclooxygenases (COX)-1 and -2 are the key enzymes in the conversion of arachidonic acid to prostaglandins. COX-2 appears to play an emerging role in inflammation and carcinogenesis. Nonsteroidal anti-inflammatory drugs (NSAIDs) are used for the treatment of numerous diseases and reduce the risk of developing colorectal cancer. Polymorphisms in the COX-2 gene could alter enzyme expression, function, and/or the response to NSAIDs. Therefore, they could modify individual risks for developing cancer and other diseases or the occurrence of side effects or sensitivity toward selective or nonselective COX inhibitors. We sequenced the COX-2 gene of 72 individuals and identified rare polymorphisms in the promoter and the coding region. A COX-2 molecular model was used to locate the coding region polymorphisms relative to functional sites in the protein, and the COX-2 V511A polymorphism was very near to the active site. This variant protein was expressed, and function was evaluated, but no difference was detected in metabolism of the COX-2 substrates, arachidonic acid, linoleic acid, and 2-arachidonyl glycerol, compared with the wild type. The Km values for arachidonic acid showed no differences between the COX-2 wild type and V511A mutant. Inhibition with selective or nonselective COX inhibitors was essentially the same for the two enzymes. The absence of functionally important polymorphisms in the COX-2 gene may suggest that there has been selective pressure against those single nucleotide polymorphisms because of the critical role of this enzyme in maintenance of homeostasis.

Free fatty acids, but not ketone bodies, protect diabetic rat hearts during low-flow ischemia

To determine whether the effects of fatty acids on the diabetic heart during ischemia involve altered glycolytic ATP and proton production, we measured energetics and intracellular pH (pH(i)) by using (31)P NMR spectroscopy plus [2-(3)H]glucose uptake in isolated rat hearts. Hearts from 7-wk streptozotocin diabetic and control rats, perfused with buffer containing 11 mM glucose, with or without 1.2 mM palmitate or the ketone bodies, 4 mM beta-hydroxybutyrate plus 1 mM acetoacetate, were subjected to 32 min of low-flow (0.3 ml x g wet wt(-1) x min(-1)) ischemia, followed by 32 min of reperfusion. In control rat hearts, neither palmitate nor ketone bodies altered the recovery of contractile function. Diabetic rat hearts perfused with glucose alone or with ketone bodies, had functional recoveries 50% lower than those of the control hearts, but palmitate restored recovery to control levels. In a parallel group with the functional recoveries, palmitate prevented the 54% faster loss of ATP in the diabetic, glucose-perfused rat hearts during ischemia, but had no effect on the rate of ATP depletion in control hearts. Palmitate decreased total glucose uptake in control rat hearts during low-flow ischemia, from 106 +/- 17 to 52 +/- 12 micromol/g wet wt, but did not alter the total glucose uptake in the diabetic rat hearts, which was 42 +/- 5 micromol/g wet wt. Recovery of contractile function was unrelated to pH(i) during ischemia; the glucose-perfused control and palmitate-perfused diabetic hearts had end-ischemic pH(i) values that were significantly different at 6.36 +/- 0.04 and 6.60 +/- 0.02, respectively, but had similar functional recoveries, whereas the glucose-perfused diabetic hearts had significantly lower functional recoveries, but their pH(i) was 6.49 +/- 0.04. We conclude that fatty acids, but not ketone bodies, protect the diabetic heart by decreasing ATP depletion, with neither having detrimental effects on the normal rat heart during low-flow ischemia.

Characterization of a human gene encoding nucleosomal binding protein NSBP1

We characterize the cDNA and genomic structure of NSBP1, and demonstrate that it is a nuclear protein and the homologue of mouse Nsbp1, which is known to encode a nucleosomal binding and transcriptional activating protein related to the HMG-14/-17 chromosomal proteins. The encoded NSBP1 protein has 86% amino acid similarity to Nsbp1, including identity in nucleosomal binding domains of the HMG-14/-17 proteins. Our radiation hybrid data localize NSBP1 and Nsbp1 to homologous regions of chromosome X, with NSBP1 in Xq13.3 between DXS983 and DXS995 and Nsbp1 in the interval DXMit65 and DXMit39. Although Nsbp1 produces one mRNA transcript, NSBP1 produces three transcripts with alternate polyadenylated sites. The 3' untranslated region (UTR) of NSPB1 mRNA also contains several AU-rich elements (AREs), which are associated with rapid mRNA turnover. Northern analysis of NSBP1/Nsbp1 shows differences in transcript abundance among adult and fetal tissues, with predominant expression in liver, kidney, trabecular bone, and bone marrow stromal cells. However, a reverse transcriptase-PCR analysis shows nearly ubiquitous expression of the three NSBP1 transcripts in all tissues examined, although the abundance of each transcript was not quantified. NSBP1 is encoded by six exons and has exon-intron boundaries identical to the HMG-14/-17 genes. The last exon and the 3' UTR of NSBP1 contain retrotransposon sequences of HAL1, HERV-H, and L1MB7, suggesting that these retrotransposons were involved in the origin of NSPB1 from an ancestral-like HMG-14/-17 gene. The similarities among NSBP1, Nsbp1, and the HMG-14/-17 proteins suggest that NSBP1 may function as a nucleosomal binding and transcriptional activating element. Further, the AREs in the 3' UTR of NSPB1 suggest that alternate poly(A) site selection may mediate the mRNA stability of this gene.

Efficacy of sperm mobility assessment in commercial flocks and the relationships of sperm mobility and insemination dose with fertility in turkeys

Our objectives were to evaluate: 1) the efficacy of the Sperm Mobility Test on commercial turkey farms, and 2) the influence of sperm mobility phenotype on fertility when insemination parameters are varied. In research flocks, differences in sperm mobility among toms are predictive of fertility. We wanted to test the efficacy of this sire selection test in practical, real-world situations, evaluating its usefulness in terms of assessing large numbers of toms, different strains of turkeys, and variable management practices. Utilizing field study results, controlled studies were then conducted to improve test parameters. For the field trials, semen from each of 405 breeder toms (11 strains or lines) was evaluated either in duplicate (n = 285) or in triplicate (n = 120). Sperm mobility was normally distributed among all toms tested, except for one strain. Because the sperm mobility indices for toms evaluated in these field trials were higher than those observed in research flocks, the Sperm Mobility Test was modified to increase the separation between high and low sperm mobility phenotypes by increasing the concentration of Accudenz. To determine the effects of sperm mobility and insemination dose on sustained fertility through time, hens from a research flock were inseminated twice before the onset of lay with sperm from toms classified as high-, average-, or low-mobility in concentrations of 25 to 400 million sperm per artificial insemination dose, and egg fertility was evaluated over a 5-wk period. Toms with the high-mobility sperm phenotype maintained higher fertility (P < 0.05) over the 5-wk period at all insemination doses compared with toms with low-mobility sperm. Toms with high-mobility sperm sired equal numbers of poults in a sperm competition study in which numbers favored low-mobility toms by 3:1. These results demonstrate that the Sperm Mobility Test can be used for on-farm evaluation of semen quality of toms in commercial flocks and that sperm mobility influences fertility and sire fitness.

High-level periplasmic expression in Escherichia coli using a eukaryotic signal peptide: importance of codon usage at the 5' end of the coding sequence

We investigated the ability of signal peptides of eukaryotic origin (human, mouse, and yeast) to efficiently direct model proteins to the Escherichia coli periplasm. These were compared against a well-characterized prokaryotic signal peptide-OmpA. Surprisingly, eukaryotic signal peptides can work very efficiently in E. coli, but require optimization of codon usage by codon-based mutagenesis of the signal peptide coding region. Analysis of the 5' of periplasmic and cytoplasmic E. coli genes shows some codon usage differences.

Differential zinc transport into testis and brain of cadmium-sensitive and -resistant murine strains

Recently, we showed that murine strain differences to the testicular toxicity of cadmium (Cd) are the result of variable transport of Cd across the blood-testis barrier. Because Cd is a nonessential trace element, it must be using the transporter for an endogenous substance. The objectives for this study were to determine the natural ligand for the transport system used by Cd to enter testis and brain, and to determine whether the transport of that natural ligand also differs among Cd-sensitive and -resistant murine strains. Because zinc (Zn) and Cd are cations of similar size and charge, and because Cd has been shown to inhibit Zn uptake in a variety of systems, we hypothesized that Cd was using Zn transporters to enter tissues. In this study we characterized Zn transport into the testis and brain of Cd-sensitive and -resistant murine strains. We found that the transport of 65Zn into testis and brain of Cd-resistant A/J mice was significantly reduced compared with that in Cd-sensitive 129/J mice. In 129/J mice, unlabeled CdCl2 significantly reduced 65Zn transport by 56% in testes and by 47% in brain. Pretreatment with Zn had no significant effect on 109Cd transport rates into testes or brain of 129/J or A/J mice, but did reduce the percentage of the injected 109Cd dose in testes of 129/J mice by 44% within 60 minutes. From these results we can conclude that Cd is using transport systems that normally function to regulate Zn levels in testes and brain. Murine strain resistance to the testicular effects of Cd is associated with a concomitant attenuation of the Zn transport system in testis.

Report of five novel and one recurrent COL2A1 mutations with analysis of genotype-phenotype correlation in patients with a lethal type II collagen disorder

Achondrogenesis II-hypochondrogenesis and severe spondyloepiphyseal dysplasia congenita (SEDC) are lethal forms of dwarfism caused by dominant mutations in the type II collagen gene (COL2A1). To identify the underlying defect in seven cases with this group of conditions, we used the combined strategy of cartilage protein analysis and COL2A1 mutation analysis. Overmodified type II collagen and the presence of type I collagen was found in the cartilage matrix of all seven cases. Five patients were heterozygous for a nucleotide change that predicted a glycine substitution in the triple helical domain (G313S, G517V, G571A, G910C, G943S). In all five cases, analysis of cartilage type II collagen suggested incorporation of the abnormal alpha1(II) chain in the extracellular collagen trimers. The G943S mutation has been reported previously in another unrelated patient with a strikingly similar phenotype, illustrating the possible specific effect of the mutation. The radiographically less severely affected patient was heterozygous for a 4 bp deletion in the splice donor site of intron 35, likely to result in aberrant splicing. One case was shown to be heterozygous for a single nucleotide change predicted to result in a T1191N substitution in the carboxy-propeptide of the proalpha1(II) collagen chain. Study of the clinical, radiographic, and morphological features of the seven cases supports evidence for a phenotypic continuum between achondrogenesis II-hypochondrogenesis and lethal SEDC and suggests a relationship between the amount of type I collagen in the cartilage and the severity of the phenotype.

Improved efficiency of site-specific copper(II) ion-catalysed protein cleavage effected by mutagenesis of cleavage site

The peptide sequence (N)DKTH(C) was previously investigated as a site for efficient, specific cleavage of a fusion protein by cupric ions using a humanized gamma1 Fab' as a model protein. Here we show that conservative mutations to three of the residues in the introduced cleavage site resulted in cleavage sites that were significantly improved. They were cleaved more efficiently by Cu(2+), such that cleavage reactions could be shorter, of lower pH or at a lower temperature. Some were even found to be measurably cleaved by Ni(2+). Use of these new cleavage sequences along with cupric ions may provide a more rapid and less harsh method for cost-effective, large-scale proteolytic cleavage of fusion proteins and peptides.

Correlation of CASA velocity and linearity parameters with sperm mobility phenotype in turkeys

Since all domestic turkeys are produced through artificial insemination, a measurable sperm characteristic that would be predictive of fertility would allow for the culling of poor males, resulting in improved reproductive efficiency. The sperm mobility test (SMT), which quantifies sperm penetration into an Accudenz solution, has been shown to correlate highly with fertilization potential of individual turkeys. Since this sire-selection test is based on the differences in sperm mobility between whole ejaculates from individual males, the objective of this study was to determine whether specific sperm velocity parameters would correlate with the SMT and to determine whether these characteristics could account for phenotypic differences in sperm mobility observed between males. The SMT was used to rank males within a flock (n = 110) in triplicate and to classify them into high, average, and low sperm mobility phenotypes on the basis of the sperm mobility index. Several sperm velocity parameters were evaluated for each male by a computer-aided sperm analysis (CASA) system, the Hobson Sperm Tracker. The types of measurements taken of 200 sperm tracks/ejaculate included the following: curvilinear velocity (VCL), average path velocity (VAP), straight-line velocity (VSL), linearity (LIN), beat-cross frequency (BCF), and mean angular displacement (MAD). Significant positive correlations were found between VSL, LIN, BCF, and sperm mobility, and a significant negative correlation was seen between MAD and sperm mobility. Subpopulations of sperm that had penetrated the Accudenz solution were isolated from each mobility phenotype and were analyzed by CASA, and significant correlations were again observed between VSL, LIN, BCF, and sperm mobility. We conclude that sperm velocity and linearity contribute to overall sperm mobility phenotype and are important characteristics of turkey sperm function. Key words: Motility, computer, spermatozoa.

Conservation of the Caenorhabditis elegans timing gene clk-1 from yeast to human: a gene required for ubiquinone biosynthesis with potential implications for aging

Mutations in the Caenorhabditis elegans gene clk-1 have a major effect on slowing development and increasing life span. The Saccharomyces cerevisiae homolog COQ7 encodes a mitochondrial protein involved in ubiquinone biosynthesis and, hence, is required for respiration and gluconeogenesis. In this study, RT-PCR and 5' RACE were used to isolate both human and mouse clk-1/COQ7 homologs. Human CLK-1 was mapped to Chr 16(p12-13.1) by Radiation Hybrid (RH) and fluorescence in situ hybridization (FISH) methods. The number and location of human CLK1 introns were determined, and the location of introns II and IV are the same as in C. elegans. Northern blot analysis showed that three different isoforms of CLK-1 mRNA are present in several tissues and that the isoforms differ in the amount of expression. The functional equivalence of human CLK-1 to the yeast COQ7 homolog was tested by introducing either a single or multicopy plasmid containing human CLK-1 cDNA into yeast coq7 deletion strains and assaying for growth on a nonfermentable carbon source. The human CLK-1 gene was able to functionally complement yeast coq7 deletion mutants. The protein similarities and the conservation of function of the CLK-1/clk-1/COQ7 gene products suggest a potential link between the production of ubiquinone and aging.

Turkey sperm mobility influences paternity in the context of competitive fertilization

We have devised a novel means of investigating competitive fertilization in turkeys, using microsatellite genotyping to identify male parentage. Our results demonstrate that sperm mobility is a mechanism responsible in part for paternity efficiency in turkeys. Sperm mobility is composed of several parameters in which sperm motility is a component. Differences between ejaculates in the number of sperm penetrating into a dense, insert, nontoxic solution were measured and used to classify males into high, average, or low sperm mobility phenotypes. Microsatellite genotyping was used to determine parentage of poults after equal numbers of sperm from 10 males (either high or average phenotype, n = 5, mixed with low phenotype, n = 5) were inseminated simultaneously. In a separate study, the numbers of sperm hydrolyzing the perivitelline layer of eggs were compared between hens inseminated with sperm from high-, average-, or low-phenotype males. Overall, heterospermic inseminations resulted in consistently fewer offspring produced by low-mobility phenotype males. This correlated with physiological data in which semen from the low-mobility males had reduced numbers of sperm at the fertilization site as determined by sperm hole counts in the perivitelline layer of eggs. This is the first illustration of a measurable sperm trait predictive of paternity success in a competitive fertilization trial in turkeys, a species that is predominately reproduced by artificial insemination of multiple-sire pools.

Isolation of sperm storage tubules from the uterovaginal junction mucosa of the turkey

This study was performed to determine whether intact sperm storage tubules (SST) could be successfully isolated from the uterovaginal junction (UVJ) mucosa of the turkey. Large White BUTA hens were inseminated and euthanatized 24 to 48 h later. Oviducts were excised, UVJ tissue removed, and SST were procured by enzymatic digestion. Recovered SST were intact and contained motile sperm. The sperm were oriented with their acrosomes pointed towards the distal end of the SST, and their long axes in parallel with the long axis of the tubule's lumen. This method for the isolation of intact SST can be readily applied for in vitro culture studies as well as for the extraction of DNA and RNA from the SST epithelium.