Copy Number Variation Analysis Facilitates Identification of Genetic Causation in Patients with Congenital Anomalies of the Kidney and Urinary Tract

Background

Congenital anomalies of the kidneys and urinary tract (CAKUT) are the most common cause of chronic kidney disease among children and adults younger than 30 yr. In our previous study, whole-exome sequencing (WES) identified a known monogenic cause of isolated or syndromic CAKUT in 13% of families with CAKUT. However, WES has limitations and detection of copy number variations (CNV) is technically challenging, and CNVs causative of CAKUT have previously been detected in up to 16% of cases.

Objective

To detect CNVs causing CAKUT in this WES cohort and increase the diagnostic yield.

Design setting and participants

We performed a genome-wide single nucleotide polymorphism (SNP)-based CNV analysis on the same CAKUT cohort for whom WES was previously conducted.

Outcome measurements and statistical analysis

We evaluated and classified the CNVs using previously published predefined criteria.

Results and limitations

In a cohort of 170 CAKUT families, we detected a pathogenic CNV known to cause CAKUT in nine families (5.29%, 9/170). There were no competing variants on genome-wide CNV analysis or WES analysis. In addition, we identified novel likely pathogenic CNVs that may cause a CAKUT phenotype in three of the 170 families (1.76%).

Conclusions

CNV analysis in this cohort of 170 CAKUT families previously examined via WES increased the rate of diagnosis of genetic causes of CAKUT from 13% on WES to 18% on WES + CNV analysis combined. We also identified three candidate loci that may potentially cause CAKUT.

Patient summary

We conducted a genetics study on families with congenital anomalies of the kidney and urinary tract (CAKUT). We identified gene mutations that can explain CAKUT symptoms in 5.29% of the families, which increased the percentage of genetic causes of CAKUT to 18% from a previous study, so roughly one in five of our patients with CAKUT had a genetic cause. These analyses can help patients with CAKUT and their families in identifying a possible genetic cause.

Reverse phenotyping facilitates disease allele calling in exome sequencing of patients with CAKUT

PURPOSE

Congenital anomalies of the kidneys and urinary tract (CAKUT) constitute the leading cause of chronic kidney disease in children. In total, 174 monogenic causes of isolated or syndromic CAKUT are known. However, syndromic features may be overlooked when the initial clinical diagnosis of CAKUT is made. We hypothesized that the yield of a molecular genetic diagnosis by exome sequencing (ES) can be increased by applying reverse phenotyping, by re-examining the case for signs/symptoms of the suspected clinical syndrome that results from the genetic variant detected by ES.

METHODS

We conducted ES in an international cohort of 731 unrelated families with CAKUT. We evaluated ES data for variants in 174 genes, in which variants are known to cause isolated or syndromic CAKUT. In cases in which ES suggested a previously unreported syndromic phenotype, we conducted reverse phenotyping.

RESULTS

In 83 of 731 (11.4%) families, we detected a likely CAKUT-causing genetic variant consistent with an isolated or syndromic CAKUT phenotype. In 19 of these 83 families (22.9%), reverse phenotyping yielded syndromic clinical findings, thereby strengthening the genotype-phenotype correlation.

CONCLUSION

We conclude that employing reverse phenotyping in the evaluation of syndromic CAKUT genes by ES provides an important tool to facilitate molecular genetic diagnostics in CAKUT.

Whole exome sequencing identifies potential candidate genes for spina bifida derived from mouse models

Spina bifida (SB) is the second most common nonlethal congenital malformation. The existence of monogenic SB mouse models and human monogenic syndromes with SB features indicate that human SB may be caused by monogenic genes. We hypothesized that whole exome sequencing (WES) allows identification of potential candidate genes by (i) generating a list of 136 candidate genes for SB, and (ii) by unbiased exome-wide analysis. We generated a list of 136 potential candidate genes from three categories and evaluated WES data of 50 unrelated SB cases for likely deleterious variants in 136 potential candidate genes, and for potential SB candidate genes exome-wide. We identified 6 likely deleterious variants in 6 of the 136 potential SB candidate genes in 6 of the 50 SB cases, whereof 4 genes were derived from mouse models, 1 gene was derived from human nonsyndromic SB, and 1 gene was derived from candidate genes known to cause human syndromic SB. In addition, by unbiased exome-wide analysis, we identified 12 genes as potential candidates for SB. Identification of these 18 potential candidate genes in larger SB cohorts will help decide which ones can be considered as novel monogenic causes of human SB.

Bi-allelic variants in SPATA5L1 lead to intellectual disability, spastic-dystonic cerebral palsy, epilepsy, and hearing loss

Spermatogenesis-associated 5 like 1 (SPATA5L1) represents an orphan gene encoding a protein of unknown function. We report 28 bi-allelic variants in SPATA5L1 associated with sensorineural hearing loss in 47 individuals from 28 (26 unrelated) families. In addition, 25/47 affected individuals (53%) presented with microcephaly, developmental delay/intellectual disability, cerebral palsy, and/or epilepsy. Modeling indicated damaging effect of variants on the protein, largely via destabilizing effects on protein domains. Brain imaging revealed diminished cerebral volume, thin corpus callosum, and periventricular leukomalacia, and quantitative volumetry demonstrated significantly diminished white matter volumes in several individuals. Immunofluorescent imaging in rat hippocampal neurons revealed localization of Spata5l1 in neuronal and glial cell nuclei and more prominent expression in neurons. In the rodent inner ear, Spata5l1 is expressed in the neurosensory hair cells and inner ear supporting cells. Transcriptomic analysis performed with fibroblasts from affected individuals was able to distinguish affected from controls by principal components. Analysis of differentially expressed genes and networks suggested a role for SPATA5L1 in cell surface adhesion receptor function, intracellular focal adhesions, and DNA replication and mitosis. Collectively, our results indicate that bi-allelic SPATA5L1 variants lead to a human disease characterized by sensorineural hearing loss (SNHL) with or without a nonprogressive mixed neurodevelopmental phenotype.

Mutations in PRDM15 Are a Novel Cause of Galloway-Mowat Syndrome

BACKGROUND

Galloway-Mowat syndrome (GAMOS) is characterized by neurodevelopmental defects and a progressive nephropathy, which typically manifests as steroid-resistant nephrotic syndrome. The prognosis of GAMOS is poor, and the majority of children progress to renal failure. The discovery of monogenic causes of GAMOS has uncovered molecular pathways involved in the pathogenesis of disease.

METHODS

Homozygosity mapping, whole-exome sequencing, and linkage analysis were used to identify mutations in four families with a GAMOS-like phenotype, and high-throughput PCR technology was applied to 91 individuals with GAMOS and 816 individuals with isolated nephrotic syndrome. In vitro and in vivo studies determined the functional significance of the mutations identified.

RESULTS

Three biallelic variants of the transcriptional regulator PRDM15 were detected in six families with proteinuric kidney disease. Four families with a variant in the protein's zinc-finger (ZNF) domain have additional GAMOS-like features, including brain anomalies, cardiac defects, and skeletal defects. All variants destabilize the PRDM15 protein, and the ZNF variant additionally interferes with transcriptional activation. Morpholino oligonucleotide-mediated knockdown of Prdm15 in Xenopus embryos disrupted pronephric development. Human wild-type PRDM15 RNA rescued the disruption, but the three PRDM15 variants did not. Finally, CRISPR-mediated knockout of PRDM15 in human podocytes led to dysregulation of several renal developmental genes.

CONCLUSIONS

Variants in PRDM15 can cause either isolated nephrotic syndrome or a GAMOS-type syndrome on an allelic basis. PRDM15 regulates multiple developmental kidney genes, and is likely to play an essential role in renal development in humans.

Mutations of the Transcriptional Corepressor ZMYM2 Cause Syndromic Urinary Tract Malformations

Congenital anomalies of the kidney and urinary tract (CAKUT) constitute one of the most frequent birth defects and represent the most common cause of chronic kidney disease in the first three decades of life. Despite the discovery of dozens of monogenic causes of CAKUT, most pathogenic pathways remain elusive. We performed whole-exome sequencing (WES) in 551 individuals with CAKUT and identified a heterozygous de novo stop-gain variant in ZMYM2 in two different families with CAKUT. Through collaboration, we identified in total 14 different heterozygous loss-of-function mutations in ZMYM2 in 15 unrelated families. Most mutations occurred de novo, indicating possible interference with reproductive function. Human disease features are replicated in X. tropicalis larvae with morpholino knockdowns, in which expression of truncated ZMYM2 proteins, based on individual mutations, failed to rescue renal and craniofacial defects. Moreover, heterozygous Zmym2-deficient mice recapitulated features of CAKUT with high penetrance. The ZMYM2 protein is a component of a transcriptional corepressor complex recently linked to the silencing of developmentally regulated endogenous retrovirus elements. Using protein-protein interaction assays, we show that ZMYM2 interacts with additional epigenetic silencing complexes, as well as confirming that it binds to FOXP1, a transcription factor that has also been linked to CAKUT. In summary, our findings establish that loss-of-function mutations of ZMYM2, and potentially that of other proteins in its interactome, as causes of human CAKUT, offering new routes for studying the pathogenesis of the disorder.

CAKUT and Autonomic Dysfunction Caused by Acetylcholine Receptor Mutations

Congenital anomalies of the kidney and urinary tract (CAKUT) are the most common cause of chronic kidney disease in the first three decades of life, and in utero obstruction to urine flow is a frequent cause of secondary upper urinary tract malformations. Here, using whole-exome sequencing, we identified three different biallelic mutations in CHRNA3, which encodes the α3 subunit of the nicotinic acetylcholine receptor, in five affected individuals from three unrelated families with functional lower urinary tract obstruction and secondary CAKUT. Four individuals from two families have additional dysautonomic features, including impaired pupillary light reflexes. Functional studies in vitro demonstrated that the mutant nicotinic acetylcholine receptors were unable to generate current following stimulation with acetylcholine. Moreover, the truncating mutations p.Thr337Asnfs^∗81 and p.Ser340^∗ led to impaired plasma membrane localization of CHRNA3. Although the importance of acetylcholine signaling in normal bladder function has been recognized, we demonstrate for the first time that mutations in CHRNA3 can cause bladder dysfunction, urinary tract malformations, and dysautonomia. These data point to a pathophysiologic sequence by which monogenic mutations in genes that regulate bladder innervation may secondarily cause CAKUT.

Novel homozygous ENPP1 mutation causes generalized arterial calcifications of infancy, thrombocytopenia, and cardiovascular and central nervous system syndrome

Generalized arterial calcifications of infancy (GACI) is caused by mutations in ENPP1. Other ENPP1-related phenotypes include pseudoxanthoma elasticum, hypophosphatemic rickets, and Cole disease. We studied four children from two Bedouin consanguineous families who presented with severe clinical phenotype including thrombocytopenia, hypoglycemia, hepatic, and neurologic manifestations. Initial working diagnosis included congenital infection; however, patients remained without a definitive diagnosis despite extensive workup. Consequently, we investigated a potential genetic etiology. Whole exome sequencing (WES) was performed for affected children and their parents. Following the identification of a novel mutation in the ENPP1 gene, we characterized this novel multisystemic presentation and revised relevant imaging studies. Using WES, we identified a novel homozygous mutation (c.556G > C; p.Gly186Arg) in ENPP1 which affects a highly conserved protein domain (somatomedin B2). ENPP1-associated genetic diseases exhibit phenotypic heterogeneity depending on mutation type and location. Follow-up clinical characterization of these families allowed us to revise and detect new features of systemic calcifications, which established the diagnosis of GACI, expanding the phenotypic spectrum associated with ENPP1 mutations. Our findings demonstrate that this novel ENPP1 founder mutation can cause a fatal multisystemic phenotype, mimicking severe congenital infection. This also represents the first reported mutation affecting the SMB2 domain, associated with GACI.

Whole-Exome Sequencing Enables a Precision Medicine Approach for Kidney Transplant Recipients

BACKGROUND

Whole-exome sequencing (WES) finds a CKD-related mutation in approximately 20% of patients presenting with CKD before 25 years of age. Although provision of a molecular diagnosis could have important implications for clinical management, evidence is lacking on the diagnostic yield and clinical utility of WES for pediatric renal transplant recipients.

METHODS

To determine the diagnostic yield of WES in pediatric kidney transplant recipients, we recruited 104 patients who had received a transplant at Boston Children's Hospital from 2007 through 2017, performed WES, and analyzed results for likely deleterious variants in approximately 400 genes known to cause CKD.

RESULTS

By WES, we identified a genetic cause of CKD in 34 out of 104 (32.7%) transplant recipients. The likelihood of detecting a molecular genetic diagnosis was highest for patients with urinary stone disease (three out of three individuals), followed by renal cystic ciliopathies (seven out of nine individuals), steroid-resistant nephrotic syndrome (nine out of 21 individuals), congenital anomalies of the kidney and urinary tract (ten out of 55 individuals), and chronic glomerulonephritis (one out of seven individuals). WES also yielded a molecular diagnosis for four out of nine individuals with ESRD of unknown etiology. The WES-related molecular genetic diagnosis had implications for clinical care for five patients.

CONCLUSIONS

Nearly one third of pediatric renal transplant recipients had a genetic cause of their kidney disease identified by WES. Knowledge of this genetic information can help guide management of both transplant patients and potential living related donors.

Clonal evolution analysis of paired anaplastic and well-differentiated thyroid carcinomas reveals shared common ancestor

Foci of papillary or follicular thyroid carcinoma are frequently noted in thyroidectomy specimens of anaplastic thyroid carcinoma (ATC). However, whether ATCs evolve from these co-existing well-differentiated thyroid carcinomas (WDTCs) has not been well-understood. To investigate the progression of ATC in patients with co-existing WDTCs, five ATC tumors with co-existing WDTCs and matching normal tissues were whole-exome sequenced. After mapping the somatic alteration landscape, evolutionary lineages were constructed by sub-clone analysis. Though each tumor harbored at least some unique private mutations, all five ATCs demonstrated numerous overlapping mutations with matched WDTCs. Clonal analysis further demonstrated that each ATC/WDTC pair shared a common ancestor, with some pairs diverging early in their evolution and others in which the ATC seems to arise directly from a sub-clone of the WDTC. Though the precise lineal relationship remains ambiguous, based on the genetic relationship, our study clearly suggests a shared origin of ATC and WDTC.

Whole-Exome Sequencing Identifies Causative Mutations in Families with Congenital Anomalies of the Kidney and Urinary Tract

BACKGROUND

Congenital anomalies of the kidney and urinary tract (CAKUT) are the most prevalent cause of kidney disease in the first three decades of life. Previous gene panel studies showed monogenic causation in up to 12% of patients with CAKUT.

METHODS

We applied whole-exome sequencing to analyze the genotypes of individuals from 232 families with CAKUT, evaluating for mutations in single genes known to cause human CAKUT and genes known to cause CAKUT in mice. In consanguineous or multiplex families, we additionally performed a search for novel monogenic causes of CAKUT.

RESULTS

In 29 families (13%), we detected a causative mutation in a known gene for isolated or syndromic CAKUT that sufficiently explained the patient's CAKUT phenotype. In three families (1%), we detected a mutation in a gene reported to cause a phenocopy of CAKUT. In 15 of 155 families with isolated CAKUT, we detected deleterious mutations in syndromic CAKUT genes. Our additional search for novel monogenic causes of CAKUT in consanguineous and multiplex families revealed a potential single, novel monogenic CAKUT gene in 19 of 232 families (8%).

CONCLUSIONS

We identified monogenic mutations in a known human CAKUT gene or CAKUT phenocopy gene as the cause of disease in 14% of the CAKUT families in this study. Whole-exome sequencing provides an etiologic diagnosis in a high fraction of patients with CAKUT and will provide a new basis for the mechanistic understanding of CAKUT.

Molecular and cellular reorganization of neural circuits in the human lineage

To better understand the molecular and cellular differences in brain organization between human and nonhuman primates, we performed transcriptome sequencing of 16 regions of adult human, chimpanzee, and macaque brains. Integration with human single-cell transcriptomic data revealed global, regional, and cell-type-specific species expression differences in genes representing distinct functional categories. We validated and further characterized the human specificity of genes enriched in distinct cell types through histological and functional analyses, including rare subpallial-derived interneurons expressing dopamine biosynthesis genes enriched in the human striatum and absent in the nonhuman African ape neocortex. Our integrated analysis of the generated data revealed diverse molecular and cellular features of the phylogenetic reorganization of the human brain across multiple levels, with relevance for brain function and disease.

Whole Exome Sequencing Reveals a Monogenic Cause of Disease in ≈43% of 35 Families With Midaortic Syndrome

Midaortic syndrome (MAS) is a rare cause of severe childhood hypertension characterized by narrowing of the abdominal aorta in children and is associated with extensive vascular disease. It may occur as part of a genetic syndrome, such as neurofibromatosis, or as consequence of a pathological inflammatory disease. However, most cases are considered idiopathic. We hypothesized that in a high percentage of these patients, a monogenic cause of disease may be detected by evaluating whole exome sequencing data for mutations in 1 of 38 candidate genes previously described to cause vasculopathy. We studied a cohort of 36 individuals from 35 different families with MAS by exome sequencing. In 15 of 35 families (42.9%), we detected likely causal dominant mutations. In 15 of 35 (42.9%) families with MAS, whole exome sequencing revealed a mutation in one of the genes previously associated with vascular disease (NF1, JAG1, ELN, GATA6, and RNF213). Ten of the 15 mutations have not previously been reported. This is the first report of ELN, RNF213, or GATA6 mutations in individuals with MAS. Mutations were detected in NF1 (6/15 families), JAG1 (4/15 families), ELN (3/15 families), and one family each for GATA6 and RNF213 Eight individuals had syndromic disease and 7 individuals had isolated MAS. Whole exome sequencing can provide conclusive molecular genetic diagnosis in a high fraction of individuals with syndromic or isolated MAS. Establishing an etiologic diagnosis may reveal genotype/phenotype correlations for MAS in the future and should, therefore, be performed routinely in MAS.

A homozygous missense variant in VWA2, encoding an interactor of the Fraser-complex, in a patient with vesicoureteral reflux

Congenital anomalies of the kidney and urinary tract (CAKUT) are the most common cause (40-50%) of chronic kidney disease (CKD) in children. About 40 monogenic causes of CAKUT have so far been discovered. To date less than 20% of CAKUT cases can be explained by mutations in these 40 genes. To identify additional monogenic causes of CAKUT, we performed whole exome sequencing (WES) and homozygosity mapping (HM) in a patient with CAKUT from Indian origin and consanguineous descent. We identified a homozygous missense mutation (c.1336C>T, p.Arg446Cys) in the gene Von Willebrand factor A domain containing 2 (VWA2). With immunohistochemistry studies on kidneys of newborn (P1) mice, we show that Vwa2 and Fraser extracellular matrix complex subunit 1 (Fras1) co-localize in the nephrogenic zone of the renal cortex. We identified a pronounced expression of Vwa2 in the basement membrane of the ureteric bud (UB) and derivatives of the metanephric mesenchyme (MM). By applying in vitro assays, we demonstrate that the Arg446Cys mutation decreases translocation of monomeric VWA2 protein and increases translocation of aggregated VWA2 protein into the extracellular space. This is potentially due to the additional, unpaired cysteine residue in the mutated protein that is used for intermolecular disulfide bond formation. VWA2 is a known, direct interactor of FRAS1 of the Fraser-Complex (FC). FC-encoding genes and interacting proteins have previously been implicated in the pathogenesis of syndromic and/or isolated CAKUT phenotypes in humans. VWA2 therefore constitutes a very strong candidate in the search for novel CAKUT-causing genes. Our results from in vitro experiments indicate a dose-dependent neomorphic effect of the Arg446Cys homozygous mutation in VWA2.

Early Assessment of Lung Cancer Immunotherapy Response via Circulating Tumor DNA

Purpose: Decisions to continue or suspend therapy with immune checkpoint inhibitors are commonly guided by tumor dynamics seen on serial imaging. However, immunotherapy responses are uniquely challenging to interpret because tumors often shrink slowly or can appear transiently enlarged due to inflammation. We hypothesized that monitoring tumor cell death in real time by quantifying changes in circulating tumor DNA (ctDNA) levels could enable early assessment of immunotherapy efficacy.Experimental Design: We compared longitudinal changes in ctDNA levels with changes in radiographic tumor size and with survival outcomes in 28 patients with metastatic non-small cell lung cancer (NSCLC) receiving immune checkpoint inhibitor therapy. CtDNA was quantified by determining the allele fraction of cancer-associated somatic mutations in plasma using a multigene next-generation sequencing assay. We defined a ctDNA response as a >50% decrease in mutant allele fraction from baseline, with a second confirmatory measurement.Results: Strong agreement was observed between ctDNA response and radiographic response (Cohen's kappa, 0.753). Median time to initial response among patients who achieved responses in both categories was 24.5 days by ctDNA versus 72.5 days by imaging. Time on treatment was significantly longer for ctDNA responders versus nonresponders (median, 205.5 vs. 69 days; P < 0.001). A ctDNA response was associated with superior progression-free survival [hazard ratio (HR), 0.29; 95% CI, 0.09-0.89; P = 0.03], and superior overall survival (HR, 0.17; 95% CI, 0.05-0.62; P = 0.007).Conclusions: A drop in ctDNA level is an early marker of therapeutic efficacy and predicts prolonged survival in patients treated with immune checkpoint inhibitors for NSCLC. Clin Cancer Res; 24(8); 1872-80. ©2018 AACR.

Contribution of rare inherited and de novo variants in 2,871 congenital heart disease probands

Congenital heart disease (CHD) is the leading cause of mortality from birth defects. Exome sequencing of a single cohort of 2,871 CHD probands including 2,645 parent-offspring trios implicated rare inherited mutations in 1.8%, including a recessive founder mutation in GDF1 accounting for ~5% of severe CHD in Ashkenazim, recessive genotypes in MYH6 accounting for ~11% of Shone complex, and dominant FLT4 mutations accounting for 2.3% of Tetralogy of Fallot. De novo mutations (DNMs) accounted for 8% of cases, including ~3% of isolated CHD patients and ~28% with both neurodevelopmental and extra-cardiac congenital anomalies. Seven genes surpassed thresholds for genome-wide significance and 12 genes not previously implicated in CHD had > 70% probability of being disease-related; DNMs in ~440 genes are inferred to contribute to CHD. There was striking overlap between genes with damaging DNMs in probands with CHD and autism.

A genome-wide association study of autism using the Simons Simplex Collection: Does reducing phenotypic heterogeneity in autism increase genetic homogeneity?

BACKGROUND

Phenotypic heterogeneity in autism has long been conjectured to be a major hindrance to the discovery of genetic risk factors, leading to numerous attempts to stratify children based on phenotype to increase power of discovery studies. This approach, however, is based on the hypothesis that phenotypic heterogeneity closely maps to genetic variation, which has not been tested. Our study examines the impact of subphenotyping of a well-characterized autism spectrum disorder (ASD) sample on genetic homogeneity and the ability to discover common genetic variants conferring liability to ASD.

METHODS

Genome-wide genotypic data of 2576 families from the Simons Simplex Collection were analyzed in the overall sample and phenotypic subgroups defined on the basis of diagnosis, IQ, and symptom profiles. We conducted a family-based association study, as well as estimating heritability and evaluating allele scores for each phenotypic subgroup.

RESULTS

Association analyses revealed no genome-wide significant association signal. Subphenotyping did not increase power substantially. Moreover, allele scores built from the most associated single nucleotide polymorphisms, based on the odds ratio in the full sample, predicted case status in subsets of the sample equally well and heritability estimates were very similar for all subgroups.

CONCLUSIONS

In genome-wide association analysis of the Simons Simplex Collection sample, reducing phenotypic heterogeneity had at most a modest impact on genetic homogeneity. Our results are based on a relatively small sample, one with greater homogeneity than the entire population; if they apply more broadly, they imply that analysis of subphenotypes is not a productive path forward for discovering genetic risk variants in ASD.

Mutations in KATNB1 cause complex cerebral malformations by disrupting asymmetrically dividing neural progenitors

Exome sequencing analysis of over 2,000 children with complex malformations of cortical development identified five independent (four homozygous and one compound heterozygous) deleterious mutations in KATNB1, encoding the regulatory subunit of the microtubule-severing enzyme Katanin. Mitotic spindle formation is defective in patient-derived fibroblasts, a consequence of disrupted interactions of mutant KATNB1 with KATNA1, the catalytic subunit of Katanin, and other microtubule-associated proteins. Loss of KATNB1 orthologs in zebrafish (katnb1) and flies (kat80) results in microcephaly, recapitulating the human phenotype. In the developing Drosophila optic lobe, kat80 loss specifically affects the asymmetrically dividing neuroblasts, which display supernumerary centrosomes and spindle abnormalities during mitosis, leading to cell cycle progression delays and reduced cell numbers. Furthermore, kat80 depletion results in dendritic arborization defects in sensory and motor neurons, affecting neural architecture. Taken together, we provide insight into the mechanisms by which KATNB1 mutations cause human cerebral cortical malformations, demonstrating its fundamental role during brain development.

Whole-exome sequencing characterizes the landscape of somatic mutations and copy number alterations in adrenocortical carcinoma

CONTEXT

Adrenocortical carcinoma (ACC) is a rare and lethal malignancy with a poorly defined etiology, and the molecular genetics of ACC are incompletely understood.

OBJECTIVE

To utilize whole-exome sequencing for genetic characterization of the underlying somatic mutations and copy number alterations present in ACC.

DESIGN

Screening for somatic mutation events and copy number alterations (CNAs) was performed by comparative analysis of tumors and matched normal samples from 41 patients with ACC.

RESULTS

In total, 966 nonsynonymous somatic mutations were detected, including 40 tumors with a mean of 16 mutations per sample and one tumor with 314 mutations. Somatic mutations in ACC-associated genes included TP53 (8/41 tumors, 19.5%) and CTNNB1 (4/41, 9.8%). Genes with potential disease-causing mutations included GNAS, NF2, and RB1, and recurrently mutated genes with unknown roles in tumorigenesis comprised CDC27, SCN7A, and SDK1. Recurrent CNAs included amplification at 5p15.33 including TERT (6/41, 14.6%) and homozygous deletion at 22q12.1 including the Wnt repressors ZNRF3 and KREMEN1 (4/41 9.8% and 3/41, 7.3%, respectively). Somatic mutations in ACC-established genes and recurrent ZNRF3 and TERT loci CNAs were mutually exclusive in the majority of cases. Moreover, gene ontology identified Wnt signaling as the most frequently mutated pathway in ACCs.

CONCLUSIONS

These findings highlight the importance of Wnt pathway dysregulation in ACC and corroborate the finding of homozygous deletion of Wnt repressors ZNRF3 and KREMEN1. Overall, mutations in either TP53 or CTNNB1 as well as focal CNAs at the ZNRF3 or TERT loci denote mutually exclusive events, suggesting separate mechanisms underlying the development of these tumors.

The contribution of de novo coding mutations to autism spectrum disorder

Whole exome sequencing has proven to be a powerful tool for understanding the genetic architecture of human disease. Here we apply it to more than 2,500 simplex families, each having a child with an autistic spectrum disorder. By comparing affected to unaffected siblings, we show that 13% of de novo missense mutations and 43% of de novo likely gene-disrupting (LGD) mutations contribute to 12% and 9% of diagnoses, respectively. Including copy number variants, coding de novo mutations contribute to about 30% of all simplex and 45% of female diagnoses. Almost all LGD mutations occur opposite wild-type alleles. LGD targets in affected females significantly overlap the targets in males of lower intelligence quotient (IQ), but neither overlaps significantly with targets in males of higher IQ. We estimate that LGD mutation in about 400 genes can contribute to the joint class of affected females and males of lower IQ, with an overlapping and similar number of genes vulnerable to contributory missense mutation. LGD targets in the joint class overlap with published targets for intellectual disability and schizophrenia, and are enriched for chromatin modifiers, FMRP-associated genes and embryonically expressed genes. Most of the significance for the latter comes from affected females.

De novo insertions and deletions of predominantly paternal origin are associated with autism spectrum disorder

Whole-exome sequencing (WES) studies have demonstrated the contribution of de novo loss-of-function single-nucleotide variants (SNVs) to autism spectrum disorder (ASD). However, challenges in the reliable detection of de novo insertions and deletions (indels) have limited inclusion of these variants in prior analyses. By applying a robust indel detection method to WES data from 787 ASD families (2,963 individuals), we demonstrate that de novo frameshift indels contribute to ASD risk (OR = 1.6; 95% CI = 1.0-2.7; p = 0.03), are more common in female probands (p = 0.02), are enriched among genes encoding FMRP targets (p = 6 × 10(-9)), and arise predominantly on the paternal chromosome (p < 0.001). On the basis of mutation rates in probands versus unaffected siblings, we conclude that de novo frameshift indels contribute to risk in approximately 3% of individuals with ASD. Finally, by observing clustering of mutations in unrelated probands, we uncover two ASD-associated genes: KMT2E (MLL5), a chromatin regulator, and RIMS1, a regulator of synaptic vesicle release.

Modest impact on risk for autism spectrum disorder of rare copy number variants at 15q11.2, specifically breakpoints 1 to 2

The proximal region of chromosome 15 is one of the genomic hotspots for copy number variants (CNVs). Among the rearrangements observed in this region, CNVs from the interval between the common breakpoints 1 and 2 (BP1 and BP2) have been reported cosegregating with autism spectrum disorder (ASD). Although evidence supporting an association between BP1-BP2 CNVs and autism accumulates, the magnitude of the effect of BP1-BP2 CNVs remains elusive, posing a great challenge to recurrence-risk counseling. To gain further insight into their pathogenicity for ASD, we estimated the penetrance of the BP1-BP2 CNVs for ASD as well as their effects on ASD-related phenotypes in a well-characterized ASD sample (n = 2525 families). Transmission disequilibrium test revealed significant preferential transmission only for the duplicated chromosome in probands (20T:9NT). The penetrance of the BP1-BP2 CNVs for ASD was low, conferring additional risks of 0.3% (deletion) and 0.8% (duplication). Stepwise regression analyses suggest a greater effect of the CNVs on ASD-related phenotype in males and when maternally inherited. Taken together, the results are consistent with BP1-BP2 CNVs as risk factors for autism. However, their effect is modest, more akin to that seen for common variants. To be consistent with the current American College of Medical Genetics guidelines for interpretation of postnatal CNV, the BP1-BP2 deletion and duplication CNVs would probably best be classified as variants of uncertain significance (VOUS): they appear to have an impact on risk, but one so modest that these CNVs do not merit pathogenic status.

Coexpression networks implicate human midfetal deep cortical projection neurons in the pathogenesis of autism

Autism spectrum disorder (ASD) is a complex developmental syndrome of unknown etiology. Recent studies employing exome- and genome-wide sequencing have identified nine high-confidence ASD (hcASD) genes. Working from the hypothesis that ASD-associated mutations in these biologically pleiotropic genes will disrupt intersecting developmental processes to contribute to a common phenotype, we have attempted to identify time periods, brain regions, and cell types in which these genes converge. We have constructed coexpression networks based on the hcASD "seed" genes, leveraging a rich expression data set encompassing multiple human brain regions across human development and into adulthood. By assessing enrichment of an independent set of probable ASD (pASD) genes, derived from the same sequencing studies, we demonstrate a key point of convergence in midfetal layer 5/6 cortical projection neurons. This approach informs when, where, and in what cell types mutations in these specific genes may be productively studied to clarify ASD pathophysiology.

Adjusting head circumference for covariates in autism: clinical correlates of a highly heritable continuous trait

BACKGROUND

Brain development follows a different trajectory in children with autism spectrum disorders (ASD) than in typically developing children. A proxy for neurodevelopment could be head circumference (HC), but studies assessing HC and its clinical correlates in ASD have been inconsistent. This study investigates HC and clinical correlates in the Simons Simplex Collection cohort.

METHODS

We used a mixed linear model to estimate effects of covariates and the deviation from the expected HC given parental HC (genetic deviation). After excluding individuals with incomplete data, 7225 individuals in 1891 families remained for analysis. We examined the relationship between HC/genetic deviation of HC and clinical parameters.

RESULTS

Gender, age, height, weight, genetic ancestry, and ASD status were significant predictors of HC (estimate of the ASD effect = .2 cm). HC was approximately normally distributed in probands and unaffected relatives, with only a few outliers. Genetic deviation of HC was also normally distributed, consistent with a random sampling of parental genes. Whereas larger HC than expected was associated with ASD symptom severity and regression, IQ decreased with the absolute value of the genetic deviation of HC.

CONCLUSIONS

Measured against expected values derived from covariates of ASD subjects, statistical outliers for HC were uncommon. HC is a strongly heritable trait, and population norms for HC would be far more accurate if covariates including genetic ancestry, height, and age were taken into account. The association of diminishing IQ with absolute deviation from predicted HC values suggests HC could reflect subtle underlying brain development and warrants further investigation.

Multilineage somatic activating mutations in HRAS and NRAS cause mosaic cutaneous and skeletal lesions, elevated FGF23 and hypophosphatemia

Pathologically elevated serum levels of fibroblast growth factor-23 (FGF23), a bone-derived hormone that regulates phosphorus homeostasis, result in renal phosphate wasting and lead to rickets or osteomalacia. Rarely, elevated serum FGF23 levels are found in association with mosaic cutaneous disorders that affect large proportions of the skin and appear in patterns corresponding to the migration of ectodermal progenitors. The cause and source of elevated serum FGF23 is unknown. In those conditions, such as epidermal and large congenital melanocytic nevi, skin lesions are variably associated with other abnormalities in the eye, brain and vasculature. The wide distribution of involved tissues and the appearance of multiple segmental skin and bone lesions suggest that these conditions result from early embryonic somatic mutations. We report five such cases with elevated serum FGF23 and bone lesions, four with large epidermal nevi and one with a giant congenital melanocytic nevus. Exome sequencing of blood and affected skin tissue identified somatic activating mutations of HRAS or NRAS in each case without recurrent secondary mutation, and we further found that the same mutation is present in dysplastic bone. Our finding of somatic activating RAS mutation in bone, the endogenous source of FGF23, provides the first evidence that elevated serum FGF23 levels, hypophosphatemia and osteomalacia are associated with pathologic Ras activation and may provide insight in the heretofore limited understanding of the regulation of FGF23.

De novo mutations in histone-modifying genes in congenital heart disease

Congenital heart disease (CHD) is the most frequent birth defect, affecting 0.8% of live births. Many cases occur sporadically and impair reproductive fitness, suggesting a role for de novo mutations. Here we compare the incidence of de novo mutations in 362 severe CHD cases and 264 controls by analysing exome sequencing of parent-offspring trios. CHD cases show a significant excess of protein-altering de novo mutations in genes expressed in the developing heart, with an odds ratio of 7.5 for damaging (premature termination, frameshift, splice site) mutations. Similar odds ratios are seen across the main classes of severe CHD. We find a marked excess of de novo mutations in genes involved in the production, removal or reading of histone 3 lysine 4 (H3K4) methylation, or ubiquitination of H2BK120, which is required for H3K4 methylation. There are also two de novo mutations in SMAD2, which regulates H3K27 methylation in the embryonic left-right organizer. The combination of both activating (H3K4 methylation) and inactivating (H3K27 methylation) chromatin marks characterizes 'poised' promoters and enhancers, which regulate expression of key developmental genes. These findings implicate de novo point mutations in several hundreds of genes that collectively contribute to approximately 10% of severe CHD.

Recessive mutations in DGKE cause atypical hemolytic-uremic syndrome

Pathologic thrombosis is a major cause of mortality. Hemolytic-uremic syndrome (HUS) features episodes of small-vessel thrombosis resulting in microangiopathic hemolytic anemia, thrombocytopenia and renal failure. Atypical HUS (aHUS) can result from genetic or autoimmune factors that lead to pathologic complement cascade activation. Using exome sequencing, we identified recessive mutations in DGKE (encoding diacylglycerol kinase ɛ) that co-segregated with aHUS in nine unrelated kindreds, defining a distinctive Mendelian disease. Affected individuals present with aHUS before age 1 year, have persistent hypertension, hematuria and proteinuria (sometimes in the nephrotic range), and develop chronic kidney disease with age. DGKE is found in endothelium, platelets and podocytes. Arachidonic acid-containing diacylglycerols (DAG) activate protein kinase C (PKC), which promotes thrombosis, and DGKE normally inactivates DAG signaling. We infer that loss of DGKE function results in a prothrombotic state. These findings identify a new mechanism of pathologic thrombosis and kidney failure and have immediate implications for treating individuals with aHUS.

Common genetic variants, acting additively, are a major source of risk for autism

Background

Autism spectrum disorders (ASD) are early onset neurodevelopmental syndromes typified by impairments in reciprocal social interaction and communication, accompanied by restricted and repetitive behaviors. While rare and especially de novo genetic variation are known to affect liability, whether common genetic polymorphism plays a substantial role is an open question and the relative contribution of genes and environment is contentious. It is probable that the relative contributions of rare and common variation, as well as environment, differs between ASD families having only a single affected individual (simplex) versus multiplex families who have two or more affected individuals.

Methods

By using quantitative genetics techniques and the contrast of ASD subjects to controls, we estimate what portion of liability can be explained by additive genetic effects, known as narrow-sense heritability. We evaluate relatives of ASD subjects using the same methods to evaluate the assumptions of the additive model and partition families by simplex/multiplex status to determine how heritability changes with status.

Results

By analyzing common variation throughout the genome, we show that common genetic polymorphism exerts substantial additive genetic effects on ASD liability and that simplex/multiplex family status has an impact on the identified composition of that risk. As a fraction of the total variation in liability, the estimated narrow-sense heritability exceeds 60% for ASD individuals from multiplex families and is approximately 40% for simplex families. By analyzing parents, unaffected siblings and alleles not transmitted from parents to their affected children, we conclude that the data for simplex ASD families follow the expectation for additive models closely. The data from multiplex families deviate somewhat from an additive model, possibly due to parental assortative mating.

Conclusions

Our results, when viewed in the context of results from genome-wide association studies, demonstrate that a myriad of common variants of very small effect impacts ASD liability.

De novo mutations revealed by whole-exome sequencing are strongly associated with autism

Multiple studies have confirmed the contribution of rare de novo copy number variations to the risk for autism spectrum disorders. But whereas de novo single nucleotide variants have been identified in affected individuals, their contribution to risk has yet to be clarified. Specifically, the frequency and distribution of these mutations have not been well characterized in matched unaffected controls, and such data are vital to the interpretation of de novo coding mutations observed in probands. Here we show, using whole-exome sequencing of 928 individuals, including 200 phenotypically discordant sibling pairs, that highly disruptive (nonsense and splice-site) de novo mutations in brain-expressed genes are associated with autism spectrum disorders and carry large effects. On the basis of mutation rates in unaffected individuals, we demonstrate that multiple independent de novo single nucleotide variants in the same gene among unrelated probands reliably identifies risk alleles, providing a clear path forward for gene discovery. Among a total of 279 identified de novo coding mutations, there is a single instance in probands, and none in siblings, in which two independent nonsense variants disrupt the same gene, SCN2A (sodium channel, voltage-gated, type II, α subunit), a result that is highly unlikely by chance.

Rare copy number variants in tourette syndrome disrupt genes in histaminergic pathways and overlap with autism

BACKGROUND

Studies of copy number variation (CNV) have characterized loci and molecular pathways in a range of neuropsychiatric conditions. We analyzed rare CNVs in Tourette syndrome (TS) to identify novel risk regions and relevant pathways, to evaluate burden of structural variation in cases versus controls, and to assess overlap of identified variations with those in other neuropsychiatric syndromes.

METHODS

We conducted a case-control study of 460 individuals with TS, including 148 parent-child trios and 1131 controls. CNV analysis was undertaken using 370 K to 1 M probe arrays, and genotyping data were used to match cases and controls for ancestry. CNVs present in < 1% of the population were evaluated.

RESULTS

While there was no significant increase in the number of de novo or transmitted rare CNVs in cases versus controls, pathway analysis using multiple algorithms showed enrichment of genes within histamine receptor (subtypes 1 and 2) signaling pathways (p = 5.8 × 10(-4) - 1.6 × 10(-2)), as well as axon guidance, cell adhesion, nervous system development, and synaptic structure and function processes. Genes mapping within rare CNVs in TS showed significant overlap with those previously identified in autism spectrum disorders but not intellectual disability or schizophrenia. Three large, likely pathogenic, de novo events were identified, including one disrupting multiple gamma-aminobutyric acid receptor genes.

CONCLUSIONS

We identify further evidence supporting recent findings regarding the involvement of histaminergic and gamma-aminobutyric acidergic mechanisms in the etiology of TS and show an overlap of rare CNVs in TS and autism spectrum disorders.

Parallel declines in cognition, motivation, and locomotion in aging mice: association with immune gene upregulation in the medial prefrontal cortex

Aging in humans is associated with parallel changes in cognition, motivation, and motoric performance. Based on the human aging literature, we hypothesized that this constellation of age-related changes is mediated by the medial prefrontal cortex and that it would be observed in aging mice. Toward this end, we performed detailed assessments of cognition, motivation, and motoric behavior in aging mice. We assessed behavioral and cognitive performance in C57Bl/6 mice aged 6, 18, and 24 months, and followed this with microarray analysis of tissue from the medial prefrontal cortex and analysis of serum cytokine levels. Multivariate modeling of these data suggested that the age-related changes in cognition, motivation, motor performance, and prefrontal immune gene expression were highly correlated. Peripheral cytokine levels were also correlated with these variables, but less strongly than measures of prefrontal immune gene upregulation. To determine whether the observed immune gene expression changes were due to prefrontal microglial cells, we isolated CD11b-positive cells from the prefrontal cortex and subject them to next-generation RNA sequencing. Many of the immune changes present in whole medial prefrontal cortex were enriched in this cell population. These data suggest that, as in humans, cognition, motivation, and motoric performance in the mouse change together with age and are strongly associated with CNS immune gene upregulation.

Functional genomic and proteomic analysis reveals disruption of myelin-related genes and translation in a mouse model of early life neglect

Early life neglect is an important public health problem which can lead to lasting psychological dysfunction. Good animal models are necessary to understand the mechanisms responsible for the behavioral and anatomical pathology that results. We recently described a novel model of early life neglect, maternal separation with early weaning (MSEW), that produces behavioral changes in the mouse that persist into adulthood. To begin to understand the mechanism by which MSEW leads to these changes we applied cDNA microarray, next-generation RNA-sequencing (RNA-seq), label-free proteomics, multiple reaction monitoring (MRM) proteomics, and methylation analysis to tissue samples obtained from medial prefrontal cortex to determine the molecular changes induced by MSEW that persist into adulthood. The results show that MSEW leads to dysregulation of markers of mature oligodendrocytes and genes involved in protein translation and other categories, an apparent downward biasing of translation, and methylation changes in the promoter regions of selected dysregulated genes. These findings are likely to prove useful in understanding the mechanism by which early life neglect affects brain structure, cognition, and behavior.

A novel heterozygous deletion within the 3' region of the PAX6 gene causing isolated aniridia in a large family group

Paired box gene 6 (PAX6) is the causative gene of aniridia. It is a dominantly inherited eye abnormality characterized by partial or complete absence of the iris. The PAX6 gene is located on chromosome 11p13 and contains 14 exons. It is expressed mainly in the developing eye and central nervous system. Submicroscopic copy number variations are common in the human genome. Submicroscopic deletions may cause several human diseases, either by disrupting coding sequences or by eliminating regulatory elements essential for expression of the gene in question. Over the past several years, array-based comparative genomic hybridization has become an increasingly useful tool for both identifying normal cytogenetic variations and characterizing chromosomal abnormalities associated with developmental delays and cancer. Our results support the notion that assessing copy number variation of the PAX6 gene itself and also of flanking regions, may contribute to the molecular diagnosis of aniridia.

Predicting neuroendocrine tumor (carcinoid) neoplasia using gene expression profiling and supervised machine learning

BACKGROUND

A more accurate taxonomy of small intestinal (SI) neuroendocrine tumors (NETs) is necessary to accurately predict tumor behavior and prognosis and to define therapeutic strategy. In this study, the authors identified a panel of such markers that have been implicated in tumorigenicity, metastasis, and hormone production and hypothesized that transcript levels of the genes melanoma antigen family D2 (MAGE-D2), metastasis-associated 1 (MTA1), nucleosome assembly protein 1-like (NAP1L1), Ki-67 (a marker of proliferation), survivin, frizzled homolog 7 (FZD7), the Kiss1 metastasis suppressor (Kiss1), neuropilin 2 (NRP2), and chromogranin A (CgA) could be used to define primary SI NETs and to predict the development of metastases.

METHODS

Seventy-three clinically and World Health Organization pathologically classified NET samples (primary tumor, n = 44 samples; liver metastases, n = 29 samples) and 30 normal human enterochromaffin (EC) cell preparations were analyzed using real-time polymerase chain reaction. Transcript levels were normalized to 3 NET housekeeping genes (asparagine-linked glycosylation 9 or ALG9, transcription factor CP2 or TFCP2, and zinc finger protein 410 or ZNF410) using geNorm analysis. A predictive gene-based model was constructed using supervised learning algorithms from the transcript expression levels.

RESULTS

Primary SI NETs could be differentiated from normal human EC cell preparations with 100% specificity and 92% sensitivity. Well differentiated NETs (WDNETs), well differentiated neuroendocrine carcinomas, and poorly differentiated NETs (PDNETs) were classified with a specificity of 78%, 78%, and 71%, respectively; whereas poorly differentiated neuroendocrine carcinomas were misclassified as either WDNETs or PDNETs. Metastases were predicted in all cases with 100% sensitivity and specificity.

CONCLUSIONS

The current results indicated that gene expression profiling and supervised machine learning can be used to classify SI NET subtypes and accurately predict metastasis. The authors believe that the application of this technique will facilitate accurate molecular pathologic delineation of NET disease, better define its extent, facilitate the assessment of prognosis, and provide a guide for the identification of appropriate strategies for individualized patient treatment.

Differential effects of diabetes induced by streptozotocin and that develops spontaneously on prostate growth in Bio Breeding (BB) rats

We investigated molecular changes in the response to insulin in prostates of spontaneously developed (Bio Breeding) and streptozotocin (STZ)-induced diabetic rats that received sufficient amounts (euglycemic group), or suboptimal doses (hyperglycemic group) of insulin for 32 weeks, using Affymetrix GeneChip analysis of gene expression. Alterations in gene expression levels identified by microarray analysis, having potential biological relevance to prostate growth, were verified by real-time reverse transcription polymerase chain reaction (RT-PCR). A significant decrease in the weight of ventral prostate was observed in the hyperglycemic STZ-induced but not spontaneously developed diabetic group. Microarray analysis revealed that gene expression profiles were distinctly different in each region of the prostate, and that hyperglycemic diabetes in spontaneously developed and STZ-diabetic rats was associated with differential changes in the prostatic expression levels of 856 genes, of which 35 were related to cell growth, proliferation and death. RT-PCR data verified significant differences in the mRNA expression levels of Igfbp6, Tieg, and Clu between euglycemic and hyperglycemic groups, whereas expression levels of these genes in control and euglycemic diabetic groups were not significantly different. In ventral prostate, the mRNA expression levels of Igfbp6 and Tieg were significantly higher in the hyperglycemic STZ-induced diabetic than in the hyperglycemic spontaneously diabetic BBDP/Wor rats. Our data demonstrate that the diabetes induced by STZ in the BBDR/Wor rats affects prostate growth and the molecular response to insulin differently than that observed in BBDP/Wor rats that develop diabetes spontaneously.

RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays

Ultra-high-throughput sequencing is emerging as an attractive alternative to microarrays for genotyping, analysis of methylation patterns, and identification of transcription factor binding sites. Here, we describe an application of the Illumina sequencing (formerly Solexa sequencing) platform to study mRNA expression levels. Our goals were to estimate technical variance associated with Illumina sequencing in this context and to compare its ability to identify differentially expressed genes with existing array technologies. To do so, we estimated gene expression differences between liver and kidney RNA samples using multiple sequencing replicates, and compared the sequencing data to results obtained from Affymetrix arrays using the same RNA samples. We find that the Illumina sequencing data are highly replicable, with relatively little technical variation, and thus, for many purposes, it may suffice to sequence each mRNA sample only once (i.e., using one lane). The information in a single lane of Illumina sequencing data appears comparable to that in a single array in enabling identification of differentially expressed genes, while allowing for additional analyses such as detection of low-expressed genes, alternative splice variants, and novel transcripts. Based on our observations, we propose an empirical protocol and a statistical framework for the analysis of gene expression using ultra-high-throughput sequencing technology.

Rapid identification of disease-causing mutations using copy number analysis within linkage intervals

SNP and comparative genome hybridization arrays (aCGH) are powerful techniques for identifying genome rearrangements, deletions, and duplications. We hypothesized that current array-based detection of copy number variation (CNV) could complement parametric linkage analysis and allow the rapid identification of functional mutations in families with inherited disorders. Herein, we demonstrate the utility of this technique by rapidly identifying a disease causing microdeletion within the PARK2 gene in a family with autosomal recessive Parkinsonism.

Small bowel carcinoid (enterochromaffin cell) neoplasia exhibits transforming growth factor-beta1-mediated regulatory abnormalities including up-regulation of C-Myc and MTA1

BACKGROUND

Although it is known that small intestinal carcinoids are derived from enterochromaffin (EC) cells, these cells remain poorly characterized and little is known of the growth regulatory mechanisms of these neuroendocrine cells. Down-regulation or loss of the transforming growth factor-beta1 (TGFbeta1) cytostatic program and activation of TGFbeta-mediated transcriptional networks is associated with uncontrolled growth and metastasis in other neural tumors, glioblastomas. Whether this phenomenon is common to small intestinal carcinoid tumors was investigated.

METHODS

The effects of TGFbeta1 on cultured normal EC cells (isolated by FACS sorting) and the neoplastic EC cell line, KRJ-I, was assessed using the MTT assay. The TGFbetaRII transcript and protein were identified in tumor cells and the effects of TGFbeta1 on SMAD2 phosphorylation and nuclear translocation quantified. The time-dependent response of SMAD4, SMAD7, c-Myc, and P21(WAF1/CIP1) protein expression and c-Myc and p21(WAF1/CIP1) transcript was measured in response to TGFbeta1 and the transcript expression of candidate downstream targets, MTA1 and E-cadherin, were assessed.

RESULTS

TGFbeta1 inhibited normal EC cell proliferation (IC(50) = 17 pM) but stimulated neoplastic EC cell proliferation (EC(50) = 22 pM). In tumor cells, significantly decreased transcript (P < .01) of TGFbetaRII was identified, but no receptor mutations were identified and protein expression was evident. TGFbeta1 (1 ng/mL) resulted in SMAD2 phosphorylation and <7% nuclear expression compared with 93% in normal EC cells. In neoplastic cells, TGFbeta1 (1 ng/mL) caused a decrease in SMAD4 (>16%, P < .05), whereas SMAD7 and c-Myc transcript and protein were respectively increased >21% (P < .05) and approximately 40% (P < .002). TGFbeta1 (1 ng/mL) also decreased p21(WAF1/CIP1) transcript by 60% (P < .001) and protein that was undetectable at 24 hours. Expression of the downstream targets of the c-Myc pathway, MTA1, was increased (20%) and E-cadherin decreased (30%).

CONCLUSIONS

The neoplastic EC cell is characterized by loss of TGFbeta-1-mediated growth inhibition and, similar to glioblastomas, utilizes the TGFbeta system to induce gene responses associated with growth promotion (c-Myc and the ERK pathway), invasion (E-cadherin), and metastasis (MTA1).

Role of CCN2/CTGF in the proliferation of Mastomys enterochromaffin-like cells and gastric carcinoid development

Mastomys enterochromaffin-like (ECL) cell proliferation is initially gastrin driven, but once neoplasia develops, cells become gastrin autonomous. We hypothesized that CCN2 (CTGF), a mitogenic growth factor, may regulate ECL cell proliferation. A Mastomys GeneChip database was examined (dCHIP) to identify CCN2 expression levels. CCN2 in normal and tumor ECL cell preparations obtained using FACS (100 nM acridine orange) was examined by real-time PCR. CCN2 protein was identified in mucosal and ECL cell preparations by immunohistochemistry. Short-term cultured cells were stimulated with either CCN2 or CCN2 + EGF, and proliferation was measured (MTT assay). The ERK1/2 inhibitor PD-98059 (0.1-100 microM) was assessed in terms of CCN2 (1 ng/ml)-mediated proliferation and ERK1/2 phosphorylation. CCN2 transcript and protein was then examined in clinical gastric carcinoids. The ccn2 transcript was upregulated in tumor samples compared with the normal mucosa (+2.36-fold, P < 0.01). PCR demonstrated that ccn2 was not expressed in FACS-prepared (>98% pure) normal ECL cells but was elevated in tumor ECL cell fractions (41.3 +/- 10.7-fold). Immunostaining of the Mastomys gastric mucosa and FACS preparations confirmed that CCN2 protein was present in ECL tumors but not in normal ECL cells. Neither CCN2 nor CCN2 + EGF stimulated normal ECL cell proliferation. CCN2 stimulated tumor proliferation (EC50 approximately 0.01 ng/ml); EGF significantly augmented (P < 0.01) CCN2-induced tumor cell proliferation (EC50 = 20 pg/ml). PD-98059 inhibited CCN2-induced proliferation (-12 +/- 3%, P < 0.05) and ERK1/2 phosphorylation (-34 +/- 5%, P < 0.05) in tumor cells. In clinical samples, both CCN2 transcript and protein were elevated in gastrin-autonomous carcinoids (P < 0.02) compared with the normal mucosa. In conclusion, CCN2 may be a proliferative regulator of Mastomys ECL neoplastic proliferation once these cells become autonomous of gastrin regulation. Identification of CCN2 in gastric carcinoid tissue may be useful both as an indicator of ECL cell transformation and may define gastrin autonomy, a criteria of gastric carcinoid malignancy.

Genetic differentiation of appendiceal tumor malignancy: a guide for the perplexed

OBJECTIVE

To use differential gene expression of candidate markers to discriminate benign appendiceal carcinoids (APCs) from malignant and mixed cell APCs.

SUMMARY BACKGROUND DATA

Controversy exists in regard to the appropriate surgical management of APCs since it is sometimes difficult to predict tumor behavior using traditional pathologic criteria. We have identified 5 differentially expressed genes (a mitosis-regulatory gene NAP1L1, an adhesin MAGE-D2, an estrogen-antagonist, the metastasis marker MTA1, the apoptotic marker NALP, and chromogranin A) that define gut neuroendocrine cell behavior.

METHODS

Total RNA was isolated using TRIzol reagent from 42 appendiceal samples, including appendiceal carcinoids identified at exploration for appendicitis (no evidence of metastasis; n = 16), appendicitis specimens (n = 11), malignant appendiceal tumors (> 1.5 cm, evidence of metastatic invasion; n = 7), and mixed (goblet) cell appendiceal adenocarcinoids (n = 3), normal appendiceal tissue (n = 5), and 5 colorectal cancers. Gene expression (CgA, NAP1L1, MAGE-D2, MTA1, and NALP1) was examined by Q-RT PCR (Applied Biosystems) and quantified against GAPDH.

RESULTS

CgA message was elevated (> 1000-fold, P < 0.05) in all tumor types. NAP1L1 was elevated (> 10-fold, P < 0.03) in both malignant and goblet cell adenocarcinoids compared with normal and incidental lesions (P < 0.006). MAGE-D2 and MTA1 message were significantly elevated (> 10-fold, P < 0.01) in the malignant and goblet cell adenocarcinoid tumors but not in the appendicitis-associated carcinoids or normal mucosa. The apoptotic marker, NALP1, was overexpressed (> 50-fold, P < 0.05) in the appendicitis-associated and malignant appendiceal carcinoids but was significantly decreased (> 10-fold, P < 0.05) in the goblet cell adenocarcinoids. Elevated CgA transcript and protein levels indicative of a carcinoid tumor were identified in one acute appendicitis sample with no histologic evidence of a tumor.

CONCLUSIONS

These data demonstrate that malignant APCs and goblet cell adenocarcinoids have elevated expression of NAP1L1, MAGE-D2, and MTA1 compared with appendiceal carcinoids identified at surgery for appendicitis. This and the differences in NALP1 gene expression (decreased in goblet cell adenocarcinoids) provide a series of molecular signatures that differentiate carcinoids of the appendix. CgA identified all appendiceal tumors as well as covert lesions, which may be more prevalent than previously recognized. The molecular delineation of malignant appendiceal tumor potential provides a scientific basis to define the appropriate surgical management as opposed to morphologic assessment alone.

Utility of molecular genetic signatures in the delineation of gastric neoplasia

BACKGROUND

Current techniques to define gastric neoplasia are limited but molecular genetic signatures can categorize tumors and provide biological rationale for predicting clinical behavior. We identified three gene signatures: Chromogranin A (CgA), MAGE-D2 (adhesion), and MTA1 (metastasis) that define gastrointestinal (GI) carcinoids and hypothesize that their expression can delineate gastric neoplasia. This strategy provides a molecular basis to define neuroendocrine gastric carcinoids (GCs), neuronal stromal tumors (GISTs), or epithelial cell (gastric adenocarcinomas [GCAs])-derived tumors.

METHODS

Total RNA was isolated from 38 GCs: Type I/II (n = 7), Type III/IV (n = 6), GISTs (n = 12), GCAs (n = 13), and normal mucosa (n = 12). Quantitative reverse transcriptase polymerase chain reaction (Q RT-PCR) gene expression was quantified against glyseraldehyde-3-phosphate dehydrogenase (GAPDH) and CgA and MTA1 protein expression levels were analyzed by immunohistochemical analyses of a gastric neoplasia microarray.

RESULTS

CgA was elevated in Type I/II (10-fold; P < .01) and Type III/IV (100-fold, P < .005), decreased in GISTs (100-fold, P < .03), and unchanged in GCAs. MAGE-D2 was 5-10-fold elevated (P < .05) in Type III/IV, GISTs, and GCAs but not in Type I/II tumors. MTA1 (> 5-fold, P < .01) was elevated in GCs (Type III/IV>I/II, P < .05), in GISTs (> 4-fold, P < .05), and GCAs. CgA protein levels were elevated in GCs (P < .005) but not in GISTs and GCAs. MTA1 levels were elevated in all tumors (P < .02) compared with normal, and especially with tumor invasion (P < .05).

CONCLUSION

CgA discriminates GCs from other gastric neoplasms; overexpression of MAGE-D2 and MTA1 differentiate Type III/IV from Type I/II GCs. GISTs share similar expression patterns with Type III/IV GCs but have decreased CgA. MTA1 is a marker of tumor invasion.

Q RT-PCR detection of chromogranin A: a new standard in the identification of neuroendocrine tumor disease

OBJECTIVE

Message and protein expression of CgA was examined to evaluate the sensitivity of a PCR-based approach in the detection of covert neuroendocrine (NE) tissue.

SUMMARY BACKGROUND DATA

Immunohistochemical (IHC) measurement of chromogranin A (CgA) discriminates gastrointestinal (GI) carcinoids from epithelial tumors. IHC is, however, an insensitive technique to identify micrometastases or delineate subpopulations of NE cells.

METHODS

CgA gene expression was examined by Q-RT PCR in GI carcinoids (small intestinal and metastases, n=17, gastric, n=5), appendiceal tumors (n=10), and adenocarcinomas (gastric, n=5, colorectal, n=6). CgA protein expression levels were quantitatively analyzed following IHC by automated quantitative analysis (AQUA) in 2 tissue microarrays (GI carcinoid and GI adenocarcinoma).

RESULTS

CgA gene was overexpressed (P<0.001) in GI carcinoids compared with GI adenocarcinomas and normal mucosa. Elevated levels (P<0.00001) were also identified in carcinoid liver and lymph node (LN) metastases. CgA levels were higher (approximately 2-4-fold) in NE appendiceal carcinoids than in adenocarcinoids, but in GI adenocarcinomas were identical to normal mucosa. Histologically normal lymph nodes expressed detectable CgA message in 30% of cases. CgA protein levels were highest in primary GI carcinoids and in liver metastases and significantly elevated (P<0.005) compared with nonmetastatic lesions. Expression in liver and LN metastases was significantly elevated (P<0.000001) compared with normal. Analysis of mRNA by Q-RT PCR was >200-fold more sensitive than by IHC.

CONCLUSIONS

Overexpression of CgA mRNA and protein in GI carcinoids can identify metastatic cells; thus, PCR for CgA can be used to identify micrometastases not evident by light microscopy or IHC as well as define tumors of ambivalent morphologic phenotype. The use of this sensitive strategy to assess NETs and apparently normal LNs and liver may be of future utility in defining therapeutic strategy.

The role of genetic markers--NAP1L1, MAGE-D2, and MTA1--in defining small-intestinal carcinoid neoplasia

BACKGROUND

Standard clinical and immunohistochemical methods cannot reliably determine whether a small intestinal carcinoid (SIC) is indolent or aggressive. We hypothesized that carcinoid malignancy could be defined by using quantitative reverse transcriptase-polymerase chain reaction (QRT-PCR) and immunohistochemical approaches that evaluate potential marker genes.

METHODS

Candidate marker gene expression (nucleosome assembly protein 1-like 1 [NAP1L1], melanoma antigen D2 [MAGE-D2], and metastasis-associated protein 1 [MTA1]) identified by Affymetrix transcriptional profiling was examined by QRT-PCR in SIC, liver, and lymph node (LN) metastases, colorectal carcinomas, and healthy tissues. Immunohistochemical expression levels of MTA1 were analyzed quantitatively by a novel automated quantitative analysis in a tissue microarray of 102 gastrointestinal carcinoids and in a breast/prostate carcinoma array.

RESULTS

Affymetrix transcriptional profiling identified three potentially useful malignancy-marker genes (out of 1709 significantly altered genes). By QRT-PCR, NAP1L1 was significantly (P < .03) overexpressed in SIC compared with colorectal carcinomas and healthy tissue. Increased levels (P < .05) were identified in both liver and LN metastases. Levels in colorectal carcinomas were the same as in healthy mucosa. MAGE-D2 and MTA1 were increased (P < .05) in primary tumors and metastases and overexpressed in carcinomas. Automated quantitative analysis demonstrated the highest levels of MTA1 immunostaining in malignant primary SICs and in metastases to the liver and LN. These were significantly increased (P < .02) compared with nonmetastatic primary tumors. MTA1 was overexpressed in breast and prostate carcinomas (P < .05).

CONCLUSIONS

SICs overexpress the neoplasia-related genes NAP1L1 (mitotic regulation), MAGE-D2 (adhesion), and MTA1 (estrogen antagonism). The ability to determine the malignant potential of these tumors and their propensity to metastasize provides a biological rationale for the management of carcinoids and may have prognostic utility.

Complement factor H polymorphism in age-related macular degeneration

Age-related macular degeneration (AMD) is a major cause of blindness in the elderly. We report a genome-wide screen of 96 cases and 50 controls for polymorphisms associated with AMD. Among 116,204 single-nucleotide polymorphisms genotyped, an intronic and common variant in the complement factor H gene (CFH) is strongly associated with AMD (nominal P value <10(-7)). In individuals homozygous for the risk allele, the likelihood of AMD is increased by a factor of 7.4 (95% confidence interval 2.9 to 19). Resequencing revealed a polymorphism in linkage disequilibrium with the risk allele representing a tyrosine-histidine change at amino acid 402. This polymorphism is in a region of CFH that binds heparin and C-reactive protein. The CFH gene is located on chromosome 1 in a region repeatedly linked to AMD in family-based studies.

Microsatellite instability and gene mutations in transforming growth factor-beta type II receptor are absent in small bowel carcinoid tumors

BACKGROUND

Microsatellite instability (MSI) with concomitant mutations in the coding region of transforming growth factor-beta type II receptor (TGFbetaRII) results in an aberrant growth-regulatory phenotype in colorectal carcinomas. The authors postulated that a similar mechanism occurred during the malignant evolution of small bowel carcinoid tumors.

METHODS

Mutational analysis of two coding regions in the TGFbetaRII gene associated with MSI and BAT-26 within intron 5 of the mismatch repair gene, hMSH2, was undertaken in small bowel carcinoids (n = 14), lymph node metasasis (n = 1) and liver metastases (n = 5). Quantitative PCR analysis [TAQMAN, Applied Biosystems, Foster City, CA] was then undertaken to examine gene alterations in mismatch repair genes (hMLH1 and hMSH2) in small bowel carcinoids (n = 7) and matched normal mucosa (n = 5). Staining was then analyzed using quantitative tissue array profiling (AQUA analysis) in a small bowel EC carcinoid tissue microarray (n = 55 tumors) with immunostaining against TGFbetaRII and MSH2.

RESULTS

Mutational examination of the TFGbetaRII gene and BAT-26 demonstrated that MSI was not present in any carcinoid material. Q RT-PCR analysis demonstrated statistically significant increased message levels of hMSH2 but not hMLH1 in carcinoid tumors. Quantitative analysis of membrane TGFbetaRII immunostaining using AQUA demonstrated that TGFbetaRII expression was down-regulated (P < 0.0002) in thirty-three primary small bowel carcinoids that exhibited lymph node and liver metastases compared to normal mucosa. AQUA analysis of nuclear MSH2 immunostaining demonstrated no differences for MSH2 between normal tissue and carcinoid tumor metastasis. Small bowel carcinoids characterized by variable expression of TGFbetaRII, did not exhibit MSI and had no differences in MSH2 expression.

CONCLUSIONS

The molecular events leading to the formation of carcinoid tumors in the small bowel were different from those resulting in epithelial carcinomas. The usually slow-growing and relatively nonaggressive carcinoid tumors had variable expression of TGFbetaRII but were associated with the retention of mismatch repair protein function and a microsatellite-stable phenotype.

Regulation of the Erk2-Elk1 signaling pathway and megakaryocytic differentiation of Bcr-Abl(+) K562 leukemic cells by Gab2

In the blast crisis phase of chronic myelogenous leukemia (CML), Bcr-Abl(+) myeloblasts fail to undergo terminal maturation. The extracellular signal-regulated kinase (Erk) mitogen-activated protein (MAP) kinase has been shown to mediate terminal differentiation of myeloid cells. Interestingly, Bcr-Abl(+) CML cell lines established from blast crisis were found to have low Erk MAP kinase activity. In this study, we analyzed the role of the Gab2 docking protein in regulation of the Erk MAP kinase in Bcr-Abl(+) K562 human CML cells. Overexpression of Gab2 in K562 cells resulted in transcriptional activation of the c-fos serum response element (SRE) promoter, whereas overexpression of SHP2, Grb2, and CrkL had no effect. Activation of the c-fos SRE transcriptional activity by Gab2 required tyrosine 604, which is a SHP2 docking site on Gab2, and the SHP2 tyrosine phosphatase activity. Elk1, c-Jun, and CHOP trans-reporting assays indicated that overexpression of Gab2 selectively activated the Erk2-Elk1 signaling pathway. To determine cellular consequences of elevating the Gab2 level in K562 cells, stable cell lines for doxycycline-inducible expression of the wild-type Gab2 (Gab2WT) and an SHP2-binding defective Gab2 (Gab2Tyr604Phe) were established. Analysis of these cell lines indicated that induction of Gab2WT expression, but not Gab2Tyr604Phe expression, led to Erk activation, growth arrest, cell spreading, and enlargement; expression of megakaryocyte/platelet lineage-specific integrins alphaIIb/beta3 (CD41/CD61); and upregulation of RNA for megakaryocyte/platelet proteins. All of these changes are characteristics of megakaryocytic differentiation. Together, these results reveal Gab2 as a limiting signaling component for Erk MAP kinase activation and terminal differentiation of K562 CML cells.