Search

Background Large β-globin gene cluster deletions (hereditary persistence of fetal hemoglobin [Hb] or β-, δβ-, γδβ-, and ϵγδβ-thalassemia), are associated with widely disparate phenotypes, including variable degrees of microcytic anemia and Hb F levels. When present, increased Hb A2 is used as a surrogate marker for β-thalassemia. Notably, ϵγδβ-thalassemias lack the essential regulatory locus control region (LCR) and cause severe transient perinatal anemia but normal newborn screen (NBS) results and Hb A2 levels. Herein, we report a novel deletion of the ϵ, Aγ, Gγ, and ψβ loci with intact LCR, δ-, and β-regions in 2 women and newborn twins. Methods Capillary electrophoresis (CE), high-performance liquid chromatography (HPLC), DNA sequencing, multiplex ligation-dependent probe amplification (MLPA), gap-polymerase chain reaction (gap-PCR), and long-read sequencing (LRS) were performed. Results NBS showed an Hb A > Hb F pattern for both twins. At 20 months, Hb A2 was increased similarly to that in the mother and an unrelated woman. Unexplained microcytosis was absent and the twins lacked severe neonatal anemia. MLPA, LRS, and gap-PCR confirmed a 32 599 base pair deletion of ϵ (HBE1) through ψβ (HBBP1) loci. Conclusions This deletion represents a hemoglobinopathy category with a distinct phenotype that has not been previously described, an ϵγ-thalassemia. Both the NBS Hb A > F pattern and the subsequent increased Hb A2 without microcytosis are unusual. A similar deletion should be considered when this pattern is encountered and appropriate test methods selected for detection. Knowledge of the clinical impact of this new category will improve genetic counselling, with distinction from the severe transient anemia associated with ϵγδβ-thalassemia.

Purpose: Androgen receptor (AR) variant AR-V7 is a ligand-independent transcription factor that promotes prostate cancer resistance to AR-targeted therapies. Accordingly, efforts are under way to develop strategies for monitoring and inhibiting AR-V7 in castration-resistant prostate cancer (CRPC). The purpose of this study was to understand whether other AR variants may be coexpressed with AR-V7 and promote resistance to AR-targeted therapies.Experimental Design: We utilized complementary short- and long-read sequencing of intact AR mRNA isoforms to characterize AR expression in CRPC models. Coexpression of AR-V7 and AR-V9 mRNA in CRPC metastases and circulating tumor cells was assessed by RNA-seq and RT-PCR, respectively. Expression of AR-V9 protein in CRPC models was evaluated with polyclonal antisera. Multivariate analysis was performed to test whether AR variant mRNA expression in metastatic tissues was associated with a 12-week progression-free survival endpoint in a prospective clinical trial of 78 CRPC-stage patients initiating therapy with the androgen synthesis inhibitor, abiraterone acetate.Results: AR-V9 was frequently coexpressed with AR-V7. Both AR variant species were found to share a common 3′ terminal cryptic exon, which rendered AR-V9 susceptible to experimental manipulations that were previously thought to target AR-V7 uniquely. AR-V9 promoted ligand-independent growth of prostate cancer cells. High AR-V9 mRNA expression in CRPC metastases was predictive of primary resistance to abiraterone acetate (HR = 4.0; 95% confidence interval, 1.31–12.2; P = 0.02).Conclusions: AR-V9 may be an important component of therapeutic resistance in CRPC. Clin Cancer Res; 23(16); 4704–15. ©2017 AACR.

Supplementary Table 1: Baseline Characteristics of the Patients; Supplementary Table 2: Baseline RNA-seq Expression Levels for AR/ARVs (spliced reads per million); Supplementary Table 3: TaqMan Gene Expression Assays with Circulating Tumor Cells; Supplementary Figure 1: Actinomycin D responses of AR-V9 and AR-V7 mRNA in 22Rv1 cells; Supplementary Figure 2: 3’ rapid amplification of cDNA ends (RACE) for single molecule real-time (SMRT) sequencing; Supplementary Figure 3: Purified AR-V9 polyclonal antibodies detect overexpressed AR-V9 in western blot with overexpressed AR-V9; Supplementary Figure 4: Relationship between tumor RNA p13 purity, AR-V7, and AR-V9 expression levels in CRPC metastases.

Supplementary Tables S1-S6 from Glutathione S-Transferase T1 and M1: Gene Sequence Variation and Functional Genomics

Purpose: The glutathione S-transferases (GSTs) catalyze the glutathione conjugation of reactive electrophiles, including carcinogens and many antineoplastic drugs. GSTT1 and GSTM1 are polymorphically deleted, but the full range of genetic variation in these two genes has not yet been explored. We set out to systematically identify common polymorphisms in GSTT1 and GSTM1, followed by functional genomic studies.Experimental Design: First, multiplex PCR was used to determine GSTT1 and GSTM1 copy number in 400 DNA samples (100 each from 4 ethnic groups). Exons, splice junctions, and 5′-flanking regions (5′-FR) were then resequenced using DNA samples that contained at least one copy of GSTT1 or GSTM1.Results: Gene deletion frequencies among ethnic groups were from 33.5% to 73.5% for GSTT1 and from 50.5% to 78.0% for GSTM1. GSTT1 deletion data correlated with the results of mRNA microarray expression studies. The 18 single nucleotide polymorphisms (SNP) observed in GSTT1 included three nonsynonymous coding SNPs (cSNPs) and one single-nucleotide deletion, whereas the 51 GSTM1 SNPs included two nonsynonymous cSNPs. Two of the GSTT1 nonsynonymous cSNPs resulted in decreases in levels of immunoreactive protein to 56% and 12% of wild type (WT), whereas those in GSTM1 resulted in modest increases in protein levels. Reporter gene assays showed that one GSTT1 5′-FR haplotype, with a frequency of 32% in African-American subjects, resulted in an increase in transcription in JEG-3 cells to 351% of that for the WT sequence, and one GSTM1 5′-FR haplotype resulted in an increase in transcription in JEG-3 cells to 129% of WT.Conclusions: These observations suggest that functionally significant pharmacogenomic variation beyond GSTT1 and GSTM1 gene deletion may contribute to carcinogenesis or individual variation in antineoplastic drug therapy response.

Supplementary Tables 1-4, Figures 1-2 from Functional Genetic Polymorphisms in the Aromatase Gene CYP19 Vary the Response of Breast Cancer Patients to Neoadjuvant Therapy with Aromatase Inhibitors

Supplementary Materials from Human Aromatase: Gene Resequencing and Functional Genomics

Glutathione S-transferase P1 (GSTP1) is of importance for cancer research because of its role in detoxifying carcinogens, activating antineoplastic prodrugs, metabolizing chemotherapeutic agents, and its involvement in cell cycle and apoptosis regulation. Two common GSTP1 genetic polymorphisms have been studied extensively. However, the full range of GSTP1 genetic variation has not been systematically characterized in the absence of disease pathology. We set out to identify common GSTP1 polymorphisms in four ethnic groups, followed by functional genomic studies. All exons, splice junctions, and the 5′-flanking region of GSTP1 were resequenced using 60 DNA samples each from four ethnic groups. The 35 single-nucleotide polymorphisms (SNP) identified included six nonsynonymous SNPs and 17 previously unreported polymorphisms. GSTP1 variant allozymes were then expressed in COS-1 cells, and five displayed significantly altered levels of enzyme activity. One decreased to 22% of the wild-type (WT) activity. Four variant allozymes had Km values that differed significantly from that of the WT, and five showed altered levels of immunoreactive protein compared with WT, with a significant correlation (r = 0.79, P < 0.007) between levels of immunoreactive protein and enzyme activity in these samples. In the Mexican American population, five linked SNPs were significantly associated with GSTP1 mRNA expression, one of which was found by electrophoretic mobility shift assay to alter protein binding. These studies have identified functionally significant genetic variation, in addition to the two frequently studied GSTP1 nonsynonymous SNPs, that may influence GSTP1's contribution to carcinogen and drug metabolism, and possibly disease pathogenesis and/or drug response. [Cancer Res 2008;68(12):4791–801]

Supplementary Data from Glutathione S-Transferase P1: Gene Sequence Variation and Functional Genomic Studies

Aromatase (CYP19) is a critical enzyme in estrogen biosynthesis and aromatase inhibitors (AI) are employed widely for endocrine therapy in postmenopausal women with breast cancer. We hypothesized that single nucleotide polymorphisms (SNPs) in the CYP19 gene may alter the effectiveness of AI therapy in the neoadjuvant setting. Genomic DNA was obtained for sequencing from 52 women pre-AI and post-AI treatment in this setting. Additionally, genomic DNA obtained from 82 samples of breast cancer and 19 samples of normal breast tissue was subjected to resequencing. No differences in CYP19 sequence were observed between tumor and germ-line DNA in the same patient. A total of 48 SNPs were identified including 4 novel SNPs when compared with previous resequencing data. For genotype-phenotype association studies, we determined the levels of aromatase activity, estrone, estradiol, and tumor size in patients pre-AI and post-AI treatment. We defined two tightly linked SNPs (rs6493497 and rs7176005 in the 5'-flanking region of CYP19 exon 1.1) that were significantly associated with a greater change in aromatase activity after AI treatment. In a follow-up study of 200 women with early-stage breast cancer who were treated with adjuvant anastrozole, these same two SNPs were also associated with higher plasma estradiol levels in patients pre-AI and post-AI treatment. Electrophoretic mobility shift and reporter gene assays confirmed likely functional effects of these two SNPs on transcription of CYP19. Our findings indicate that two common genetic polymorphisms in the aromatase gene CYP19 vary the response of breast cancer patients to aromatase inhibitors. Cancer Res; 70(1); 319–28

Coronavirus disease 2019 (COVID-19) is associated with significant morbidity and mortality, albeit with considerable heterogeneity among affected individuals. It remains unclear which factors determine severity of illness and long-term outcomes. Emerging evidence points towards an important role of preexisting host factors, such as a deregulated inflammatory response at the time of infection. Here, we demonstrate the negative impact of clonal hematopoiesis, a prevalent clonal disorder of ageing individuals, on COVID-19-related cytokine release severity and mortality. We show that mutations in the gene coding for the methylcytosine dioxygenase TET2 promotes amplification of classical and intermediate monocyte subsets. Using single cell multiomic sequencing approaches, we identify cell-specific gene expression changes associated with the loss of TET2 and significant epigenomic deregulation affecting enhancer accessibility of a subset of transcription factors involved in monocyte differentiation. We further identify EGR1 down-regulation secondary to TET2-mediated hypermethylation, resulting in overexpression of MALAT1, a lncRNA that plays a role in hematopoietic stem cell differentiation and monocyte lineage commitment. Together, these data provide a mechanistic insight to the poor prognostic value of clonal hematopoiesis in patients infected with Sars-COV2.

SARS-CoV-2 has had an unprecedented impact on human health and highlights the need for genomic epidemiology studies to increase our understanding of the evolution and spread of pathogens and to inform policy decisions. Most efforts have focused on international or country-wide transmission, which are unable to highlight state-wide trends. We sequenced virus genomes from over 22,000 patients tested at Mayo Clinic Laboratories between 2020-2022 and leveraged detailed patient metadata to describe county-to-county spread in Minnesota. Our findings indicate that spread in the state was mostly dominated by viruses from Hennepin County, which contains the largest metropolis. For many counties, we found that state government restrictions eventually led to a decrease in the diversity of circulating viruses from other counties and that their complete removal in May of 2021 saw a drastic revert to levels at or greater than those observed during the months before. We also linked over 14,000 genomes with patient risk characteristics and infection-related phenotypes from the Mayo Clinic electronic health record. We found that the genetic relationship of Omicron viruses was structured by clinical outcomes when stratifying by patient risk factor and variant of concern. However, we were unable to identify nucleotide variants that drove this association.

SARS-CoV-2 has had an unprecedented impact on human health and highlights the need for genomic epidemiology studies to increase our understanding of virus evolution and spread, and to inform policy decisions. We sequenced viral genomes from over 22,000 patient samples tested at Mayo Clinic Laboratories between 2020-2022 and use Bayesian phylodynamics to describe county and regional spread in Minnesota.The earliest introduction into Minnesota was to Hennepin County from a domestic source around January 22, 2020; six weeks before the first confirmed case in the state. This led to the virus spreading to Northern Minnesota, and eventually, the rest of the state. International introductions were most abundant in Hennepin (home to the Minneapolis/St. Paul International (MSP) airport) totaling 45 (out of 107) over the two-year period. Southern Minnesota counties were most common for domestic introductions with 19 (out of 64), potentially driven by bordering states such as Iowa and Wisconsin as well as Illinois which is nearby. Hennepin also was, by far, the most dominant source of in-state transmissions to other Minnesota locations (n=772) over the two-year period.We also analyzed the diversity of the location source of SARS-CoV-2 viruses in each county and noted the timing of state-wide policies as well as trends in clinical cases. Neither the number of clinical cases or major policy decisions, such as the end of the lockdown period in 2020 or the end of all restrictions in 2021, appeared to have impact on virus diversity across each individual county.ImportanceWe analyzed over 22,000 SARS-CoV-2 genomes of patient samples tested at Mayo Clinic Laboratories during a two-year period in the COVID-19 pandemic that included Alpha, Delta, and Omicron VoCs to examine the roles and relationships of Minnesota virus transmission.We found that Hennepin County, the most populous county, drove the transmission of SARS-CoV-2 viruses in the state after including the formation of earlier clades including 20A, 20C, and 20G, as well as variants of concern Alpha and Delta. We also found that Hennepin County was the source for most of the county-to-county introductions after its initial introduction with the virus in early 2020 from an international source, while other counties acted as transmission “sinks”. In addition, major policies such as the end of the lockdown period in 2020 or the end of all restrictions in 2021, did not appear to have an impact on virus diversity across individual counties.

DNA junctions (DNAJs) frequently impact clinically relevant genes in tumors and are important for diagnostic and therapeutic purposes. Although routinely screened through fluorescence in situ hybridization assays, such testing only allows the interrogation of single-gene regions or known fusion partners. Comprehensive assessment of DNAJs present across the entire genome can only be determined from whole-genome sequencing. Structural variance analysis from whole-genome paired-end sequencing data is, however, frequently restricted to copy number changes without DNAJ detection. Through optimized whole-genome sequencing and specialized bioinformatics algorithms, complete structural variance analysis is reported, including DNAJs, from formalin-fixed DNA. Selective library assembly from larger fragments (500 bp) and economical sequencing depths (300 to 400 million reads) provide representative genomic coverage profiles and increased allelic coverage to levels compatible with DNAJ calling (40× to 60×). Although consistently fragmented, more recently formalin-fixed, specimens (2 years' storage) revealed consistent populations of larger DNA fragments. Optimized bioinformatics efficiently detected90% of DNAJs in two prostate tumors (approximately 60% tumor) previously analyzed by mate-pair sequencing on fresh frozen tissue, with evidence of at least one spanning-read in 99% of DNAJs. Rigorous masking with data from unrelated formalin-fixed tissue progressively eliminated many false-positive DNAJs, without loss of true positives, resulting in low numbers of false-positive passing current filters. This methodology enables more comprehensive clinical genomics testing on formalin-fixed clinical specimens.

To examine the length of a hexanucleotide expansion in C9orf72, which represents the most frequent genetic cause of frontotemporal lobar degeneration and motor neuron disease, we employed a targeted amplification-free long-read sequencing technology: No-Amp sequencing. In our cross-sectional study, we assessed cerebellar tissue from 28 well-characterized C9orf72 expansion carriers. We obtained 3507 on-target circular consensus sequencing reads, of which 814 bridged the C9orf72 repeat expansion (23%). Importantly, we observed a significant correlation between expansion sizes obtained using No-Amp sequencing and Southern blotting (P = 5.0 × 10−4). Interestingly, we also detected a significant survival advantage for individuals with smaller expansions (P = 0.004). Additionally, we uncovered that smaller expansions were significantly associated with higher levels of C9orf72 transcripts containing intron 1b (P = 0.003), poly(GP) proteins (P = 1.3 × 10− 5), and poly(GA) proteins (P = 0.005). Thorough examination of the composition of the expansion revealed that its GC content was extremely high (median: 100%) and that it was mainly composed of GGGGCC repeats (median: 96%), suggesting that expanded C9orf72 repeats are quite pure. Taken together, our findings demonstrate that No-Amp sequencing is a powerful tool that enables the discovery of relevant clinicopathological associations, highlighting the important role played by the cerebellar size of the expanded repeat in C9orf72-linked diseases., DeJesus-Hernandez et al. employ an innovative long-read sequencing method to examine the length of the expanded C9orf72 repeat that represents the most common genetic cause of frontotemporal dementia and amyotrophic lateral sclerosis. Smaller expansion size confers a survival benefit for patients.

Background: COVID-19 causes significant morbidity and mortality, albeit with considerable heterogeneity among affected individuals. It remains unclear which host factors determine disease severity and survival. Given the propensity of clonal hematopoiesis (CH) to promote inflammation in healthy individuals, we investigated its effect on COVID-19 outcomes. Methods: We performed a multi-omics interrogation of the genome, epigenome, transcriptome, and proteome of peripheral blood mononuclear cells from COVID-19 patients (n=227). We obtained clinical data, laboratory studies, and survival outcomes. We determined CH status and TET2-related DNA methylation. We performed single-cell proteogenomics to understand clonal composition in relation to cell phenotype. We interrogated single-cell gene expression in isolation and in conjunction with DNA accessibility. We integrated these multi-omics data to understand the effect of CH on clonal composition, gene expression, methylation of cis-regulatory elements, and lineage commitment in COVID-19 patients. We performed shRNA knockdowns to validate the effect of one candidate transcription factor in myeloid cell lines. Results: The presence of CH was strongly associated with COVID-19 severity and all-cause mortality, independent of age (HR 3.48, 95% CI 1.45-8.36, p=0.005). Differential methylation of promoters and enhancers was prevalent in TET2-mutant, but not DNMT3A-mutant CH. TET2-mutant CH was associated with enhanced classical/intermediate monocytosis and single-cell proteogenomics confirmed an enrichment of TET2 mutations in these cell types. We identified cell-type specific gene expression changes associated with TET2 mutations in 102,072 single cells (n=34). Single-cell RNA-seq confirmed the skewing of hematopoiesis towards classical and intermediate monocytes and demonstrated the downregulation of EGR1 (a transcription factor important for monocyte differentiation) along with up-regulation of the lncRNA MALAT1 in monocytes. Combined scRNA-/scATAC-seq in 43,160 single cells (n=18) confirmed the skewing of hematopoiesis and up-regulation of MALAT1 in monocytes along with decreased accessibility of EGR1 motifs in known cis-regulatory elements. Using myeloid cell lines for functional validation, shRNA knockdowns of EGR1 confirmed the up-regulation of MALAT1 (in comparison to wildtype controls). Conclusions: CH is an independent prognostic factor in COVID-19 and skews hematopoiesis towards monocytosis. TET2-mutant CH is characterized by differential methylation and accessibility of enhancers binding myeloid transcriptions factors including EGR1. The ensuing loss of EGR1 expression in monocytes causes MALAT1 overexpression, a factor known to promote monocyte differentiation and inflammation. These data provide a mechanistic insight to the adverse prognostic impact of CH in COVID-19. Citation Format: Moritz Binder, Terra L. Lasho, Wazim Mohammed Ismail, Nana A. Ben-Crentsil, Jenna A. Fernandez, Minsuk Kim, Susan M. Geyer, Amelia Mazzone, Christy M. Finke, Abhishek A. Mangaonkar, Jeong-Heon Lee, Kwan Hyun Kim, Vernadette A. Simon, Fariborz Rakhshan Rohakthar, Amik Munankarmy, Susan M. Schwager, Jonathan J. Harrington, Melissa R. Snyder, Nathalie M. Droin, Eric Solary, Keith D. Robertson, Eric D. Wieben, Eric Padron, Nicholas Chia, Alexandre Gaspar-Maia, Mrinal M. Patnaik. Enhancer deregulation inTET2-mutant clonal hematopoiesis is associated with increased COVID-19 severity and mortality [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 806.

Expansion of CTG trinucleotide repeats (TNR) in the transcription factor 4 (TCF4) gene is highly associated with Fuchs Endothelial Corneal Dystrophy (FECD). Due to limitations in the availability of DNA from diseased corneal endothelium, sizing of CTG repeats in FECD patients has typically been determined using DNA samples isolated from peripheral blood leukocytes. However, it is non-feasible to extract enough DNA from surgically isolated FECD corneal endothelial tissue to determine repeat length based on current technology. To circumvent this issue, total RNA was isolated from FECD corneal endothelium and sequenced using long-read sequencing. Southern blotting of DNA samples isolated from primary cultures of corneal endothelium from these same affected individuals was also assessed. Both long read sequencing and Southern blot analysis showed significantly longer CTG TNR expansion (>1000 repeats) in the corneal endothelium from FECD patients than those characterized in leukocytes from the same individuals (TCF4 CTG repeat expansions in the FECD corneal endothelium are much longer than those found in leukocytes.

Purpose: CTG trinucleotide repeat (TNR) expansion is frequently found in transcription factor 4 (TCF4) in Fuchs' endothelial corneal dystrophy (FECD), though the effect of TNR expansion on FECD pathophysiology remains unclear. The purpose of this study was to evaluate the effect of TNR expansion on TCF4 expression in corneal endothelium of patients with FECD. \ud\udMethods: Peripheral blood DNA and Descemet membrane with corneal endothelium were obtained from 203 German patients with FECD. The CTG TNR repeat length in TCF4 was determined by short tandem repeat (STR) assays and Southern blotting using genomic DNA. Genotyping of rs613872 in TCF4 was performed by PCR. TCF4 mRNA levels in corneal endothelium were evaluated by quantitative PCR using three different probes. Control corneal endothelial samples were obtained from 35 non-FECD subjects. \ud\udResults: The STR assay and Southern blotting showed that 162 of the 203 patients with FECD (80%) harbored CTG trinucleotide repeat lengths larger than 50. Quantitative PCR using all three probes demonstrated that TCF4 mRNA is significantly upregulated in the corneal endothelium of patients with FECD, regardless of the presence of TNR expansion. However, the length of the TNR tended to show a positive correlation with TCF4 expression level. No correlation was shown between the genotype of TCF4 SNP, rs613872, and the level of TCF4 expression. \ud\udConclusions: Our findings showed that TCF4 mRNA is upregulated in the corneal endothelium of patients with FECD. Further studies on the effects of TCF4 upregulation on corneal endothelial cell function will aid in understanding the pathophysiology of FECD.

PURPOSE To investigate the association of body mass index (BMI) with Fuchs endothelial corneal dystrophy (FECD) severity and TCF4 CTG18.1 expansion. METHODS A total of 343 patients with FECD were enrolled from the Mayo Clinic. FECD severity was graded by slit-lamp biomicroscopy. BMI values were obtained from the electronic medical records. DNA extracted from leukocytes was analyzed for CTG18.1 expansion length, with ≥40 repeats considered expanded. Wilcoxon signed-rank tests were used to compare FECD grade and CTG18.1 expansion length in patients by BMI (

Purpose: To investigate candidate genes associated with familial strabismus and propose a theory of their interaction in familial strabismus associated with early neurodevelopment. Methods: Eighteen families, including 53 patients diagnosed with strabismus and 34 unaffected family members, were analyzed. All patients with strabismus and available unaffected family members were evaluated using whole exome sequencing. The primary outcome was to identify rare occurring variants among affected individuals and investigate the evidence of their genetic heterogeneity. These results were compared with exome sequencing analysis to build a comprehensive genetic profile of the study families. Results: We observed 60 variants from 58 genes in 53 patients diagnosed with strabismus. We prioritized the most credible risk variants, which showed clear segregation in family members affected by strabismus. As a result, we found risk variants in four genes (FAT3, KCNH2, CELSR1, and TTYH1) in five families, suggesting their role in development of familial strabismus. In other families, there were several rare genetic variants in affected cases, but we did not find clear segregation pattern across family members. Conclusion: Genomic sequencing holds great promise in elucidating the genetic causes of strabismus, further research with larger cohorts or other related approaches are warranted.

Purpose To characterize inheritance, penetrance, and trinucleotide repeat expansion stability in Fuchs endothelial corneal dystrophy (FECD). Methods One thousand unrelated and related subjects with and without FECD were prospectively recruited. CTG18.1 repeat length (CTG18.1L) was determined via short tandem repeat assay and Southern blotting of leukocyte DNA. Multivariable logistic regression and generalized estimating equation models were employed. Results There were 546 unrelated FECD cases (67.6% female; 70 ± 10 years) and 235 controls (63.8% female; 73 ± 8 years; all ≥ 50 years). CTG18.1 expansion (CTG18.1exp+) was observed in 424 (77.7%) cases and 18 (7.7%) controls (P = 2.48 × 10-44). CTG18.1 expansion was associated with FECD severity (P = 5.62 × 10-7). The family arm of the study included 331 members from 112 FECD-affected families; 87 families were CTG18.1exp+. Autosomal dominant inheritance with variable expression of FECD was observed, regardless of expansion status. FECD penetrance of CTG18.1 expansion increased with age, ranging from 44.4% in the youngest (19-46 years) to 86.2% in the oldest (64-91 years) age quartiles. Among 62 parent-offspring transmissions of CTG18.1exp+, 48 (77.4%) had a change in CTG18.1L ≤ 10 repeats, and eight (12.9%) were ≥50 repeats, including five large expansions (∼1000-2000 repeats) that contracted. Among 44 offspring who did not inherit the CTG18.1exp+ allele, eight (18.2%) exhibited FECD. Conclusions CTG18.1 expansion was highly associated with FECD but demonstrated incomplete penetrance. CTG18.1L instability occurred in a minority of parent-offspring transmissions, with large expansions exhibiting contraction. The observation of FECD without CTG18.1 expansion among family members in CTG18.1exp+ families highlights the complexity of the relationship between the FECD phenotype and CTG18.1 expansion.

Background Bone allotransplant viability can be maintained long-term by implanting arteriovenous (AV) bundles and creating an autogenous neoangiogenic circulation. Only short-term immunosuppression is required. This study investigates the origin of viable osteocytes observed in areas of active bone remodeling in orthotopically transplanted tibiae in a Yucatan mini-pig model. Methods Segmental tibial defects created in female Yucatan minipigs (N = 14) were reconstructed with a matched vascularized composite allotransplant from a male donor. The circulation was microsurgically restored, with simultaneous autogenous AV-bundle implantation in group 1 (N = 7). A ligated AV-bundle was implanted as a no-angiogenesis control in group 2 (N = 7). After 20-weeks, repopulation of the allotransplant was assessed by real-time qPCR measurement of relative copy numbers of a Y chromosome-specific gene (SRY) and an autosomal housekeeping gene, ribosomal protein L4 (RPL4). A lower SRY/RPL4 ratio demonstrates replacement of male allogeneic cells with female, autogenous cells in the sample. Genomic DNA was extracted from cross-sections of the allotransplant, liver and spleen. Additionally, areas of new bone formation within the allotransplant were sampled by laser capture microdissection. A comparison was made between groups as well as male control samples. RNA was extracted from bone as well, as a measure of metabolically active cells. Results Laser-captured areas of new bone formation in animals with both normal and ligated AV-bundles were found to have significantly lower relative copy numbers of SRY (p = 0.03) than control specimens from male bone, indicating replacement by female (autogenous) bone-forming cells. Analysis of an entire segment of the allotransplant from Group 1 was similarly reduced (p = 0.04), unlike that from Group 2. RNA expression of SRY was observed in both groups. No chimerism could be found in non-bone tissues (liver and spleen). Conclusion We observed a significant level of transplant chimerism in areas of new bone formation sampled by laser capture microdissection. The migration of autogenous cells including osteocytes was seen in both groups. Survival of some allogeneic (male) cells was also demonstrable. No microchimerism was found in liver and spleen.

We report a female infant identified by newborn screening for severe combined immunodeficiencies (NBS SCID) with T cell lymphopenia (TCL). The patient had persistently elevated alpha-fetoprotein (AFP) with IgA deficiency, and elevated IgM. Gene sequencing for a SCID panel was uninformative. We sought to determine the cause of the immunodeficiency in this infant. We performed whole-exome sequencing (WES) on the patient and parents to identify a genetic diagnosis. Based on the WES result, we developed a novel flow cytometric panel for rapid assessment of DNA repair defects using blood samples. We also performed whole transcriptome sequencing (WTS) on fibroblast RNA from the patient and father for abnormal transcript analysis. WES revealed a pathogenic paternally inherited indel in ATM. We used the flow panel to assess several proteins in the DNA repair pathway in lymphocyte subsets. The patient had absent phosphorylation of ATM, resulting in absent or aberrant phosphorylation of downstream proteins, including γH2AX. However, ataxia-telangiectasia (AT) is an autosomal recessive condition, and the abnormal functional data did not correspond with a single ATM variant. WTS revealed in-frame reciprocal fusion transcripts involving ATM and SLC35F2 indicating a chromosome 11 inversion within 11q22.3, of maternal origin. Inversion breakpoints were identified within ATM intron 16 and SLC35F2 intron 7. We identified a novel ATM-breaking chromosome 11 inversion in trans with a pathogenic indel (compound heterozygote) resulting in non-functional ATM protein, consistent with a diagnosis of AT. Utilization of several molecular and functional assays allowed successful resolution of this case.

Background Genomic aberrations have been identified in metastatic castration-resistant prostate cancer (mCRPC), but molecular predictors of resistance to abiraterone acetate/prednisone (AA/P) treatment are not known. Patients and methods In a prospective clinical trial, mCRPC patients underwent whole-exome sequencing (n = 82) and RNA sequencing (n = 75) of metastatic biopsies before initiating AA/P with the objective of identifying genomic alterations associated with resistance to AA/P. Primary resistance was determined at 12 weeks of treatment using criteria for progression that included serum prostate-specific antigen measurement, bone and computerized tomography imaging and symptom assessments. Acquired resistance was determined using the end point of time to treatment change (TTTC), defined as time from enrollment until change in treatment from progressive disease. Associations of genomic and transcriptomic alterations with primary resistance were determined using logistic regression, Fisher‘s exact test, single and multivariate analyses. Cox regression models were utilized for determining association of genomic and transcriptomic alterations with TTTC. Results At 12 weeks, 32 patients in the cohort had progressed (nonresponders). Median study follow-up was 32.1 months by which time 58 patients had switched treatments due to progression. Median TTTC was 10.1 months (interquartile range: 4.4–24.1). Genes in the Wnt/s-catenin pathway were more frequently mutated and negative regulators of Wnt/s-catenin signaling were more frequently deleted or displayed reduced mRNA expression in nonresponders. Additionally, mRNA expression of cell cycle regulatory genes was increased in nonresponders. In multivariate models, increased cell cycle proliferation scores (≥50) were associated with shorter TTTC (hazard ratio = 2.11, 95% confidence interval: 1.17–3.80; P = 0.01). Conclusions Wnt/s-catenin pathway activation and increased cell cycle progression scores can serve as molecular markers for predicting resistance to AA/P therapy.

The εγδβ thalassemias (EGDBT) are infrequently occurring deletions of the β globin gene cluster on chromosome 11. EGDBTs can be an unsuspected cause of severe neonatal anemia that resolves during the first year of life into a normal Hb A2 thalassemic phenotype, presumed to be due to an imbalance of β-like chains during development. First described in 1972, they are categorized into two groups, both with loss of the five DNase I hypersensitivity sites upstream of the β gene (HBB) called the locus control region (LCR). In group 1 deletions, all or most of the β globin gene cluster including HBB is deleted; in group 2 deletions, removal of the LCR silences an intact HBB. Large deletions that do not involve LCR but remove the embryonic (HBE) and fetal genes (HBG2/HBG1) [LCR(ε Gγ Aγ) 0δβ] have not been reported. It is not known whether a deletion with intact LCR as well as δ and β genes would similarly manifest as neonatal anemia that resolves into a normocytic phenotype with normal hemoglobin fraction percentages after infancy. Four individuals with an LCR(ε Gγ Aγ) 0δβ deletion have been evaluated. These include monochorionic, diamniotic twin girls whose newborn screen (NBS) returned significant for Hb A percentage greater than Hb F for both twins. NBS was repeated two weeks later with the same result. The newborns were induced late pre-term (36 weeks and 6 days) after an uncomplicated pregnancy, including no features of twin to twin transfusion. Both infants had no anemia and bilirubin levels remained in the low to low-intermediate risk zones. Neither twin received transfusions or phototherapy during admission. At a follow up visit at 20 months of age a hemoglobin electrophoresis was performed which showed borderline Hb A2 and no Hb F in both twins (Table 1). Ferritin levels were low-normal in Twin A (14 µg/L) and low in Twin B (9 µg/L). Both infants were found to have a heterozygous large deletion affecting HBE, HBG2, HBG1, and HBBP1 loci by Multiplex Ligation-dependent Probe Amplification (MLPA), a method that uses multiple DNA probe pairs to determine the size of the DNA segment deleted. MLPA deletion/duplication testing on the α genes was negative. No β globin variants were identified by DNA sequencing. An evaluation was performed on the 25 year-old healthy mother and showed similar findings including the heterozygous large deletion. Mother and infants continue to be asymptomatic with normal red cell indices for age (Table 1). The fourth was an unrelated 40-year-old female with mild anemia screened by hemoglobin electrophoresis for obstetric evaluation who was found to have borderline Hb A2 (see table 1). MLPA showed a heterozygous deletion involving the same loci. Again, no β globin variants were identified by DNA sequencing. All four cases showed the same MLPA pattern except one non-specific probe commonly altered in common gamma gene conversion events. Long range sequencing performed on one case confirmed a single contiguous 32,599 bp deletion that matched the MLPA data (g.5262276-5294875). The detected deletion has not been reported in literature and differs from EGDBT mutations in that the LCR is intact. As the LCR controls expression of the entire gene cluster, this mutation is expected to display different phenotypic features than classic EGDBT. HBD and HBB are not deleted and therefore adequate transcription of Hb A2 and Hb A are expected with no imbalance of α/β chains after the neonatal period. It is postulated this finding explains the decreased Hb F levels in the NBS test result. The Hb A2 increase is likely a product of this deletion and does not indicate β-thalassemia. The heterozygous absence of HBE and HBG2/HBG1 may result in upregulation of HBD and HBB expression through loss of FKLF and FKLF-2 binding. In addition, HBB-EKLF binding would have less competition as the gamma promoter elements have been heterozygously deleted. In summary, in contrast to classic EGDBTs which cause transient severe neonatal anemia and resolve to normal Hb A2 thalassemic indices, this novel LCR(ε Gγ Aγ) 0δβ deletion was associated with normal CBC values with an absence of severe neonatal anemia, inverted Hb F/Hb A percentages at birth and borderline Hb A2 levels. Figure 1 Figure 1. Disclosures No relevant conflicts of interest to declare.

Primary cultures of human corneal endothelial cells (HCECs) are an important model system for studying the pathophysiology of corneal endothelium. The purpose of this study was to identify and validate an optimal primary culture model of normal and Fuchs endothelial corneal dystrophy (FECD) endothelial cells by comparing cell morphology and marker expression under different media conditions to in vivo donor tissues. Primary and immortalized HCECs, isolated from normal and FECD donors, were cultured in proliferation media (Joyce, M4, Bartakova) alone or sequentially with maturation media (F99, Stabilization 1, M5). CD56, CD73 and CD166 expressions were quantified in confluent and matured cell lines by flow cytometry. HCECs that were allowed to proliferate in Joyce’s medium followed by maturation in low-mitogen containing media yielded cells with similar morphology to corneal endothelial tissues. Elevated expression of CD56 and CD166 and low expression of CD73 correlated with regular, hexagonal-like HCEC morphology. CD56:CD73 > 2.5 was most consistent with normal HCEC morphology and mimicked corneal endothelial tissue. Immortalization of normal HCECs by hTERT transduction showed morphology and CD56:CD73 ratios similar to parental cell lines. HCECs established from FECD donors showed reduced CD56:CD73 ratios compared to normal HCECs which coincided with reduced uniformity and regularity of cell monolayers. Overall, a dual media system with Joyce’s medium for proliferation and a low-mitogen media for maturation, provided normal cultures with regular, hexagonal-like cell morphologies consistent with corneal endothelial cells in vivo. CD56:CD73 expression ratio >2.5 was predictive of in vivo-like cellular morphology.

Background Congenital myopathies due to ryanodine receptor (RYR1) mutations are increasingly identified and correlate with a wide range of phenotypes, most commonly that of malignant hyperthermia susceptibility and central cores on muscle biopsy with rare reports of distal muscle weakness, but in the setting of early onset global weakness. Methods We report a case of a patient presenting with childhood onset hand stiffness and adult onset progressive hand weakness and jaw contractures discovered to have two variants in the RYR1 gene. Results The patient manifested with distal upper limb weakness which progressed to involve the distal lower limb, proximal upper limb, as well as the face in addition to limited jaw opening. Creatine kinase was mildly elevated with EMG findings supporting a myopathy. Muscle biopsy showed features consistent with centronuclear myopathy. Whole exome sequencing revealed a novel heterozygous pathogenic variant in RYR1 (c.12315_12328delAGAAATCCAGTTCC, p.Glu4106Alafs*8), and a heterozygous missense variant (c.10648C>T, p.Arg3550Trp) of unknown significance in compound heterozygous state. Conclusion We expand the spectrum of RYR1-related myopathy with the description of a novel phenotype in an adult patient presenting with hand weakness and suggest considering RYR1 analysis in the diagnosis of distal myopathies.

// Xueju Wang 1, 2 , Surendra Dasari 3 , Grzegorz S. Nowakowski 4 , Konstantinos N. Lazaridis 5, 6 , Eric D. Wieben 5, 7 , Marshall E. Kadin 8 , Andrew L. Feldman 1, * , Rebecca L. Boddicker 1, * 1 Department of Laboratory Medicine and Pathology, Mayo Clinic, Rochester, Minnesota, United States of America 2 Department of Pathology, China-Japan Union Hospital of Jilin Province, Changchun, Jilin Province, China 3 Division of Biomedical Statistics and Informatics, Department of Health Sciences Research, Mayo Clinic, Rochester, Minnesota, United States of America 4 Division of Hematology, Department of Medicine, Mayo Clinic, Rochester, Minnesota, United States of America 5 Center for Individualized Medicine, Mayo Clinic, Rochester, Minnesota, United States of America 6 Division of Gastroenterology and Hepatology, Department of Medicine, Mayo Clinic, Rochester, Minnesota, United States of America 7 Department of Biochemistry and Molecular Biology, Mayo Clinic, Rochester, Minnesota, United States of America 8 Department of Pathology and Laboratory Medicine, Rhode Island Hospital and Department of Dermatology, Roger Williams Medical Center, Providence, Rhode Island, United States of America * These authors have contributed equally to this work Correspondence to: Andrew L. Feldman, email: feldman.andrew@mayo.edu Rebecca L. Boddicker, email: boddicker.rebecca@mayo.edu Keywords: T-cell lymphoma, retinoids, retinoic acid receptor alpha, all-trans retinoic acid, cell cycle, individualized medicine Received: October 04, 2016 Accepted: February 06, 2017 Published: February 17, 2017 ABSTRACT Peripheral T-cell lymphomas (PTCLs) are aggressive non-Hodgkin lymphomas with generally poor outcomes following standard therapy. Few candidate therapeutic targets have been identified to date. Retinoic acid receptor alpha (RARA) is a transcription factor that modulates cell growth and differentiation in response to retinoids. While retinoids have been used to treat some cutaneous T-cell lymphomas (CTCLs), their mechanism of action and the role of RARA in CTCL and other mature T-cell lymphomas remain poorly understood. After identifying a PTCL with a RARA R394Q mutation, we sought to characterize the role of RARA in T-cell lymphoma cells. Overexpressing wild-type RARA or RARA R394Q significantly increased cell growth in RARA low cell lines, while RARA knockdown induced G1 arrest and decreased expression of cyclin-dependent kinases CDK2/4/6 in RARA high cells. The retinoids, AM80 (tamibarotene) and all- trans retinoic acid, caused dose-dependent growth inhibition, G 1 arrest, and CDK2/4/6 down-regulation. Genes down-regulated in transcriptome data were enriched for cell cycle and G1-S transition. Finally, RARA overexpression augmented chemosensitivity to retinoids. In conclusion, RARA drives cyclin-dependent kinase expression, G 1 -S transition, and cell growth in T-cell lymphoma. Synthetic retinoids inhibit these functions in a dose-dependent fashion and are most effective in cells with high RARA expression, indicating RARA may represent a therapeutic target in some PTCLs.

Purpose To identify RNA missplicing events in human corneal endothelial tissue isolated from Fuchs' endothelial corneal dystrophy (FECD). Methods Total RNA was isolated and sequenced from corneal endothelial tissue obtained during keratoplasty from 12 patients with FECD and 4 patients undergoing keratoplasty or enucleation for other indications. The length of the trinucleotide repeat (TNR) CTG in the transcription factor 4 (TCF4) gene was determined using leukocyte-derived DNA analyzed by a combination of Southern blotting and Genescan analysis. Commercial statistical software was used to quantify expression of alternatively spliced genes. Validation of selected alternative splicing events was performed by using RT-PCR. Gene sets identified were analyzed for overrepresentation using Web-based analysis system. Results Corneal endothelial tissue from FECD patients containing a CTG TNR expansion sequence in the TCF4 gene revealed widespread changes in mRNA splicing, including a novel splicing event involving FGFR2. Differential splicing of NUMA1, PPFIBP1, MBNL1, and MBNL2 transcripts were identified in all FECD samples containing a TNR expansion. The differentially spliced genes were enriched for products that localize to the cell cortex and bind cytoskeletal and cell adhesion proteins. Conclusions Corneal endothelium from FECD patients harbors a unique signature of mis-splicing events due to CTG TNR expansion in the TCF4 gene, consistent with the hypothesis that RNA toxicity contributes to the pathogenesis of FECD. Changes to the endothelial barrier function, a known event in the development of FECD, was identified as a key biological process influenced by the missplicing events.

Following publication of the original article [1], the author explained that Table 2 is displayed incorrectly. The correct Table 2 is given below. The original article has been corrected.

PURPOSE: To investigate single nucleotide polymorphisms (SNPs) and trinucleotide repeat (TNR) expansion in the transcription factor 4 (TCF4) gene in a large cohort of German patients with Fuchs endothelial corneal dystrophy (FECD). METHODS: Genomic DNA was obtained from 398 patients with FECD and from 58 non-FECD controls. Thirty-seven previously reported SNPs were evaluated by genotyping. The 398 FECD samples were analyzed for TNR expansions by short tandem repeat (STR) assays and Southern blotting. The possible associations between TNR length and clinical parameters (age, sex, visual acuity, and central corneal thickness) were analyzed in 132 patients. RESULT: The SNPs in COL8A2, TCF8, LOXHD1, and AGBL1 showed no heterogeneity in 36 cases, although SLCA411 showed 3 nonsense mutations. SNPs were detected for TCF4 (rs613872, rs2123392, rs17089887, rs1452787, and rs1348047), but only rs613872 showed a significant association with FECD (P=9.93×10(−12)). Overall, 315/398 (79%) patients harbored TNR lengths >50, while no non-FECD control harbored TNR lengths >50. The TCF4 SNP rs613872 genotype was TT:39 (67%), TG:18 (31%), and GG:1 (2%) in non-FECD controls; TT:39 (47%), TG: 38 (46%), and GG:6 (7%) in FECD cases harboring TNR 50 (P=2.93×10(−25)). No significant association was detected between clinical parameters and TNR length. CONCLUSIONS: Our large German cohort demonstrated a significant association between the risk allele G in rs613872 and FECD, irrespective of TNR expansion, although this risk allele was more frequent in FECD cases with TNR expansion than without.

Amplification of a CAG trinucleotide motif (CTG18.1) within the TCF4 gene has been strongly associated with Fuchs Endothelial Corneal Dystrophy (FECD). Nevertheless, a small minority of clinically unaffected elderly patients who have expanded CTG18.1 sequences have been identified. To test the hypothesis that the CAG expansions in these patients are protected from FECD because they have interruptions within the CAG repeats, we utilized a combination of an amplification-free, long-read sequencing method and a new target-enrichment sequence analysis tool developed by Pacific Biosciences to interrogate the sequence structure of expanded repeats. The sequencing was successful in identifying a previously described interruption within an unexpanded allele and provided sequence data on expanded alleles greater than 2000 bases in length. The data revealed considerable heterogeneity in the size distribution of expanded repeats within each patient. Detailed analysis of the long sequence reads did not reveal any instances of interruptions to the expanded CAG repeats, but did reveal novel variants within the AGG repeats that flank the CAG repeats in two of the five samples from clinically unaffected patients with expansions. This first examination of the sequence structure of CAG repeats in CTG18.1 suggests that factors other than interruptions to the repeat structure account for the absence of disease in some elderly patients with repeat expansions in the TCF4 gene.

Background Use of the Genome Analysis Toolkit (GATK) continues to be the standard practice in genomic variant calling in both research and the clinic. Recently the toolkit has been rapidly evolving. Significant computational performance improvements have been introduced in GATK3.8 through collaboration with Intel in 2017. The first release of GATK4 in early 2018 revealed rewrites in the code base, as the stepping stone toward a Spark implementation. As the software continues to be a moving target for optimal deployment in highly productive environments, we present a detailed analysis of these improvements, to help the community stay abreast with changes in performance. Results We re-evaluated multiple options, such as threading, parallel garbage collection, I/O options and data-level parallelization. Additionally, we considered the trade-offs of using GATK3.8 and GATK4. We found optimized parameter values that reduce the time of executing the best practices variant calling procedure by 29.3% for GATK3.8 and 16.9% for GATK4. Further speedups can be accomplished by splitting data for parallel analysis, resulting in run time of only a few hours on whole human genome sequenced to the depth of 20X, for both versions of GATK. Nonetheless, GATK4 is already much more cost-effective than GATK3.8. Thanks to significant rewrites of the algorithms, the same analysis can be run largely in a single-threaded fashion, allowing users to process multiple samples on the same CPU. Conclusions In time-sensitive situations, when a patient has a critical or rapidly developing condition, it is useful to minimize the time to process a single sample. In such cases we recommend using GATK3.8 by splitting the sample into chunks and computing across multiple nodes. The resultant walltime will be nnn.4 hours at the cost of $41.60 on 4 c5.18xlarge instances of Amazon Cloud. For cost-effectiveness of routine analyses or for large population studies, it is useful to maximize the number of samples processed per unit time. Thus we recommend GATK4, running multiple samples on one node. The total walltime will be ∼34.1 hours on 40 samples, with 1.18 samples processed per hour at the cost of $2.60 per sample on c5.18xlarge instance of Amazon Cloud.

Fuchs endothelial corneal dystrophy (FECD) is a common cause for heritable visual loss in the elderly. Since the first description of an association between FECD and common polymorphisms situated within the transcription factor 4 (TCF4) gene, genetic and molecular studies have implicated an intronic CTG trinucleotide repeat (CTG18.1) expansion as a causal variant in the majority of FECD patients. To date, several non-mutually exclusive mechanisms have been proposed that drive and/or exacerbate the onset of disease. These mechanisms include (i) TCF4 dysregulation; (ii) toxic gain-of-function from TCF4 repeat-containing RNA; (iii) toxic gain-of-function from repeat-associated non-AUG dependent (RAN) translation; and (iv) somatic instability of CTG18.1. However, the relative contribution of these proposed mechanisms in disease pathogenesis is currently unknown. In this review, we summarise research implicating the repeat expansion in disease pathogenesis, define the phenotype-genotype correlations between FECD and CTG18.1 expansion, and provide an update on research tools that are available to study FECD as a trinucleotide repeat expansion disease. Furthermore, ongoing international research efforts to develop novel CTG18.1 expansion-mediated FECD therapeutics are highlighted and we provide a forward-thinking perspective on key unanswered questions that remain in the field., Highlights • FECD is a common, age-related corneal dystrophy. • The majority of cases are associated with expansion of a CTG repeat (CTG18.1). • FECD is the most common trinucleotide repeat expansion disease in humans. • Evidence supports multiple molecular mechanisms underlying the pathophysiology. • Novel CTG18.1-targeted therapeutics are in development.

Objective To describe the experience and outcome of performing whole-exome sequencing (WES) for resolution of patients on a diagnostic odyssey in the first 18 months of an individualized medicine clinic (IMC). Patients and Methods The IMC offered WES to physicians of Mayo Clinic practice for patients with suspected genetic disease. DNA specimens of the proband and relatives were submitted to WES laboratories. We developed the Genomic Odyssey Board with multidisciplinary expertise to determine the appropriateness for IMC services, review WES reports, and make the final decision about whether the exome findings explain the disease. This study took place from September 30, 2012, to March 30, 2014. Results In the first 18 consecutive months, the IMC received 82 consultation requests for patients on a diagnostic odyssey. The Genomic Odyssey Board deferred 7 cases and approved 75 cases to proceed with WES. Seventy-one patients met with an IMC genomic counselor. Fifty-one patients submitted specimens for WES testing, and the results have been received for all. There were 15 cases in which a diagnosis was made on the basis of WES findings; thus, the positive diagnostic yield of this practice was 29%. The mean cost per patient for this service was approximately $8000. Medicaid supported 27% of the patients, and 38% of patients received complete or partial insurance coverage. Conclusion The significant diagnostic yield, moderate cost, and notable health marketplace acceptance for WES compared with conventional genetic testing make the former method a rational diagnostic approach for patients on a diagnostic odyssey.

Background: Both HR deficiency and BRCA mutation status predict response to platinum-based therapy and BRCA mutation status predicts docetaxel resistance. However, the association of either biomarker with response to the individual elements of either AC or taxanes (T) is unknown since T is commonly given concomitantly with or after anthracyclines (A). We evaluated the association of HRD and BRCA mutation status with response to neoadjuvant weekly T followed by AC or (F)EC in high-risk breast cancer. Methods: We studied 140 high risk Stage I-III breast cancer patients (pts), enrolled in the breast cancer genome guided therapy study (BEAUTY), obtaining biopsies for DNA/RNA sequencing and MRI imaging to assess response to neoadjuvant weekly T (+trastuzumab+/-pertuzumab for HER2+ disease) followed by AC or (F)EC. Germline BRCA status and HR status of tumor samples (Myriad laboratories) were obtained. HR deficient tumor was defined as HRD score ≥42 or BRCA mutation. MRI response by changes in tumor size after 12 weeks of T was classified by WHO criteria. pCR was defined as ypT0/Tis ypN0. Both MRI response after T and pCR (after T and AC) were examined in terms of germline BRCA mutation (gBRCAmut vs. gBRCAwt) and tumor HR deficiency. Results: Of 140 pts enrolled, 8 withdrew consent and 2 carboplatin treated pts were excluded. Germline data were available for 124/130 pts. 12 patients had BRCA deleterious germline mutations (4 BRCA1, 8 BRCA2). MRI partial (PR)/complete response (CR) rate to T was 47.3% (95% CI: 37.8-57.0%) in the BRCAwt group and 66.7% (95% CI: 34.9-90.1%) in the BRCAmut group. No MRI CR's were observed in BRCA1 mut pts. In contrast, pCR rate was 50% in the 12 gBRCAmut pts (95% CI: 21.1-78.9%) and 31.3% in the 112 gBRCAwt pts (95% CI: 22.8-40.7%). HR deficiency status has thus far been determined for 74 pts: 26 pts have HD deficient tumors: 18 TNBC, 5 Luminal B, 2 ER-/HER2+; and 1 ER+/HER2+. Determination of HR deficiency is ongoing and will be reported for the full cohort in terms of 12 week MRI response to T and pCR to T+AC. HR deficientMolecular Subtypeyes (%)no (%)TBD (%)Luminal A0/112/11 (18.2)9/11 (81.8)Luminal B5/37 (13.5)13/37 (35.1)19/37 (51.3)Luminal NOS0/21/2 (50)1/2 (50)ER+/Her2+1/17 (5.8)14/17 (82.4)2/17 (11.8)ER-/Her2+2/20 (10)11/20 (55)7/20 (35)Triple Negative18/43 (41.9)6/43 (18.6)17/43 (39.5)germline BRCA statusMRI partial response after T (%)MRI complete response after T (%)pCR after T&AC (%)BRCA11/4 (25)0/42/4 (50)BRCA25/8 (62.5)2/8 (25)4/8 (50)BRCAwt35/112 (31.3)18/112 (16.1)35/112 (31.3) Conclusion: In the setting of neoadjuvant weekly T followed by AC, pCR rates were non-significantly higher in pts with BRCA1 mutations. While we observed no overall association between BRCA mutation status and response rates to taxanes; nearly all MRI responses to taxanes (partial and complete) were observed in the BRCA2 group. Prospective studies are needed to validate these findings and to determine whether BRCA status can be used to select therapy. HR deficiency is uncommon in luminal A and HER2+, frequent in TNBC, and the association of HRD with both MRI response to taxanes and pCR will be reported at the meeting. Citation Format: Boughey JC, Kalari KR, Suman VJ, McLaughlin SA, Moreno Aspitia A, Moyer AM, Northfelt DW, Gray RJ, Vedell PT, Tang X, Dockter TJ, Jones KN, Felten SJ, Conners AL, Hart SN, Visscher DW, Wieben ED, Ingle JN, Hartman A-R, Timms K, Elkin E, Jones J, Wang L, Weinshilboum RW, Goetz MP. Role of germline BRCA status and tumor homologous recombination (HR) deficiency in response to neoadjuvant weekly paclitaxel followed by anthracycline-based chemotherapy. [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 P3-07-29.

Background: Tumorigenesis involves both genetic and epigenetic changes. Epigenetic alterations are reversible and are promising cancer therapeutic targets. Decitabine (5-aza-2'-deoxycytidine), a DNA methyltransferase inhibitor, is FDA approved for hematological malignancies. However, the effect of decitabine in breast cancer is not completely understood. Previous reports indicated that one decitabine mechanism involves regulation of protein levels for DNMT1, the major DNA methyltransferase that methylates hemimethylated CpG di-nucleotides in DNA. However, the E3 ligase involved in this process has not been identified. Whether decitabine also regulates DNMT3A and 3B in a similar fashion remains unclear. Therefore, our goals were to 1) understand mechanisms underlying decitabine action, 2) test the antitumor activity of decitabine in breast cancer models and 3) identify biomarkers associated with response to decitabine. Methods and Results: Western blots of breast cancer cell lines showed that DNMT1, DNMT3A, and DNMT3B protein levels decreased following decitabine treatment without a reduction in mRNA levels. Bioinformatic analysis of DNA methyltransferase sequences revealed a potential TRAF6 binding motif, and the interaction with TRAF6 (TNF receptor-associated factor 6) was confirmed by IP. TRAF6 functions as an E3 ligase. To determine whether TRAF6 might be the E3 ligase responsible for the degradation of DNMTs after decitabine treatment, we knocked down TRAF6 by RNA interference or knocked out the TRAF6 gene by CRISPR/Cas9. Down regulation of TRAF6 attenuated DNMT ubiquitination and increased DNMT protein levels, suggesting that TRAF6 might mediate proteasome-dependent degradation of all three DNMTs. This was further confirmed by reconstituting the knockout cells with WT and a TRAF6-C70A mutant, followed by assessing DNMT protein levels. Global DNA methylation was also increased after TRAF6 depletion and was confirmed in TRAF6 knock out cells in which DNMT levels were unaffected by decitabine. Cell cytotoxicity and colony forming assays showed that TRAF6 knockout cells were resistant to decitabine, suggesting that a major decitabine mechanism of action is through the regulation of TRAF6 which, in turn, degrades DNMTs, leading to decreased global methylation. Finally, decitabine significantly induced TRAF6 at both mRNA and protein levels, a process that might create positive feedback leading to increased degradation of DNMT proteins upon decitabine treatment. Based on these results, we further hypothesized that levels of the three DNMTs might influence decitabine response. Using 18 breast cancer patient derived xenograft (PDX) models, we found a wide range of DNMT protein levels regardless of ER/HER2 status. DNMT levels in the PDX models were directly associated with sensitivity to decitabine treatment, confirming our hypothesis. Conclusion: Our data showed that decitabine might be an effective agent for treating breast cancer and revealed a novel mechanism underlying decitabine treatment. Baseline DNMT protein levels may serve as a biomarker for predicting decitabine drug response. Citation Format: Yu J, Qin B, Boughey JC, Moyer AM, Visscher DW, Sinnwell JP, Yin P, Thompson KJ, Docter TJ, Kalari KR, Suman VJ, Wieben ED, Felten SJ, Conners AL, Jones KN, McLaughlin SA, Copland JA III, Moreno Aspitia A, Northfelt DW, Gray RJ, Ingle JN, Lou Z, Weinshilboum R, Goetz MP, Wang L. Regulation of DNA methyltransferases via TRAF6 determines breast cancer response to decitabine. [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 P3-07-51.

Background A 4-year-old boy born at 37 weeks' gestation with intrauterine growth retardation presented with developmental delay with pronounced language and gross motor delay, axial hypotonia, and dynamic hypertonia of the extremities. Investigations including the Minnesota Newborn Screen, thyroid stimulating hormone/thyroxin, and inborn errors of metabolism screening were negative. Cerebral magnetic resonance imaging and spectroscopy were normal. Genetic testing was negative for coagulopathy, Smith-Lemli-Opitz, fragile X, and Prader-Willi/Angelman syndromes. Whole genome array analysis was unremarkable. Methods Whole exome sequencing was performed through a commercial testing laboratory to elucidate the underlying etiology for the child's presentation. A de novo mutation was hypothesized. In attempt to establish pathogenicity of our candidate variant, cellular electrophysiologic functional analysis of the putative de novo mutation was performed using patch-clamp technology. Results Whole exome sequencing revealed a p.P1353L variant in the CACNA1A gene, which encodes for the α1-subunit of the brain-specific P/Q-type calcium channel (Ca V 2.1). This presynaptic high-voltage-gated channel couples neuronal excitation to the vesicular release of neurotransmitter and is implicated in several neurologic disorders. DNA Sanger sequencing confirmed that the de novo mutation was absent in both parents and present in the child only. Electrophysiologic analysis of P1353L-CACNA1A demonstrated near complete loss of function, with a 95% reduction in peak current density. Conclusions Whole exome sequencing coupled with cellular electrophysiologic functional analysis of a de novo CACNA1A missense mutation has elucidated the probable underlying pathophysiologic mechanism responsible for the child's phenotype. Genetic testing of CACNA1A in patients with congenital hypotonia and developmental delay may be warranted.

As reliable, efficient genome sequencing becomes more ubiquitous, the need for similarly reliable and efficient variant calling becomes increasingly important. The Genome Analysis Toolkit (GATK), maintained by the Broad Institute, is currently the widely accepted standard for variant calling software. However, alternative solutions may provide faster variant calling without sacrificing accuracy. One such alternative is Sentieon DNASeq, a toolkit analogous to GATK but built on a highly optimized backend. We evaluated the DNASeq single-sample variant calling pipeline in comparison to that of GATK. Our results confirm the near-identical accuracy of the two software packages, showcase perfect scalability and great speed from Sentieon, and describe computational performance considerations for the deployment of Sentieon DNASeq.

Introduction: Globin gene sequencing has long been impeded by repetitive regions, homology and structural variation. Currently, long-read sequencing technology has been limited by high cost and error rate. These challenges can be mitigated by sample multiplexing of targeted capture sequencing and the use of consensus reads, respectively. This ongoing study analyzes the utility of long range sequencing technology in identifying historically problematic beta globin cluster mutations. Materials and methods: Long range sequencing was performed on complex beta globin gene cluster cases utilizing the Sequel platform (Pacific Biosciences, Menlo Park, California, United States). Sample input was 1 microgram of DNA extracted from whole blood. Circular Consensus Sequence (CCS) FASTQ sequences files were mapped to reference genome (GRCh37/hg19) using NGMLR, followed by annotation of variants using a custom bioinformatics workflow. Comparison to multiplex ligation-dependent probe amplification (MLPA), Sanger sequencing, array CGH, short-range NGS (300 bp), hemoglobin protein, and clinical data was performed. Results: CCS read length varied from 1 to 10 kb, average 3.5. Analysis of the beta globin cluster region showed impressive superiority and comprehensiveness to Sanger, MLPA, array CGH and short range sequencing technology: 1) single nucleotide variants (SNVs) were identified to a sensitivity and specificity of 99.5%; 2) large structural variants (SVs) such as large deletions, duplications, insertions, crossovers and fusions spanning many kilobases were characterized in the heterozygous, homozygous and compound heterozygous state to a precise genomic coordinate (table 1); 3) phasing SNVs identified Hb S haplotypes (Central African region (CAR), Benin, Senegal, Arab-Indian and Cameroon). Conclusion: 1) Long range sequencing holds significant promise in genotyping hemoglobin disorders. This method improves efficiency with potential for a single test offering with increased resolution, particularly in large deletion analysis. Complex cases greatly benefited from improved range and resolution. 2) The decreased requirement of blood volume for extensive testing is of particular interest for the pediatric population. 3) Genotyping of Hb S patients detected by positive newborn screening will help predict need for prophylaxis or therapy by assessing variant zygosity, phasing and haplotypes. In conclusion, we believe that long range sequencing has immense potential to contribute extensively to thalassemia and hemoglobinopathy diagnostics and is far superior to the currently available technologies for beta globin cluster analysis. Disclosures Jen: Celgene: Employment.

Background: Pathogenic germline mutations in the CDKN2A tumor suppressor gene are rare and associated with highly penetrant familial melanoma and pancreatic cancer in non-Hispanic whites (NHW). To date, the prevalence and impact of CDKN2A rare coding variants (RCV) in racial minority groups remain poorly characterized. We examined the role of CDKN2A RCVs on the risk of pancreatic cancer among minority subjects. Methods: We sequenced CDKN2A in 220 African American (AA) pancreatic cancer cases, 900 noncancer AA controls, and 183 Nigerian controls. RCV frequencies were determined for each group and compared with that of 1,537 NHW patients with pancreatic cancer. Odds ratios (OR) and 95% confidence intervals (CI) were calculated for both a case–case comparison of RCV frequencies in AAs versus NHWs, and case–control comparison between AA cases versus noncancer AA controls plus Nigerian controls. Smaller sets of Hispanic and Native American cases and controls also were sequenced. Results: One novel missense RCV and one novel frameshift RCV were found among AA patients: 400G>A and 258_278del. RCV carrier status was associated with increased risk of pancreatic cancer among AA cases (11/220; OR, 3.3; 95% CI, 1.5–7.1; P = 0.004) compared with AA and Nigerian controls (17/1,083). Further, AA cases had higher frequency of RCVs: 5.0% (OR, 13.4; 95% CI, 4.9–36.7; P < 0.001) compared with NHW cases (0.4%). Conclusions: CDKN2A RCVs are more common in AA than in NHW patients with pancreatic cancer and associated with moderately increased pancreatic cancer risk among AAs. Impact: RCVs in CDKN2A are frequent in AAs and are associated with risk for pancreatic cancer. Cancer Epidemiol Biomarkers Prev; 27(11); 1364–70. ©2018 AACR.

Use of the Genome Analysis Toolkit (GATK) continues to be the standard practice in genomic variant calling in both research and the clinic. Recently the toolkit has been rapidly evolving. Significant computational performance improvements have been introduced in GATK3.8 through collaboration with Intel in 2017. The first release of GATK4 in early 2018 revealed significant rewrites in the code base, as the stepping stone toward a Spark implementation. As the software continues to be a moving target for optimal deployment in highly productive environments, we present a detailed analysis of these improvements, to help the community stay abreast with changes in performance. We re-evaluated the options previously identified as advantageous, such as threading, parallel garbage collection, I/O options and data-level parallelization. Based on our results, we consider the performance and cost trade-offs of using GATK3.8 and GATK4 for different types of analyses.

Objective We identified a pedigree over five generations with 49 members, some of whom had chronic megacolon presenting in adolescence or adulthood. We aimed to assess the genetic cause of chronic megacolon through clinical and DNA studies. Design After ethical approval and informed consent, family members provided answers to standard bowel disease questionnaires, radiological or surgical records, and DNA (buccal mucosal scraping). Exome DNA sequencing of colon tissue or blood DNA from seven family members with colon or duodenal dilatation, or no megacolon (n = 1) was carried out. Sanger sequencing was performed in 22 additional family members to further evaluate candidate variants. The study focused on genes of potential relevance to enteric nerve (ENS) maturation and Hirschsprung's disease or megacolon, based on the literature (GFRA1, NKX2-1, KIF26A, TPM3, ACTG2, SCN10A, and C17orf107 [CHRNE]) and other genetic variants that co-segregated with megacolon in the six affected family members. Results Information was available in all except five members alive at time of study; among 30 members who provided DNA, six had definite megacolon, one megaduodenum, seven significant constipation without bowel dilatation, and 16 normal bowel function by questionnaire. Among genes studied, SEMA3F (g.3:50225360A>G; c1873A>G) was found in 6/6 family members with megacolon. The SEMA3F gene variant was assessed as potentially pathogenic, based on M-CAP in silico prediction. SEMA3F function is associated with genes (KIT and PDGFRB) that impact intestinal pacemaker function. Conclusion Familial chronic megacolon appears to be associated with SEMA3F, which is associated with genes impacting enteric nerve or pacemaker function.