Genome sequencing reveals loci under artificial selection that underlie disease phenotypes in the laboratory rat.

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.

Clinical heterogeneity of mitochondrial NAD kinase deficiency caused by a NADK2 start loss variant.

Mitochondrial NAD kinase deficiency (NADK2D, OMIM #615787) is a rare autosomal recessive disorder of NADPH biosynthesis that can cause hyperlysinemia and dienoyl-CoA reductase deficiency (DECRD, OMIM #616034). NADK2 deficiency has been reported in only three unrelated patients. Two had severe, unremitting disease; one died at 4 months and the other at 5 years of age. The third was a 10 year old female with CNS anomalies, ataxia, and incoordination. In two cases mutations in NADK2 have been demonstrated. Here, we report the fourth known case, a 15 year old female with normal intelligence and a mild clinical and biochemical phenotype presumably without DECRD. Her clinical symptoms, which are now stable, became evident at the age of 9 with the onset of decreased visual acuity, bilateral optic atrophy, nystagmus, episodic lower extremity weakness, peripheral neuropathy, and gait abnormalities. Plasma amino acid levels were within normal limits except for mean lysine and proline levels that were 3.7 and 2.5 times the upper limits of normal. Whole exome sequencing (WES) revealed homozygosity for a g.36241900 A>G p. Met1Val start loss mutation in the primary NADK2 transcript (NM_001085411.1) encoding the 442 amino acid isoform. This presumed hypomorphic mutation has not been previously reported and is absent from the v1000GP, EVS, and ExAC databases. Our patient's normal intelligence and stable disease expands the clinical heterogeneity and the prognosis associated with NADK2 deficiency. Our findings also clarify the mechanism underlying NADK2 deficiency and suggest that this disease should be ruled out in cases of hyperlysinemia, especially those with visual loss, and neurological phenotypes.

Identification of somatic mutations in non-small cell lung carcinomas using whole-exome sequencing.

Lung cancer is the leading cause of cancer-related death, with non-small cell lung cancer (NSCLC) being the predominant form of the disease. Most lung cancer is caused by the accumulation of genomic alterations due to tobacco exposure. To uncover its mutational landscape, we performed whole-exome sequencing in 31 NSCLCs and their matched normal tissue samples. We identified both common and unique mutation spectra and pathway activation in lung adenocarcinomas and squamous cell carcinomas, two major histologies in NSCLC. In addition to identifying previously known lung cancer genes (TP53, KRAS, EGFR, CDKN2A and RB1), the analysis revealed many genes not previously implicated in this malignancy. Notably, a novel gene CSMD3 was identified as the second most frequently mutated gene (next to TP53) in lung cancer. We further demonstrated that loss of CSMD3 results in increased proliferation of airway epithelial cells. The study provides unprecedented insights into mutational processes, cellular pathways and gene networks associated with lung cancer. Of potential immediate clinical relevance, several highly mutated genes identified in our study are promising druggable targets in cancer therapy including ALK, CTNNA3, DCC, MLL3, PCDHIIX, PIK3C2B, PIK3CG and ROCK2.

Making a definitive diagnosis: successful clinical application of whole exome sequencing in a child with intractable inflammatory bowel disease.

PURPOSE: We report a male child who presented at 15 months with perianal abscesses and proctitis, progressing to transmural pancolitis with colocutaneous fistulae, consistent with a Crohn disease-like illness. The age and severity of the presentation suggested an underlying immune defect; however, despite comprehensive clinical evaluation, we were unable to arrive at a definitive diagnosis, thereby restricting clinical management. METHODS: We sought to identify the causative mutation(s) through exome sequencing to provide the necessary additional information required for clinical management. RESULTS: After sequencing, we identified 16,124 variants. Subsequent analysis identified a novel, hemizygous missense mutation in the X-linked inhibitor of apoptosis gene, substituting a tyrosine for a highly conserved and functionally important cysteine. X-linked inhibitor of apoptosis was not previously associated with Crohn disease but has a central role in the proinflammatory response and bacterial sensing through the NOD signaling pathway. The mutation was confirmed by Sanger sequencing in a licensed clinical laboratory. Functional assays demonstrated an increased susceptibility to activation-induced cell death and defective responsiveness to NOD2 ligands, consistent with loss of normal X-linked inhibitor of apoptosis protein function in apoptosis and NOD2 signaling. CONCLUSIONS: Based on this medical history, genetic and functional data, the child was diagnosed as having an X-linked inhibitor of apoptosis deficiency. Based on this finding, an allogeneic hematopoietic progenitor cell transplant was performed to prevent the development of life-threatening hemophagocytic lymphohistiocytosis, in concordance with the recommended treatment for X-linked inhibitor of apoptosis deficiency. At >42 days posttransplant, the child was able to eat and drink, and there has been no recurrence of gastrointestinal disease, suggesting this mutation also drove the gastrointestinal disease. This report describes the identification of a novel cause of inflammatory bowel disease. Equally importantly, it demonstrates the power of exome sequencing to render a molecular diagnosis in an individual patient in the setting of a novel disease, after all standard diagnoses were exhausted, and illustrates how this technology can be used in a clinical setting.

Genetically hypertensive Brown Norway congenic rat strains suggest intermediate traits underlying genetic hypertension.

AIM: To determine the independent and combined effects of three quantitative trait loci (QTL) for blood pressure in the Genetically Hypertensive (GH/Omr) rat by generating and characterizing single and combined congenic strains that have QTL on rat chromosomes (RNO) 2, 6, and 18 from the GH rat introduced into a hypertension resistant Brown Norway (BN) background. METHODS: Linkage analysis and QTL identification (genome wide QTL scan) were performed with MapMaker/EXP to build the genetic maps and MapMaker/QTL for linking the phenotypes to the genetic map. The congenic strains were derived using marker-assisted selection strategy from a single male F1 offspring of an intercross between the male GH/Omr and female BN/Elh, followed by 10 generations of selective backcrossing to the female BN progenitor strain. Single congenic strains generated were BN.GH-(D2Rat22-D2Mgh11)/Mcwi (BN.GH2); BN.GH-(D6Mit12-D6Rat15)/Mcwi (BN.GH6); and BN.GH-(D18Rat41-D18Mgh4)/Mcwi (BN.GH18). Blood pressure measurements were obtained either via a catheter placed in the femoral artery or by radiotelemetry. Responses to angiotensin II (ANGII), norepinephrine (NE), and baroreceptor sensitivity were measured in the single congenics. RESULTS: Transferring one or more QTL from the hypertensive GH into normotensive BN strain was not sufficient to cause hypertension in any of the developed congenic strains. There were no differences between the parental and congenic strains in their response to NE. However, BN.GH18 rats revealed significantly lower baroreceptor sensitivity (beta=-1.25-/+0.17), whereas BN.GH2 (beta=0.66-/+0.09) and BN.GH18 (beta=0.71-/+0.07) had significantly decreased responses to ANGII from those observed in the BN (beta=0.88-/+0.08). CONCLUSION: The failure to alter blood pressure levels by introducing the hypertensive QTL from the GH into the hypertension resistant BN background suggests that the QTL effects are genome background-dependent in the GH rat. BN.GH2 and BN.GH18 rats reveal significant differences in response to ANGII and impaired baroreflex sensitivity, suggesting that we may have captured a locus responsible for the genetic control of baroreceptor sensitivity, which would be considered an intermediate phenotype of blood pressure.

De novo repeat interruptions are associated with reduced somatic instability and mild or absent clinical features in myotonic dystrophy type 1.

Myotonic dystrophy type 1 (DM1) is a multisystem disorder, caused by expansion of a CTG trinucleotide repeat in the 3'-untranslated region of the DMPK gene. The repeat expansion is somatically unstable and tends to increase in length with time, contributing to disease progression. In some individuals, the repeat array is interrupted by variant repeats such as CCG and CGG, stabilising the expansion and often leading to milder symptoms. We have characterised three families, each including one person with variant repeats that had arisen de novo on paternal transmission of the repeat expansion. Two individuals were identified for screening due to an unusual result in the laboratory diagnostic test, and the third due to exceptionally mild symptoms. The presence of variant repeats in all three expanded alleles was confirmed by restriction digestion of small pool PCR products, and allele structures were determined by PacBio sequencing. Each was different, but all contained CCG repeats close to the 3'-end of the repeat expansion. All other family members had inherited pure CTG repeats. The variant repeat-containing alleles were more stable in the blood than pure alleles of similar length, which may in part account for the mild symptoms observed in all three individuals. This emphasises the importance of somatic instability as a disease mechanism in DM1. Further, since patients with variant repeats may have unusually mild symptoms, identification of these individuals has important implications for genetic counselling and for patient stratification in DM1 clinical trials.

Plasma exosomal miRNAs-based prognosis in metastatic kidney cancer.

Plasma exosomal miRNAs were evaluated for prognosis in an initial set of 44 metastatic renal cell cancer (mRCC) patients by RNA sequencing. Among ∼3.49 million mappable reads per patient, miRNAs accounted for 93.1% of the mapped RNAs. 227 miRNAs with high abundance were selected for survival analysis. Cox regression analysis identified association of 6 miRNAs with overall survival (OS) (P<0.01, False discovery rate (FDR) < 0.3). Five of the associated miRNAs were quantified in an independent follow-up cohort of 65 mRCC patients by TaqMan-based miRNA assays. Kaplan-Meier analysis confirmed the significant OS association of three miRs; miR-let-7i-5p (P=0.018, HR=0.49, 95% CI=0.21-0.84), miR-26a-1-3p (P=0.025, HR=0.43, 95% CI=0.10-0.84) and miR-615-3p (P=0.0007, HR=0.36, 95% CI=0.11-0.54). A multivariate analysis of miR-let-7i-5p with the clinical factor-based Memorial Sloan-Kettering Cancer Center (MSKCC) score improved survival prediction from an area under the curve (AUC) of 0.58 for MSKCC score to an average AUC of 0.64 across 48-month follow-up time. The multivariate model was able to define a high-risk group with median survival of 14 months and low risk group of 39 months (P=0.0002, HR=3.43, 95%CI, 2.73-24.15). Further validation of miRNA-based prognostic biomarkers are needed to improve current clinic-pathologic based prognostic models in patients with mRCC.

Plasma genetic and genomic abnormalities predict treatment response and clinical outcome in advanced prostate cancer.

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.

Utilization of a labeled tracking oligonucleotide for visualization and quality control of spotted 70-mer arrays.

BACKGROUND: Spotted 70-mer oligonucleotide arrays offer potentially greater specificity and an alternative to expensive cDNA library maintenance and amplification. Since microarray fabrication is a considerable source of data variance, we previously directly tagged cDNA probes with a third fluorophore for prehybridization quality control. Fluorescently modifying oligonucleotide sets is cost prohibitive, therefore, a co-spotted Staphylococcus aureus-specific fluorescein-labeled "tracking" oligonucleotide is described to monitor fabrication variables of a Mycobacterium tuberculosis oligonucleotide microarray. RESULTS: Significantly (p < 0.01) improved DNA retention was achieved printing in 15% DMSO/1.5 M betaine compared to the vendor recommended buffers. Introduction of tracking oligonucleotide did not effect hybridization efficiency or introduce ratio measurement bias in hybridizations between M. tuberculosis H37Rv and M. tuberculosis mprA. Linearity between the mean log Cy3/Cy5 ratios of genes differentially expressed from arrays either possessing or lacking the tracking oligonucleotide was observed (R2 = 0.90, p < 0.05) and there were no significant differences in Pearson's correlation coefficients of ratio data between replicates possessing (0.72 +/- 0.07), replicates lacking (0.74 +/- 0.10), or replicates with and without (0.70 +/- 0.04) the tracking oligonucleotide. ANOVA analysis confirmed the tracking oligonucleotide introduced no bias. Titrating target-specific oligonucleotide (40 microM to 0.78 microM) in the presence of 0.5 microM tracking oligonucleotide, revealed a fluorescein fluorescence inversely related to target-specific oligonucleotide molarity, making tracking oligonucleotide signal useful for quality control measurements and differentiating false negatives (synthesis failures and mechanical misses) from true negatives (no gene expression). CONCLUSIONS: This novel approach enables prehybridization array visualization for spotted oligonucleotide arrays and sets the stage for more sophisticated slide qualification and data filtering applications.

Genomics in clinical practice: lessons from the front lines.

The price of whole-genome and -exome sequencing has fallen to the point where these methods can be applied to clinical medicine. Here, we outline the lessons we have learned in converting a sequencing laboratory designed for research into a fully functional clinical program.

Mapping the genetic determinants of hypertension, metabolic diseases, and related phenotypes in the lyon hypertensive rat.

The complex nature of hypertension makes identifying the pathophysiology and its genetic contributions a challenging task. One powerful approach for the genetic dissection of blood pressure regulation is studying inbred rat models of hypertension, as they provide natural allele variants but reduced heterogeneity (both genetic and etiologic). Furthermore, the detailed physiologic studies to which the rat is amenable allow for the determination of intermediate phenotypes. We have performed a total genome scan in offspring of an F2 intercross between the Lyon hypertensive (LH) and Lyon normotensive rat strains to identify linkage of anthropometric, blood pressure, renal, metabolic, and endocrine phenotypes. Quantitative trait locus (QTL) regions involved in blood pressure regulation, end-stage organ damage, body and organ weight, and lipid metabolism in the LH rat were identified on chromosomes 1, 2, 3, 5, 7, 10, 13, and 17, with 2 phenotypes associated with the metabolic syndrome identified on chromosomes 1 and 17. Regions on chromosomes 2, 13, and 17 were revealed to be important for blood pressure regulation. Regions on chromosome 17 were found to significantly contribute to both metabolic homeostasis and blood pressure regulation; 2 aggregates of a total of 23 QTLs were identified, including several "intermediate phenotypes." These intermediate phenotypes may be used as closer surrogates to the mechanisms leading to hypertension and metabolic dysfunction in the LH rat.