Search

The development of immune checkpoint–based immunotherapies has been a major advancement in the treatment of cancer, with a subset of patients exhibiting durable clinical responses. A predictive biomarker for immunotherapy response is the preexisting T-cell infiltration in the tumor immune microenvironment (TIME). Bulk transcriptomics–based approaches can quantify the degree of T-cell infiltration using deconvolution methods and identify additional markers of inflamed/cold cancers at the bulk level. However, bulk techniques are unable to identify biomarkers of individual cell types. Although single-cell RNA sequencing (scRNA-seq) assays are now being used to profile the TIME, to our knowledge there is no method of identifying patients with a T-cell inflamed TIME from scRNA-seq data. Here, we describe a method, iBRIDGE, which integrates reference bulk RNA-seq data with the malignant subset of scRNA-seq datasets to identify patients with a T-cell inflamed TIME. Using two datasets with matched bulk data, we show iBRIDGE results correlated highly with bulk assessments (0.85 and 0.9 correlation coefficients). Using iBRIDGE, we identified markers of inflamed phenotypes in malignant cells, myeloid cells, and fibroblasts, establishing type I and type II interferon pathways as dominant signals, especially in malignant and myeloid cells, and finding the TGFβ-driven mesenchymal phenotype not only in fibroblasts but also in malignant cells. Besides relative classification, per-patient average iBRIDGE scores and independent RNAScope quantifications were used for threshold-based absolute classification. Moreover, iBRIDGE can be applied to in vitro grown cancer cell lines and can identify the cell lines that are adapted from inflamed/cold patient tumors.

Supplementary Figures S1-S6

Turan et. al Supplementary Methods File S1

The development of immune checkpoint–based immunotherapies has been a major advancement in the treatment of cancer, with a subset of patients exhibiting durable clinical responses. A predictive biomarker for immunotherapy response is the preexisting T-cell infiltration in the tumor immune microenvironment (TIME). Bulk transcriptomics–based approaches can quantify the degree of T-cell infiltration using deconvolution methods and identify additional markers of inflamed/cold cancers at the bulk level. However, bulk techniques are unable to identify biomarkers of individual cell types. Although single-cell RNA sequencing (scRNA-seq) assays are now being used to profile the TIME, to our knowledge there is no method of identifying patients with a T-cell inflamed TIME from scRNA-seq data. Here, we describe a method, iBRIDGE, which integrates reference bulk RNA-seq data with the malignant subset of scRNA-seq datasets to identify patients with a T-cell inflamed TIME. Using two datasets with matched bulk data, we show iBRIDGE results correlated highly with bulk assessments (0.85 and 0.9 correlation coefficients). Using iBRIDGE, we identified markers of inflamed phenotypes in malignant cells, myeloid cells, and fibroblasts, establishing type I and type II interferon pathways as dominant signals, especially in malignant and myeloid cells, and finding the TGFβ-driven mesenchymal phenotype not only in fibroblasts but also in malignant cells. Besides relative classification, per-patient average iBRIDGE scores and independent RNAScope quantifications were used for threshold-based absolute classification. Moreover, iBRIDGE can be applied to in vitro grown cancer cell lines and can identify the cell lines that are adapted from inflamed/cold patient tumors.

Supplementary Tables S1-S16

The majority of causative variants in familial breast cancer remain unknown. Of the known risk variants, most are tumor cell autonomous, and little attention has been paid yet to germline variants that may affect the tumor microenvironment. In this study, we developed a system called the Consomic Xenograft Model (CXM) to map germline variants that affect only the tumor microenvironment. In CXM, human breast cancer cells are orthotopically implanted into immunodeficient consomic strains and tumor metrics are quantified (e.g., growth, vasculogenesis, and metastasis). Because the strain backgrounds vary, whereas the malignant tumor cells do not, any observed changes in tumor progression are due to genetic differences in the nonmalignant microenvironment. Using CXM, we defined genetic variants on rat chromosome 3 that reduced relative tumor growth and hematogenous metastasis in the SS.BN3IL2Rγ consomic model compared with the SSIL2Rγ parental strain. Paradoxically, these effects occurred despite an increase in the density of tumor-associated blood vessels. In contrast, lymphatic vasculature and lymphogenous metastasis were unaffected by the SS.BN3IL2Rγ background. Through comparative mapping and whole-genome sequence analysis, we narrowed candidate variants on rat chromosome 3 to six genes with a priority for future analysis. Collectively, our results establish the utility of CXM to localize genetic variants affecting the tumor microenvironment that underlie differences in breast cancer risk. Cancer Res; 74(22); 6419–29. ©2014 AACR.

Expanded Materials and Methods, Supplemental Figures, and Supplemental Tables Figure S1. Flow cytometric analysis of peripheral blood from SS, SS.BN3, SSIL2Rγ, and SS.BN3IL2Rγ rats. Figure S2.SH2B3 specifically marks CD31+ blood vessels, but not LYVE1+ lymphatic vessels. Table S1. Peripheral blood analysis of SS, SS.BN3, SSIL2Rγ and SS.BN3IL2Rγ rats. Table S2. Genes located in the overlapping portions of the rat chromosome 3 cancer QTLs. Table S3. Variants that cosegregate between tumor resistant (BN and Cop) and tumor susceptible (ACI, F344, and SS) rat strains.

Background The genes and mechanisms involved in the association between diabetes or hypertension and CKD risk are unclear. Previous studies have implicated a role for γ-adducin (ADD3), a cytoskeletal protein encoded by Add3. Methods We investigated renal vascular function in vitro and in vivo and the susceptibility to CKD in rats with wild-type or mutated Add3 and in genetically modified rats with overexpression or knockout of ADD3. We also studied glomeruli and primary renal vascular smooth muscle cells isolated from these rats. Results This study identified a K572Q mutation in ADD3 in fawn-hooded hypertensive (FHH) rats-a mutation previously reported in Milan normotensive (MNS) rats that also develop kidney disease. Using molecular dynamic simulations, we found that this mutation destabilizes a critical ADD3-ACTIN binding site. A reduction of ADD3 expression in membrane fractions prepared from the kidney and renal vascular smooth muscle cells of FHH rats was associated with the disruption of the F-actin cytoskeleton. Compared with renal vascular smooth muscle cells from Add3 transgenic rats, those from FHH rats had elevated membrane expression of BKα and BK channel current. FHH and Add3 knockout rats exhibited impairments in the myogenic response of afferent arterioles and in renal blood flow autoregulation, which were rescued in Add3 transgenic rats. We confirmed these findings in a genetic complementation study that involved crossing FHH and MNS rats that share the ADD3 mutation. Add3 transgenic rats showed attenuation of proteinuria, glomerular injury, and kidney fibrosis with aging and mineralocorticoid-induced hypertension. Conclusions This is the first report that a mutation in ADD3 that alters ACTIN binding causes renal vascular dysfunction and promotes the susceptibility to kidney disease.

In precision oncology, reliable identification of tumor-specific DNA mutations requires sequencing tumor DNA and non-tumor DNA (so-called “matched normal”) from the same patient. The normal sample allows researchers to distinguish acquired (somatic) and hereditary (germline) variants. The ability to distinguish somatic and germline variants facilitates estimation of tumor mutation burden (TMB), which is a recently FDA-approved pan-cancer marker for highly successful cancer immunotherapies; in tumor-only variant calling (i.e., without a matched normal), the difficulty in discriminating germline and somatic variants results in inflated and unreliable TMB estimates. We apply machine learning to the task of somatic vs germline classification in tumor-only samples using TabNet, a recently developed attentive deep learning model for tabular data that has achieved state of the art performance in multiple classification tasks (Arik and Pfister 2019). We constructed a training set for supervised classification using features derived from tumor-only variant calling and drawing somatic and germline truth-labels from an independent pipeline incorporating the patient-matched normal samples. Our trained model achieved state-of-the-art performance on two hold-out test datasets: a TCGA dataset including sarcoma, breast adenocarcinoma, and endometrial carcinoma samples (F1-score: 88.3), and a metastatic melanoma dataset, (F1-score 79.8). Concordance between matched-normal and tumor-only TMB improves from R2 = 0.006 to 0.705 with the addition of our classifier. And importantly, this approach generalizes across tumor tissue types and capture kits and has a call rate of 100%. The interpretable feature masks of the attentive deep learning model explain the reasons for misclassified variants. We reproduce the recent finding that tumor-only TMB estimates for Black patients are extremely inflated relative to that of White patients due to the racial biases of germline databases. We show that our machine learning approach appreciably reduces this racial bias in tumor-only variant-calling.

Background The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants. Methods To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles. Results We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33). Conclusions This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas.

Contains fulltext : 202800.pdf (Publisher’s version ) (Open Access) SMARCC2 (BAF170) is one of the invariable core subunits of the ATP-dependent chromatin remodeling BAF (BRG1-associated factor) complex and plays a crucial role in embryogenesis and corticogenesis. Pathogenic variants in genes encoding other components of the BAF complex have been associated with intellectual disability syndromes. Despite its significant biological role, variants in SMARCC2 have not been directly associated with human disease previously. Using whole-exome sequencing and a web-based gene-matching program, we identified 15 individuals with variable degrees of neurodevelopmental delay and growth retardation harboring one of 13 heterozygous variants in SMARCC2, most of them novel and proven de novo. The clinical presentation overlaps with intellectual disability syndromes associated with other BAF subunits, such as Coffin-Siris and Nicolaides-Baraitser syndromes and includes prominent speech impairment, hypotonia, feeding difficulties, behavioral abnormalities, and dysmorphic features such as hypertrichosis, thick eyebrows, thin upper lip vermilion, and upturned nose. Nine out of the fifteen individuals harbor variants in the highly conserved SMARCC2 DNA-interacting domains (SANT and SWIRM) and present with a more severe phenotype. Two of these individuals present cardiac abnormalities. Transcriptomic analysis of fibroblasts from affected individuals highlights a group of differentially expressed genes with possible roles in regulation of neuronal development and function, namely H19, SCRG1, RELN, and CACNB4. Our findings suggest a novel SMARCC2-related syndrome that overlaps with neurodevelopmental disorders associated with variants in BAF-complex subunits.

The UK Biobank Exome Sequencing Consortium (UKB-ESC) is a unique private/public partnership between the UK Biobank and eight biopharma companies that will sequence the exomes of all ∼500,000 UK Biobank participants. Here we describe early results from the exome sequence data generated by this consortium for the first ∼200,000 UKB subjects and the key features of this project that enabled the UKB-ESC to come together and generate this data.Exome sequencing data from the first 200,643 UKB enrollees are now accessible to the research community. Approximately 10M variants were observed within the targeted regions, including: 8,086,176 SNPs, 370,958 indels and 1,596,984 multi-allelic variants. Of the ∼8M variants observed, 84.5% are coding variants and include 2,139,318 (25.3%) synonymous, 4,549,694 (53.8%) missense, 453,733 (5.4%) predicted loss-of-function (LOF) variants (initiation codon loss, premature stop codons, stop codon loss, splicing and frameshift variants) affecting at least one coding transcript. This open access data provides a rich resource of coding variants for rare variant genetic studies, and is particularly valuable for drug discovery efforts that utilize rare, functionally consequential variants.Over the past decade, the biopharma industry has increasingly leveraged human genetics as part of their drug discovery and development strategies. This shift was motivated by technical advances that enabled cost-effective human genetics research at scale, the emergence of electronic health records and biobanks, and a maturing understanding of how human genetics can increase the probability of successful drug development. Recognizing the need for large-scale human genetics data to drive drug discovery, and the unique value of the open data access policies and contribution terms of the UK Biobank, the UKB-ESC was formed. This precompetitive collaboration has further strengthened the ties between academia and industry and provided teams an unprecedented opportunity to interact with and learn from the wider research community.

Response to a drug often differs widely among individual patients. This variability is frequently observed not only with respect to effective responses but also with adverse drug reactions. Matching patients to the drugs that are most likely to be effective and least likely to cause harm is the goal of effective therapeutics. Pharmacogenomics (PGx) holds the promise of precision medicine through elucidating the genetic determinants responsible for pharmacological outcomes and using them to guide drug selection and dosing. Here, we survey the US landscape of research programs in PGx implementation, review current advances and clinical applications of PGx, summarize the obstacles that have hindered PGx implementation, and identify the critical knowledge gaps and possible studies needed to help to address them.

Genome sequencing has revolutionized the diagnosis of genetic diseases. Close collaborations between basic scientists and clinical genomicists are now needed to link genetic variants with disease causation. To facilitate such collaborations, we recommend prioritizing clinically relevant genes for functional studies, developing reference variant-phenotype databases, adopting phenotype description standards, and promoting data sharing.

Diagnosis at the edges of our knowledge calls upon clinicians to be data driven, cross-disciplinary, and collaborative in unprecedented ways. Exact disease recognition, an element of the concept of precision in medicine, requires new infrastructure that spans geography, institutional boundaries, and the divide between clinical care and research. The National Institutes of Health (NIH) Common Fund supports the Undiagnosed Diseases Network (UDN) as an exemplar of this model of precise diagnosis. Its goals are to forge a strategy to accelerate the diagnosis of rare or previously unrecognized diseases, to improve recommendations for clinical management, and to advance research, especially into disease mechanisms. The network will achieve these objectives by evaluating patients with undiagnosed diseases, fostering a breadth of expert collaborations, determining best practices for translating the strategy into medical centers nationwide, and sharing findings, data, specimens, and approaches with the scientific and medical communities. Building the UDN has already brought insights to human and medical geneticists. The initial focus has been on data sharing, establishing common protocols for institutional review boards and data sharing, creating protocols for referring and evaluating patients, and providing DNA sequencing, metabolomic analysis, and functional studies in model organisms. By extending this precision diagnostic model nationally, we strive to meld clinical and research objectives, improve patient outcomes, and contribute to medical science.

Whole-exome sequencing (WES) has increasingly enabled new pathogenic gene variant identification for undiagnosed neurodevelopmental disorders and provided insights into both gene function and disease biology. Here, we describe seven children with a neurodevelopmental disorder characterized by microcephaly, profound developmental delays and/or intellectual disability, cataracts, severe epilepsy including infantile spasms, irritability, failure to thrive, and stereotypic hand movements. Brain imaging in these individuals reveals delay in myelination and cerebral atrophy. We observe an identical recurrent de novo heterozygous c.892C>T (p.Arg298Trp) variant in the nucleus accumbens associated 1 ( NACC1 ) gene in seven affected individuals. One of the seven individuals is mosaic for this variant. NACC1 encodes a transcriptional repressor implicated in gene expression and has not previously been associated with germline disorders. The probability of finding the same missense NACC1 variant by chance in 7 out of 17,228 individuals who underwent WES for diagnoses of neurodevelopmental phenotypes is extremely small and achieves genome-wide significance (p = 1.25 × 10 −14 ). Selective constraint against missense variants in NACC1 makes this excess of an identical missense variant in all seven individuals more remarkable. Our findings are consistent with a germline recurrent mutational hotspot associated with an allele-specific neurodevelopmental phenotype in NACC1 .

Early B cell factor 3 (EBF3) is a member of the highly evolutionarily conserved Collier/Olf/EBF (COE) family of transcription factors. Prior studies on invertebrate and vertebrate animals have shown that EBF3 homologs are essential for survival and that loss-of-function mutations are associated with a range of nervous system developmental defects, including perturbation of neuronal development and migration. Interestingly, aristaless-related homeobox (ARX), a homeobox-containing transcription factor critical for the regulation of nervous system development, transcriptionally represses EBF3 expression. However, human neurodevelopmental disorders related to EBF3 have not been reported. Here, we describe three individuals who are affected by global developmental delay, intellectual disability, and expressive speech disorder and carry de novo variants in EBF3. Associated features seen in these individuals include congenital hypotonia, structural CNS malformations, ataxia, and genitourinary abnormalities. The de novo variants affect a single conserved residue in a zinc finger motif crucial for DNA binding and are deleterious in a fly model. Our findings indicate that mutations in EBF3 cause a genetic neurodevelopmental syndrome and suggest that loss of EBF3 function might mediate a subset of neurologic phenotypes shared by ARX-related disorders, including intellectual disability, abnormal genitalia, and structural CNS malformations.

Although cystic fibrosis (CF) is attributed to dysfunction of a single gene, the relationships between the abnormal gene product and the development of inflammation and progression of lung disease are not fully understood, which limits our ability to predict an individual patient’s clinical course and treatment response. To better understand CF progression, we characterized the molecular signatures of CF disease status with plasma-based functional genomics. Peripheral blood mononuclear cells (PBMCs) from healthy donors were cultured with plasma samples from CF patients ( n = 103) and unrelated, healthy controls ( n = 31). Gene expression levels were measured with an Affymetrix microarray (GeneChip Human Genome U133 Plus 2.0). Peripheral blood samples from a subset of the CF patients ( n = 40) were immunophenotyped by flow cytometry, and the data were compared with historical data for age-matched healthy controls ( n = 351). Plasma samples from another subset of CF patients ( n = 56) and healthy controls ( n = 16) were analyzed by multiplex enzyme-linked immunosorbent assay (ELISA) for numerous cytokines and chemokines. Principal component analysis and hierarchical clustering of induced transcriptional data revealed disease-specific plasma-induced PBMC profiles. Among 1,094 differentially expressed probe sets, 51 genes were associated with pancreatic sufficient status, and 224 genes were associated with infection with Pseudomonas aeruginosa. The flow cytometry and ELISA data confirmed that various immune modulators are relevant contributors to the CF molecular signature. This study provides strong evidence for distinct molecular signatures among CF patients. An understanding of these molecular signatures may lead to unique molecular markers that will enable more personalized prognoses, individualized treatment plans, and rapid monitoring of treatment response.

The balance between immune-stimulatory and immune-suppressive mechanisms in the tumour microenvironment is associated with tumour rejection and can predict the efficacy of immune checkpoint-inhibition therapies. We consider the observed differences between the transcriptional programmes associated with cancer types where the levels of immune infiltration predict a favourable prognosis versus those in which the immune infiltration predicts an unfavourable prognosis and defined a score named Mediators of Immune Response Against Cancer in soLid microEnvironments (MIRACLE). MIRACLE deconvolves T cell infiltration, from inhibitory mechanisms, such as TGFβ, EMT and PI3Kγ signatures. Our score outperforms current state-of-the-art immune signatures as a predictive marker of survival in TCGA (n = 9305, HR: 0.043, p value: 6.7 × 10−36). In a validation cohort (n = 7623), MIRACLE predicts better survival compared to other immune metrics (HR: 0.1985, p value: 2.73 × 10−38). MIRACLE also predicts response to checkpoint-inhibitor therapies (n = 333). The tumour-intrinsic factors inversely associated with the reported score such as EGFR, PRKAR1A and MAP3K1 are frequently associated with immune-suppressive phenotypes. The association of cancer outcome with the level of infiltrating immune cells is mediated by the balance of activatory and suppressive factors. MIRACLE accounts for this balance and predicts favourable cancer outcomes.

A pilot program was initiated using whole genome sequencing (WGS) to diagnose suspected genetic disorders in the Genetics Clinic at Children's Hospital of Wisconsin. Twenty-two patients underwent WGS between 2010 and 2013. Initially, we obtained a 14% (3/22) diagnosis rate over 2 years; with subsequent reanalysis, this increased to 36% (8/22). Disease causing variants were identified in SKIV2L, CECR1, DGKE, PYCR2, RYR1, PDGFRB, EFTUD2, and BCS1L. In 75% (6/8) of diagnosed cases, the diagnosis affected treatment and/or medical surveillance. Additionally, one case demonstrated a homozygous A18V variant in VLDLR that appears to be associated with a previously undescribed phenotype.

Renal preglomerular arterioles regulate vascular tone to ensure a large pressure gradient over short distances, a function that is extremely important for maintaining renal microcirculation. Regulation of renal microvascular tone is impaired in salt-sensitive (SS) hypertension-induced nephropathy, but the molecular mechanisms contributing to this impairment remain elusive. Here, we assessed the contribution of the SH2 adaptor protein p66Shc (encoded by Shc1) in regulating renal vascular tone and the development of renal vascular dysfunction associated with hypertension-induced nephropathy. We generated a panel of mutant rat strains in which specific modifications of Shc1 were introduced into the Dahl SS rats. In SS rats, overexpression of p66Shc was linked to increased renal damage. Conversely, deletion of p66Shc from these rats restored the myogenic responsiveness of renal preglomerular arterioles ex vivo and promoted cellular contraction in primary vascular smooth muscle cells (SMCs) that were isolated from renal vessels. In primary SMCs, p66Shc restricted the activation of transient receptor potential cation channels to attenuate cytosolic Ca2+ influx, implicating a mechanism by which overexpression of p66Shc impairs renal vascular reactivity. These results establish the adaptor protein p66Shc as a regulator of renal vascular tone and a driver of impaired renal vascular function in hypertension-induced nephropathy.

Pronucleus-like vesicle formation following premature chromosome condensation (PCC) of the donor cell nucleus is the key event for successful generation of cloned rodents by nuclear transplantation (NT). However in rat cloning, this change is difficult to induce in enucleated recipient oocytes because of their inability to maintain maturation-promoting factor levels. In this study, intact oocytes retrieved from nuclear-visualized H2B-tdTomato knock-in rats were injected with Venus-labeled cell nuclei. Because the incidence of PCC under MG-132 treatment significantly increased with the culture period (0%, 10.8%, 36.8%, and 87.5% at 0, 0.5, 1, and 2 h postinjection, respectively), the metaphase plate of the oocyte was removed 1-2 h after the nuclear injection. The NT-derived rat zygotes (n = 748) were activated with ionomycin/cycloheximide and transferred into temporal host mothers, resulting in the harvest of three blastocysts (0.4%) with Venus fluorescence. Two blastocysts were examined for their potential to commit to NT-derived embryonic stem cells (ntESCs). One ntESC line was established successfully and found to be competent in terms of karyotype, stem cell marker expression, and pluripotency. In conclusion, time-lagged enucleation of visualized oocyte nuclei allows the PCC incidence of donor nuclei and generation of NT blastocysts, and the blastocysts can commit to germline-competent ntESCs.

A 1.37 Mbp region of chromosome 13 previously identified by exclusion mapping was consistently associated with a reduction of salt-induced hypertension in the Dahl salt-sensitive (SS) rat. This region contained five genes that were introgressed from the salt-insensitive Brown Norway (BN) rat. The goal of the present study was to further narrow that region to identify the gene(s) most likely to protect from salt-induced hypertension. The studies yielded a subcongenic SS rat strain containing a 0.71 Mbp insert from BN (26-P strain) in which salt-induced hypertension was reduced by 24 mmHg. The region contained two protein-coding genes ( Astn1 and Pappa2) and a microRNA ( miR-488). Pappa2 mRNA in the renal cortex of the protected 26-P was 6- to 10-fold greater than in SS fed a 0.4% NaCl diet but was reduced to levels observed in SS when fed 8.0% NaCl diet for 7 days. Compared with brain nuclei (NTS, RVLM, CVLM) and the adrenal gland, Pappa2 in the renal cortex was the only gene found to be differentially expressed between SS and 26-P and that responded to changes of salt diet. Immunohistochemistry studies found Pappa2 localized in the cytosol of the epithelial cells of the cortical thick ascending limbs. In more distal segments of the renal tubules, it was observed within tubular lumens and most notably bound to the apical membranes of the intercalated cells of collecting ducts. We conclude that we have identified a variant form of Pappa2 that can protect against salt-induced hypertension in the Dahl S rat.

BACKGROUND: Many patients remain without a diagnosis despite extensive medical evaluation. The Undiagnosed Diseases Network (UDN) was established to apply a multidisciplinary model in the evaluation of the most challenging cases and to identify the biologic characteristics of newly discovered diseases. The UDN, which is funded by the National Institutes of Health, was formed in 2014 as a network of seven clinical sites, two sequencing cores, and a coordinating center. Later, a central biorepository, a metabolomics core, and a model organisms screening center were added. METHODS: We evaluated patients who were referred to the UDN over a period of 20 months. The patients were required to have an undiagnosed condition despite thorough evaluation by a health care provider. We determined the rate of diagnosis among patients who subsequently had a complete evaluation, and we observed the effect of diagnosis on medical care. RESULTS: A total of 1519 patients (53% female) were referred to the UDN, of whom 601 (40%) were accepted for evaluation. Of the accepted patients, 192 (32%) had previously undergone exome sequencing. Symptoms were neurologic in 40% of the applicants, musculoskeletal in 10%, immunologic in 7%, gastrointestinal in 7%, and rheumatologic in 6%. Of the 382 patients who had a complete evaluation, 132 received a diagnosis, yielding a rate of diagnosis of 35%. A total of 15 diagnoses (11%) were made by clinical review alone, and 98 (74%) were made by exome or genome sequencing. Of the diagnoses, 21% led to recommendations regarding changes in therapy, 37% led to changes in diagnostic testing, and 36% led to variant-specific genetic counseling. We defined 31 new syndromes. CONCLUSIONS: The UDN established a diagnosis in 132 of the 382 patients who had a complete evaluation, yielding a rate of diagnosis of 35%. (Funded by the National Institutes of Health Common Fund.)

Interferon regulatory factor 2 binding protein-like (IRF2BPL) encodes a member of the IRF2BP family of transcriptional regulators. Currently the biological function of this gene is obscure, and the gene has not been associated with a Mendelian disease. Here we describe seven individuals who carry damaging heterozygous variants in IRF2BPL and are affected with neurological symptoms. Five individuals who carry IRF2BPL nonsense variants resulting in a premature stop codon display severe neurodevelopmental regression, hypotonia, progressive ataxia, seizures, and a lack of coordination. Two additional individuals, both with missense variants, display global developmental delay and seizures and a relatively milder phenotype than those with nonsense alleles. The IRF2BPL bioinformatics signature based on population genomics is consistent with a gene that is intolerant to variation. We show that the fruit-fly IRF2BPL ortholog, called pits (protein interacting with Ttk69 and Sin3A), is broadly detected, including in the nervous system. Complete loss of pits is lethal early in development, whereas partial knockdown with RNA interference in neurons leads to neurodegeneration, revealing a requirement for this gene in proper neuronal function and maintenance. The identified IRF2BPL nonsense variants behave as severe loss-of-function alleles in this model organism, and ectopic expression of the missense variants leads to a range of phenotypes. Taken together, our results show that IRF2BPL and pits are required in the nervous system in humans and flies, and their loss leads to a range of neurological phenotypes in both species.

Purpose Sixty-75% of individuals with rare and undiagnosed phenotypes remain undiagnosed after whole exome sequencing (ES). With standard ES reanalysis resolving 10–15% of the ES negatives, further approaches are necessary to maximize diagnoses in these individuals. Methods In 38 ES negative patients an individualized genomic-phenotypic approach was employed utilizing: A) Phenotyping; B) Reanalyses of FASTQ files, with innovative bioinformatics; C) Targeted molecular testing; D) Whole genome sequencing (WGS) and E) Conferring of clinical diagnoses when pathognomonic clinical findings occurred. Results Certain and Highly Likely diagnoses were made in 18/38 (47%) individuals, including identifying two new developmental disorders. The majority of diagnoses (>70%) were due to our bioinformatics, phenotyping and targeted testing identifying variants that were undetected or not prioritized on prior ES. WGS diagnosed 3/18 individuals, with structural variants not amenable to ES. Additionally, Tentative diagnoses were made in three (8%) and in five individuals (13%) candidate genes were identified. Overall, diagnoses/potential leads were identified in 26/38 (68%). Conclusions Our comprehensive approach to ES negatives maximizes the ES and clinical data for both diagnoses and candidate gene identification, without WGS in the majority. This iterative approach is cost-effective and is pertinent to the current conundrum of ES negatives.

Background Interpreting genetic variants is one of the greatest challenges impeding analysis of rapidly increasing volumes of genomic data from patients. For example, SHROOM3 is an associated risk gene for CKD, yet causative mechanism(s) of SHROOM3 allele(s) are unknown.Methods We used our analytic pipeline that integrates genetic, computational, biochemical, CRISPR/Cas9 editing, molecular, and physiologic data to characterize coding and noncoding variants to study the human SHROOM3 risk locus for CKD.Results We identified a novel SHROOM3 transcriptional start site, which results in a shorter isoform lacking the PDZ domain and is regulated by a common noncoding sequence variant associated with CKD (rs17319721, allele frequency: 0.35). This variant disrupted allele binding to the transcription factor TCF7L2 in podocyte cell nuclear extracts and altered transcription levels of SHROOM3 in cultured cells, potentially through the loss of repressive looping between rs17319721 and the novel start site. Although common variant mechanisms are of high utility, sequencing is beginning to identify rare variants involved in disease; therefore, we used our biophysical tools to analyze an average of 112,849 individual human genome sequences for rare SHROOM3 missense variants, revealing 35 high-effect variants. The high-effect alleles include a coding variant (P1244L) previously associated with CKD (P=0.01, odds ratio=7.95; 95% CI, 1.53 to 41.46) that we find to be present in East Asian individuals at an allele frequency of 0.0027. We determined that P1244L attenuates the interaction of SHROOM3 with 14-3-3, suggesting alterations to the Hippo pathway, a known mediator of CKD.Conclusions These data demonstrate multiple new SHROOM3-dependent genetic/molecular mechanisms that likely affect CKD.

Human genome-wide association studies have linked SH2B adaptor protein 3 ( SH2B3 , LNK ) to hypertension and renal disease, although little experimental investigation has been performed to verify a role for SH2B3 in these pathologies. SH2B3, a member of the SH2B adaptor protein family, is an intracellular adaptor protein that functions as a negative regulator in many signaling pathways, including inflammatory signaling processes. To explore a mechanistic link between SH2B3 and hypertension, we targeted the SH2B3 gene for mutation on the Dahl salt-sensitive (SS) rat genetic background with zinc-finger nucleases. The resulting mutation was a 6-bp, in-frame deletion within a highly conserved region of the Src homology 2 (SH2) domain of SH2B3 . This mutation significantly attenuated Dahl SS hypertension and renal disease. Also, infiltration of leukocytes into the kidneys, a key mediator of Dahl SS pathology, was significantly blunted in the Sh2b3 em1Mcwi mutant rats. To determine whether this was because of differences in immune signaling, bone marrow transplant studies were performed in which Dahl SS and Sh2b3 em1Mcwi mutants underwent total body irradiation and were then transplanted with Dahl SS or Sh2b3 em1Mcwi mutant bone marrow. Rats that received Sh2b3 em1Mcwi mutant bone marrow had a significant reduction in mean arterial pressure and kidney injury when placed on a high salt diet (4% NaCl). These data further support a role for the immune system as a modulator of disease severity in the pathogenesis of hypertension and provide insight into inflammatory mechanisms at play in human hypertension and renal disease.

Liquid biopsies, examinations of tumor components in body fluids, have shown promise for predicting clinical outcomes. To evaluate tumor-associated genomic and genetic variations in plasma cell-free DNA (cfDNA) and their associations with treatment response and overall survival, we applied whole genome and targeted sequencing to examine the plasma cfDNAs derived from 20 patients with advanced prostate cancer. Sequencing-based genomic abnormality analysis revealed locus-specific gains or losses that were common in prostate cancer, such as 8q gains, AR amplifications, PTEN losses and TMPRSS2-ERG fusions. To estimate tumor burden in cfDNA, we developed a Plasma Genomic Abnormality (PGA) score by summing the most significant copy number variations. Cox regression analysis showed that PGA scores were significantly associated with overall survival (p < 0.04). After androgen deprivation therapy or chemotherapy, targeted sequencing showed significant mutational profile changes in genes involved in androgen biosynthesis, AR activation, DNA repair, and chemotherapy resistance. These changes may reflect the dynamic evolution of heterozygous tumor populations in response to these treatments. These results strongly support the feasibility of using non-invasive liquid biopsies as potential tools to study biological mechanisms underlying therapy-specific resistance and to predict disease progression in advanced prostate cancer.

Background— Genome-wide association studies are powerful tools for nominating pathogenic variants, but offer little insight as to how candidate genes affect disease outcome. Such is the case for SH2B adaptor protein 3 ( SH2B3 ), which is a negative regulator of multiple cytokine signaling pathways and is associated with increased risk of myocardial infarction (MI), but its role in post-MI inflammation and fibrosis is completely unknown. Methods and Results— Using an experimental model of MI (left anterior descending artery occlusion/reperfusion injury) in wild-type and Sh2b3 knockout rats (Sh2b3 em2Mcwi ), we assessed the role of Sh2b3 in post-MI fibrosis, leukocyte infiltration, angiogenesis, left ventricle contractility, and inflammatory gene expression. Compared with wild-type, Sh2b3 em2Mcwi rats had significantly increased fibrosis (2.2-fold; P 2-fold; P P em2Mcwi rats (1.7-fold; P em2Mcwi rats. In total, 12 genes were significantly elevated in the post left anterior descending artery occluded/reperfused hearts of Sh2b3 em2Mcwi rats relative to wild-type, of which 3 ( NLRP12 , CCR2 , and IFNγ ) were significantly elevated in the left ventricle of heart failure patients carrying the MI-associated rs3184504 [T] SH2B3 risk allele. Conclusions— These data demonstrate for the first time that SH2B3 is a crucial mediator of post-MI inflammation and fibrosis.

We have reported that a reduction in renal production of 20-HETE contributes to development of hypertension in Dahl salt-sensitive (SS) rats. The present study examined whether 20-HETE production is also reduced in the cerebral vasculature of SS rats and whether this impairs the myogenic response and autoregulation of cerebral blood flow (CBF). The production of 20-HETE, the myogenic response of middle cerebral arteries (MCA), and autoregulation of CBF were compared in SS, SS-5BN rats and a newly generated CYP4A1 transgenic rat. 20-HETE production was 6-fold higher in cerebral arteries of CYP4A1 and SS-5BN than in SS rats. The diameter of the MCA decreased to 70 ± 3% to 65 ± 6% in CYP4A1 and SS-5BN rats when pressure was increased from 40 to 140 mmHg. In contrast, the myogenic response of MCA isolated from SS rats did not constrict. Administration of a 20-HETE synthesis inhibitor, HET0016, abolished the myogenic response of MCA in CYP4A1 and SS-5BN rats but had no effect in SS rats. Autoregulation of CBF was impaired in SS rats compared with CYP4A1 and SS-5BN rats. Blood-brain barrier leakage was 5-fold higher in the brain of SS rats than in SS-5BN and SS.CYP4A1 rats. These findings indicate that a genetic deficiency in the formation of 20-HETE contributes to an impaired myogenic response in MCA and autoregulation of CBF in SS rats and this may contribute to vascular remodeling and cerebral injury following the onset of hypertension.

Supplementing diets with high potassium helps reduce hypertension in humans. Inwardly rectifying K+ channels Kir4.1 (Kcnj10) and Kir5.1 (Kcnj16) are highly expressed in the basolateral membrane of distal renal tubules and contribute to Na+ reabsorption and K+ secretion through the direct control of transepithelial voltage. To define the importance of Kir5.1 in blood pressure control under conditions of salt-induced hypertension, we generated a Kcnj16 knockout in Dahl salt-sensitive (SS) rats (SSKcnj16-/-). SSKcnj16-/- rats exhibited hypokalemia and reduced blood pressure, and when fed a high-salt diet (4% NaCl), experienced 100% mortality within a few days triggered by salt wasting and severe hypokalemia. Electrophysiological recordings of basolateral K+ channels in the collecting ducts isolated from SSKcnj16-/- rats revealed activity of only homomeric Kir4.1 channels. Kir4.1 expression was upregulated in SSKcnj16-/- rats, but the protein was predominantly localized in the cytosol in SSKcnj16-/- rats. Benzamil, but not hydrochlorothiazide or furosemide, rescued this phenotype from mortality on a high-salt diet. Supplementation of high-salt diet with increased potassium (2% KCl) prevented mortality in SSKcnj16-/- rats and prevented or mitigated hypertension in SSKcnj16-/- or control SS rats, respectively. Our results demonstrate that Kir5.1 channels are key regulators of renal salt handling in SS hypertension.

PurposeTo describe examples of missed pathogenic variants on whole-exome sequencing (WES) and the importance of deep phenotyping for further diagnostic testing.MethodsGuided by phenotypic information, three children with negative WES underwent targeted single-gene testing.ResultsIndividual 1 had a clinical diagnosis consistent with infantile systemic hyalinosis, although WES and a next-generation sequencing (NGS)-based ANTXR2 test were negative. Sanger sequencing of ANTXR2 revealed a homozygous single base pair insertion, previously missed by the WES variant caller software. Individual 2 had neurodevelopmental regression and cerebellar atrophy, with no diagnosis on WES. New clinical findings prompted Sanger sequencing and copy number testing of PLA2G6. A novel homozygous deletion of the noncoding exon 1 (not included in the WES capture kit) was detected, with extension into the promoter, confirming the clinical suspicion of infantile neuroaxonal dystrophy. Individual 3 had progressive ataxia, spasticity, and magnetic resonance image changes of vanishing white matter leukoencephalopathy. An NGS leukodystrophy gene panel and WES showed a heterozygous pathogenic variant in EIF2B5; no deletions/duplications were detected. Sanger sequencing of EIF2B5 showed a frameshift indel, probably missed owing to failure of alignment.ConclusionThese cases illustrate potential pitfalls of WES/NGS testing and the importance of phenotype-guided molecular testing in yielding diagnoses.

One major challenge encountered with interpreting human genetic variants is the limited understanding of the functional impact of genetic alterations on biological processes. Furthermore, there remains an unmet demand for an efficient survey of the wealth of information on human homologs in model organisms across numerous databases. To efficiently assess the large volume of publically available information, it is important to provide a concise summary of the most relevant information in a rapid user-friendly format. To this end, we created MARRVEL (model organism aggregated resources for rare variant exploration). MARRVEL is a publicly available website that integrates information from six human genetic databases and seven model organism databases. For any given variant or gene, MARRVEL displays information from OMIM, ExAC, ClinVar, Geno2MP, DGV, and DECIPHER. Importantly, it curates model organism-specific databases to concurrently display a concise summary regarding the human gene homologs in budding and fission yeast, worm, fly, fish, mouse, and rat on a single webpage. Experiment-based information on tissue expression, protein subcellular localization, biological process, and molecular function for the human gene and homologs in the seven model organisms are arranged into a concise output. Hence, rather than visiting multiple separate databases for variant and gene analysis, users can obtain important information by searching once through MARRVEL. Altogether, MARRVEL dramatically improves efficiency and accessibility to data collection and facilitates analysis of human genes and variants by cross-disciplinary integration of 18 million records available in public databases to facilitate clinical diagnosis and basic research.

Protein modeling and molecular dynamics hold a unique toolset to aide in the characterization of clinical variants that may result in disease. Not only do these techniques offer the ability to study under characterized proteins, but they do this with the speed that is needed for time-sensitive clinical cases. In this paper we retrospectively study a clinical variant in the XIAP protein, C203Y, while addressing additional variants seen in patients with similar gastrointestinal phenotypes as the C203Y mutation. In agreement with the clinical tests performed on the C203Y patient, protein modeling and molecular dynamics suggest that direct interactions with RIPK2 and Caspase3 are altered by the C203Y mutation and subsequent loss of Zn coordination in the second BIR domain of XIAP. Interestingly, the variant does not appear to alter interactions with SMAC, resulting in further damage to the caspase and NOD2 pathways. To expand the computational strategy designed when studying XIAP, we have applied the molecular modeling tools to a list of 140 variants seen in CFTR associated with cystic fibrosis, and a list of undiagnosed variants in 17 different genes. This paper shows the exciting applications of molecular modeling in the classification and characterization of genetic variants identified in next generation sequencing. Graphical abstract XIAP in Caspase 3 and NOD2 signaling pathways.

The majority of causative variants in familial breast cancer remain unknown. Of the known risk variants, most are tumor cell autonomous, and little attention has been paid yet to germline variants that may affect the tumor microenvironment. In this study, we developed a system called the Consomic Xenograft Model (CXM) to map germline variants that affect only the tumor microenvironment. In CXM, human breast cancer cells are orthotopically implanted into immunodeficient consomic strains and tumor metrics are quantified (e.g., growth, vasculogenesis, and metastasis). Because the strain backgrounds vary, whereas the malignant tumor cells do not, any observed changes in tumor progression are due to genetic differences in the nonmalignant microenvironment. Using CXM, we defined genetic variants on rat chromosome 3 that reduced relative tumor growth and hematogenous metastasis in the SS.BN3IL2Rγ consomic model compared with the SSIL2Rγ parental strain. Paradoxically, these effects occurred despite an increase in the density of tumor-associated blood vessels. In contrast, lymphatic vasculature and lymphogenous metastasis were unaffected by the SS.BN3IL2Rγ background. Through comparative mapping and whole-genome sequence analysis, we narrowed candidate variants on rat chromosome 3 to six genes with a priority for future analysis. Collectively, our results establish the utility of CXM to localize genetic variants affecting the tumor microenvironment that underlie differences in breast cancer risk. Cancer Res; 74(22); 6419–29. ©2014 AACR.

Genome-wide association studies (GWAS) identify regions of the genome correlated with disease risk but are restricted in their ability to identify the underlying causative mechanism(s). Thus, GWAS are useful “roadmaps” that require functional analysis to establish the genetic and mechanistic structure of a particular locus. Unfortunately, direct functional testing in humans is limited, demonstrating the need for complementary approaches. Here we used an integrated approach combining zebrafish, rat, and human data to interrogate the function of an established GWAS locus (SHROOM3) lacking prior functional support for chronic kidney disease (CKD). Congenic mapping and sequence analysis in rats suggested Shroom3 was a strong positional candidate gene. Transferring a 6.1-Mb region containing the wild-type Shroom3 gene significantly improved the kidney glomerular function in FHH (fawn-hooded hypertensive) rat. The wild-type Shroom3 allele, but not the FHH Shroom3 allele, rescued glomerular defects induced by knockdown of endogenous shroom3 in zebrafish, suggesting that the FHH Shroom3 allele is defective and likely contributes to renal injury in the FHH rat. We also show for the first time that variants disrupting the actin-binding domain of SHROOM3 may cause podocyte effacement and impairment of the glomerular filtration barrier.

Previously, we found that transferring 6.1 Mb of salt-sensitive (SS) chromosome 12 (13.4–19.5 Mb) onto the consomic SS-12 BN background significantly elevated mean arterial pressure in response to an 8% NaCl diet (178±7 versus 144±2 mm Hg; P BN consomic, the 830-kb blood pressure locus was associated with a ∆+15 mm Hg ( P P P P P Chst12 (carbohydrate chondroitin 4 sulfotransferase 12) was most consistently downregulated by 2.6- to 4.5-fold ( P Chst12 that are unique to the hypertensive SS strain when compared with the normotensive Brown Norway, Dahl salt-resistant, and Wistar-Kyoto strains. Finally, the 830-kb interval is syntenic to a region on human chromosome 7 that has been genetically linked to blood pressure, suggesting that insight gained from our SS-12 BN congenic strain may be translated to a better understanding of human hypertension.

Motivation: RNA-Seq (also called whole-transcriptome sequencing) is an emerging technology that uses the capabilities of next-generation sequencing to detect and quantify entire transcripts. One of its important applications is the improvement of existing genome annotations. RNA-Seq provides rapid, comprehensive and cost-effective tools for the discovery of novel genes and transcripts compared with expressed sequence tag (EST), which is instrumental in gene discovery and gene sequence determination. The rat is widely used as a laboratory disease model, but has a less well-annotated genome as compared with humans and mice. In this study, we incorporated deep RNA-Seq data from three rat tissues—bone marrow, brain and kidney—with EST data to improve the annotation of the rat genome. Results: Our analysis identified 32 197 transcripts, including 13 461 known transcripts, 13 934 novel isoforms and 4802 new genes, which almost doubled the numbers of transcripts in the current public rat genome database (rn5). Comparisons of our predicted protein-coding gene sets with those in public datasets suggest that RNA-Seq significantly improves genome annotation and identifies novel genes and isoforms in the rat. Importantly, the large majority of novel genes and isoforms are supported by direct evidence of RNA-Seq experiments. These predicted genes were integrated into the Rat Genome Database (RGD) and can serve as an important resource for functional studies in the research community. Availability and implementation: The predicted genes are available at http://rgd.mcw.edu . Contact: hmxu@zju.edu.cn or pliu@mcw.edu or yanlu76@zju.edu.cn Supplementary information: Supplementary data are available at Bioinformatics online.

Chromosome 8q24 locus contains regulatory variants that modulate genetic risk to various cancers including prostate cancer (PC). However, the biological mechanism underlying this regulation is not well understood. Here, we developed a chromosome conformation capture (3C)-based multi-target sequencing technology and systematically examined three PC risk regions at the 8q24 locus and their potential regulatory targets across human genome in six cell lines. We observed frequent physical contacts of this risk locus with multiple genomic regions, in particular, inter-chromosomal interaction with CD96 at 3q13 and intra-chromosomal interaction with MYC at 8q24. We identified at least five interaction hot spots within the predicted functional regulatory elements at the 8q24 risk locus. We also found intra-chromosomal interaction genes PVT1, FAM84B and GSDMC and inter-chromosomal interaction gene CXorf36 in most of the six cell lines. Other gene regions appeared to be cell line-specific, such as RRP12 in LNCaP, USP14 in DU-145 and SMIN3 in lymphoblastoid cell line. We further found that the 8q24 functional domains more likely interacted with genomic regions containing genes enriched in critical pathways such as Wnt signaling and promoter motifs such as E2F1 and TCF3. This result suggests that the risk locus may function as a regulatory hub by physical interactions with multiple genes important for prostate carcinogenesis. Further understanding genetic effect and biological mechanism of these chromatin interactions will shed light on the newly discovered regulatory role of the risk locus in PC etiology and progression.

Endothelial progenitor cells (EPCs) promote angiogenesis, and clinical trials suggest autologous EPC-based therapy may be effective in treatment of vascular diseases. Albeit promising, variability in the efficacy of EPCs associated with underlying disease states has hindered the realization of EPC-based therapy. Here we first identify and characterize EPC dysfunction in a rodent model of vascular disease (SS/Mcwi rat) that exhibits impaired angiogenesis. To identify molecular candidates that mediate the angiogenic potential of these cells, we performed a broad analysis of cell surface protein expression using chemical labeling combined with mass spectrometry. Analysis revealed EPCs derived from SS/Mcwi rats express significantly more type 2 low-affinity immunoglobulin Fc-gamma (FCGR2) and natural killer 2B4 (CD244) receptors compared with controls. Genome-wide sequencing (RNA-seq) and qt-PCR confirmed isoforms of CD244 and FCGR2a transcripts were increased in SS/Mcwi EPCs. EPCs with elevated expression of FCGR2a and CD244 receptors are predicted to increase the probability of SS/Mcwi EPCs being targeted for death, providing a mechanistic explanation for their reduced angiogenic efficacy in vivo. Pathway analysis supported this contention, as “key” molecules annotated to cell death paths were differentially expressed in the SS/Mcwi EPCs. We speculate that screening and neutralization of cell surface proteins that “tag” and impair EPC function may provide an alternative approach to utilizing incompetent EPCs in greater numbers, as circulating EPCs are depleted in patients with vascular disease. Overall, novel methods to identify putative targets for repair of EPCs using discovery-based technologies will likely provide a major advance in the field of regenerative medicine.

The rat has been widely used as a disease model in a laboratory setting, resulting in an abundance of genetic and phenotype data from a wide variety of studies. These data can be found at the Rat Genome Database (RGD, http://rgd.mcw.edu/ ), which provides a platform for researchers interested in linking genomic variations to phenotypes. Quantitative trait loci (QTLs) form one of the earliest and core datasets, allowing researchers to identify loci harboring genes associated with disease. These QTLs are not only important for those using the rat to identify genes and regions associated with disease, but also for cross-organism analyses of syntenic regions on the mouse and the human genomes to identify potential regions for study in these organisms. Currently, RGD has data on >1,900 rat QTLs that include details about the methods and animals used to determine the respective QTL along with the genomic positions and markers that define the region. RGD also curates human QTLs (>1,900) and houses >4,000 mouse QTLs (imported from Mouse Genome Informatics). Multiple ontologies are used to standardize traits, phenotypes, diseases, and experimental methods to facilitate queries, analyses, and cross-organism comparisons. QTLs are visualized in tools such as GBrowse and GViewer, with additional tools for analysis of gene sets within QTL regions. The QTL data at RGD provide valuable information for the study of mapped phenotypes and identification of candidate genes for disease associations.

Genome-wide association studies (GWAS) are useful for nominating candidate genes, but typically are unable to establish disease causality or differentiate between the effects of variants in linkage disequilibrium (LD). Additionally, some GWAS loci might contain multiple causative variants or genes that contribute to the overall disease susceptibility at a single locus. However, the majority of current GWAS lack the statistical power to test whether multiple causative genes underlie the same locus, prompting us to adopt an alternative approach to testing multiple GWAS genes empirically. We used gene targeting in a disease-susceptible rat model of genetic hypertension to test all six genes at the Agtrap-Plod1 locus (Agtrap, Mthfr, Clcn6, Nppa, Nppb, and Plod1) for blood pressure (BP) and renal phenotypes. This revealed that the majority of genes at this locus (five out of six) can impact hypertension by modifying BP and renal phenotypes. Mutations of Nppa, Plod1, and Mthfr increased disease susceptibility, whereas Agtrap and Clcn6 mutations decreased hypertension risk. Reanalysis of the human AGTRAP-PLOD1 locus also implied that disease-associated haplotype blocks with polygenic effects were not only possible, but rather were highly plausible. Combined, these data demonstrate for the first time that multiple modifiers of hypertension can cosegregate at a single GWAS locus.

Summary Large numbers of inbred laboratory rat strains have been developed for a range of complex disease phenotypes. To gain insights into the evolutionary pressures underlying selection for these phenotypes, we sequenced the genomes of 27 rat strains, including 11 models of hypertension, diabetes, and insulin resistance, along with their respective control strains. Altogether, we identified more than 13 million single-nucleotide variants, indels, and structural variants across these rat strains. Analysis of strain-specific selective sweeps and gene clusters implicated genes and pathways involved in cation transport, angiotensin production, and regulators of oxidative stress in the development of cardiovascular disease phenotypes in rats. Many of the rat loci that we identified overlap with previously mapped loci for related traits in humans, indicating the presence of shared pathways underlying these phenotypes in rats and humans. These data represent a step change in resources available for evolutionary analysis of complex traits in disease models. PaperClip, Graphical Abstract, Highlights • Genomes of 27 rat strains were sequenced; >13 million sequence variants identified • Selective sweeps and coevolved gene clusters were detected in 11 disease models • Previously identified and new disease genes and pathways were identified • This is first evolutionary analysis of artificial selection for disease phenotypes, Evolution analysis of artificial selection for disease phenotypes, such as hypertension and diabetes, in 27 rat strains reveals disease-related variants and loci.

Renin was the first blood pressure (BP) quantitative trait locus mapped by linkage analysis in the rat. Subsequent BP linkage and congenic studies capturing different portions of the renin region have returned conflicting results, suggesting that multiple interdependent BP loci may be residing in the chromosome 13 BP quantitative trait locus that includes Renin . We used SS-13 BN congenic strains to map 2 BP loci in the Renin region (chr13: 45.2–49.0 Mb). We identified a 1.1-Mb protective Brown Norway region around Renin (chr13: 46.1–47.2 Mb) that significantly decreased BP by 32 mm Hg. The Renin protective BP locus was offset by an adjacent hypertensive locus (chr13: 47.2–49.0 Mb) that significantly increased BP by 29 mm Hg. Sequence analysis of the protective and hypertensive BP loci revealed 1433 and 2063 variants between Dahl salt-sensitive/Mcwi and Brown Norway rats, respectively. To further reduce the list of candidate variants, we regenotyped an overlapping SS-13 SR congenic strain (S/ren rr ) with a previously reported BP phenotype. Sequence comparison among Dahl salt-sensitive, Dahl R, and Brown Norway reduced the number of candidate variants in the 2 BP loci by 42% for further study. Combined with previous studies, these data suggest that at least 4 BP loci reside within the 30-cM chromosome 13 BP quantitative trait locus that includes Renin .

Hypertension and renal damage in Dahl SS rats are associated with increased infiltrating immune cells in the kidney. To examine the role of infiltrating immune cells in this disease process, a zinc finger nuclease targeting bases 672–706 of recombination-activating gene 1 (Rag1) was injected into the pronucleus of Dahl SS (SS/JrHsdMcwi) strain embryos and implanted in pseudopregnant females. This strategy yielded a rat strain with a 13-base frame-shift mutation in the target region of Rag1 and a deletion of immunoreactive Rag1 protein in the thymus. Flow cytometry demonstrated that the Rag1-null mutant rats have a significant reduction in T and B lymphocytes in the circulation and spleen. Studies were performed on SS and Rag1-null rats fed a 4.0% NaCl diet for 3 wk. The infiltration of T cells into the kidney following high-salt intake was significantly blunted in the Rag1-null rats (1.7 ± 0.6 × 105 cells/kidney) compared with the Dahl SS (5.6 ± 0.9 × 105 cells/kidney). Accompanying the reduction in infiltration of immune cells in the kidney, mean arterial blood pressure and urinary albumin excretion rate were significantly lower in Rag1-null mutants (158 ± 3 mmHg and 60 ± 16 mg/day, respectively) than in SS rats (180 ± 11 mmHg and 251 ± 37 mg/day). Finally, a histological analysis revealed that the glomerular and tubular damage in the kidneys of the SS rats fed a high-salt diet was also attenuated in the Rag1 mutants. These studies demonstrate the importance of renal infiltration of immune cells in the pathogenesis of hypertension and renal damage in Dahl SS rats.