Common mitochondrial deletions in RNA-Seq: evaluation of bulk, single-cell, and spatial transcriptomic datasets

Common mitochondrial DNA (mtDNA) deletions are large structural variants in the mitochondrial genome that accumulate in metabolically active tissues with age and have been investigated in various diseases. We applied the Splice-Break2 pipeline (designed for high-throughput quantification of mtDNA deletions) to human RNA-Seq datasets and describe the methodological considerations for evaluating common deletions in bulk, single-cell, and spatial transcriptomics datasets. A robust evaluation of 1570 samples from 14 RNA-Seq studies showed: (i) the abundance of some common deletions detected in PCR-amplified mtDNA correlates with levels observed in RNA-Seq data; (ii) RNA-Seq library preparation method has a strong effect on deletion detection; (iii) deletions had a significant, positive correlation with age in brain and muscle; (iv) deletions were enriched in cortical grey matter, specifically in layers 3 and 5; and (v) brain regions with dopaminergic neurons (i.e., substantia nigra, ventral tegmental area, and caudate nucleus) had remarkable enrichment of common mtDNA deletions.

Construction of high coverage whole-genome sequencing libraries from single colon crypts without DNA extraction or whole-genome amplification

OBJECTIVE

Comprehensive and reliable genome-wide variant analysis of a small number of cells has been challenging due to genome coverage bias, PCR over-cycling, and the requirement of expensive technologies. To comprehensively identify genome alterations in single colon crypts that reflect genome heterogeneity of stem cells, we developed a method to construct whole-genome sequencing libraries from single colon crypts without DNA extraction, whole-genome amplification, or increased PCR enrichment cycles.

RESULTS

We present post-alignment statistics of 81 single-crypts (each contains four- to eight-fold less DNA than the requirement of conventional methods) and 16 bulk-tissue libraries to demonstrate the consistent success in obtaining reliable coverage, both in depth (≥ 30X) and breadth (≥ 92% of the genome covered at ≥ 10X depth), of the human genome. These single-crypt libraries are of comparable quality as libraries generated with the conventional method using high quality and quantities of purified DNA. Conceivably, our method can be applied to small biopsy samples from many tissues and can be combined with single cell targeted sequencing to comprehensively profile cancer genomes and their evolution. The broad potential application of this method offers expanded possibilities in cost-effectively examining genome heterogeneity in small numbers of cells at high resolution.

Impact of Genomic and Clinical Factors on Outcome of Children ≥18 Months of Age with Stage 3 Neuroblastoma with Unfavorable Histology and without MYCN Amplification: A Children's Oncology Group (COG) Report

PURPOSE

Patients ≥18 months of age with International Neuroblastoma Staging System (INSS) stage 3 unfavorable histology (UH), MYCN-nonamplified (MYCN-NA) tumors have favorable survival rates compared with other high-risk neuroblastoma populations. The impact of select clinical and biological factors on overall survival (OS) and event-free survival (EFS) were evaluated.

EXPERIMENTAL DESIGN

Patients enrolled on Children's Oncology Group (COG) A3973 (n = 34), ANBL0532 (n = 27), and/or biology protocol ANBL00B1 (n = 72) were analyzed. Tumors with available DNA (n = 65) and RNA (n = 42) were subjected to whole-exome sequencing (WES) and RNA sequencing. WES analyses and gene expression profiling were evaluated for their impact on survival. Multivariate analyses of EFS/OS using significant factors from univariate analyses were performed.

RESULTS

5-year EFS/OS for patients treated with high-risk therapy on A3973 and ANBL0532 were 73.0% ± 8.1%/87.9% ± 5.9% and 61.4% ± 10.2%/73.0% ± 9.2%, respectively (P = 0.1286 and P = 0.2180). In the A3973/ANBL0532 cohort, patients with less than partial response (PR; n = 5) at end-induction had poor outcomes (5-year EFS/OS: 0%/20.0% ± 17.9%. Univariate analyses of WES data revealed that subjects whose tumors had chromosome 1p or 11q loss/LOH and chromosome 5 or 9 segmental chromosomal aberrations had inferior EFS compared with those with tumors without these aberrations. Multivariate analysis revealed that 11q loss/LOH was an independent predictor of inferior OS [HR, 3.116 (95% confidence interval, 1.034-9.389), P = 0.0435].

CONCLUSIONS

Patients ≥18 months of age at diagnosis who had tumors with UH and MYCN-NA INSS stage 3 neuroblastoma assigned to high-risk therapy had an 81.6% ± 5.3% 5-year OS. Less than PR to induction therapy and chromosome 11q loss/LOH are independent predictors of inferior outcome and identify patients who should be eligible for future high-risk clinical trials.

Evolutionary conservation of the fidelity of transcription

Accurate transcription is required for the faithful expression of genetic information. However, relatively little is known about the molecular mechanisms that control the fidelity of transcription, or the conservation of these mechanisms across the tree of life. To address these issues, we measured the error rate of transcription in five organisms of increasing complexity and found that the error rate of RNA polymerase II ranges from 2.9 × 10-6 ± 1.9 × 10-7/bp in yeast to 4.0 × 10-6 ± 5.2 × 10-7/bp in worms, 5.69 × 10-6 ± 8.2 × 10-7/bp in flies, 4.9 × 10-6 ± 3.6 × 10-7/bp in mouse cells and 4.7 × 10-6 ± 9.9 × 10-8/bp in human cells. These error rates were modified by various factors including aging, mutagen treatment and gene modifications. For example, the deletion or modification of several related genes increased the error rate substantially in both yeast and human cells. This research highlights the evolutionary conservation of factors that control the fidelity of transcription. Additionally, these experiments provide a reasonable estimate of the error rate of transcription in human cells and identify disease alleles in a subunit of RNA polymerase II that display error-prone transcription. Finally, we provide evidence suggesting that the error rate and spectrum of transcription co-evolved with our genetic code.

Analysis of the genomic landscapes of Barbadian and Nigerian women with triple negative breast cancer

PURPOSE

Triple negative breast cancer (TNBC) is an aggressive breast cancer subtype that disproportionately affects women of African ancestry (WAA) and is often associated with poor survival. Although there is a high prevalence of TNBC across West Africa and in women of the African diaspora, there has been no comprehensive genomics study to investigate the mutational profile of ancestrally related women across the Caribbean and West Africa.

METHODS

This multisite cross-sectional study used 31 formalin-fixed paraffin-embedded (FFPE) samples from Barbadian and Nigerian TNBC participants. High-resolution whole exome sequencing (WES) was performed on the Barbadian and Nigerian TNBC samples to identify their mutational profiles and comparisons were made to African American, European American and Asian American sequencing data obtained from The Cancer Genome Atlas (TCGA). Whole exome sequencing was conducted on tumors with an average of 382 × coverage and 4335 × coverage for pooled germline non-tumor samples.

RESULTS

Variants detected at high frequency in our WAA cohorts were found in the following genes NBPF12, PLIN4, TP53 and BRCA1. In the TCGA TNBC cases, these genes had a lower mutation rate, except for TP53 (32% in our cohort; 63% in TCGA-African American; 67% in TCGA-European American; 63% in TCGA-Asian). For all altered genes, there were no differences in frequency of mutations between WAA TNBC groups including the TCGA-African American cohort. For copy number variants, high frequency alterations were observed in PIK3CA, TP53, FGFR2 and HIF1AN genes.

CONCLUSION

This study provides novel insights into the underlying genomic alterations in WAA TNBC samples and shines light on the importance of inclusion of under-represented populations in cancer genomics and biomarker studies.

High prevalence of deleterious mutations in concomitant nonsyndromic cleft and outflow tract heart defects

Our previous work demonstrating enrichment of outflow tract (OFT) congenital heart disease (CHD) in children with cleft lip and/or palate (CL/P) suggests derangements in common underlying developmental pathways. The current pilot study examines the underlying genetics of concomitant nonsyndromic CL/P and OFT CHD phenotype. Of 575 patients who underwent CL/P surgery at Children's Hospital Los Angeles, seven with OFT CHD, negative chromosomal microarray analysis, and no recognizable syndromic association were recruited with their parents (as available). Whole genome sequencing of blood samples paired with whole‐blood‐based RNA sequencing for probands was performed. A pathogenic or potentially pathogenic variant was identified in 6/7 (85.7%) probands. A total of seven candidate genes were mutated (CHD7, SMARCA4, MED12, APOB, RNF213, SETX, and JAG1). Gene ontology analysis of variants predicted involvement in binding (100%), regulation of transcription (42.9%), and helicase activity (42.9%). Four patients (57.1%) expressed gene variants (CHD7, SMARCA4, MED12, and RNF213) previously involved in the Wnt signaling pathway. Our pilot analysis of a small cohort of patients with combined CL/P and OFT CHD phenotype suggests a potentially significant prevalence of deleterious mutations. In our cohort, an overrepresentation of mutations in molecules associated with Wnt‐signaling was found. These variants may represent an expanded phenotypic heterogeneity within known monogenic disease genes or provide novel evidence of shared developmental pathways. The mechanistic implications of these mutations and subsequent developmental derangements resulting in the CL/P and OFT CHD phenotype require further analysis in a larger cohort of patients.

Abstract PO-142: Analysis of the genomic landscapes of Barbadian and Nigerian women with triple negative breast cancer

Women of African ancestry (WAA) are disproportionately affected by the aggressive triple negative breast cancer (TNBC) subtype that is often associated with high recurrence rates and metastasis. Although there is a high prevalence of TNBC across West Africa and in women of the African diaspora, there has been no comprehensive genomics study to investigate the mutational profile of ancestrally related women across the Caribbean and West Africa. To shed more light on this phenomenon, whole exome sequencing (WES) was performed on 31 formalin-fixed paraffin-embedded TNBC tissues from ancestrally related Barbadian and Nigerian women. We compared these genomics profiles with data from The Cancer Genome Atlas (TCGA) for African American (TCGA-AA), European American (TCGA-EA) women with TNBC. With an average coverage of 382x for tumour samples (n= 31) and 4335x for pooled germline (n=22) non-tumor samples, the most mutated genes in our cohorts include NBPF12, PLIN4, TP53 and BRCA1. For TCGA TNBC cases, these genes had a lower mutation rate, except for TP53 (32% in our cohort; 63% in TCGA-AA; 67% in TCGA-EA). For all altered genes, there were no differences in frequency of mutations between WAA TNBC groups including the TCGA-AA cohort. Additionally, we observed a high frequency of copy number variant alterations in PIK3CA, TP53, FGFR2 and HIF1AN genes. This study provides in-depth insights into the underlying genomic alterations in WAA-TNBC samples and shines light on the importance of inclusion of non-European populations in cancer genomics and biomarker studies.

Citation Format: Shawn M. Hercules, Xiyu Liu, Blessing I. Bassey-Archibong, Desiree H.A. Skeete, Suzanne Smith Connell, Adetola Daramola, Adekunbiola A.F. Banjo, Godwin Ebughe, Thomas Agan, Ima-Obong Ekanem, Joe E. Udosen, Christopher Obiorah, Aaron C. Ojule, Michael A. Misauno, Ayuba M. Dauda, Ejike C. Egbujo, Jevon C. Hercules, Amna Ansari, Ian Brain, Christine MacColl, Yili Xu, Yuxin Jin, Sharon Chang, John D. Carpten, André Bédard, Gregory R. Pond, Kim R.M. Blenman, Zarko Manojlovic, Juliet M. Daniel. Analysis of the genomic landscapes of Barbadian and Nigerian women with triple negative breast cancer [abstract]. In: Proceedings of the AACR Virtual Conference: 14th AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2021 Oct 6-8. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2022;31(1 Suppl):Abstract nr PO-142.

Wild-Type KRAS Allele Effects on Druggable Targets in KRAS Mutant Lung Adenocarcinomas

KRAS mutations are one of the most common oncogenic drivers in non-small cell lung cancer (NSCLC) and in lung adenocarcinomas in particular. Development of therapeutics targeting KRAS has been incredibly challenging, prompting indirect inhibition of downstream targets such as MEK and ERK. Such inhibitors, unfortunately, come with limited clinical efficacy, and therefore the demand for developing novel therapeutic strategies remains an urgent need for these patients. Exploring the influence of wild-type (WT) KRAS on druggable targets can uncover new vulnerabilities for the treatment of KRAS mutant lung adenocarcinomas. Using commercially available KRAS mutant lung adenocarcinoma cell lines, we explored the influence of WT KRAS on signaling networks and druggable targets. Expression and/or activation of 183 signaling proteins, most of which are targets of FDA-approved drugs, were captured by reverse-phase protein microarray (RPPA). Selected findings were validated on a cohort of 23 surgical biospecimens using the RPPA. Kinase-driven signatures associated with the presence of the KRAS WT allele were detected along the MAPK and AKT/mTOR signaling pathway and alterations of cell cycle regulators. FoxM1 emerged as a potential vulnerability of tumors retaining the KRAS WT allele both in cell lines and in the clinical samples. Our findings suggest that loss of WT KRAS impacts on signaling events and druggable targets in KRAS mutant lung adenocarcinomas.

Role of Ki-67, MRE11, and PD-L1 as Predictive Biomarkers for Recurrence Pattern in Muscle-invasive Bladder Cancer

BACKGROUND/AIM

Muscle invasive bladder cancer (MIBC) is an aggressive disease with high rates of local recurrence following radical cystectomy (RC). Currently, there are no clinically validated biomarkers to predict local only recurrence (LOR) and guide adjuvant treatment decisions. This pilot study evaluated the role of Ki-67, MRE11 and PD-L1 as predictive biomarkers for recurrence patterns in patients undergoing RC for MIBC.

PATIENTS AND METHODS

Our institutional cystectomy database containing cases from 1992-2014 was queried for patients with local only recurrence (LOR), and case-matched to patients with distant recurrence (DR) and no recurrence (NR). Clinicopathological data were collected and a tissue microarray was analyzed for presence of Ki-67, MRE11, and PD-L1 using immunofluorescence and immunohistochemistry.

RESULTS

Pathologic specimens from 42 patients (18 NR, 16 LOR, and 8 DR) were reviewed. Compared to normal bladder tissue, tumors had increased expression of Ki-67 (p<0.01) and PD-L1 (p<0.05). High Ki-67 was associated with recurrence pattern (local vs. distant) on univariate analysis (p<0.05). Ki-67 cell density varied by recurrence type: LOR (1354 cells/mm²), DR (557 cells/mm²) and NR (1111 cells/mm²) (p=0.034).

CONCLUSION

Our selected biomarkers could distinguish MIBC from normal bladder tissue but could not classify samples by recurrence pattern.

Abstract 862: Transcriptomics analysis reveals unique immunoinflammatory signatures of nigerian prostate cancer

Background: Prostate cancer (PCa) is the leading cause of cancer death among Nigerian men, with most men presenting late and with high grade disease. Furthermore, improved PCa diagnosis and treatment could rely on genomics-driven precision approaches. Mounting evidence supports potential population enriched molecular features of cancer, however the extent of which is hampered due to a dearth of genomic data in diverse cohorts. This study, therefore, aims at utilizing transcriptomics analysis to deduce transcriptional states in PCa from Nigerian men.

Methodology: RNA was extracted from macro-dissected prostate tumor and adjacent normal tissue from 35 archival formalin fixed paraffin-embedded (FFPE) blocks from Gleason 8-9 cancers, of which 25 had cribriform morphologic architecture. High resolution hybrid-capture mRNA-seq libraries were sequenced using Illumina NovaSeq 6000 system and generated an average of 124 million read pairs per sample. Raw reads were processed using a custom pipeline based upon industry standard algorithms along with the deduction of genetic ancestry. Molecular annotation was also performed to investigate tumor and immune cell heterogeneity.

Result: The Nigerian data were compared with high grade PCa data of Caucasian obtained from TCGA. The Nigerian PCa tumors where characterized with uniquely enriched immune-inflammatory microenvironment which include crosstalk between dendritic cells and natural killer cells, and natural killer cell signaling. While the Nigerian tumors share more features of distal metastatic tumor, their enrichment in NKT Cells is a unique feature.

Conclusion: This represents the first report of RNA-seq analysis of PCa in Sub-Saharan Africa and the unique transcriptomics features of Nigerian PCa presents opportunities for precision intervention in addressing its burden in Nigeria.

Citation Format: Solomon Oladapo Rotimi, Abidemi E. Omonisi, Faruk Mohammed, Paul Jibrin, Babatunde Duduyemi, Daniel B. Olusina, Yili Xu, Yuxin Jin, Xiyu Liu, Guang Q. Xiao, Folakemi Odedina, John Carpten, Renee Reams, Zarko Manojlovic. Transcriptomics analysis reveals unique immunoinflammatory signatures of nigerian prostate cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 862.

Imaging of type I procollagen biosynthesis in cells reveals biogenesis in highly organized bodies; Collagenosomes

Mechanistic aspects of type I procollagen biosynthesis in cells are poorly understood. To provide more insight into this process we designed a system to directly image type I procollagen biogenesis by co-expression of fluorescently labeled full size procollagen α1(I) and one α2(I) polypeptides. High resolution images show that collagen α1(I) and α2(I) polypeptides are produced in coordination in discrete structures on the ER membrane, which we termed the collagenosomes. Collagenosomes are disk shaped bodies, 0.5-1 μM in diameter and 200-400 nm thick, in the core of which folding of procollagen takes place. Collagenosomes are intimately associated with the ER membrane and their formation requires intact translational machinery, suggesting that they are the sites of nascent procollagen biogenesis. Collagenosomes show little co-localization with the COPII transport vesicles, which export type I procollagen from the ER, suggesting that these two structures are distinct. LARP6 is the protein which regulates translation of type I collagen mRNAs. The characteristic organization of collagenosomes depends on binding of LARP6 to collagen mRNAs. Without LARP6 regulation, collagenosomes are poorly organized and the folding of α1(I) and α2(I) polypeptides into procollagen in their cores is diminished. This indicates that formation of collagenosomes is dependent on regulated translation of collagen mRNAs. In live cells the size, number and shape of collagenosomes show little change within several hours, suggesting that they are stable structures of type I procollagen biogenesis. This is the first report of structural organization of type I collagen biogenesis in collagenosomes, while the fluorescent reporter system based on simultaneous imaging of both type I collagen polypeptides will enable the detailed elucidation of their structure and function.

Abstract SS1-07: Gene expression profiles of Ghanaian and Ethiopian triple-negative breast tumors

Background: Triple-negative breast cancer (TNBC) remains the most aggressive molecular subtype of breast cancer, with worse survival outcomes compared to other breast cancer subtypes. TNBC prevalence is highest among women of African descent worldwide, and through our previous work we have established a connection between West African ancestry (WAa) and higher rates of TNBC. Specifically, we have shown that prevalence rates of TNBC among West African and African American women are similar and higher than that of East African and White American women. We have also shown that quantified African ancestry is higher among TNBC cases compared to non-TNBC cases. To determine the influence African ancestry on the TNBC tumor environment, we sought to determine any differences in gene expression profiles of Ghanaian (West African) compared to Ethiopian (East African) women. Methods: RNA was extracted and sequenced from a pilot cohort of archival FFPE tumor tissue among Ghanaian (n = 19) and Ethiopian (n = 20) women. RNAseq reads were aligned, and quality of alignments were assessed, where de-duplicated samples with counts above 10M reads were included in the final analysis. Genetic ancestry was quantified by obtaining SNVs called from the RNAseq alignments, using GATK best practices. Differentially expressed genes lists were determined comparing Ghanaian vs. Ethiopian TNBC tumors, and also by identifying genes that were associated with increasing African ancestry. These gene lists, and log-fold change between comparison groups, were used as input for Ingenuity Pathway Analysis (IPA), to identify canonical pathways and de novo networks that are specific to Ghanaian or Ethiopian TNBC tumors. Results: Using 1KG populations as our reference to quantify genetic ancestry, we show that Ghanaian samples have &gt;94% AFR ancestry, specifically matching population groups representative of WAa. The Ethiopian samples showed between 37-48% AFR ancestry, primarily represented by East African groups. Interestingly, there seems to be a significant proportion of EUR ancestry among the Ethiopians samples (30-49%), primarily represented by Italian ancestry. We have conducted the differential gene expression analysis in two ways. First, we have compared gene expression profiles between Ghanaian and Ethiopian tumors. In our preliminary analysis, we identified &gt;600 genes (p &lt; 0.01) that were differentially expressed between Ghanaian and Ethiopian TNBC tumors. Second, we used AFR ancestry as a continuous variable, where we conducted a linear regression analysis to identify genes associated with AFR ancestry. We identified &gt;900 genes associated with AFR ancestry (p &lt; 0.01), and this gene signature distinguished Ghanaian from Ethiopian tumors in an unsupervised hierarchical clustering. In comparing the differentially expressed gene lists from these two approaches, approximately 200 genes were shared, indicating the distinct value of both analyses. Using these gene lists as input for IPA analysis, we have begun to identify canonical pathways that have been altered by our differentially expressed genes, alongside de novo networks that differ between our Ghanaian and Ethiopian tumors. In our overlapping gene list, we see predicted differences in functions such as quantity of T lymphocytes, where genes downregulated in Ethiopian tumors may indicated reduced presence of these immune cells. Using CIBERSORT and xCell deconvolution methods, validation of these findings are ongoing. Conclusions and Ongoing work: This work highlights how ancestry-specific gene regulation can delineate differences in the tumor microenvironment among a cohort of African tumors. We are currently evaluating distribution of TNBC subtypes and estimation of immune cell populations in these tumors, to determine ancestry-specific differences in tumor heterogeneity and immune response.

Citation Format: Rachel Martini, Endale Gebregzabher, Princesca Dorsaint, Timothy Chu, Kanika Arora, Lee Gibbs, Zarko Manojlovic, Nicolas Robine, Andrea Sboner, Olivier Elemento, John Carpten, Lisa Newman, Melissa Davis. Gene expression profiles of Ghanaian and Ethiopian triple-negative breast tumors [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr SS1-07.

Abstract PO-187: Identification of ancestry-specific gene expression profiles of African triple- negative breast tumors

Women with shared West African ancestry suffer from disproportionately higher breast cancer (BC) mortality and are more frequently diagnosed with the most aggressive form of BC, triple negative BC (TNBC). We and others have shown a higher prevalence of TNBC in West African (Ghanaian) and AA women, compared to lower levels of TNBC in East African (Ethiopian) and White/European American (EA) women, and that higher levels of quantified West African ancestry are found among women with TNBC compared to non-TNBC breast tumors. To further understand the role of West African ancestry in the TNBC tumor microenvironment, we conducted RNA sequencing on a pilot cohort of archival FFPE tumor tissue obtained from Ghanaian (n = 19) and Ethiopian (n = 20) women with TNBC disease to study differences in gene expression profiles. SNV calls were obtained from the RNAseq data, and 1K Genome subpopulations were used as references to quantify ancestry of each patient. We observed that our Ghanaian samples have over 94% African ancestry (corresponding to West African populations), where our Ethiopian samples showed between 37-48% African ancestry (corresponding to East African populations). Quantified ancestry among our Ethiopian samples also revealed a significant proportion of European ancestry (30-49%), corresponding to Italian ancestry. Differentially expressed genes (DEGs) were identified by comparing self- reported ethnicity and correlations with genetic ancestry. We have identified over 600 DEGs between Ghanaian vs Ethiopian TNBC samples (p &lt; 0.01). We found over 900 genes that were associated with quantified African ancestry (p &lt; 0.01), using African ancestry as a continuous variable with linear regression methods.

Unsupervised hierarchical clustering of this gene set was able to separate Ghanaian and Ethiopian samples. We compared the output of both approaches, and found that only ~200 genes were shared, showing the distinct value of both approaches. We have begun to investigate the role of our differential gene lists in canonical pathways and gene networks using Ingenuity Pathway analysis methods. Using our overlapping differential gene list from both approaches, we observed differences such as quantity of T lymphocytes, where downregulation of genes may indicate differing levels this immune cell population in the TNBC tumor environment. In our ongoing analysis, we are using deconvolution methods, such as CIBERSORT and xCell, to validate and further characterize the presence of various immune cell populations in our TNBC tumors. This study shows that quantified genetic ancestry can be used to identify ancestry-specific gene regulation in African TNBC tumors. Specifically, with our ongoing analysis of the immune cell representation in our samples, we can define differences in tumor heterogeneity and immune response in an ancestry- specific manor.

Citation Format: Rachel Martini, Endale Gebregzabher, Princesca Dorsaint, Timothy Chu, Kanika Arora, Lee Gibbs, Zarko Manojlovic, Nicolas Rovine, Andrea Sboner, Olivier Elemento, John Carpten, Lisa Newman, Melissa Davis. Identification of ancestry-specific gene expression profiles of African triple- negative breast tumors [abstract]. In: Proceedings of the AACR Virtual Conference: Thirteenth AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2020 Oct 2-4. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2020;29(12 Suppl):Abstract nr PO-187.

Abstract B067: Whole-exome sequencing of matched longitudinal primary and metastatic prostate cancers from African American men

The incidence and associated mortality and morbidity from prostate cancer (PCa) are disproportionately higher among African Americans (AA). Most deaths from PCa are attributable to the metastatic phase of the disease; the molecular underpinnings of which remain poorly understood. To identify somatic genomic variations associated with PCa metastases, we performed whole-exome sequencing (WES) on matched normal, primary, and metastatic prostate cancer tissue trios from AA patients. Microdissected formalin fixed paraffin embedded archival tissue was the source of the samples. The analysis to identify significant driver mutations and somatic copy number events was deduced using publicly available tools and custom bioinformatics processes. The cohort was comprised of matched tissue trios from 18 AA prostate cancer cases. Mean WES coverage for normal, primary and metastatic tumor were 79X, 76X and 72.5X, respectively. Percent African ancestry ranged from 0.755 to 0.997 based upon admixture informative markers from WES data. Somatic mutation counts varied among individual primary (26 to1,457) and metastatic (22 to 323) tumors. The extent of mutation overlaps between primary and metastatic PCa genomes varied substantially from 0 (0%) to 53 (14%). Recurrent mutations in the TECPR2 and PARP4 genes were found unique to metastatic cancers. Three missense mutations in the TECPR2 gene were found in the metastatic lesions of three individual patients: G344R (novel), S514N (novel) and R1066C (rs776056002). In silico functional annotation of these mutations showed to have a potentially deleterious (SIFT) and possibly/probably damaging (PolyPhen) effect. Two potentially deleterious missense mutations in the PARP4 gene namely E216V (novel) and L1080F (rs201405094) were found in the metastatic lesions in two patients. We observed differentially altered copy number events between primary and metastatic tumors. Significantly higher recurrent amplifications of chromosomal regions harboring JAK3 and ANK2 genes were observed in 26% of metastatic cases (p&lt;0.05). Copy number losses (p&lt;0.05) specific to metastatic tumors included regions harboring TPO, AFF3, and CYP7B1 genes. Here we report identification of recurrent potentially deleterious missense mutations in the TECPR2 and PARP4 genes in a subset of metastatic PCa that were not present in matched primary tumors. Several copy number alterations including amplifications of the JAK3 and ANK2 genes were also found unique to metastases. Further examination of the biological role of the identified novel mutations in the etiopathogenesis of PCa metastases is warranted, and is currently underway.

Citation Format: Zarko Manojlovic, Babak Shokrani, John D Carpten, Mezbah U Faruque. Whole-exome sequencing of matched longitudinal primary and metastatic prostate cancers from African American men [abstract]. In: Proceedings of the Twelfth AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2019 Sep 20-23; San Francisco, CA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2020;29(6 Suppl_2):Abstract nr B067.

Abstract IA26: Mechanistic differences in early-onset colorectal cancer

African Americans (AAs) have higher incidence and mortality rates of colorectal cancer (CRC) compared to other US populations. AAs present with more right-sided, microsatellite stable disease, and they are diagnosed at earlier ages compared to non-Hispanic whites (NHWs). Data from the Chicago Colorectal Cancer Consortium (CCCC) suggest that the right-sided and early-onset CRCs are distinct diseases, because early-onset CRC arises more often in the distal colon (1). In addition, analysis of stage at presentation data from the SEER database suggests that early-onset CRC is on average a more rapidly developing disease that is less likely to be prevented by colonoscopy (2). To better understand these trends, we conducted exome sequencing (n=45), copy number (n=33), and methylation analysis (n=11) of microsatellite stable AA CRCs. Results were compared to data from The Cancer Genome Atlas (TCGA). In the 43 non-hypermutable tumors, only 27 (63%) contained loss-of-function mutations in APC as compared to 80% of TCGA NHW CRCs. Importantly, APC mutation-negative CRCs were associated with an earlier age of onset of CRC (p=0.01). In the TCGA, APC mutation-negative CRCs were also associated with an earlier age of onset of CRC (p=10−5). In CCCC CRCs, APC mutation-negative CRCs were also associated with previous cancer, lower overall mutation burden, and fewer copy number variants. We conducted an analysis of DNA methylation patterns and found an epigenetic signature that was distinct from the CpG island methylator phenotype characterized in microsatellite unstable disease, referred to as CIMP. Included in the list of genes that were differentially hypermethylated in APC mutation-negative CRCs were genes that regulate the WNT signaling pathway, such as SOX9, GATA6, TET1, GLIS1, and FAT1. Using the most variable differentially methylated regions from the CCCC data, we found that these regions similarly clustered in APC mutation-negative CRCs from the TCGA. These data strongly suggest that a novel epigenetic mechanism accounts for cancer development in early-onset CRCs. Contrary to the mechanism that predominates in later-onset CRC, the early-onset mechanism does not depend on mutation in the APC gene but is associated with differential methylation of WNT pathway regulating genes instead. Our data support the claim that early-onset CRC is driven by a distinct subtype of CRC that is associated with lack of APC mutation, microsatellite and chromosome stability, lower mutation burden, and distinctive DNA methylation changes. CRC driven by epigenetic changes is consistent with the epidemiologic data suggesting that early-onset CRC develops as a more rapidly advancing disease. A deeper understanding is needed of the pathways affected by the epigenetic changes and the exposures that drive those changes in order to develop therapeutic approaches to early-onset CRC.

References

1. Xicola et al. Clin Cancer Res 2014;20:4962-70.

2. Augustus et al. PLOS One 2018 13:e0200462.

Citation Format: Rosa Xicola, Zarko Manojlovic, Gaius Augustus, Sonia Kupfer, Rajyasree Emmadi, Victoria Alagiozian-Angelova, Tim Triche, Salhia Bodour, John Carpten, Xavier Llor, Nathan Ellis. Mechanistic differences in early-onset colorectal cancer [abstract]. In: Proceedings of the Eleventh AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2018 Nov 2-5; New Orleans, LA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2020;29(6 Suppl):Abstract nr IA26.

Cancer transcriptomic profiling from rapidly enriched circulating tumor cells

Transcriptomic profiling of metastatic cancer can illuminate mechanisms of progression and lead to new therapies, but standard biopsy is invasive and reflects only a single metastatic site. In contrast, circulating tumor cell (CTC) profiling is noninvasive and repeatable, reflecting the dynamic and systemic nature of advanced disease. To date, transcriptomic profiling of CTCs has not delivered on its full potential, because white blood cells (WBCs) vastly outnumber CTCs. Current profiling strategies either lack cancer sensitivity and specificity or require specialized CTC capture protocols that are not readily scalable to large patient cohorts. Here, we describe a new strategy for rapid CTC enrichment and transcriptomic profiling using commercially available WBC depletion, microfluidic enrichment and RNA sequencing. When applied to blood samples from patients with advanced prostate cancer (PC), transcriptomes from enriched samples cluster with cancer positive controls and previously undetectable prostate-specific transcripts become readily measurable. Gene set enrichment analysis reveals multiple significantly enriched signaling pathways associated with PC, as well as novel pathways that merit further study. This accessible and scalable approach yields cancer-specific transcriptomic data and can be applied repeatedly and noninvasively in large cancer patient cohorts to discover new therapeutic targets in advanced disease.

Lack of APC somatic mutation is associated with early-onset colorectal cancer in African Americans

African Americans (AAs) have higher incidence and mortality rates of colorectal cancer (CRC) compared with other US populations. They present with more right-sided, microsatellite stable disease and are diagnosed at earlier ages compared with non-Hispanic Whites (NHWs). To gain insight into these trends, we conducted exome sequencing (n = 45), copy number (n = 33) and methylation analysis (n = 11) of microsatellite stable AA CRCs. Results were compared with data from The Cancer Genome Atlas (TCGA). Two of the 45 tumors contained POLE mutations. In the remaining 43 tumors, only 27 (63%) contained loss-of-function mutations in APC compared with 80% of TCGA NHW CRCs. APC-mutation-negative CRCs were associated with an earlier onset of CRC (P = 0.01). They were also associated with lower overall mutation burden, fewer copy number variants and a DNA methylation signature that was distinct from the CpG island methylator phenotype characterized in microsatellite unstable disease. Three of the APC-mutation-negative CRCs had loss-of-function mutations in BCL9L. Mutations in driver genes identified by TCGA exome analysis were less frequent in AA CRC cases than TCGA NHWs. Genes that regulate the WNT signaling pathway, including SOX9, GATA6, TET1, GLIS1 and FAT1, were differentially hypermethylated in APC-mutation-negative CRCs, suggesting a novel mechanism for cancer development in these tumors. In summary, we have identified a subtype of CRC that is associated with younger age of diagnosis, lack of APC mutation, microsatellite and chromosome stability, lower mutation burden and distinctive methylation changes.

Simple and rapid enrichment of circulating tumor cells (CTCs) for RNAseq in metastatic castrate resistant prostate cancer (mCRPC)

e16587

Background: CTCs have the potential to reflect not only genomic alterations but also cancer-relevant transcriptomic phenotypes. However, CTC gene expression has been hampered by signal-to-noise: rare CTC-derived transcripts are drowned out by abundant leukocyte-derived RNA. To date, a few specialized labs have achieved CTC RNAseq by capturing and analyzing single cells, a laborious and expensive approach not suitable for routine analysis of numerous samples. To address this need, we developed and validated a simple, rapid method for enrichment of live CTCs for RNAseq. Methods: Blood was drawn with informed consent under an IRB-approved protocol. Prostate cancer cell line spike-in samples were used to optimize live CTC enrichment by sequential leukocyte depletion (RosetteSep, Stem Cell Technologies) and size-based enrichment (Parsortix, Angle). Cancer-specific gene expression was first measured by multiplexed prostate specific qRT-PCR and subsequently by whole transcriptome amplification (WTA, SMARTer V2, Clontech) and RNAseq. Four patient samples were similarly analyzed by enrichment and RNAseq, along with spike-in positive controls and matched unenriched buffy coat negative controls. Results: Processing “from patient to RNA” took < 3 hrs. and achieved mean CTC recovery of 30% (range 28-33%) and mean leukocyte background of 100 (range 47-179), a 100,000-fold enrichment. Prostate specific genes (AR, PSA, PSMA) were consistently detected by qRT-PCR from enriched samples but not from unenriched samples. When analyzed by RNAseq, patient samples clustered with spike-in positive controls and away from matched buffy coat controls by principle component analysis and by unsupervised hierarchical clustering. Differential gene expression (enriched vs. matched buffy coat) identified prostate cancer-relevant transcripts. Conclusions: We developed a simple and efficient method for live CTC enrichment and expression profiling, applicable to large numbers of patient samples. This approach can be used serially over time to detect known cancer-specific transcripts and to discover new gene expression signatures that reflect tumor biology and inform disease management.

Discovery and evaluation of inhibitor of LARP6 as specific antifibrotic compound

Fibrosis is characterized by excessive production of type I collagen. Biosynthesis of type I collagen in fibrosis is augmented by binding of protein LARP6 to the 5' stem-loop structure (5'SL), which is found exclusively in type I collagen mRNAs. A high throughput screen was performed to discover inhibitors of LARP6 binding to 5'SL, as potential antifibrotic drugs. The screen yielded one compound (C9) which was able to dissociate LARP6 from 5' SL RNA in vitro and to inactivate the binding of endogenous LARP6 in cells. Treatment of hepatic stellate cells (liver cells responsible for fibrosis) with nM concentrations of C9 reduced secretion of type I collagen. In precision cut liver slices, as an ex vivo model of hepatic fibrosis, C9 attenuated the profibrotic response at 1 μM. In prophylactic and therapeutic animal models of hepatic fibrosis C9 prevented development of fibrosis or hindered the progression of ongoing fibrosis when administered at 1 mg/kg. Toxicogenetics analysis revealed that only 42 liver genes changed expression after administration of C9 for 4 weeks, suggesting minimal off target effects. Based on these results, C9 represents the first LARP6 inhibitor with significant antifibrotic activity.

Comprehensive molecular profiling of 718 Multiple Myelomas reveals significant differences in mutation frequencies between African and European descent cases

Multiple Myeloma (MM) is a plasma cell malignancy with significantly greater incidence and mortality rates among African Americans (AA) compared to Caucasians (CA). The overall goal of this study is to elucidate differences in molecular alterations in MM as a function of self-reported race and genetic ancestry. Our study utilized somatic whole exome, RNA-sequencing, and correlated clinical data from 718 MM patients from the Multiple Myeloma Research Foundation CoMMpass study Interim Analysis 9. Somatic mutational analyses based upon self-reported race corrected for ancestry revealed significant differences in mutation frequency between groups. Of interest, BCL7A, BRWD3, and AUTS2 demonstrate significantly higher mutation frequencies among AA cases. These genes are all involved in translocations in B-cell malignancies. Moreover, we detected a significant difference in mutation frequency of TP53 and IRF4 with frequencies higher among CA cases. Our study provides rationale for interrogating diverse tumor cohorts to best understand tumor genomics across populations.

This study represents the largest comprehensive molecular analysis of ethnically defined newly diagnosed treatment naïve Multiple Myeloma (MM). We revealed significant differences in mutation frequencies for important cancer genes between AA and CA MM. This study provides support for interrogating diverse tumor cohorts to best understand tumor genomics across populations.

Abstract A32: Comprehensive analysis of molecular pathogenesis of multiple myeloma by genetic ancestry

Introduction: Multiple Myeloma (MM) is a complex malignancy of plasma cells triggered by immunoglobulin gene rearrangements and well-described hyperdyploidy, accounting for slightly more than 10% of all hematologic cancers. MM is the most common hematologic malignancy in African American population with conspicuous racial disparities in both, mortality rates and incidence in cancer compared to European American. This observation evolved into a central hypothesis that MM has distinct biological differences across different ethnicities with yet unidentified race specific markers of tumor heterogeneity. Clear understanding of these molecular differences among ethnic minorities with MM will fulfill a major unmet medical need and eliminate racial disparity.

Methods: We acquired data from the Multiple Myeloma Research Foundation (MMRF) initiated comprehensive longitudinal study (CoMMpass) with an overall goal to profile 1,000 multiple myeloma patients at diagnosis, with multiple follow-up points throughout the course of the disease. To generate population subgroups based on genetic ancestry, we used a population stratification principle component analysis (PCA) and STRUCTURE to stratify myeloma patients by Ancestry Informative Markers. These well-established methods have allowed us to avoid confounders associated with self-reporting, and thus stratify the myeloma samples by genetic ancestry mapped along with anchor populations developed by 1000 genome project. We then assessed mutational frequencies as a function of PCA for each ancestry group using complex bioinformatics algorithms.

Results: We confirmed known commonly mutated genes in MM including KRAS, NRAS, FAM46C, and DIS3. Among the most striking and novel observations in our preliminary analysis of CoMMpass data using genetic ancestry and PCA was a significant difference in the frequency of nonsilent mutations in TP53, with a frequency of 7.1% (33/464) in patients clustering within the European ancestry compared to none (0/142) in African ancestry populations. Further analysis of enrichment of differentially mutated key factors within the TP53 pathway showed ATM as another gene with a significantly (p= 0.019) higher mutation frequency in EA PCA 4.7% (22/464) compared to AA PCA 0.7% (1/142). Analysis of clinical outcomes data showed poorer overall survival in patients harboring TP53 alterations. Furthermore, a comprehensive mutation analysis across samples identified a novel candidate PTCHD3 (p = 7.07E-06) with a significantly higher mutation occurrence in patients of African ancestry. Moreover, the frequency of copy number alterations known to be associated with poor prognosis revealed notable, but not significant (p=0.259) lower frequency of 1q gain in tumors from African compared to European descent. Lastly, we also observed a significant (p=0.0157) two-fold increase in early age of onset of MM in patients of African descent compared to those of European descent.

Conclusion: CoMMpass has constructed a fruitful discovery environment at nexus of high-resolution next generation deep sequencing with detailed clinical data allowing to elucidate potential ancestral drivers of MM paving the way to personalized treatments. Ultimately, these data may help us further delineate the influence of percent admixture on biological factors that drive differences in incidence and outcomes among multi-ethnic MM patients.

Citation Format: Zarko Manojlovic, Austin Christofferson, Gil Speyer, Seungchan Kim, Winnie Liang, Mary Derome, Daniel Auclair, David Craig, Jonathan Keats, John Carpten. Comprehensive analysis of molecular pathogenesis of multiple myeloma by genetic ancestry [abstract]. In: Proceedings of the AACR International Conference: New Frontiers in Cancer Research; 2017 Jan 18-22; Cape Town, South Africa. Philadelphia (PA): AACR; Cancer Res 2017;77(22 Suppl):Abstract nr A32.

A method to reduce ancestry related germline false positives in tumor only somatic variant calling

BACKGROUND

Significant clinical and research applications are driving large scale adoption of individualized tumor sequencing in cancer in order to identify tumors-specific mutations. When a matched germline sample is available, somatic mutations may be identified using comparative callers. However, matched germline samples are frequently not available such as with archival tissues, which makes it difficult to distinguish somatic from germline variants. While population databases may be used to filter out known germline variants, recent studies have shown private germline variants result in an inflated false positive rate in unmatched tumor samples, and the number germline false positives in an individual may be related to ancestry.

METHODS

First, we examined the relationship between the germline false positives and ancestry. Then we developed and implemented a tumor only caller (LumosVar) that leverages differences in allelic frequency between somatic and germline variants in impure tumors. We used simulated data to systematically examine how copy number alterations, tumor purity, and sequencing depth should affect the sensitivity of our caller. Finally, we evaluated the caller on real data.

RESULTS

We find the germline false-positive rate is significantly higher for individuals of non-European Ancestry largely due to the limited diversity in public polymorphism databases and due to population-specific characteristics such as admixture or recent expansions. Our Bayesian tumor only caller (LumosVar) is able to greatly reduce false positives from private germline variants, and our sensitivity is similar to predictions based on simulated data.

CONCLUSIONS

Taken together, our results suggest that studies of individuals of non-European ancestry would most benefit from our approach. However, high sensitivity requires sufficiently impure tumors and adequate sequencing depth. Even in impure tumors, there are copy number alterations that result in germline and somatic variants having similar allele frequencies, limiting the sensitivity of the approach. We believe our approach could greatly improve the analysis of archival samples in a research setting where the normal is not available.

Abstract 3389: Integrated genomic analyses reveal frequent TERT aberrations in acral melanoma

Genomic analyses of cutaneous melanoma (CM) have yielded biological and therapeutic insights, but understanding of non-ultraviolet (UV)-derived CMs remains limited. Deeper analysis of acral lentiginous melanoma (ALM), a rare sun-shielded melanoma subtype associated with worse survival than CM, is needed to delineate non-UV oncogenic mechanisms. We thus performed comprehensive genomic and transcriptomic analysis of 34 ALM patients. Unlike CM, somatic alterations were dominated by structural variation and absence of UV-derived mutation signatures. Only 38% of patients demonstrated driver BRAF/NRAS/NF1 mutations. Contrasting with CM, we observed PAK1 copy gains in 15% of patients, and somatic TERT translocations, copy gains, and missense and promoter mutations, or germline events, in 41% of patients. We further show that in vitro TERT inhibition has cytotoxic effects on primary ALM cells. These findings provide insight into the role of TERT in ALM tumorigenesis, and reveal preliminary evidence that TERT inhibition represents a potential therapeutic strategy in ALM.

Citation Format: Winnie S. Liang, William Hendricks, Jeffrey Kiefer, Jessica Schmidt, Shobana Sekar, John Carpten, David W. Craig, Jonathan Adkins, Lori Cuyugan, Zarko Manojlovic, Rebecca F. Halperin, Adrienne Helland, Sara Nasser, Christophe Legendre, Laurence H. Hurley, Karthigayini Sivaprakasam, Douglas B. Johnson, Holly Crandall, Klaus J. Busam, Victoria Zismann, Valerie De Luca, Jeeyun Lee, Aleksandar Sekulic, Charlotte E. Ariyan, Jeffrey Sosman, Jeffrey Trent. Integrated genomic analyses reveal frequent TERT aberrations in acral melanoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 3389. doi:10.1158/1538-7445.AM2017-3389

Integrated genomic analyses reveal frequent TERT aberrations in acral melanoma

Genomic analyses of cutaneous melanoma (CM) have yielded biological and therapeutic insights, but understanding of non-ultraviolet (UV)-derived CMs remains limited. Deeper analysis of acral lentiginous melanoma (ALM), a rare sun-shielded melanoma subtype associated with worse survival than CM, is needed to delineate non-UV oncogenic mechanisms. We thus performed comprehensive genomic and transcriptomic analysis of 34 ALM patients. Unlike CM, somatic alterations were dominated by structural variation and absence of UV-derived mutation signatures. Only 38% of patients demonstrated driver BRAF/NRAS/NF1 mutations. In contrast with CM, we observed PAK1 copy gains in 15% of patients, and somatic TERT translocations, copy gains, and missense and promoter mutations, or germline events, in 41% of patients. We further show that in vitro TERT inhibition has cytotoxic effects on primary ALM cells. These findings provide insight into the role of TERT in ALM tumorigenesis and reveal preliminary evidence that TERT inhibition represents a potential therapeutic strategy in ALM.

La-related protein 6 controls ciliated cell differentiation

Background

La-related protein 6 (LARP6) is an evolutionally conserved RNA-binding protein. Vertebrate LARP6 binds the 5′ stem-loop found in mRNAs encoding type I collagen to regulate their translation, but other target mRNAs and additional functions for LARP6 are unknown. The aim of this study was to elucidate an additional function of LARP6 and to evaluate the importance of its function during development.

Methods

To uncover the role of LARP6 in development, we utilized Morpholino Oligos to deplete LARP6 protein in Xenopus embryos. Then, embryonic phenotypes and ciliary structures of LAPR6 morphants were examined. To identify the molecular mechanism underlying ciliogenesis regulated by LARP6, we tested the expression level of cilia-related genes, which play important roles in ciliogenesis, by RT-PCR or whole mount in situ hybridization (WISH).

Results

We knocked down LARP6 in Xenopus embryos and found neural tube closure defects. LARP6 mutant, which compromises the collagen synthesis, could rescue these defects. Neural tube closure defects are coincident with lack of cilia, antenna-like cellular organelles with motility- or sensory-related functions, in the neural tube. The absence of cilia at the epidermis was also observed in LARP6 morphants, and this defect was due to the absence of basal bodies which are formed from centrioles and required for ciliary assembly. In the process of multi-ciliated cell (MCC) differentiation, mcidas, which activates the transcription of genes required for centriole formation during ciliogenesis, could partially restore MCCs in LARP6 morphants. In addition, LARP6 likely controls the expression of mcidas in a Notch-independent manner.

Conclusions

La-related protein 6 is involved in ciliated cell differentiation during development by controlling the expression of cilia-related genes including mcidas. This LARP6 function involves a mechanism that is distinct from its established role in binding to collagen mRNAs and regulating their translation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13630-017-0047-7) contains supplementary material, which is available to authorized users.

Molecular Genetic Profiling of Adolescent Glassy Cell Carcinoma of the Cervix Reveals Targetable EGFR Amplification with Potential Therapeutic Implications

Glassy cell carcinoma of the cervix (GCCC) is a very rare and aggressive form of cervical cancer. An adolescent female with advanced metastatic disease was enrolled in our genomic profiling research protocol. We identified high-level amplification of epidermal growth factor receptor (EGFR) and Yes-associated protein-1 (YAP1), which led to the addition of EGFR inhibitors to the chemotherapy regimen. Here, we report the first genetically profiled case of GCCC with potential therapeutic implications.

RFX7 is required for the formation of cilia in the neural tube

Regulatory Factor X (RFX) transcription factors are important for development and are likely involved in the pathogenesis of serious human diseases including ciliopathies. While seven RFX genes have been identified in vertebrates and several RFX transcription factors have been reported to be regulators of ciliogenesis, the role of RFX7 in development including ciliogenesis is not known. Here we show that RFX7 in Xenopus laevis is expressed in the neural tube, eye, otic vesicles, and somites. Knockdown of RFX7 in Xenopus embryos resulted in a defect of ciliogenesis in the neural tube and failure of neural tube closure. RFX7 controlled the formation of cilia by regulating the expression of RFX4 gene, which has been reported to be required for ciliogenesis in the neural tube. Moreover, ectopic expression of Foxj1, which is a master regulator of motile cilia formation, suppressed the expression of RFX4 but not RFX7. Taken together, RFX7 plays an important role in the process of neural tube closure at the top of the molecular cascade which controls ciliogenesis in the neural tube.

Serine-threonine kinase receptor-associated protein (STRAP) regulates translation of type I collagen mRNAs

Type I collagen is the most abundant protein in the human body and is composed of two α1(I) and one α2(I) polypeptides which assemble into a triple helix. For the proper assembly of the collagen triple helix, the individual polypeptides must be translated in coordination. Here, we show that serine-threonine kinase receptor-associated protein (STRAP) is tethered to collagen mRNAs by interaction with LARP6. LARP6 is a protein which directly binds the 5' stem-loop (5'SL) present in collagen α1(I) and α2(I) mRNAs, but it interacts with STRAP with its C-terminal domain, which is not involved in binding 5'SL. Being tethered to collagen mRNAs, STRAP prevents unrestricted translation, primarily that of collagen α2(I) mRNAs, by interacting with eukaryotic translation initiation factor 4A (eIF4A). In the absence of STRAP, more collagen α2(I) mRNA can be pulled down with eIF4A, and collagen α2(I) mRNA is unrestrictedly loaded onto the polysomes. This results in an imbalance of synthesis of α1(I) and α2(I) polypeptides, in hypermodifications of α1(I) polypeptide, and in inefficient assembly of the polypeptides into a collagen trimer and their secretion as monomers. These defects can be partially restored by supplementing STRAP. Thus, we discovered STRAP as a novel regulator of the coordinated translation of collagen mRNAs.

Tacrolimus (FK506) prevents early stages of ethanol induced hepatic fibrosis by targeting LARP6 dependent mechanism of collagen synthesis

Tacrolimus (FK506) is a widely used immunosuppressive drug. Its effects on hepatic fibrosis have been controversial and attributed to immunosuppression. We show that in vitro FK506, inhibited synthesis of type I collagen polypeptides, without affecting expression of collagen mRNAs. In vivo, administration of FK506 at a dose of 4 mg/kg completely prevented development of alcohol/carbon tetrachloride induced liver fibrosis in rats. Activation of hepatic stellate cells (HSCs) was absent in the FK506 treated livers and expression of collagen α2(I) mRNA was at normal levels. Collagen α1(I) mRNA was increased in the FK506 treated livers, but this mRNA was not translated into α1(I) polypeptide. No significant inflammation was associated with the fibrosis model used. FK506 binding protein 3 (FKBP3) is one of cellular proteins which binds FK506 with high affinity. We discovered that FKBP3 interacts with LARP6 and LARP6 is the major regulator of translation and stability of collagen mRNAs. In the presence of FK506 the interaction between FKBP3 and LARP6 is weakened and so is the pull down of collagen mRNAs with FKBP3. We postulate that FK506 inactivates FKBP3 and that lack of interaction of LARP6 and FKBP3 results in aberrant translation of collagen mRNAs and prevention of fibrosis. This is the first report of such activity of FK506 and may renew the interest in using this drug to alleviate hepatic fibrosis.

A novel role of RNA helicase A in regulation of translation of type I collagen mRNAs

Type I collagen is composed of two α1(I) polypeptides and one α2(I) polypeptide and is the most abundant protein in the human body. Expression of type I collagen is primarily controlled at the level of mRNA stability and translation. Coordinated translation of α(I) and α2(I) mRNAs is necessary for efficient folding of the corresponding peptides into the collagen heterotrimer. In the 5' untranslated region (5' UTR), collagen mRNAs have a unique 5' stem-loop structure (5' SL). La ribonucleoprotein domain family member 6 (LARP6) is the protein that binds 5' SL with high affinity and specificity and coordinates their translation. Here we show that RNA helicase A (RHA) is tethered to the 5' SL of collagen mRNAs by interaction with the C-terminal domain of LARP6. In vivo, collagen mRNAs immunoprecipitate with RHA in an LARP6-dependent manner. Knockdown of RHA prevents formation of polysomes on collagen mRNAs and dramatically reduces synthesis of collagen protein, without affecting the level of the mRNAs. A reporter mRNA with collagen 5' SL is translated three times more efficiently in the presence of RHA than the same reporter without the 5' SL, indicating that the 5' SL is the cis-acting element conferring the regulation. During activation of quiescent cells into collagen-producing cells, expression of RHA is highly up-regulated. We postulate that RHA is recruited to the 5' UTR of collagen mRNAs by LARP6 to facilitate their translation. Thus, RHA has been discovered as a critical factor for synthesis of the most abundant protein in the human body.