Gene characteristics predicting missense, nonsense and frameshift mutations in tumor samples.

BACKGROUND: Because driver mutations provide selective advantage to the mutant clone, they tend to occur at a higher frequency in tumor samples compared to selectively neutral (passenger) mutations. However, mutation frequency alone is insufficient to identify cancer genes because mutability is influenced by many gene characteristics, such as size, nucleotide composition, etc. The goal of this study was to identify gene characteristics associated with the frequency of somatic mutations in the gene in tumor samples. RESULTS: We used data on somatic mutations detected by genome wide screens from the Catalog of Somatic Mutations in Cancer (COSMIC). Gene size, nucleotide composition, expression level of the gene, relative replication time in the cell cycle, level of evolutionary conservation and other gene characteristics (totaling 11) were used as predictors of the number of somatic mutations. We applied stepwise multiple linear regression to predict the number of mutations per gene. Because missense, nonsense, and frameshift mutations are associated with different sets of gene characteristics, they were modeled separately. Gene characteristics explain 88% of the variation in the number of missense, 40% of nonsense, and 23% of frameshift mutations. Comparisons of the observed and expected numbers of mutations identified genes with a higher than expected number of mutations- positive outliers. Many of these are known driver genes. A number of novel candidate driver genes was also identified. CONCLUSIONS: By comparing the observed and predicted number of mutations in a gene, we have identified known cancer-associated genes as well as 111 novel cancer associated genes. We also showed that adding the number of silent mutations per gene reported by genome/exome wide screens across all cancer type (COSMIC data) as a predictor substantially exceeds predicting accuracy of the most popular cancer gene predicting tool - MutsigCV.

A recurrent mutation in PARK2 is associated with familial lung cancer.

PARK2, a gene associated with Parkinson disease, is a tumor suppressor in human malignancies. Here, we show that c.823C>T (p.Arg275Trp), a germline mutation in PARK2, is present in a family with eight cases of lung cancer. The resulting amino acid change, p.Arg275Trp, is located in the highly conserved RING finger 1 domain of PARK2, which encodes an E3 ubiquitin ligase. Upon further analysis, the c.823C>T mutation was detected in three additional families affected by lung cancer. The effect size for PARK2 c.823C>T (odds ratio = 5.24) in white individuals was larger than those reported for variants from lung cancer genome-wide association studies. These data implicate this PARK2 germline mutation as a genetic susceptibility factor for lung cancer. Our results provide a rationale for further investigations of this specific mutation and gene for evaluation of the possibility of developing targeted therapies against lung cancer in individuals with PARK2 variants by compensating for the loss-of-function effect caused by the associated variation.

Drosophila pheromone-sensing neurons expressing the ppk25 ion channel subunit stimulate male courtship and female receptivity.

As in many species, gustatory pheromones regulate the mating behavior of Drosophila. Recently, several ppk genes, encoding ion channel subunits of the DEG/ENaC family, have been implicated in this process, leading to the identification of gustatory neurons that detect specific pheromones. In a subset of taste hairs on the legs of Drosophila, there are two ppk23-expressing, pheromone-sensing neurons with complementary response profiles; one neuron detects female pheromones that stimulate male courtship, the other detects male pheromones that inhibit male-male courtship. In contrast to ppk23, ppk25, is only expressed in a single gustatory neuron per taste hair, and males with impaired ppk25 function court females at reduced rates but do not display abnormal courtship of other males. These findings raised the possibility that ppk25 expression defines a subset of pheromone-sensing neurons. Here we show that ppk25 is expressed and functions in neurons that detect female-specific pheromones and mediates their stimulatory effect on male courtship. Furthermore, the role of ppk25 and ppk25-expressing neurons is not restricted to responses to female-specific pheromones. ppk25 is also required in the same subset of neurons for stimulation of male courtship by young males, males of the Tai2 strain, and by synthetic 7-pentacosene (7-P), a hydrocarbon normally found at low levels in both males and females. Finally, we unexpectedly find that, in females, ppk25 and ppk25-expressing cells regulate receptivity to mating. In the absence of the third antennal segment, which has both olfactory and auditory functions, mutations in ppk25 or silencing of ppk25-expressing neurons block female receptivity to males. Together these results indicate that ppk25 identifies a functionally specialized subset of pheromone-sensing neurons. While ppk25 neurons are required for the responses to multiple pheromones, in both males and females these neurons are specifically involved in stimulating courtship and mating.

Two Drosophila DEG/ENaC channel subunits have distinct functions in gustatory neurons that activate male courtship.

Trimeric sodium channels of the DEG/ENaC family have important roles in neurons, but the specific functions of different subunits present in heteromeric channels are poorly understood. We previously reported that the Drosophila DEG/ENaC subunit Ppk25 is essential in a small subset of gustatory neurons for activation of male courtship behavior, likely through detection of female pheromones. Here we show that, like mutations in ppk25, mutations in another Drosophila DEG/ENaC subunit gene, nope, specifically impair male courtship of females. nope regulatory sequences drive reporter gene expression in gustatory neurons of the labellum wings, and legs, including all gustatory neurons in which ppk25 function is required for male courtship of females. In addition, gustatory-specific knockdown of nope impairs male courtship. Further, the impaired courtship response of nope mutant males to females is rescued by targeted expression of nope in the subset of gustatory neurons in which ppk25 functions. However, nope and ppk25 have nonredundant functions, as targeted expression of ppk25 does not compensate for the lack of nope and vice versa. Moreover, Nope and Ppk25 form specific complexes when coexpressed in cultured cells. Together, these data indicate that the Nope and Ppk25 polypeptides have specific, nonredundant functions in a subset of gustatory neurons required for activation of male courtship in response to females, and suggest the hypothesis that Nope and Ppk25 function as subunits of a heteromeric DEG/ENaC channel required for gustatory detection of female pheromones.

A Drosophila DEG/ENaC subunit functions specifically in gustatory neurons required for male courtship behavior.

Detection of specific female pheromones stimulates courtship behavior in Drosophila melanogaster males, but the chemosensory molecules, cells, and mechanisms involved remain poorly understood. Here we show that ppk25, a DEG/ENaC ion channel subunit required for normal male response to females, is expressed at highest levels in a single sexually dimorphic gustatory neuron of most taste hairs on legs and wings, but not in neurons that detect courtship-inhibiting pheromones or food. Synaptic inactivation of ppk25-expressing neurons, or knockdown of ppk25 expression in all gustatory neurons, significantly impairs male response to females, whereas gustatory expression of ppk25 rescues the courtship behavior of ppk25 mutant males. Remarkably, the only other detectable albeit significantly weaker expression of ppk25 occurs in olfactory neurons implicated in modulation of courtship behavior. However, expression of ppk25 in olfactory neurons is not required for male courtship under our experimental conditions. These data show that ppk25 functions specifically in peripheral taste neurons involved in activation of courtship behavior, an unexpected function for this type of channel. Furthermore, our work identifies a small subset of gustatory neurons with an essential role in activation of male courtship behavior, most likely in response to female pheromones.

Genetic Variation and Recurrent Haplotypes on Chromosome 6q23-25 Risk Locus in Familial Lung Cancer.

Although lung cancer is known to be caused by environmental factors, it has also been shown to have genetic components, and the genetic etiology of lung cancer remains understudied. We previously identified a lung cancer risk locus on 6q23-25 using microsatellite data in families with a history of lung cancer. To further elucidate that signal, we performed targeted sequencing on nine of our most strongly linked families. Two-point linkage analysis of the sequencing data revealed that the signal was heterogeneous and that different families likely had different risk variants. Three specific haplotypes were shared by some of the families: 6q25.3-26 in families 42 and 44, 6q25.2-25.3 in families 47 and 59, and 6q24.2-25.1 in families 30, 33, and 35. Region-based logarithm of the odds scores and expression data identified the likely candidate genes for each haplotype overlap: ARID1B at 6q25.3, MAP3K4 at 6q26, and UTRN (6q24.1) and PHACTR2 (6q24.2). Further annotation was used to zero in on potential risk variants in those genes. All four genes are good candidate genes for lung cancer risk, having been linked to either lung cancer specifically or other cancers. However, this is the first time any of these genes has been implicated in germline risk. Functional analysis of these four genes is planned for future work. SIGNIFICANCE: This study identifies four genes associated with lung cancer risk, which could help guide future lung cancer prevention and treatment approaches.

Rare deleterious germline variants and risk of lung cancer.

Recent studies suggest that rare variants exhibit stronger effect sizes and might play a crucial role in the etiology of lung cancers (LC). Whole exome plus targeted sequencing of germline DNA was performed on 1045 LC cases and 885 controls in the discovery set. To unveil the inherited causal variants, we focused on rare and predicted deleterious variants and small indels enriched in cases or controls. Promising candidates were further validated in a series of 26,803 LCs and 555,107 controls. During discovery, we identified 25 rare deleterious variants associated with LC susceptibility, including 13 reported in ClinVar. Of the five validated candidates, we discovered two pathogenic variants in known LC susceptibility loci, ATM p.V2716A (Odds Ratio [OR] 19.55, 95%CI 5.04-75.6) and MPZL2 p.I24M frameshift deletion (OR 3.88, 95%CI 1.71-8.8); and three in novel LC susceptibility genes, POMC c.*28delT at 3' UTR (OR 4.33, 95%CI 2.03-9.24), STAU2 p.N364M frameshift deletion (OR 4.48, 95%CI 1.73-11.55), and MLNR p.Q334V frameshift deletion (OR 2.69, 95%CI 1.33-5.43). The potential cancer-promoting role of selected candidate genes and variants was further supported by endogenous DNA damage assays. Our analyses led to the identification of new rare deleterious variants with LC susceptibility. However, in-depth mechanistic studies are still needed to evaluate the pathogenic effects of these specific alleles.

Whole Exome Sequencing of Highly Aggregated Lung Cancer Families Reveals Linked Loci for Increased Cancer Risk on Chromosomes 12q, 7p, and 4q.

BACKGROUND: Lung cancer kills more people than any other cancer in the United States. In addition to environmental factors, lung cancer has genetic risk factors as well, though the genetic etiology is still not well understood. We have performed whole exome sequencing on 262 individuals from 28 extended families with a family history of lung cancer. METHODS: Parametric genetic linkage analysis was performed on these samples using two distinct analyses-the lung cancer only (LCO) analysis, where only patients with lung cancer were coded as affected, and the all aggregated cancers (AAC) analysis, where other cancers seen in the pedigree were coded as affected. RESULTS: The AAC analysis yielded a genome-wide significant result at rs61943670 in POLR3B at 12q23.3. POLR3B has been implicated somatically in lung cancer, but this germline finding is novel and is a significant expression quantitative trait locus in lung tissue. Interesting genome-wide suggestive haplotypes were also found within individual families, particularly near SSPO at 7p36.1 in one family and a large linked haplotype spanning 4q21.3-28.3 in a different family. The 4q haplotype contains potential causal rare variants in DSPP at 4q22.1 and PTPN13 at 4q21.3. CONCLUSIONS: Regions on 12q, 7p, and 4q are linked to increased cancer risk in highly aggregated lung cancer families, 12q across families and 7p and 4q within a single family. POLR3B, SSPO, DSPP, and PTPN13 are currently the best candidate genes. IMPACT: Functional work on these genes is planned for future studies and if confirmed would lead to potential biomarkers for risk in cancer.

A Drosophila protein specific to pheromone-sensing gustatory hairs delays males' copulation attempts.

In insects, increasing evidence suggests that small secreted pheromone binding proteins (PBPs) and odorant binding proteins (OBPs) are important for normal olfactory detection of airborne pheromones and odorants far from their source. In contrast, it is unknown whether extracellular ligand binding proteins participate in perception of less volatile chemicals, including many pheromones, that are detected by direct contact with chemosensory organs. CheB42a, a small Drosophila melanogaster protein unrelated to known PBPs or OBPs, is expressed and likely secreted in only a small subset of gustatory sensilla on males' front legs, the site of gustatory perception of contact pheromones. Here we show that CheB42a is expressed specifically in the sheath cells surrounding the taste neurons expressing Gr68a, a putative gustatory pheromone receptor for female cuticular hydrocarbons that stimulate male courtship. Surprisingly, however, CheB42a mutant males attempt to copulate with females earlier and more frequently than control males. Furthermore, CheB42a mutant males also attempt to copulate more frequently with other males that secrete female-specific cuticular hydrocarbon pheromones, but not with females lacking cuticular hydrocarbons. Together, these data indicate that CheB42a is required for a normal gustatory response to female cuticular hydrocarbon pheromones that modulate male courtship.

Multi-Omics Analysis Reveals a HIF Network and Hub Gene EPAS1 Associated with Lung Adenocarcinoma.

Recent technological advancements have permitted high-throughput measurement of the human genome, epigenome, metabolome, transcriptome, and proteome at the population level. We hypothesized that subsets of genes identified from omic studies might have closely related biological functions and thus might interact directly at the network level. Therefore, we conducted an integrative analysis of multi-omic datasets of non-small cell lung cancer (NSCLC) to search for association patterns beyond the genome and transcriptome. A large, complex, and robust gene network containing well-known lung cancer-related genes, including EGFR and TERT, was identified from combined gene lists for lung adenocarcinoma. Members of the hypoxia-inducible factor (HIF) gene family were at the center of this network. Subsequent sequencing of network hub genes within a subset of samples from the Transdisciplinary Research in Cancer of the Lung-International Lung Cancer Consortium (TRICL-ILCCO) consortium revealed a SNP (rs12614710) in EPAS1 associated with NSCLC that reached genome-wide significance (OR = 1.50; 95% CI: 1.31-1.72; p = 7.75 × 10-9). Using imputed data, we found that this SNP remained significant in the entire TRICL-ILCCO consortium (p = .03). Additional functional studies are warranted to better understand interrelationships among genetic polymorphisms, DNA methylation status, and EPAS1 expression.

Rare Variants in Known Susceptibility Loci and Their Contribution to Risk of Lung Cancer.

BACKGROUND: Genome-wide association studies are widely used to map genomic regions contributing to lung cancer (LC) susceptibility, but they typically do not identify the precise disease-causing genes/variants. To unveil the inherited genetic variants that cause LC, we performed focused exome-sequencing analyses on genes located in 121 genome-wide association study-identified loci previously implicated in the risk of LC, chronic obstructive pulmonary disease, pulmonary function level, and smoking behavior. METHODS: Germline DNA from 260 case patients with LC and 318 controls were sequenced by utilizing VCRome 2.1 exome capture. Filtering was based on enrichment of rare and potential deleterious variants in cases (risk alleles) or controls (protective alleles). Allelic association analyses of single-variant and gene-based burden tests of multiple variants were performed. Promising candidates were tested in two independent validation studies with a total of 1773 case patients and 1123 controls. RESULTS: We identified 48 rare variants with deleterious effects in the discovery analysis and validated 12 of the 43 candidates that were covered in the validation platforms. The top validated candidates included one well-established truncating variant, namely, BRCA2, DNA repair associated gene (BRCA2) K3326X (OR = 2.36, 95% confidence interval [CI]: 1.38-3.99), and three newly identified variations, namely, lymphotoxin beta gene (LTB) p.Leu87Phe (OR = 7.52, 95% CI: 1.01-16.56), prolyl 3-hydroxylase 2 gene (P3H2) p.Gln185His (OR = 5.39, 95% CI: 0.75-15.43), and dishevelled associated activator of morphogenesis 2 gene (DAAM2) p.Asp762Gly (OR = 0.25, 95% CI: 0.10-0.79). Burden tests revealed strong associations between zinc finger protein 93 gene (ZNF93), DAAM2, bromodomain containing 9 gene (BRD9), and the gene LTB and LC susceptibility. CONCLUSION: Our results extend the catalogue of regions associated with LC and highlight the importance of germline rare coding variants in LC susceptibility.

Wnt pathway activation by ADP-ribosylation.

Wnt/ß-catenin signalling directs fundamental processes during metazoan development and can be aberrantly activated in cancer. Wnt stimulation induces the recruitment of the scaffold protein Axin from an inhibitory destruction complex to a stimulatory signalosome. Here we analyse the early effects of Wnt on Axin and find that the ADP-ribose polymerase Tankyrase (Tnks)--known to target Axin for proteolysis-regulates Axin's rapid transition following Wnt stimulation. We demonstrate that the pool of ADP-ribosylated Axin, which is degraded under basal conditions, increases immediately following Wnt stimulation in both Drosophila and human cells. ADP-ribosylation of Axin enhances its interaction with the Wnt co-receptor LRP6, an essential step in signalosome assembly. We suggest that in addition to controlling Axin levels, Tnks-dependent ADP-ribosylation promotes the reprogramming of Axin following Wnt stimulation; and propose that Tnks inhibition blocks Wnt signalling not only by increasing destruction complex activity, but also by impeding signalosome assembly.

Focused Analysis of Exome Sequencing Data for Rare Germline Mutations in Familial and Sporadic Lung Cancer.

INTRODUCTION: The association between smoking-induced chronic obstructive pulmonary disease (COPD) and lung cancer (LC) is well documented. Recent genome-wide association studies (GWAS) have identified 28 susceptibility loci for LC, 10 for COPD, 32 for smoking behavior, and 63 for pulmonary function, totaling 107 nonoverlapping loci. Given that common variants have been found to be associated with LC in genome-wide association studies, exome sequencing of these high-priority regions has great potential to identify novel rare causal variants. METHODS: To search for disease-causing rare germline mutations, we used a variation of the extreme phenotype approach to select 48 patients with sporadic LC who reported histories of heavy smoking-37 of whom also exhibited carefully documented severe COPD (in whom smoking is considered the overwhelming determinant)-and 54 unique familial LC cases from families with at least three first-degree relatives with LC (who are likely enriched for genomic effects). RESULTS: By focusing on exome profiles of the 107 target loci, we identified two key rare mutations. A heterozygous p.Arg696Cys variant in the coiled-coil domain containing 147 (CCDC147) gene at 10q25.1 was identified in one sporadic and two familial cases. The minor allele frequency (MAF) of this variant in the 1000 Genomes database is 0.0026. The p.Val26Met variant in the dopamine ß-hydroxylase (DBH) gene at 9q34.2 was identified in two sporadic cases; the minor allele frequency of this mutation is 0.0034 according to the 1000 Genomes database. We also observed three suggestive rare mutations on 15q25.1: iron-responsive element binding protein neuronal 2 (IREB2); cholinergic receptor, nicotinic, alpha 5 (neuronal) (CHRNA5); and cholinergic receptor, nicotinic, beta 4 (CHRNB4). CONCLUSIONS: Our results demonstrated highly disruptive risk-conferring CCDC147 and DBH mutations.

Sexy DEG/ENaC channels involved in gustatory detection of fruit fly pheromones.

Hydrocarbon pheromones on the cuticle of Drosophila melanogaster modulate the complex courtship behavior of males. Recently, three members of the degenerin/epithelial Na+ channel (DEG/ENaC) family of sodium channel subunits, Ppk25, Ppk23, and Ppk29 (also known as Nope), have been shown to function in gustatory perception of courtship-modulating contact pheromones. All three proteins are required for the activation of male courtship by female pheromones. Specific interactions between two of them have been demonstrated in cultured cells, suggesting that, in a subset of cells where they are coexpressed, these three subunits function within a common heterotrimeric DEG/ENaC channel. Such a DEG/ENaC channel may be gated by pheromones, either directly or indirectly, or alternatively may control the excitability of pheromone-sensing cells. In addition, these studies identify taste neurons that respond specifically to courtship-modulating pheromones and mediate their effects on male behavior. Two types of pheromone-sensing taste neurons, F and M cells, have been defined on the basis of their specific response to either female or male pheromones. These reports set the stage for the dissection of the molecular and cellular mechanisms that mediate gustatory detection of contact pheromones.

Roles of the Drosophila SK channel (dSK) in courtship memory.

A role for SK channels in synaptic plasticity has been very well-characterized. However, in the absence of simple genetic animal models, their role in behavioral memory remains elusive. Here, we take advantage of Drosophila melanogaster with its single SK gene (dSK) and well-established courtship memory assay to investigate the contribution of this channel to memory. Using two independent dSK alleles, a null mutation and a dominant negative subunit, we show that while dSK negatively regulates the acquisition of short-term memory 30 min after a short training session, it is required for normal long-term memory 24 h after extended training. These findings highlight important functions for dSK in courtship memory and suggest that SK channels can mediate multiple forms of behavioral plasticity.

Drosophila CheB proteins involved in gustatory detection of pheromones are related to a human neurodegeneration factor.

The Drosophila CheBs proteins are expressed in a variety of sexually dimorphic subsets of taste hairs, some of which have been directly implicated in pheromone detection. Their remarkable collection of expression patterns suggests that CheBs have specialized roles in gustatory detection of pheromones. Indeed, mutations in the CheB42a gene specifically alter male response to female-specific cuticular hydrocarbons. Furthermore, CheBs belong to the large ML (MD-2-like) superfamily of lipid-binding proteins and share amino acids with an essential role in the function of human GM2-activator protein (GM2-AP), a protein whose absence results in neurodegeneration and death. As GM2-AP binds specifically to the GM2 ganglioside, we have proposed that CheB42a and other CheBs function by interacting directly with the lipid-like cuticular hydrocarbons of Drosophila melanogaster and modulating their detection by transmembrane receptors. Here I review the current knowledge of the CheB family and discuss possible models for their function.

A Drosophila protein family implicated in pheromone perception is related to Tay-Sachs GM2-activator protein.

Low volatility, lipid-like cuticular hydrocarbon pheromones produced by Drosophila melanogaster females play an essential role in triggering and modulating mating behavior, but the chemosensory mechanisms involved remain poorly understood. Recently, we showed that the CheB42a protein, which is expressed in only 10 pheromone-sensing taste hairs on the front legs of males, modulates progression to late stages of male courtship behavior in response to female-specific cuticular hydrocarbons. Here we report that expression of all 12 genes in the CheB gene family is predominantly or exclusively gustatory-specific, and occurs in many different, often non-overlapping patterns. Only the Gr family of gustatory receptor genes displays a comparable variety of gustatory-specific expression patterns. Unlike Grs, however, expression of all but one CheB gene is sexually dimorphic. Like CheB42a, other CheBs may therefore function specifically in gustatory perception of pheromones. We also show that CheBs belong to the ML superfamily of lipid-binding proteins, and are most similar to human GM2-activator protein (GM2-AP). In particular, GM2-AP residues involved in ligand binding are conserved in CheBs but not in other ML proteins. Finally, CheB42a is specifically secreted into the inner lumen of pheromone-sensing taste hairs, where pheromones interact with membrane-bound receptors. We propose that CheB proteins interact directly with lipid-like Drosophila pheromones and modulate their detection by the gustatory signal transduction machinery. Furthermore, as loss of GM2-AP in Tay-Sachs disease prevents degradation of GM2 gangliosides and results in neurodegeneration, the function of CheBs in pheromone response may involve biochemical mechanisms critical for lipid metabolism in human neurons.

A Drosophila DEG/ENaC channel subunit is required for male response to female pheromones.

Odorants and pheromones as well as sweet- and bitter-tasting small molecules are perceived through activation of G protein-coupled chemosensory receptors. In contrast, gustatory detection of salty and sour tastes may involve direct gating of sodium channels of the DEG/ENaC family by sodium and hydrogen ions, respectively. We have found that ppk25, a Drosophila melanogaster gene encoding a DEG/ENaC channel subunit, is expressed at highest levels in the male appendages responsible for gustatory and olfactory detection of female pheromones: the legs, wings, and antennae. Mutations in the ppk25 gene reduce or even abolish male courtship response to females in the dark, conditions under which detection of female pheromones is an essential courtship-activating sensory input. In contrast, the same mutations have no effect on other behaviors tested. Importantly, ppk25 mutant males that show no response to females in the dark execute all of the normal steps of courtship behavior in the presence of visible light, suggesting that ppk25 is required for activation of courtship behavior by chemosensory perception of female pheromones. Finally, a ppk25 mutant allele predicted to encode a truncated protein has dominant-negative properties, suggesting that the normal Ppk25 protein acts as part of a multiprotein complex. Together, these results indicate that ppk25 is necessary for response to female pheromones by D. melanogaster males, and suggest that members of the DEG/ENaC family of genes play a wider role in chemical senses than previously suspected.