Rare variants create synthetic genome-wide associations.

Genome-wide association studies (GWAS) have now identified at least 2,000 common variants that appear associated with common diseases or related traits (http://www.genome.gov/gwastudies), hundreds of which have been convincingly replicated. It is generally thought that the associated markers reflect the effect of a nearby common (minor allele frequency >0.05) causal site, which is associated with the marker, leading to extensive resequencing efforts to find causal sites. We propose as an alternative explanation that variants much less common than the associated one may create "synthetic associations" by occurring, stochastically, more often in association with one of the alleles at the common site versus the other allele. Although synthetic associations are an obvious theoretical possibility, they have never been systematically explored as a possible explanation for GWAS findings. Here, we use simple computer simulations to show the conditions under which such synthetic associations will arise and how they may be recognized. We show that they are not only possible, but inevitable, and that under simple but reasonable genetic models, they are likely to account for or contribute to many of the recently identified signals reported in genome-wide association studies. We also illustrate the behavior of synthetic associations in real datasets by showing that rare causal mutations responsible for both hearing loss and sickle cell anemia create genome-wide significant synthetic associations, in the latter case extending over a 2.5-Mb interval encompassing scores of "blocks" of associated variants. In conclusion, uncommon or rare genetic variants can easily create synthetic associations that are credited to common variants, and this possibility requires careful consideration in the interpretation and follow up of GWAS signals.

Rare hereditary COL4A3/COL4A4 variants may be mistaken for familial focal segmental glomerulosclerosis.

Focal segmental glomerulosclerosis (FSGS) is a histological lesion with many causes, including inherited genetic defects, with significant proteinuria being the predominant clinical finding at presentation. Mutations in COL4A3 and COL4A4 are known to cause Alport syndrome (AS), thin basement membrane nephropathy, and to result in pathognomonic glomerular basement membrane (GBM) findings. Secondary FSGS is known to develop in classic AS at later stages of the disease. Here, we present seven families with rare or novel variants in COL4A3 or COL4A4 (six with single and one with two heterozygous variants) from a cohort of 70 families with a diagnosis of hereditary FSGS. The predominant clinical finding at diagnosis was proteinuria associated with hematuria. In all seven families, there were individuals with nephrotic-range proteinuria with histologic features of FSGS by light microscopy. In one family, electron microscopy showed thin GBM, but four other families had variable findings inconsistent with classical Alport nephritis. There was no recurrence of disease after kidney transplantation. Families with COL4A3 and COL4A4 variants that segregated with disease represent 10% of our cohort. Thus, COL4A3 and COL4A4 variants should be considered in the interpretation of next-generation sequencing data from such patients. Furthermore, this study illustrates the power of molecular genetic diagnostics in the clarification of renal phenotypes.

Visualizing Rare Watson-Crick-Like Tautomeric and Anionic Mismatches in DNA and RNA

The central dogma of molecular biology relies on the correct Watson-Crick (WC) geometry of canonical deoxyribonucleic acid (DNA) dG•dC and dA•dT base pairs to replicate and transcribe genetic information with speed and an astonishing level of fidelity. In addition, the Watson-Crick geometry of canonical ribonucleic acid (RNA) rG•rC and rA•rU base pairs is highly conserved to ensure that proteins are translated with high fidelity. However, numerous other potential nucleobase tautomeric and ionic configurations are possible that can give rise to entirely new pairing modes between the nucleotide bases. Very early on, James Watson and Francis Crick recognized their importance and in 1953 postulated that if bases adopted one of their less energetically disfavored tautomeric forms (and later ionic forms) during replication it could lead to the formation of a mismatch with a Watson-Crick-like geometry and could give rise to “natural mutations.”

Since this time numerous studies have provided evidence in support of this hypothesis and have expanded upon it; computational studies have addressed the energetic feasibilities of different nucleobases’ tautomeric and ionic forms in siico; crystallographic studies have trapped different mismatches with WC-like geometries in polymerase or ribosome active sites. However, no direct evidence has been given for (i) the direct existence of these WC-like mismatches in canonical DNA duplex, RNA duplexes, or non-coding RNAs; (ii) which, if any, tautomeric or ionic form stabilizes the WC-like geometry. This thesis utilizes nuclear magnetic resonance (NMR) spectroscopy and rotating frame relaxation dispersion (R1ρ RD) in combination with density functional theory (DFT), biochemical assays, and targeted chemical perturbations to show that (i) dG•dT mismatches in DNA duplexes, as well as rG•rU mismatches RNA duplexes and non-coding RNAs, transiently adopt a WC-like geometry that is stabilized by (ii) an interconnected network of rapidly interconverting rare tautomers and anionic bases. These results support Watson and Crick’s tautomer hypothesis, but additionally support subsequent hypotheses invoking anionic mismatches and ultimately tie them together. This dissertation shows that a common mismatch can adopt a Watson-Crick-like geometry globally, in both DNA and RNA, and whose geometry is stabilized by a kinetically linked network of rare tautomeric and anionic bases. The studies herein also provide compelling evidence for their involvement in spontaneous replication and translation errors.

VCP/p97 is essential for maturation of ubiquitin-containing autophagosomes and this function is impaired by mutations that cause IBMPFD.

VCP (VCP/p97) is a ubiquitously expressed member of the AAA(+)-ATPase family of chaperone-like proteins that regulates numerous cellular processes including chromatin decondensation, homotypic membrane fusion and ubiquitin-dependent protein degradation by the proteasome. Mutations in VCP cause a multisystem degenerative disease consisting of inclusion body myopathy, Paget disease of bone, and frontotemporal dementia (IBMPFD). Here we show that VCP is essential for autophagosome maturation. We generated cells stably expressing dual-tagged LC3 (mCherry-EGFP-LC3) which permit monitoring of autophagosome maturation. We determined that VCP deficiency by RNAi-mediated knockdown or overexpression of dominant-negative VCP results in significant accumulation of immature autophagic vesicles, some of which are abnormally large, acidified and exhibit cathepsin B activity. Furthermore, expression of disease-associated VCP mutants (R155H and A232E) also causes this autophagy defect. VCP was found to be essential to autophagosome maturation under basal conditions and in cells challenged by proteasome inhibition, but not in cells challenged by starvation, suggesting that VCP might be selectively required for autophagic degradation of ubiquitinated substrates. Indeed, a high percentage of the accumulated autophagic vesicles contain ubiquitin-positive contents, a feature that is not observed in autophagic vesicles that accumulate following starvation or treatment with Bafilomycin A. Finally, we show accumulation of numerous, large LAMP-1 and LAMP-2-positive vacuoles and accumulation of LC3-II in myoblasts derived from patients with IBMPFD. We conclude that VCP is essential for maturation of ubiquitin-containing autophagosomes and that defect in this function may contribute to IBMPFD pathogenesis.

Individuals with mutations in XPNPEP3, which encodes a mitochondrial protein, develop a nephronophthisis-like nephropathy.

The autosomal recessive kidney disease nephronophthisis (NPHP) constitutes the most frequent genetic cause of terminal renal failure in the first 3 decades of life. Ten causative genes (NPHP1-NPHP9 and NPHP11), whose products localize to the primary cilia-centrosome complex, support the unifying concept that cystic kidney diseases are "ciliopathies". Using genome-wide homozygosity mapping, we report here what we believe to be a new locus (NPHP-like 1 [NPHPL1]) for an NPHP-like nephropathy. In 2 families with an NPHP-like phenotype, we detected homozygous frameshift and splice-site mutations, respectively, in the X-prolyl aminopeptidase 3 (XPNPEP3) gene. In contrast to all known NPHP proteins, XPNPEP3 localizes to mitochondria of renal cells. However, in vivo analyses also revealed a likely cilia-related function; suppression of zebrafish xpnpep3 phenocopied the developmental phenotypes of ciliopathy morphants, and this effect was rescued by human XPNPEP3 that was devoid of a mitochondrial localization signal. Consistent with a role for XPNPEP3 in ciliary function, several ciliary cystogenic proteins were found to be XPNPEP3 substrates, for which resistance to N-terminal proline cleavage resulted in attenuated protein function in vivo in zebrafish. Our data highlight an emerging link between mitochondria and ciliary dysfunction, and suggest that further understanding the enzymatic activity and substrates of XPNPEP3 will illuminate novel cystogenic pathways.

Expression signatures of TP53 mutations in serous ovarian cancers.

BACKGROUND: Mutations in the TP53 gene are extremely common and occur very early in the progression of serous ovarian cancers. Gene expression patterns that relate to mutational status may provide insight into the etiology and biology of the disease. METHODS: The TP53 coding region was sequenced in 89 frozen serous ovarian cancers, 40 early stage (I/II) and 49 advanced stage (III/IV). Affymetrix U133A expression data was used to define gene expression patterns by mutation, type of mutation, and cancer stage. RESULTS: Missense or chain terminating (null) mutations in TP53 were found in 59/89 (66%) ovarian cancers. Early stage cancers had a significantly higher rate of null mutations than late stage disease (38% vs. 8%, p < 0.03). In advanced stage cases, mutations were more prevalent in short term survivors than long term survivors (81% vs. 30%, p = 0.0004). Gene expression patterns had a robust ability to predict TP53 status within training data. By using early versus late stage disease for out of sample predictions, the signature derived from early stage cancers could accurately (86%) predict mutation status of late stage cancers. CONCLUSIONS: This represents the first attempt to define a genomic signature of TP53 mutation in ovarian cancer. Patterns of gene expression characteristic of TP53 mutation could be discerned and included several genes that are known p53 targets or have been described in the context of expression signatures of TP53 mutation in breast cancer.

Altered ultrasonic vocalization and impaired learning and memory in Angelman syndrome mouse model with a large maternal deletion from Ube3a to Gabrb3.

Angelman syndrome (AS) is a neurobehavioral disorder associated with mental retardation, absence of language development, characteristic electroencephalography (EEG) abnormalities and epilepsy, happy disposition, movement or balance disorders, and autistic behaviors. The molecular defects underlying AS are heterogeneous, including large maternal deletions of chromosome 15q11-q13 (70%), paternal uniparental disomy (UPD) of chromosome 15 (5%), imprinting mutations (rare), and mutations in the E6-AP ubiquitin ligase gene UBE3A (15%). Although patients with UBE3A mutations have a wide spectrum of neurological phenotypes, their features are usually milder than AS patients with deletions of 15q11-q13. Using a chromosomal engineering strategy, we generated mutant mice with a 1.6-Mb chromosomal deletion from Ube3a to Gabrb3, which inactivated the Ube3a and Gabrb3 genes and deleted the Atp10a gene. Homozygous deletion mutant mice died in the perinatal period due to a cleft palate resulting from the null mutation in Gabrb3 gene. Mice with a maternal deletion (m-/p+) were viable and did not have any obvious developmental defects. Expression analysis of the maternal and paternal deletion mice confirmed that the Ube3a gene is maternally expressed in brain, and showed that the Atp10a and Gabrb3 genes are biallelically expressed in all brain sub-regions studied. Maternal (m-/p+), but not paternal (m+/p-), deletion mice had increased spontaneous seizure activity and abnormal EEG. Extensive behavioral analyses revealed significant impairment in motor function, learning and memory tasks, and anxiety-related measures assayed in the light-dark box in maternal deletion but not paternal deletion mice. Ultrasonic vocalization (USV) recording in newborns revealed that maternal deletion pups emitted significantly more USVs than wild-type littermates. The increased USV in maternal deletion mice suggests abnormal signaling behavior between mothers and pups that may reflect abnormal communication behaviors in human AS patients. Thus, mutant mice with a maternal deletion from Ube3a to Gabrb3 provide an AS mouse model that is molecularly more similar to the contiguous gene deletion form of AS in humans than mice with Ube3a mutation alone. These mice will be valuable for future comparative studies to mice with maternal deficiency of Ube3a alone.

Rare targets are rarely missed in correctable search.

Failing to find a tumor in an x-ray scan or a gun in an airport baggage screening can have dire consequences, making it fundamentally important to elucidate the mechanisms that hinder performance in such visual searches. Recent laboratory work has indicated that low target prevalence can lead to disturbingly high miss rates in visual search. Here, however, we demonstrate that misses in low-prevalence searches can be readily abated. When targets are rarely present, observers adapt by responding more quickly, and miss rates are high. Critically, though, these misses are often due to response-execution errors, not perceptual or identification errors: Observers know a target was present, but just respond too quickly. When provided an opportunity to correct their last response, observers can catch their mistakes. Thus, low target prevalence may not be a generalizable cause of high miss rates in visual search.

Mutations in NGLY1 cause an inherited disorder of the endoplasmic reticulum-associated degradation pathway.

PURPOSE: The endoplasmic reticulum-associated degradation pathway is responsible for the translocation of misfolded proteins across the endoplasmic reticulum membrane into the cytosol for subsequent degradation by the proteasome. To define the phenotype associated with a novel inherited disorder of cytosolic endoplasmic reticulum-associated degradation pathway dysfunction, we studied a series of eight patients with deficiency of N-glycanase 1. METHODS: Whole-genome, whole-exome, or standard Sanger sequencing techniques were employed. Retrospective chart reviews were performed in order to obtain clinical data. RESULTS: All patients had global developmental delay, a movement disorder, and hypotonia. Other common findings included hypolacrima or alacrima (7/8), elevated liver transaminases (6/7), microcephaly (6/8), diminished reflexes (6/8), hepatocyte cytoplasmic storage material or vacuolization (5/6), and seizures (4/8). The nonsense mutation c.1201A>T (p.R401X) was the most common deleterious allele. CONCLUSION: NGLY1 deficiency is a novel autosomal recessive disorder of the endoplasmic reticulum-associated degradation pathway associated with neurological dysfunction, abnormal tear production, and liver disease. The majority of patients detected to date carry a specific nonsense mutation that appears to be associated with severe disease. The phenotypic spectrum is likely to enlarge as cases with a broader range of mutations are detected.

Evaluation of Altered Kras Codon Bias and NOS Inhibition During Lung Tumorigenesis

The small GTPases HRAS, NRAS and KRAS are mutated in approximately one-third of all human cancers, rendering the proteins constitutively active and oncogenic. Lung cancer is the leading cause of cancer deaths worldwide, and more than 20% of human lung cancers harbor mutations in RAS, with 98% of those occurring in the KRAS isoform. While there have been many advances in the understanding of KRAS–driven lung tumorigenesis, it remains a therapeutic challenge. To further this understanding and assess novel approaches for treatment, I have investigated two aspects of Kras–driven tumorigenesis in the lung:

(I) Despite nearly identical protein sequences, the three RAS proto-oncogenes exhibit divergent codon usage. Of the three isoforms, KRAS contains the most rare codons resulting in lower levels of KRAS protein expression relative to HRAS and NRAS. To determine the consequences of rare codon bias during de novo tumorigenesis, we created a knock-in Krasex3op mouse in which synonymous mutations in exon 3 converted codons from rare to common. These mice had reduced tumor burden and fewer oncogenic mutations in the Krasex3op allele following carcinogen exposure. The reduction in tumorigenesis appeared to be a product of rare codons affecting both the oncogenic and non–oncogenic alleles. Converting rare codons to common codons yielded a more potent oncogenic allele that promoted growth arrest and enhanced tumor suppression by the non-oncogenic allele. Thus, rare codons play an integral role in Kras tumorigenesis.

(II) Lung cancer patients exhale higher levels of NO and iNOS-/- mice are resistant to chemically induced lung tumorigenesis. I hypothesize that NO promotes Kras–driven lung adenocarcinoma, and NOS inhibition may decrease Kras–driven lung tumorigenesis. To test this hypothesis, I assessed efficacy of the NOS inhibitor L–NAME in a genetically engineered mouse model of Kras-driven lung adenocarcinoma. Adenoviral Cre recombinase was delivered into the lungs intranasally, resulting in expression of oncogenic KrasG12D and dominant-negative Trp53R172H in lung epithelial cells. L–NAME treatment was provided in the water and continued until survival endpoints. In this model, L–NAME treatment decreased tumor growth and prolonged survival. These data establish a potential clinical role for NOS inhibition in lung cancer treatment.

Mitochondrial mutations in adenoid cystic carcinoma of the salivary glands.

BACKGROUND: The MitoChip v2.0 resequencing array is an array-based technique allowing for accurate and complete sequencing of the mitochondrial genome. No studies have investigated mitochondrial mutation in salivary gland adenoid cystic carcinomas. METHODOLOGY: The entire mitochondrial genome of 22 salivary gland adenoid cystic carcinomas (ACC) of salivary glands and matched leukocyte DNA was sequenced to determine the frequency and distribution of mitochondrial mutations in ACC tumors. PRINCIPAL FINDINGS: Seventeen of 22 ACCs (77%) carried mitochondrial mutations, ranging in number from 1 to 37 mutations. A disproportionate number of mutations occurred in the D-loop. Twelve of 17 tumors (70.6%) carried mutations resulting in amino acid changes of translated proteins. Nine of 17 tumors (52.9%) with a mutation carried an amino acid changing mutation in the nicotinamide adenine dinucleotide dehydrogenase (NADH) complex. CONCLUSIONS/SIGNIFICANCE: Mitochondrial mutation is frequent in salivary ACCs. The high incidence of amino acid changing mutations implicates alterations in aerobic respiration in ACC carcinogenesis. D-loop mutations are of unclear significance, but may be associated with alterations in transcription or replication.

Dominant Splice Site Mutations in PIK3R1 Cause Hyper IgM Syndrome, Lymphadenopathy and Short Stature.

The purpose of this research was to use next generation sequencing to identify mutations in patients with primary immunodeficiency diseases whose pathogenic gene mutations had not been identified. Remarkably, four unrelated patients were found by next generation sequencing to have the same heterozygous mutation in an essential donor splice site of PIK3R1 (NM_181523.2:c.1425 + 1G > A) found in three prior reports. All four had the Hyper IgM syndrome, lymphadenopathy and short stature, and one also had SHORT syndrome. They were investigated with in vitro immune studies, RT-PCR, and immunoblotting studies of the mutation's effect on mTOR pathway signaling. All patients had very low percentages of memory B cells and class-switched memory B cells and reduced numbers of naïve CD4+ and CD8+ T cells. RT-PCR confirmed the presence of both an abnormal 273 base-pair (bp) size and a normal 399 bp size band in the patient and only the normal band was present in the parents. Following anti-CD40 stimulation, patient's EBV-B cells displayed higher levels of S6 phosphorylation (mTOR complex 1 dependent event), Akt phosphorylation at serine 473 (mTOR complex 2 dependent event), and Akt phosphorylation at threonine 308 (PI3K/PDK1 dependent event) than controls, suggesting elevated mTOR signaling downstream of CD40. These observations suggest that amino acids 435-474 in PIK3R1 are important for its stability and also its ability to restrain PI3K activity. Deletion of Exon 11 leads to constitutive activation of PI3K signaling. This is the first report of this mutation and immunologic abnormalities in SHORT syndrome.