Exome Capture Sequencing Research Articles

Genome-Wide Association Study in Bread Wheat Identifies Genomic Regions Associated with Grain Yield and Quality under Contrasting Water Availability.

The objectives of this study were to identify genetic loci in the bread wheat genome that would influence yield stability and quality under water stress, and to identify accessions that can be recommended for cultivation in dry and hot regions. We performed a genome-wide association study (GWAS) using a panel of 232 wheat accessions spanning diverse ecogeographic regions. Plants were evaluated in the Israeli Northern Negev, under two environments: water-limited (D; 250 mm) and well-watered (W; 450 mm) conditions; they were genotyped with ~71,500 SNPs derived from exome capture sequencing. Of the 14 phenotypic traits evaluated, 12 had significantly lower values under D compared to W conditions, while the values for two traits were higher under D. High heritability (H2 = 0.5–0.9) was observed for grain yield, spike weight, number of grains per spike, peduncle length, and plant height. Days to heading and grain yield could be partitioned based on accession origins. GWAS identified 154 marker-trait associations (MTAs) for yield and quality-related traits, 82 under D and 72 under W, and identified potential candidate genes. We identified 24 accessions showing high and/or stable yields under D conditions that can be recommended for cultivation in regions under the threat of global climate change.

Identification of the powdery mildew resistance gene in wheat breeding line Yannong 99102-06188 via bulked segregant exome capture sequencing.

Powdery mildew of wheat (Triticum aestivum), caused by Blumeria graminis f.sp. tritici (Bgt), is a destructive disease that seriously threatens the yield and quality of its host. Identifying resistance genes is the most attractive and effective strategy for developing disease-resistant cultivars and controlling this disease. In this study, a wheat breeding line Yannong 99102-06188 (YN99102), an elite derivative line from the same breeding process as the famous wheat cultivar Yannong 999, showed high resistance to powdery mildew at the whole growth stages. Genetic analysis was carried out using Bgt isolate E09 and a population of YN99102 crossed with a susceptible parent Jinhe 13–205 (JH13–205). The result indicated that a single recessive gene, tentatively designated pmYN99102, conferred seedling resistance to the Bgt isolate E09. Using bulked segregant exome capture sequencing (BSE-Seq), pmYN99102 was physically located to a ~33.7 Mb (691.0–724.7 Mb) interval on the chromosome arm 2BL, and this interval was further locked in a 1.5 cM genetic interval using molecular markers, which was aligned to a 9.0 Mb physical interval (699.2–708.2 Mb). Based on the analysis of physical location, origin, resistant spectrum, and inherited pattern, pmYN99102 differed from those of the reported powdery mildew (Pm) resistance genes on 2BL, suggesting pmYN99102 is most likely a new Pm gene/allele in the targeted interval. To transfer pmYN99102 to different genetic backgrounds using marker-assisted selection (MAS), 18 closely linked markers were tested for their availability in different genetic backgrounds for MAS, and all markers expect for YTU103-97 can be used in MAS for tracking pmYN99102 when it transferred into those susceptible cultivars.

Genetic Analysis of Adult Plant Resistance to Stripe Rust in Common Wheat Cultivar "Pascal".

Wheat stripe rust is an important foliar disease that affects the wheat yield globally. Breeding for resistant wheat varieties is one of the most economically and environmentally effective ways to control this disease. The common wheat (Triticum aestivum L.) cultivar “Pascal” exhibited susceptibility to stripe rust at the seedling stage but it showed high resistance to stripe rust at the adult plant stage over 20 years in Gansu, a hotspot of the disease in northwestern China. To understand the genetic mechanism of stripe rust resistance in this cultivar, a 55K SNP array was used to analyze the two parents and the 220 recombinant inbred lines (RILs) derived from the cross of “Huixianhong” × “Pascal.” We detected three new stripe rust adult plant resistance (APR) quantitative trait locus (QTL) contributed by Pascal, viz. QYr.gaas-1AL, QYr.gaas-3DL, and QYr.gaas-5AS, using the inclusive composite interval mapping method. They were flanked by SNP markers AX-111218361—AX-110577861, AX-111460455—AX-108798599, and AX-111523523—AX-110028503, respectively, and explained the phenotypic variation ranging from 11.0 to 23.1%. Bulked segregant exome capture sequencing (BSE-Seq) was used for fine mapping of QYr.gaas-1AL and selection of candidate genes, TraesCS1A02G313700, TraesCS1A02G313800, and TraesCS1A02G314900 for QYr.gaas-1AL. KASP markers BSE-1A-12 and HXPA-3D for QYr.gaas-1AL and QYr.gaas-3DL were developed for breeders to develop durable stripe rust-resistant wheat varieties.

A Multi-Omics Approach for Rapid Identification of Large Genomic Lesions at the Wheat Dense Spike (wds) Locus.

Optimal spike architecture provides a favorable structure for grain development and yield improvement. However, the number of genes cloned to underlie wheat spike architecture is extremely limited. Here, we obtained a wheat dense spike mutant (wds) induced by 60Co treatment of a common wheat landrace Huangfangzhu that exhibited significantly reduced spike and grain lengths. The shortened spike length was caused by longitudinal reduction in number and length of rachis cells. We adopted a multi-omics approach to identify the genomic locus underlying the wds mutant. We performed Exome Capture Sequencing (ECS) and identified two large deletion segments, named 6BL.1 at 334.8∼424.3 Mb and 6BL.2, 579.4∼717.8 Mb in the wds mutant. RNA-seq analysis confirmed that genes located in these regions lost their RNA expression. We then found that the 6BL.2 locus was overlapping with a known spike length QTL, qSL6B.2. Totally, 499 genes were located within the deleted region and two of them were found to be positively correlated with long spike accessions but not the ones with short spike. One of them, TraesCS6B01G334600, a well-matched homolog of the rice OsBUL1 gene that works in the Brassinosteroids (BR) pathway, was identified to be involved in cell size and number regulation. Further transcriptome analysis of young spikes showed that hormone-related genes were enriched among differentially expressed genes, supporting TraesCS6B01G334600 as a candidate gene. Our work provides a strategy to rapid locate genetic loci with large genomic lesions in wheat and useful resources for future wheat study.

Exome Sequencing from Bulked Segregant Analysis Identifies a Gene for All-Stage Resistance to Stripe Rust on Chromosome 1AL in Chinese Wheat Landrace 'Xiaohemai'.

Stripe rust caused by Puccinia striiformis f. sp. tritici is one of the most destructive diseases of wheat. Identifying novel resistance genes applicable for developing disease-resistant cultivars is important for the sustainable control of wheat stripe rust. Chinese wheat landrace 'Xiaohemai' ('XHM') is an elite germplasm line with all-stage resistance (ASR) effective against predominant Chinese P. striiformis f. sp. tritici races. In this study, we performed a bulked segregant analysis coupled with exome capture sequencing (BSE-seq) to identify a candidate genomic region strongly associated with stripe rust resistance on chromosome 1AL in 173 F2:3 lines derived from the cross 'XHM' × 'Avocet S'. The gene, designated as YrXH-1AL, was validated by a conventional quantitative trait locus analysis using newly developed Kompetitive allele-specific PCR (KASP) markers, explaining up to 48.50% of the phenotypic variance. By testing a secondary mapping population comprising 144 lines from the same cross at the seedling stage with prevalent P. striiformis f. sp. tritici race CYR34, YrXH-1AL was identified as a single Mendelian factor in a 1.5-cM interval flanked by KASP markers KP1A_484.33 and KP1A_490.09. This region corresponded to a 5.76-Mb genomic interval on 'Chinese Spring' chromosome 1AL. Furthermore, two cosegregating KASP markers showed high polymorphisms among 130 Chinese wheat cultivars and could be used for marker-assisted selection. Because no other Yr genes for ASR that originated from common wheat have been detected on chromosome 1AL, YrXH-1AL is likely a novel gene that can be incorporated into modern breeding materials to develop wheat cultivars with enhanced stripe rust resistance.

CCCH Zinc finger genes in Barley: genome-wide identification, evolution, expression and haplotype analysis

BackgroundCCCH transcription factors are important zinc finger transcription factors involved in the response to biotic and abiotic stress and physiological and developmental processes. Barley (Hordeum vulgare) is an agriculturally important cereal crop with multiple uses, such as brewing production, animal feed, and human food. The identification and assessment of new functional genes are important for the molecular breeding of barley.ResultsIn this study, a total of 53 protein-encoding CCCH genes unevenly dispersed on seven different chromosomes were identified in barley. Phylogenetic analysis categorized the barley CCCH genes (HvC3Hs) into eleven subfamilies according to their distinct features, and this classification was supported by intron–exon structure and conserved motif analysis. Both segmental and tandem duplication contributed to the expansion of CCCH gene family in barley. Genetic variation of HvC3Hs was characterized using publicly available exome-capture sequencing datasets. Clear genetic divergence was observed between wild and landrace barley populations in HvC3H genes. For most HvC3Hs, nucleotide diversity and the number of haplotype polymorphisms decreased during barley domestication. Furthermore, the HvC3H genes displayed distinct expression profiles for different developmental processes and in response to various types of stresses. The HvC3H1, HvC3H2 and HvC3H13 of arginine-rich tandem CCCH zinc finger (RR-TZF) genes were significantly induced by multiple types of abiotic stress and/or phytohormone treatment, which might make them as excellent targets for the molecular breeding of barley.ConclusionsOverall, our study provides a comprehensive characterization of barley CCCH transcription factors, their diversity, and their biological functions.

A major-effect genetic locus, ApRVII, controlling resistance against both adapted and non-adapted aphid biotypes in pea.

A genome-wide association study for pea resistance against a pea-adapted biotype and a non-adapted biotype of the aphid, Acyrthosiphon pisum, identified a genomic region conferring resistance to both biotypes. In a context of reduced insecticide use, the development of cultivars resistant to insect pests is crucial for an integrated pest management. Pea (Pisum sativum) is a crop of major importance among cultivated legumes, for the supply of dietary proteins and nitrogen in low-input cropping systems. However, yields of the pea crop have become unstable due to plant parasites. The pea aphid (Acyrthosiphon pisum) is an insect pest species forming a complex of biotypes, each one adapted to feed on one or a few related legume species. This study aimed to identify resistance to A. pisum and the underlying genetic determinism by examining a collection of 240 pea genotypes. The collection was screened against a pea-adapted biotype and a non-adapted biotype of A. pisum to characterize their resistant phenotype. Partial resistance was observed in some pea genotypes exposed to the pea-adapted biotype. Many pea genotypes were completely resistant to non-adapted biotype, but some exhibited partial susceptibility. A genome-wide association study, using pea exome-capture sequencing data, enabled the identification of the major-effect quantitative trait locus ApRVII on the chromosome 7. ApRVII includes linkage disequilibrium blocks significantly associated with resistance to one or both of the two aphid biotypes studied. Finally, we identified candidate genes underlying ApRVII that are potentially involved in plant-aphid interactions and marker haplotypes linked with aphid resistance. This study sets the ground for the functional characterization of molecular pathways involved in pea defence to the aphids but also is a step forward for breeding aphid-resistant cultivars.

Identification of QTL Associated with Regrowth Vigor Using the Nested Association Mapping Population in Switchgrass.

Switchgrass (Panicum virgatum L.) is a warm-season perennial grass species that is utilized as forage for livestock and biofuel feedstock. The stability of biomass yield and regrowth vigor under changing harvest frequency would help manage potential fluctuations in the feedstock market and would provide a continuous supply of quality forage for livestock. This study was conducted to (i) assess the genetic variation and (ii) identify the quantitative trait loci (QTL) associated with regrowth vigor after multiple cuttings in lowland switchgrass. A nested association mapping (NAM) population comprising 2000 pseudo F2 progenies was genotyped with single nucleotide polymorphism (SNP) markers derived from exome-capture sequencing and was evaluated for regrowth vigor in 2017 and 2018. The results showed significant variation among the NAM families in terms of regrowth vigor (p < 0.05). A total of 10 QTL were detected on 6 chromosomes: 1B, 5A, 5B, 6B, 7B, and 8A, explaining the phenotypic variation by up to 4.7%. The additive genetic effects of an individual QTL ranged from −0.13 to 0.26. No single QTL showed a markedly large effect, suggesting complex genetics underlying regrowth vigor in switchgrass. The homologs of candidate genes that play a variety of roles in developmental processes, including plant hormonal signal transduction, nucleotide biosynthesis, secondary metabolism, senescence, and responses to both biotic and abiotic stresses, were identified in the vicinity of QTL.

Association of Complement and MAPK Activation With SARS-CoV-2–Associated Myocardial Inflammation

Myocardial injury is a common feature of patients with SARS-CoV-2 infection. However, the cardiac inflammatory processes associated with SARS-CoV-2 infection are not completely understood. To investigate the inflammatory cardiac phenotype associated with SARS-CoV-2 infection compared with viral myocarditis, immune-mediated myocarditis, and noninflammatory cardiomyopathy by integrating histologic, transcriptomic, and proteomic profiling. This case series was a cooperative study between the Ludwig Maximilian University Hospital Munich and the Cardiopathology Referral Center at the University of Tübingen in Germany. A cohort of 19 patients with suspected myocarditis was examined; of those, 5 patients were hospitalized with SARS-CoV-2 infection between March and May 2020. Cardiac tissue specimens from those 5 patients were compared with specimens from 5 patients with immune-mediated myocarditis, 4 patients with non-SARS-CoV-2 viral myocarditis, and 5 patients with noninflammatory cardiomyopathy, collected from January to August 2019. Endomyocardial biopsy. The inflammatory cardiac phenotypes were measured by immunohistologic analysis, RNA exome capture sequencing, and mass spectrometry-based proteomic analysis of endomyocardial biopsy specimens. Among 19 participants, the median age was 58 years (range, 37-76 years), and 15 individuals (79%) were male. Data on race and ethnicity were not collected. The abundance of CD163+ macrophages was generally higher in the cardiac tissue of patients with myocarditis, whereas lymphocyte counts were lower in the tissue of patients with SARS-CoV-2 infection vs patients with non-SARS-CoV-2 virus-associated and immune-mediated myocarditis. Among those with SARS-CoV-2 infection, components of the complement cascade, including C1q subunits (transcriptomic analysis: 2.5-fold to 3.6-fold increase; proteomic analysis: 2.0-fold to 3.4-fold increase) and serine/cysteine proteinase inhibitor clade G member 1 (transcriptomic analysis: 1.7-fold increase; proteomic analysis: 2.6-fold increase), belonged to the most commonly upregulated transcripts and differentially abundant proteins. In cardiac macrophages, the abundance of C1q was highest in SARS-CoV-2 infection. Assessment of important signaling cascades identified an upregulation of the serine/threonine mitogen-activated protein kinase pathways. This case series found that the cardiac immune signature varied in inflammatory conditions with different etiologic characteristics. Future studies are needed to examine the role of these immune pathways in myocardial inflammation.

Development of the Wheat Practical Haplotype Graph database as a resource for genotyping data storage and genotype imputation.

To improve the efficiency of high-density genotype data storage and imputation in bread wheat (Triticum aestivum L.), we applied the Practical Haplotype Graph (PHG) tool. The Wheat PHG database was built using whole-exome capture sequencing data from a diverse set of 65 wheat accessions. Population haplotypes were inferred for the reference genome intervals defined by the boundaries of the high-quality gene models. Missing genotypes in the inference panels, composed of wheat cultivars or recombinant inbred lines genotyped by exome capture, genotyping-by-sequencing (GBS), or whole-genome skim-seq sequencing approaches, were imputed using the Wheat PHG database. Though imputation accuracy varied depending on the method of sequencing and coverage depth, we found 92% imputation accuracy with 0.01× sequence coverage, which was slightly lower than the accuracy obtained using the 0.5× sequence coverage (96.6%). Compared to Beagle, on average, PHG imputation was ∼3.5% (P-value < 2 × 10−14) more accurate, and showed 27% higher accuracy at imputing a rare haplotype introgressed from a wild relative into wheat. We found reduced accuracy of imputation with independent 2× GBS data (88.6%), which increases to 89.2% with the inclusion of parental haplotypes in the database. The accuracy reduction with GBS is likely associated with the small overlap between GBS markers and the exome capture dataset, which was used for constructing PHG. The highest imputation accuracy was obtained with exome capture for the wheat D genome, which also showed the highest levels of linkage disequilibrium and proportion of identity-by-descent regions among accessions in the PHG database. We demonstrate that genetic mapping based on genotypes imputed using PHG identifies SNPs with a broader range of effect sizes that together explain a higher proportion of genetic variance for heading date and meiotic crossover rate compared to previous studies.

High rate of hypomorphic variants as the cause of inherited ataxia and related diseases: study of a cohort of 366 families

PurposeDiagnosis of inherited ataxia and related diseases represents a real challenge given the tremendous heterogeneity and clinical overlap of the various causes. We evaluated the efficacy of molecular diagnosis of these diseases by sequencing a large cohort of undiagnosed families. MethodsWe analyzed 366 unrelated consecutive patients with undiagnosed ataxia or related disorders by clinical exome-capture sequencing. In silico analysis was performed with an in-house pipeline that combines variant ranking and copy-number variant (CNV) searches. Variants were interpreted according to American College of Medical Genetics and Genomics/Association for Molecular Pathology (ACMG/AMP) guidelines. ResultsWe established the molecular diagnosis in 46% of the cases. We identified 35 mildly affected patients with causative variants in genes that are classically associated with severe presentations. These cases were explained by the occurrence of hypomorphic variants, but also rarely suspected mechanisms such as C-terminal truncations and translation reinitiation. ConclusionA significant fraction of the clinical heterogeneity and phenotypic overlap is explained by hypomorphic variants that are difficult to identify and not readily predicted. The hypomorphic C-terminal truncation and translation reinitiation mechanisms that we identified may only apply to few genes, as it relies on specific domain organization and alterations. We identified PEX10 and FASTKD2 as candidates for translation reinitiation accounting for mild disease presentation.

Genetic loci associated with winter survivorship in diverse lowland switchgrass populations.

High winter mortality limits biomass yield of lowland switchgrass (Panicum virgatumL.) planted in the northern latitudes of North America. Breeding of cold tolerant switchgrass cultivars requires many years due to its perennial growth habit and the unpredictable winter selection pressure that is required to identify winter-hardy individuals. Identification of causal genetic variants for winter survivorship would accelerate the improvement of switchgrass biomass production. The objective of this study was to identify allelic variation associated with winter survivorship in lowland switchgrass populations using bulk segregant analysis (BSA). Twenty-nine lowland switchgrass populations were evaluated for winter survival at two locations in southern Wisconsin and 21 populations with differential winter survivorship were used for BSA. A maximum of 10% of the individuals (8-20) were bulked to create survivor and nonsurvivor DNA pools from each population and location. The DNA pools were evaluated using exome capture sequencing, and allele frequencies were used to conduct statistical tests. The BSA tests revealed nine quatitative trait loci (QTL) from tetraploid populations and seven QTL from octoploid populations. Many QTL were population-specific, but some were identified in multiple populations that originated across a broad geographic landscape. Four QTL (at positions 88 Mb on chromosome 2N, 115 Mb on chromosome 5K, and 1 and 100 Mb on chromosome 9N) were potentially the most useful QTL. Markers associated with winter survivorship in this study can be used to accelerate breeding cycles of lowland switchgrass populations and should lead to improvements in adaptation within USDA hardiness zones 4 and 5.

Screening for monogenic subtypes of gestational diabetes in a high prevalence island population - A whole exome sequencing study.

The reported frequency of monogenic defects of beta cell function in gestational diabetes (GDM) varies extensively. This study aimed to evaluate the frequency and molecular spectrum of variants in genes associated with monogenic/atypical diabetes in non-obese females of Maltese ethnicity with GDM. 50 non-obese females who met the International Association of the Diabetes and Pregnancy Study Groups (IADPSG) criteria for diagnosis of GDM and with a first-degree relative with non-autoimmune diabetes were included in this study. Whole exome capture and high throughput sequencing was carried out. Rare sequence variants were filtered, annotated, and prioritised according to the American College for Medical Genetics guidelines. For selected missense variants we explored effects on protein stability and structure through in-silico tools. We identified three pathogenic variants in GCK, ABCC8 and HNF1A and several variants of uncertain significance in the cohort. Genotype-phenotype correlations and post-pregnancy follow-up data are described. This study provides the first insight into an underlying monogenic aetiology in non-obese females with GDM from an island population having a high prevalence of diabetes. It suggests that monogenic variants constitute an underestimated cause of diabetes detected in pregnancy, and that careful evaluation of GDM probands to identify monogenic disease subtypes is indicated.

Abstract PS5-24: Novel genomic variants and pathways associated with baseline serum thymidine kinase 1 levels in HR-positive HER2-negative MBC patients commencing palbociclib and letrozole

Abstract Background: Cyclin dependent 4/6 kinase inhibitors (CDK4/6i) and endocrine therapy (ET) have improved progression-free survival (PFS) and overall survival in hormone-receptor (HR)-positive metastatic breast cancer (MBC), but progression of disease is inevitable. Serum thymidine kinase-1 (TK1) is a secreted marker of proliferation that is prognostic in patients (pts) with HR-positive, HER2-negative MBC and may be predictive of ET and CDK 4/6i response. PROMISE (NCT0281902) is a prospective study enrolling women with HR-positive MBC starting palbociclib (Pb) + letrozole (L) (1st line) or Pb + fulvestrant (2nd line). We undertook a comprehensive “omic” assessment of blood, tumor, urine and the fecal microbiome in order to identify novel genomic variants and pathways associated with an early decline in TK1 (measured after 2 months) and PFS. Additionally, patient derived xenografts/organoids were generated at baseline and progression to test new therapeutic approaches to overcome resistance to CDK4/6i and ET. We report the initial association between the baseline genomic landscape and baseline TK1 levels. Methods: FFPE tumor biopsies were obtained for DNA/RNA sequencing (TempusTM) and blood samples for TK1 (Divitum® assay, Biovica) were collected pretreatment (pre-Pb) and after 2 cycles of Pb + ET (post-Pb2). Both whole-exome (exome capture) sequencing (WES) and RNA-Seq used the Integrated DNA Technologies xGen Exome Research Panel v1.0 capture kit. TK1+ disease was defined as &amp;gt; 200 Du/L and TK1- disease as below limit of detection up to 200 Du/L. We tested the association between genomic and transcriptomic characteristics with baseline TK1 data in pretreatment samples where both WES and RNA-seq and TK1 was available. The data were analyzed using bioinformatics pipelines for somatic and germline mutations and copy number alterations. The current analysis focuses on baseline 1st-line pre-Pb omics data in conjunction with baseline TK1 levels. The database was locked for analysis on 5/29/2020. Results: Thirty-three pts (median age: 59 yrs.) were evaluable, with paired samples for TK1 in 32. Six pts had TK1+ disease pre-Pb and post-Pb2. Twenty-two pts had TK1- disease pre-Pb and post-Pb2. Four pts had a decrease in TK1 after 2 cycles of treatment that altered the classification from TK1+ to TK1-. Both baseline RNA seq and serum TK1 (n=16) were available for 4 TK1+ and 12 TK1- pts. In this group, 476 genes were differentially regulated (398 upregulated; 78 downregulated). Pathway analysis demonstrated enrichment in complement and coagulation cascade pathway, PPAR signaling pathway, and metabolism-related pathways related to up-regulation of CYP and UGT gene families. Further testing for the association of WES data with baseline TK1+ (n=8) and TK1- (n=16) disease demonstrated somatic copy number variations on chromosomes 6, 11, 12 and 15. CDK4 copy number gains were observed in 3/8 TK1+ pts and 0/16 TK1- pts. We also observed that somatic mutations (LOH, copy number and/or SNV/INDELs) were more prevalent in the TK1+ compared to the TK1- pre-Pb group in several cancer-associated genes (FAS [p=0.06] PTEN, PIK3CB, NAB2, SOX9 and FAT1 [p=0.08], TP53, and MAP2K4 [p=0.22]). Conversely, we also noted that 6/7 pts with GATA3 mutations had TK1- disease (p=0.23). Conclusions: Using a comprehensive “omics” approach, our data suggest that a secreted biomarker of proliferation (TK1) obtained prior to initiating CDK4/6i and ET for the first line treatment of HR+ MBC is associated with established and novel genes and pathways associated with prognosis of pts receiving ET and CDK 4/6i. Analysis of on-treatment (after 2 cycles) tumor RNA seq and its association with change in TK1 as well as data from the 2nd-line cohort will be presented at the meeting. Citation Format: Ciara C O Sullivan, Krishna R Kalari, Vera J Suman, Peter T Vedell, Ann Moyer, Abraham D Eyman Casey, Jason Sinnwell, Xiaojia Tang, Kevin Thompson, Alvaro Moreno-Aspitia, Donald W Northfelt, Minetta C Liu, Tufia C Haddad, Saranya Chumsri, Prema Peethambaram, Kathryn J Ruddy, Karthik V Giridhar, Roberto A Leon-Ferre, Adrian Nordmark, Mattias Bergqvist, Brendan P McMenomy, Richard M Weinshilboum, Liewei Wang, Matthew P Goetz. Novel genomic variants and pathways associated with baseline serum thymidine kinase 1 levels in HR-positive HER2-negative MBC patients commencing palbociclib and letrozole [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 PS5-24.

A functional genomics approach to dissect spotted alfalfa aphid resistance in Medicago truncatula

Aphids are virus-spreading insect pests affecting crops worldwide and their fast population build-up and insecticide resistance make them problematic to control. Here, we aim to understand the molecular basis of spotted alfalfa aphid (SAA) or Therioaphis trifolii f. maculata resistance in Medicago truncatula, a model organism for legume species. We compared susceptible and resistant near isogenic Medicago lines upon SAA feeding via transcriptome sequencing. Expression of genes involved in defense and stress responses, protein kinase activity and DNA binding were enriched in the resistant line. Potentially underlying some of these changes in gene expression was the finding that members of the MYB, NAC, AP2 domain and ERF transcription factor gene families were differentially expressed in the resistant versus susceptible lines. A TILLING population created in the resistant cultivar was screened using exome capture sequencing and served as a reverse genetics tool to functionally characterise genes involved in the aphid resistance response. This screening revealed three transcription factors (a NAC, AP2 domain and ERF) as important regulators in the defence response, as a premature stop-codon in the resistant background led to a delay in aphid mortality and enhanced plant susceptibility. This combined functional genomics approach will facilitate the future development of pest resistant crops by uncovering candidate target genes that can convey enhanced aphid resistance.

Combining a New Exome Capture Panel With an Effective varBScore Algorithm Accelerates BSA-Based Gene Cloning in Wheat.

The discovery of functional genes underlying agronomic traits is of great importance for wheat improvement. Here we designed a new wheat exome capture probe panel based on IWGSC RefSeq v1.0 genome sequence information and developed an effective algorithm, varBScore, that can sufficiently reduce the background noise in gene mapping and identification. An effective method, termed bulked segregant exome capture sequencing (BSE-Seq) for identifying causal mutations or candidate genes was established by combining the use of a newly designed wheat exome capture panel, sequencing of bulked segregant pools from segregating populations, and the robust algorithm varBScore. We evaluated the effectiveness of varBScore on SNP calling using the published dataset for mapping and cloning the yellow rust resistance gene Yr7 in wheat. Furthermore, using BSE-Seq, we rapidly identified a wheat yellow leaf mutant gene, ygl1, in an ethyl methanesulfonate (EMS) mutant population and found that a single mutation of G to A at 921 position in the wild type YGL1 gene encoding magnesium-chelatase subunit chlI caused the leaf yellowing phenotype. We further showed that mutation of YGL1 through CRISPR/Cas9 gene editing led to a yellow phenotype on the leaves of transgenic wheat, indicating that ygl1 is the correct causal gene responsible for the mutant phenotype. In summary, our approach is highly efficient for discovering causal mutations and gene cloning in wheat.

Abstract A14: Next-generation sequencing of Bacille Calmette-Guerin (BCG)-unresponsive tumors identifies actionable alterations for potential targeted therapy in non-muscle invasive bladder cancer

Abstract Introduction and Objective: Intravesical BCG remains standard of care for patients with intermediate- and high-risk non-muscle invasive bladder cancer (NMIBC). Although most patients initially respond to therapy, there is a group of patients who will experience recurrence and subsequent progression to invasive disease. There are limited treatment options for patients with BCG-unresponsive bladder cancer and alternative therapeutic strategies are needed. Next-generation sequencing helps identify and prioritize therapeutic opportunities for clinical trial of novel targeted agents. Methods: All patients with urothelial carcinoma of the bladder who progressed to muscle-invasive disease after at least one course of intravesical BCG were identified within our institutional genomic database and in The Cancer Genome Atlas (TCGA) dataset. DNA from BCG-treated, chemotherapy naive, secondary MIBC tumors and matched normal DNA underwent targeted exome capture sequencing or whole-exome sequencing. Potentially actionable genomic alterations were defined according to OncoKB levels of evidence for each alteration as a predictive biomarker (OncoKB.org). Results: The cohort comprised 77 secondary MIBC specimens, each from a unique patient. Sixty-six (86%) were male and the median age was 67 years (IQR 61-75). At least one potentially actionable genomic alteration was identified in 85% of the specimens. RTK/MAPK pathway alterations, most commonly FGFR3 (18%) and ERBB2 (15%), were seen in 42% of patients. PIK3CA mutations, seen in 21%, often co-occurred with RTK/MAPK alterations, potentially representing a resistance mechanism to FGFR or ERBB2-directed monotherapy. CDKN2A loss, hypothesized to be a driver of progression in FGFR3-mutated tumors, occurred in 10% of the samples. Alterations in cell cycle regulators were seen in 65% of specimens, of which 40% also had an intact RB gene and could potentially be vulnerable to CDK4/6 inhibitors. MDM2 amplification with “wild-type” TP53 was seen in 13% of specimens that might be sensitive to small-molecule antagonists. Inactivating mutations in histone acetyltransferases, CREBBP or EP300, occurred in 17% of specimens that might confer sensitivity to HDAC inhibitors. Truncating mutations in the SWI/SNF chromatin-modifying complex occurred in 33%, predominately in ARID1A, which we recently reported to be associated with worse recurrence-free survival after BCG. Truncating KDM6A mutations were seen in 19%, which along with ARID1A mutated tumors may be vulnerable to EZH2 inhibitors. Conclusion: Limited effective treatments beyond radical cystectomy are available for BCG-unresponsive NMIBC. Genomic profiling identified several actionable alterations that offer potential therapeutic strategies for management of BCG-unresponsive disease. Targeted therapy for NMIBC patients warrants further investigation as further advances are made in development of more selective systemic inhibitors and improved intravesical drug delivery. Citation Format: Aleksandra Walasek, Nima Almassi, Victor McPherson, Shawn Dason, Nicole Benfante, Tanya Klein, Michael F. Berger, Nikolaus Schultz, Guido Dalbagni, Dean F. Bajorin, Jonathan E. Rosenberg, Hikmat Al-Ahmadie, Bernard H. Bochner, David B. Solit, Gopakumar Iyer, Eugene J. Pietzak. Next-generation sequencing of Bacille Calmette-Guerin (BCG)-unresponsive tumors identifies actionable alterations for potential targeted therapy in non-muscle invasive bladder cancer [abstract]. In: Proceedings of the AACR Special Conference on Bladder Cancer: Transforming the Field; 2019 May 18-21; Denver, CO. Philadelphia (PA): AACR; Clin Cancer Res 2020;26(15_Suppl):Abstract nr A14.

USH2A variants in Chinese patients with Usher syndrome type II and non-syndromic retinitis pigmentosa

AimsTo reveal the Usher syndrome type IIA (USH2A) gene variant profile in a large cohort of Chinese patients with non-syndromic retinitis pigmentosa (RP) or Usher syndrome type II (USH2) and...

Kinome capture sequencing of high-grade serous ovarian carcinoma reveals novel mutations in the JAK3 gene.

High-grade serous ovarian carcinoma (HGSOC) remains the deadliest form of epithelial ovarian cancer and despite major efforts little improvement in overall survival has been achieved. Identification of recurring “driver” genetic lesions has the potential to enable design of novel therapies for cancer. Here, we report on a study to find such new therapeutic targets for HGSOC using exome-capture sequencing approach targeting all kinase genes in 127 patient samples. Consistent with previous reports, the most frequently mutated gene was TP53 (97% mutation frequency) followed by BRCA1 (10% mutation frequency). The average mutation frequency of the kinase genes mutated from our panel was 1.5%. Intriguingly, after BRCA1, JAK3 was the most frequently mutated gene (4% mutation frequency). We tested the transforming properties of JAK3 mutants using the Ba/F3 cell-based in vitro functional assay and identified a novel gain-of-function mutation in the kinase domain of JAK3 (p.T1022I). Importantly, p.T1022I JAK3 mutants displayed higher sensitivity to the JAK3-selective inhibitor Tofacitinib compared to controls. For independent validation, we re-sequenced the entire JAK3 coding sequence using tagged amplicon sequencing (TAm-Seq) in 463 HGSOCs resulting in an overall somatic mutation frequency of 1%. TAm-Seq screening of CDK12 in the same population revealed a 7% mutation frequency.Our data confirms that the frequency of mutations in kinase genes in HGSOC is low and provides accurate estimates for the frequency of JAK3 and CDK12 mutations in a large well characterized cohort. Although p.T1022I JAK3 mutations are rare, our functional validation shows that if detected they should be considered as potentially actionable for therapy. The observation of CDK12 mutations in 7% of HGSOC cases provides a strong rationale for routine somatic testing, although more functional and clinical characterization is required to understand which nonsynonymous mutations alterations are associated with homologous recombination deficiency.

Effects of landscapes and range expansion on population structure and local adaptation.

Understanding the origin and distribution of genetic diversity across landscapes is critical for predicting the future of organisms in changing climates. This study investigated how adaptive and demographic forces have shaped diversity and population structure in Pinus densata, a keystone species on Qinghai-Tibetan Plateau (QTP). We examined the distribution of genomic diversity across the range of P.densata using exome capture sequencing. We applied spatially explicit tests to dissect the impacts of allele surfing, geographic isolation and environmental gradients on population differentiation and forecasted how this genetic legacy may limit the persistence of P.densata in future climates. We found that allele surfing from range expansion could explain the distribution of 39% of the c. 48000 genotyped single nucleotide polymorphisms (SNPs). Uncorrected, these allele frequency clines severely confounded inferences of selection. After controlling for demographic processes, isolation-by-environment explained 9.2-19.5% of the genetic structure, with c. 4.0% of loci being affected by selection. Allele surfing and genotype-environment associations resulted in genomic mismatch under projected climate scenarios. We illustrate that significant local adaptation, when coupled with reduced diversity as a result of demographic history, constrains potential evolutionary response to climate change. The strong signal of genomic vulnerability in P.densata may be representative for other QTP endemics.