The transcriptional map of the common eliminated region 1 (C3CER1) in 3p21.3

Genomic characterization and transcription analysis of European sea bass (Dicentrarchus labrax) rtp3 genes

Fish RTP3, belonging to the receptor-transporting protein family, display several functions, including a putative antiviral role as virus-responsive gene. In this work, we have identified and characterized two different European sea bass rtp3 genes. In addition, an in vivo transcription analysis in response to LPS, poly I:C and betanodavirus infection (RGNNV genotype) has been performed. The sequence analysis showed that European sea bass displays two rtp3 genes, X1 and X2, composed of two exons and a single intron (1007-bp and 888-bp long, respectively), located within the ORF sequence. The full-length cDNA is 1969 bp for rtp3 X1, and 1491 bp for rtp3 X2. Several ATTTA motifs have been found in the intron sequence of both genes, whereas rtp3 X1 also contains this motif in both untranslated regions. The transcription analyses revealed significant level of rtp3 X2 mRNA in brain and head kidney after LPS and poly I:C inoculation; however, the induction elicited by RGNNV infection was much higher, suggesting an essential role for this protein in controlling NNV infections.

GWAS reveals genetic basis of a predisposition to severe COVID-19 through in silico modeling of the FYCO1 protein

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, is heavily reliant on its natural ability to "hack" the host's genetic and biological pathways. The genetic susceptibility of the host is a key factor underlying the severity of the disease. Polygenic risk scores are essential for risk assessment, risk stratification, and the prevention of adverse outcomes. In this study, we aimed to assess and analyze the genetic predisposition to severe COVID-19 in a large representative sample of the Russian population as well as to build a reliable but simple polygenic risk score model with a lower margin of error. Another important goal was to learn more about the pathogenesis of severe COVID-19. We examined the tertiary structure of the FYCO1 protein, the only gene with mutations in its coding region and discovered changes in the coiled-coil domain. Our findings suggest that FYCO1 may accelerate viral intracellular replication and excessive exocytosis and may contribute to an increased risk of severe COVID-19. We found significant associations between COVID-19 and LZTFL1, FYCO1, XCR1, CCR9, TMLHE-AS1, and SCYL2 at 3p21.31. Our findings further demonstrate the polymorphic nature of the severe COVID-19 phenotype.

Comparative proteomics reveals genetic mechanisms of body weight in Hu sheep and Dorper sheep

Body weight (BW) is a critical economic trait for meat production in sheep, and it is a complex trait affected by numerous elements. The aim of this study was to investigate the genetic mechanisms of sheep BW by a label-free proteomics approach. The result showed, a total of 27, 14, 61, and 65 differentially abundant proteins (DAPs) were identified in the Hu_HBW vs. Hu_LBW, DP_HBW vs. DP_LBW, Hu_HBW vs. DP_HBW, and Hu_LBW vs. DP_LBW comparisons, respectively. Five proteins (including ILK, AHCYL2, MLIP, CYB5A, and SMTNL1) related to fat synthesis and muscle development were detected in the Hu sheep group. In the Dorper sheep group, the screened DAPs strictly related to muscle development and fat synthesis were significantly enriched in MAP kinase activity (MAPK12), Arachidonic acid metabolism, and Steroid hormone biosynthesis (PGFS, LOC101107119) pathways. Several DAPs related to immune responses (SERPINA1, FGG, SERPINC1, and LOC101108131), fat deposition (APOH, GC, AHSG, SKP1, ACSL1, ACAT1, and ACADS), and muscle development (LMOD3 and LRRC39) were detected in the Hu vs. Dorper sheep comparison. These analyses indicated that the BW of sheep is regulated via a variety of pathways, and these DAPs can be further investigated as candidate markers for predicting the BW of sheep. SIGNIFICANCE: Body weight is one of the key traits in sheep and involves multiple coordinated regulatory mechanisms, but the genetic mechanism of BW is still unclear in sheep. In the current study, the label-free method was used to identify the proteins and pathways related to BW using LT muscle of Hu sheep and Dorper sheep with different BW. These findings will provide new candidate proteins and vital pathways into the molecular mechanisms involved growth traits in sheep.

FYCO1 Frameshift Deletion in Wirehaired Pointing Griffon Dogs with Juvenile Cataract

Different breed-specific inherited cataracts have been described in dogs. In this study, we investigated an inbred family of Wirehaired Pointing Griffon dogs in which three offspring were affected by juvenile cataract. The pedigree suggested monogenic autosomal recessive inheritance of the trait. Whole-genome sequencing of an affected dog revealed 12 protein-changing variants that were not present in 566 control genomes, of which two were located in functional candidate genes, FYCO1 and CRYGB. Targeted genotyping of both variants in the investigated family excluded CRYGB and revealed perfect co-segregation of the FYCO1 variant with the juvenile cataract phenotype. This variant, FYCO1:c.2024delG, represents a 1 bp frameshift deletion predicted to truncate ~50% of the open reading frame p.(Ser675Thrfs*5). FYCO1 encodes the FYVE and coiled-coil domain autophagy adaptor 1, a known regulator of lens autophagy, which is required for the normal homeostasis in the eye. In humans, at least 37 pathogenic variants in FYCO1 have been shown to cause autosomal recessive cataract. Fcyo1^−/− knockout mice also develop cataracts. Together with the current knowledge on FYCO1 variants and their functional impact in humans and mice, our data strongly suggest FYCO1:c.2024delG as a candidate causative variant for the observed juvenile cataract in Wirehaired Pointing Griffon dogs. To the best of our knowledge, this study represents the first report of a FYCO1-related cataract in domestic animals.

Identification of a novel nonsense variant in FYCO1 gene associated with infantile cataract and cortical atrophy

INTRODUCTION

Cataract is a major condition characterized by ocular lens opacification, resulting from alteration in the lens architecture, lens proteins or both. It is responsible for about one-third of infants' blindness worldwide. Variants in the FYCO1 gene have been associated with autosomal recessive infantile cataract.

MATERIAL AND METHODS

We conducted whole exome sequencing (WES) in a nine months old male patient who was referred for genetic investigation because of infantile cataract. WES analysis revealed the presence of a homozygous pathogenic variant (c.2365C>T) in exon 8 of the FYCO1 gene.

RESULTS AND DISCUSSION

This is the first report on a Lebanese infant with infantile cataract and cortical atrophy which was not previously reported, resulting from a novel homozygous FYCO1 variant; thus expanding the clinical phenotypic spectrum of FYCO1 involvement.

FYCO1 Regulates Cardiomyocyte Autophagy and Prevents Heart Failure Due to Pressure Overload In Vivo

•

FYCO1, a component of the autophagic machinery, is highly expressed in the heart and a potent inducer of cardiomyocyte autophagy.

•

Loss of FYCO1 in vivo inhibits adaptation to starvation or biomechanical stress of the heart by an abrogated increase of autophagic flux and results in contractile dysfunction.

•

Heart specific overexpression of FYCO1 improves autophagic flux and rescues contractile dysfunction following pressure overload.

Autophagy is a cellular degradation process that has been implicated in diverse disease processes. The authors provide evidence that FYCO1, a component of the autophagic machinery, is essential for adaptation to cardiac stress. Although the absence of FYCO1 does not affect basal autophagy in isolated cardiomyocytes, it abolishes induction of autophagy after glucose deprivation. Likewise, Fyco1-deficient mice subjected to starvation or pressure overload are unable to respond with induction of autophagy and develop impaired cardiac function. FYCO1 overexpression leads to induction of autophagy in isolated cardiomyocytes and transgenic mouse hearts, thereby rescuing cardiac dysfunction in response to biomechanical stress.

Epigenetics of Skeletal Muscle-Associated Genes in the ASB, LRRC, TMEM, and OSBPL Gene Families

Much remains to be discovered about the intersection of tissue-specific transcription control and the epigenetics of skeletal muscle (SkM), a very complex and dynamic organ. From four gene families, Leucine-Rich Repeat Containing (LRRC), Oxysterol Binding Protein Like (OSBPL), Ankyrin Repeat and Socs Box (ASB), and Transmembrane Protein (TMEM), we chose 21 genes that are preferentially expressed in human SkM relative to 52 other tissue types and analyzed relationships between their tissue-specific epigenetics and expression. We also compared their genetics, proteomics, and descriptions in the literature. For this study, we identified genes with little or no previous descriptions of SkM functionality (ASB4, ASB8, ASB10, ASB12, ASB16, LRRC14B, LRRC20, LRRC30, TMEM52, TMEM233, OSBPL6/ORP6, and OSBPL11/ORP11) and included genes whose SkM functions had been previously addressed (ASB2, ASB5, ASB11, ASB15, LRRC2, LRRC38, LRRC39, TMEM38A/TRIC-A, and TMEM38B/TRIC-B). Some of these genes have associations with SkM or heart disease, cancer, bone disease, or other diseases. Among the transcription-related SkM epigenetic features that we identified were: super-enhancers, promoter DNA hypomethylation, lengthening of constitutive low-methylated promoter regions, and SkM-related enhancers for one gene embedded in a neighboring gene (e.g., ASB8-PFKM, LRRC39-DBT, and LRRC14B-PLEKHG4B gene-pairs). In addition, highly or lowly co-expressed long non-coding RNA (lncRNA) genes probably regulate several of these genes. Our findings give insights into tissue-specific epigenetic patterns and functionality of related genes in a gene family and can elucidate normal and disease-related regulation of gene expression in SkM.

Characterization of a novel disease resistance gene rtp3 and its association with VNN disease resistance in Asian seabass

Asian seabass, an important food fish in Southeast Asia, has suffered from nervous necrosis virus (NNV) infection, resulting in massive mortality of Asian seabass larvae and enormous economic losses. Identification and characterization of disease resistance genes is important. Previous transcriptome analysis of Asians seabass epithelial cells after NNV infection revealed a highly inducible gene, receptor-transporting protein 3 (rtp3), indicating it could play an important role in Asian seabass - NNV interaction. To characterize this gene, we determined its expression pattern and subcellular localization. The rtp3 was highly induced in most examined tissues and organs of Asian seabass after NNV infection, and protein Rtp3 was localized in cytoplasm. Further association study in multiple families revealed that a microsatellite marker, (GT)ntt(GT)n, in the 3' UTR of rtp3 was significantly associated with VNN disease resistance in Asian seabass. Our results imply that rtp3 may be a novel disease resistance gene in Asian seabass. This data could improve our understanding of molecular interaction between Asian seabass and NNV, and has the potential to be applied in marker-assisted selection for disease resistance breeding in Asian seabass.

A nuclear-directed human pancreatic ribonuclease (PE5) targets the metabolic phenotype of cancer cells

Ribonucleases represent a new class of antitumor RNA-damaging drugs. However, many wild-type members of the vertebrate secreted ribonuclease family are not cytotoxic because they are not able to evade the cytosolic ribonuclease inhibitor. We previously engineered the human pancreatic ribonuclease to direct it to the cell nucleus where the inhibitor is not present. The best characterized variant is PE5 that kills cancer cells through apoptosis mediated by the p21^WAF1/CIP1 induction and the inactivation of JNK. Here, we have used microarray-derived transcriptional profiling to identify PE5 regulated genes on the NCI/ADR-RES ovarian cancer cell line. RT-qPCR analyses have confirmed the expression microarray findings. The results show that PE5 cause pleiotropic effects. Among them, it is remarkable the down-regulation of multiple genes that code for enzymes involved in deregulated metabolic pathways in cancer cells.

Expression of pre-selected TMEMs with predicted ER localization as potential classifiers of ccRCC tumors

Background

VHL inactivation is the most established molecular characteristic of clear cell renal cell carcinoma (ccRCC), with only a few additional genes implicated in development of this kidney tumor. In recently published ccRCC gene expression meta-analysis study we identified a number of deregulated genes with limited information available concerning their biological role, represented by gene transcripts belonging to transmembrane proteins family (TMEMs). TMEMs are predicted to be components of cellular membranes, such as mitochondrial membranes, ER, lysosomes and Golgi apparatus. Interestingly, the function of majority of TMEMs remains unclear. Here, we analyzed expression of ten TMEM genes in the context of ccRCC progression and development, and characterized these proteins bioinformatically.

Methods

The expression of ten TMEMs (RTP3, SLC35G2, TMEM30B, TMEM45A, TMEM45B, TMEM61, TMEM72, TMEM116, TMEM207 and TMEM213) was measured by qPCR. T-test, Pearson correlation, univariate and multivariate logistic and Cox regression were used in statistical analysis. The topology of studied proteins was predicted with Metaserver, together with PSORTII, Pfam and Localizome tools.

Results

We observed significant deregulation of expression of 10 analyzed TMEMs in ccRCC tumors. Cluster analysis of expression data suggested the down-regulation of all tested TMEMs to be a descriptor of the most advanced tumors. Logistic and Cox regression potentially linked TMEM expression to clinical parameters such as: metastasis, Fuhrman grade and overall survival. Topology predictions classified majority of analyzed TMEMs as type 3 and type 1 transmembrane proteins, with predicted localization mainly in ER.

Conclusions

The massive down-regulation of expression of TMEM family members suggests their importance in the pathogenesis of ccRCC and the bioinformatic analysis of TMEM topology implies a significant involvement of ER proteins in ccRCC pathology.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-015-1530-4) contains supplementary material, which is available to authorized users.

VEGFR-1 Pseudogene Expression and Regulatory Function in Human Colorectal Cancer Cells

A large number of pseudogenes have been found to be transcribed in human cancers. However, only a few pseudogenes are functionally characterized. Here, we identified a transcribed pseudogene of VEGFR1, or fms-related tyrosine kinase 1 (FLT1), in human colorectal cancer cells. Interestingly, this pseudogene (designated as FLT1P1) was found to be transcribed bidirectionally and functionally modulated cognate VEGFR1 protein expression in the cells. Mechanistically, expression of FLT1P1 antisense transcript not only inhibited the VEGFR1 expression, but also inhibited non-cognate VEGF-A expression through interaction with miR-520a. Perturbation of FLT1P1 expression by RNA interference (RNAi) markedly inhibited tumor cell proliferation and xenograft tumor growth. This study identifies FLT1P1 antisense as a critical regulator of VEGFR1 and VEGF-A expression in colorectal cancer cells, and highlights its role in regulation of the pathogenesis of colorectal cancer.

IMPLICATIONS

The VEGFR1 pseudogene, FLT1P1, is a novel and functional regulator of VEGF signaling and its targeting could be an alternative strategy to modulate its cognate/target gene expression and downstream activity in cancer.

Mutations in FYCO1 cause autosomal-recessive congenital cataracts

Congenital cataracts (CCs), responsible for about one-third of blindness in infants, are a major cause of vision loss in children worldwide. Autosomal-recessive congenital cataracts (arCC) form a clinically diverse and genetically heterogeneous group of disorders of the crystalline lens. To identify the genetic cause of arCC in consanguineous Pakistani families, we performed genome-wide linkage analysis and fine mapping and identified linkage to 3p21-p22 with a summed LOD score of 33.42. Mutations in the gene encoding FYVE and coiled-coil domain containing 1 (FYCO1), a PI(3)P-binding protein family member that is associated with the exterior of autophagosomes and mediates microtubule plus-end-directed vesicle transport, were identified in 12 Pakistani families and one Arab Israeli family in which arCC had previously been mapped to the overlapping CATC2 region. Nine different mutations were identified, including c.3755 delC (p.Ala1252AspfsX71), c.3858_3862dupGGAAT (p.Leu1288TrpfsX37), c.1045 C>T (p.Gln349X), c.2206C>T (p.Gln736X), c.2761C>T (p.Arg921X), c.2830C>T (p.Arg944X), c.3150+1 G>T, c.4127T>C (p.Leu1376Pro), and c.1546C>T (p.Gln516X). Fyco1 is expressed in the mouse embryonic and adult lens and peaks at P12d. Expressed mutant proteins p.Leu1288TrpfsX37 and p.Gln736X are truncated on immunoblots. Wild-type and p.L1376P FYCO1, the only missense mutant identified, migrate at the expected molecular mass. Both wild-type and p. Leu1376Pro FYCO1 proteins expressed in human lens epithelial cells partially colocalize to microtubules and are found adjacent to Golgi, but they primarily colocalize to autophagosomes. Thus, FYCO1 is involved in lens development and transparency in humans, and mutations in this gene are one of the most common causes of arCC in the Pakistani population.

FYCO1 is a Rab7 effector that binds to LC3 and PI3P to mediate microtubule plus end-directed vesicle transport

FYCO1 recognition of LC3 on autophagosomes facilitates microtubule-mediated cytosolic transport of this degradative organelle.

Autophagy is the main eukaryotic degradation pathway for long-lived proteins, protein aggregates, and cytosolic organelles. Although the protein machinery involved in the biogenesis of autophagic vesicles is well described, very little is known about the mechanism of cytosolic transport of autophagosomes. In this study, we have identified an adaptor protein complex, formed by the two autophagic membrane-associated proteins LC3 and Rab7 and the novel FYVE and coiled-coil (CC) domain–containing protein FYCO1, that promotes microtubule (MT) plus end–directed transport of autophagic vesicles. We have characterized the LC3-, Rab7-, and phosphatidylinositol-3-phosphate–binding domains in FYCO1 and mapped part of the CC region essential for MT plus end–directed transport. We also propose a mechanism for selective autophagosomal membrane recruitment of FYCO1.

Genome-wide association study for femoral neck bone geometry

Poor femoral neck bone geometry at the femur is an important risk factor for hip fracture. We conducted a genome-wide association study (GWAS) of femoral neck bone geometry, examining approximately 379,000 eligible single-nucleotide polymorphisms (SNPs) in 1000 Caucasians. A common genetic variant, rs7430431 in the receptor transporting protein 3 (RTP3) gene, was identified in strong association with the buckling ratio (BR, P = 1.6 x 10(-7)), an index of bone structural instability, and with femoral cortical thickness (CT, P = 1.9 x 10(-6)). The RTP3 gene is located in 3p21.31, a region that we found to be linked with CT (LOD = 2.19, P = 6.0 x 10(-4)) in 3998 individuals from 434 pedigrees. The replication analyses in 1488 independent Caucasians and 2118 Chinese confirmed the association of rs7430431 to BR and CT (combined P = 7.0 x 10(-3) for BR and P = 1.4 x 10(-2) for CT). In addition, 350 hip fracture patients and 350 healthy control individuals were genotyped to assess the association of the RTP3 gene with the risk of hip fracture. Significant association between a nearby common SNP, rs10514713 of the RTP3 gene, and hip fracture (P = 1.0 x 10(-3)) was found. Our observations suggest that RTP3 may be a novel candidate gene for femoral neck bone geometry.

The interferon-alpha responsive gene TMEM7 suppresses cell proliferation and is downregulated in human hepatocellular carcinoma

Multiple regions on the chromosome arm 3p are frequently affected by loss of heterozygosity in human cancers. A candidate tumor suppressor gene is TMEM7, at 3p21.3, which encodes a transmembrane protein. TMEM7 is expressed specifically in the liver, and the encoded protein shares substantial sequence homology with human and mouse 28-kDa interferon-alpha (IFN-alpha) responsive protein. In investigation of the possible role of TMEM7 in development of hepatocellular carcinoma (HCC), we examined TMEM7 expression in 20 primary HCC and 18 HCC cell lines and found recurrent functional alterations. Although TMEM7 mRNA was expressed in normal hepatic cells, downregulation or inactivation of the gene was detected in 85% of primary HCC and 33% of HCC cell lines. To identify the mechanisms responsible, we examined genomic deletion and mutation, and also the effect of inhibitors of DNA methyltransferase and histone deacetylase on cells with low or no endogenous TMEM7 expression. Homozygous deletion of TMEM7 was not detected in 17 pairs of human HCC and corresponding noncancerous liver tissues, nor in any of the 18 HCC cell lines. TMEM7 mutation was not detected in the 18 HCC cell lines (low or normal TMEM7 expression). Treatment of two of six cell lines exhibiting downregulation or loss of TMEM7 with 5-aza-2'-deoxycytidine and trichostatin A yielded additive increase in TMEM7 expression, implicating aberrant DNA methylation and histone deacetylation in transcriptional silencing of this gene. Ectopic expression of TMEM7 in two TMEM7-deficient HCC lines suppressed cell proliferation, colony formation, and cell migration in vitro and reduced tumor formation in nude mice. Treatment of two highly invasive HCC cell lines with IFN-alpha for 7 days significantly increased TMEM7 expression and inhibited cell migration. These findings implicate loss of TMEM7 expression in hepatocarcinogenesis and suggest that modification of TMEM7 expression by IFN-alpha may have therapeutic relevance in a subset of HCC.

Genetic abnormalities involved in t(12;21) TEL-AML1 acute lymphoblastic leukemia: analysis by means of array-based comparative genomic hybridization

The TEL (ETV6)-AML1 (RUNX1) chimeric gene fusion is the most common genetic abnormality in childhood acute lymphoblastic leukemias. Evidence suggests that this chimeric gene fusion constitutes an initiating mutation that is necessary but insufficient for the development of leukemia. In a search for additional genetic events that could be linked to the development of leukemia, we applied a genome-wide array-comparative genomic hybridization technique to 24 TEL-AML1 leukemia samples and two cell lines. It was found that at least two chromosomal imbalances were involved in all samples. Recurrent regions of chromosomal imbalance (>10% of cases) and representative involved genes were gain of chromosomes 10 (17%) and 21q (25%; RUNX1) and loss of 12p13.2 (87%; TEL), 9p21.3 (29%; p16INK4a/ARF), 9p13.2 (25%; PAX5), 12q21.3 (25%; BTG1), 3p21 (21%; LIMD1), 6q21 (17%; AIM1 and BLIMP1), 4q31.23 (17%; NR3C2), 11q22-q23 (13%; ATM) and 19q13.11-q13.12 (13%; PDCD5). Enforced expression of TEL and to a lesser extent BTG1, both single genes known to be located in their respective minimum common region of loss, inhibited proliferation of the TEL-AML1 cell line Reh. Together, these findings suggest that some of the genes identified as lost by array-comparative genomic hybridization may partly account for the development of leukemia.

Medulloblastomas Derived fromCxcr6Mutant Mice Respond to Treatment with a Smoothened Inhibitor

The sonic hedgehog (Shh) pathway is activated in approximately 30% of human medulloblastoma resulting in increased expression of downstream target genes. In about half of these cases, this has been shown to be a consequence of mutations in regulatory genes within the pathway, including Ptc1, Smo, and Sufu. However, for some tumors, no mutations have been detected in known pathway genes. This suggests that either mutations in other genes promote tumorigenesis or that epigenetic alterations increase pathway activity in these tumors. Here, we report that 3% to 4% of mice lacking either one or both functional copies of Cxcr6 develop medulloblastoma. Although CXCR6 is not known to be involved in Shh signaling, tumors derived from Cxcr6 mutant mice expressed Shh pathway target genes including Gli1, Gli2, Ptc2, and Sfrp1, indicating elevated pathway activity. Interestingly, the level of Ptc1 expression was decreased in tumor cells although two normal copies of Ptc1 were retained. This implies that reduced CXCR6 function leads to suppression of Ptc1 thereby increasing Smoothened function and promoting tumorigenesis. We used a direct transplant model to test the sensitivity of medulloblastoma arising in Cxcr6 mutant mice to a small-molecule inhibitor of Smoothened (HhAntag). We found that transplanted tumors were dramatically inhibited in mice treated for only 4 days with HhAntag. These findings suggest that HhAntag may be effective against tumors lacking mutations in known Shh pathway genes.

Differential gene expression profiles in the hippocampus of senescence-accelerated mouse

The senescence-accelerated mouse (SAM) is an animal model for studying senescence and age-associated disorders due to its inherited aging phenotype. The SAM/prone8 (SAMP8) is a useful animal model to investigate the fundamental mechanisms involved in age-related learning and memory deficits that may have relevance to age-associated AD, while SAM/resistant1 (SAMR1) shows normal. To identify genes rendering the cognitive deterioration with aging, the subtractive cDNA libraries containing 1924 clones with the positive ratio of 96.18% were generated and the microarray containing 3136 cDNA was prepared. The results of screening libraries by the microarray showed that of all 91 differentially expressed genes, 50 were over-expressed and 41 were low-expressed in SAMP8. Some of the identified genes were confirmed by the real time quantitative RT-PCR. These results indicated the profiles of gene expression in the hippocampus of SAMP8 and SAMR1 were significantly different, which may play important roles in the age-related cognitive deficit in SAMP8, suggesting those genes related to the cognitive deficient or pathology change in the brain of SAMP8 may be potential gene targets for Alzheimer's disease therapy.

Modeling non-random deletions in cancer

Chromosome deletions do abound in cancer and are detected in certain regions in a non-random manner. Although their relevance remains elusive, it is a general agreement that segmental losses provide the cell with selective growth advantage. Consequently these may contain genes and/or regulatory sequences that control normal growth and inhibit malignancy. We have developed a monochromosomal hybrid based experimental model for the generation and functional analysis of deletions, that is called "elimination test" (Et). Focused on human chromosome 3 - that was known to carry multiple 3p deletions - the Et was expected to restrict a 3p tumor suppressor region to a sufficiently small segment that permits the selection of a critically important candidate gene. Surprisingly, we detected three regions that were lost in all or majority of tumors: CER1 (3p21.3, Mb: 43.32-45.74), CER2 (3p22, Mb: 37.83-39.06) and FER (3p14.3-p21.2, Mb: 50.12-58.03). In contrast a 3q26-qter region (CRR) was regularly retained. CER1 - our main focus - contains multiple genes that may inhibit tumor growth, but 3 genes, RIS1, LF (LTF) and LIMD1 have already the necessary experimental support to be considered bona fide tumor suppressors. Tumor suppressor region borders display instability features including: (1) they break in evolution and in tumors, (2) they evolve horizontally, and (3) they are enriched with pseudogene insertions. The most remarkable features at the breakpoint cluster regions were segmental duplications that drive horizontal evolution and contribute to cancer associated instability.

Predicting Patient Survival from Microarray Data by Accelerated Failure Time Modeling Using Partial Least Squares and LASSO

We consider the problem of predicting survival times of cancer patients from the gene expression profiles of their tumor samples via linear regression modeling of log-transformed failure times. The partial least squares (PLS) and least absolute shrinkage and selection operator (LASSO) methodologies are used for this purpose where we first modify the data to account for censoring. Three approaches of handling right censored data-reweighting, mean imputation, and multiple imputation-are considered. Their performances are examined in a detailed simulation study and compared with that of full data PLS and LASSO had there been no censoring. A major objective of this article is to investigate the performances of PLS and LASSO in the context of microarray data where the number of covariates is very large and there are extremely few samples. We demonstrate that LASSO outperforms PLS in terms of prediction error when the list of covariates includes a moderate to large percentage of useless or noise variables; otherwise, PLS may outperform LASSO. For a moderate sample size (100 with 10,000 covariates), LASSO performed better than a no covariate model (or noise-based prediction). The mean imputation method appears to best track the performance of the full data PLS or LASSO. The mean imputation scheme is used on an existing data set on lung cancer. This reanalysis using the mean imputed PLS and LASSO identifies a number of genes that were known to be related to cancer or tumor activities from previous studies.

Multipoint interphase FISH analysis of chromosome 3 abnormalities in 28 childhood AML patients

We detected non-random 3p losses and 3q gains on well-determined regions in both murine and human tumors using a microcell hybrid-based model system called 'elimination test'. We suggest that these are general malignancy-associated aberrations not necessarily linked to a particular tissue of origin. To examine chromosome 3 abnormalities, in 28 childhood acute myeloid leukemia bone marrow samples, we performed interphase multipoint-fluorescence in situ hybridization using 84 chromosome 3-specific probes and detected clonal chromosome 3 aberrations in nine cases, which is of a higher frequency than the previously reported one. In 3/28 children, a chromosome 3 abnormality was detected which was not visible using conventional cytogenetic analysis. We did not detect any 3p deletion. Increased copy number of 3q was found in four cases with trisomy of whole chromosome 3 and one case with 3q tetrasomy (isodisomy). We identified rare structural rearrangements in childhood acute myeloblastic leukemia, involving 3q21 and 3q26 loci around RPN1 and MDS1/EVI1 respectively. The poor outcome in pediatric patients with 3q rearrangements appears to be quite uniform.

Evolutionarily plastic regions at human 3p21.3 coincide with tumor breakpoints identified by the "elimination test"

We have previously found with the microcell hybrid-based "elimination test" that human chromosome 3 transferred into murine or human tumor cells regularly lost certain 3p regions during tumor growth in SCID mice. The most common eliminated region, CER1, is approximately 2.4 Mb at 3p21.3. CER1 breakpoints were clustered in approximately 200-kb regions at both telomeric and centromeric borders. We have also shown, earlier, that tumor-related deletions often coincide with human/mouse synteny breakpoints on 3p12-p22. Here we describe the results of a comparative genomic analysis on the CER1 region in Caenorhabditis elegans, Drosophila melanogaster, Fugu rubripes, Gallus gallus, Mus musculus, Rattus norvegicus, and Canis familiaris. First, four independent synteny breaks were found within the CER1 telomeric breakpoint cluster region, comparing human, dog, and chicken genomes, and two independent synteny breaks within the CER1 centromeric breakpoint cluster region, comparing human, mouse, and chicken genomes, suggesting a nonrandom involvement of tumor breakpoint regions in chromosome evolution. Second, both CER1 breakpoint cluster regions show recent tandem duplications (seven Zn finger protein family genes at the telomeric and eight chemokine receptor genes at the centromeric side). Finally, all genes from these regions underwent horizontal evolution in mammals, with formation of new genes and expansion of gene families, which were displayed in the human genome as tandem gene duplications and pseudogene insertions. In contrast the CER1 middle region contained evolutionarily well-conserved solitary genes and a minimal amount of retroposed genes. The coincidence of evolutionary plasticity with CER1 breakpoints may suggest that regional structural instability is expressed in both evolutionary and cancer-associated chromosome rearrangements.

Down regulation of 3p genes,LTF, SLC38A3 andDRR1, upon growth of human chromosome 3–mouse fibrosarcoma hybrids in severe combined immunodeficiency mice

We have applied a functional test for tumour antagonizing genes based on human chromosome 3 (chr3)-mouse fibrosarcoma A9 MCHs that were studied in vitro and after growth as tumours in severe combined immunodeficiency (SCID) mice. Previously, we reported that 9 out of the 36 SCID-tumours maintained the transferred chr3 ("chr3+" tumours), but lost the expression of the known human TSG fragile histidine triad gene (FHIT) in contrast to 14 other 3p-genes examined. Here we report the results of the duplex RT-PCR analysis of 9 "chr3+" tumours and 3 parental MCHs. We have examined the expression of 34 human 3p-genes from known cancer-related regions of instability, including 13 genes from CER1 defined by us previously at 3p21.33-p21.31 and 10 genes from the LUCA region at 3p21.31. We have found that in addition to FHIT, expression of the LTF gene from CER1 at 3p21.33-p21.31 was lost in all 9 tumours analyzed. The transcript of the solute carrier family 38 member 3 gene (SLC38A3) gene from LUCA region at 3p21.31 was not found in 8 and was greatly reduced in 1 out of these 9 tumours. Expression of the down-regulated in renal cell carcinoma gene (DRR1) gene at 3p14.2 was lost in 7 and down regulated in 2 "chr3+" tumours. In the SCID-tumour derived cell lines treatment with 5-aza-2'-deoxycytidine restored the mRNA expression of LTF, indicating the integrity of DNA sequences. Notably that transcription of the LTF and 2 flanking genes, LRRC2 and TMEM7, as well as transcription of the SLC38A3 gene, were also impaired in all 5 RCC cell lines analyzed. Our data indicate these genes as putative tumour suppressor genes.

Luminal kidney and intestine SLC6 amino acid transporters of B0AT-cluster and their tissue distribution in Mus musculus

The B degrees transport system mediates the Na(+)-driven uptake of a broad range of neutral amino acids into epithelial cells of small intestine and kidney proximal tubule. A corresponding transporter was identified in 2004 (A. Broer, K. Klingel, S. Kowalczuk, J. E. Rasko, J. Cavanaugh, and S. Broer. J Biol Chem 279: 24467-24476, 2004) within the SLC6 family and named B degrees AT1 (SLC6A19). A phylogenetically related transporter known as XT3 in human (SLC6A20) and XT3s1 in mouse was shown to function as an imino acid transporter, to localize also to kidney and small intestine and renamed SIT1 or Imino(B). Besides these two transporters with known functions, there are two other gene products belonging to the same phylogenetic B degrees AT-cluster, XT2 (SLC6A18) and rodent XT3 that are still "orphans." Quantitative real-time RT-PCR showed that the mRNAs of the four B degrees AT-cluster members are abundant in kidney, whereas only those of B degrees AT1 and XT3s1/SIT1 are elevated in small intestine. In brain, the XT3s1/SIT1 mRNA is more abundant than the other B degrees AT-cluster mRNAs. We show here by immunofluorescence that all four mouse B degrees AT-cluster transporters localize, with differential axial gradients, to the brush-border membrane of proximal kidney tubule and, with the possible exception of XT3, also of intestine. Deglycosylation and Western blotting of brush-border proteins demonstrated the glycosylation and confirmed the luminal localization of B degrees AT1, XT2, and XT3. In summary, this study shows the luminal brush-border localization of the Na(+)-dependent amino and imino acid transporters B degrees AT1 and XT3s1/SIT1 in kidney and intestine. It also shows that the structurally highly similar orphan transporters XT2 and XT3 have the same luminal but a slightly differing axial localization along the kidney proximal tubule.

Absence of fatty liver in familial hypobetalipoproteinemia linked to chromosome 3p21

Our aim was to ascertain whether fatty liver may be present in the genetic form of familial hypobetalipoproteinemia (FHBL) linked to a susceptibility locus on chromosome 3p21. Three genetic forms of FHBL exist: (a) FHBL caused by truncation-specifying mutations of apolipoprotein B (apoB), (b) FHBL linked to chr3p21, and (c) FHBL not linked either to APOB or to chr3p21. Fatty liver is common in apoB-defective FHBL. Hepatic fat contents were quantified by magnetic resonance spectroscopy in 16 subjects with 3p21-linked FHBL, 32 subjects with apoB-defective FHBL, and 39 sex- and age-matched controls. Mean liver fat of 3p21 subjects was similar to controls and approximately 60% lower than apoB-defective FHBL subjects ( P = .0012). Indices of adiposity (body mass index, waist/hip ratio) and masses of abdominal subcutaneous, retroperitoneal, and intraperitoneal adipose tissue (IPAT) were quantified by MR imaging. Mean measures of adiposity were similar in the 3 groups, suggesting that adiposity per se was not responsible for differences in liver fat. Liver fat content was positively correlated with IPAT. The intercepts of regression lines of IPAT on liver fat content were similar in controls and 3p21, but higher in apoB-defective FHBL subjects. The slopes of the lines were steepest in apoB-defective, intermediate in 3p21, and flattest in controls. Lipoprotein profiles and very low density lipoprotein-apoB100 kinetics of 3p21 and apoB-defective groups also differed. Thus, 2 genetic subtypes of FHBL also differ in several phenotypic features.

LIM domains-containing protein 1 (LIMD1), a tumor suppressor encoded at chromosome 3p21.3, binds pRB and represses E2F-driven transcription

LIM domains-containing protein 1 (LIMD1) is encoded at chromosome 3p21.3, a region commonly deleted in many solid malignancies. However, the function of LIMD1 is unknown. Here we show that LIMD1 specifically interacts with retinoblastoma protein (pRB), inhibits E2F-mediated transcription, and suppresses the expression of the majority of genes with E2F1-responsive elements. LIMD1 blocks tumor growth in vitro and in vivo and is down-regulated in the majority of human lung cancer samples tested. Our data indicate that LIMD1 is a tumor-suppressor gene, the protein product of which functionally interacts with pRB and the loss of which promotes lung carcinogenesis.

Search for unknown tumor-antagonizing genes

Following the ingenious prediction of Alfred Knudson in 1971, the first tumor suppressor gene, RB1, has been isolated. Its product, the RB1 protein, was found to play a major role in the control of the cell cycle. The loss of heterozygosity (LOH) technique, introduced by Cavenee and colleagues, was an important milestone toward the confirmation of Knudson's hypothesis and the identification of the gene. Subsequently, the LOH technique has provided important clues that have led to the discovery of other tumor suppressor genes. Most of them play important roles in the regulation of the cell cycle and/or of apoptosis. Circumstantial evidence suggests that still other and perhaps many unknown genes may participate in the protection of the organism against malignant growth. The numerous genome losses in tumors, detected by LOH, comparative genomic hybridization, and by cytogenetic techniques, support this possibility. The early work of one of us (G.K.), together with Henry Harris and Francis Wiener, had shown that the malignant phenotype can be suppressed by hybridizing malignant with low- or non-tumorigenic cells. However, analysis of this phenomenon failed to assign the inhibition of tumorigenicity to any particular gene. We have pursued the search for new tumor-antagonizing genes with two unconventional approaches, focusing on human chromosomal subband 3p21.3, a region frequently targeted by cytogenetically detectable deletions. We have detected four clusters of candidate tumor suppressor genes at 3p21.3 by a combination of deletion mapping and the "elimination test." These findings raise the question whether the number and variety of genes that may contribute to the defense against uncontrolled proliferation may have been underestimated.

Interstitial deletions including chromosome 3 common eliminated region 1 (C3CER1) prevail in human solid tumors from 10 different tissues

A human chromosomal segment regularly lost during tumor formation of microcell hybrids in SCID mice has been mapped to 3p21.3. This segment, called chromosome 3 common eliminated region 1 (C3CER1, also referred to as CER1), may harbor multiple tumor-suppressor genes. Because it was found that similar regions were eliminated in an inter- and intraspecies system and in two tumor types (mouse fibrosarcoma and human renal cell carcinoma), we hypothesized that the importance of C3CER1 would transgress tissue specificity, that is, it could occur in tumors derived from multiple tissues. To evaluate the loss of C3CER1 in various human tumor types, we conducted loss of heterozygosity (LOH) analysis of 576 human solid tumors from 10 different tissues and compared the frequency of deletion in the C3CER1 area to that in two other regions on 3p: the FHIT/FRA3B region, at 3p14.2, and the VHL region, at 3p25.3. Deletions were detected in the C3CER1 region in 83% of informative tumors. Half (47%) the LOH-positive tumors showed LOH at all informative markers, indicating a large deletion. The other half (53%) had a discontinuous LOH pattern, suggesting interstitial deletions or breakpoints. The proportion of tumors with C3CER1 deletions was high in all tumor types investigated, ranging from 70% to 94%, except for the soft-tissue sarcomas (40%). In the VHL and FHIT regions, deletions were observed in 73% and 43%, respectively, of the tumors. Of the three 3p regions analyzed, the highest deletion frequency was observed in the C3CER1 region. Furthermore, we demonstrated that the interstitial deletions including C3CER1 prevail over 3p14.2-pter losses in solid tumors.

The human phosphatidylinositol phosphatase SAC1 interacts with the coatomer I complex

The Saccharomyces cerevisiae SAC1 gene encodes an integral membrane protein of the endoplasmic reticulum (ER) and the Golgi apparatus. Yeast SAC1 mutants display a wide array of phenotypes including inositol auxotrophy, cold sensitivity, secretory defects, disturbed ATP transport into the ER, or suppression of actin gene mutations. At present, it is not clear how these phenotypes relate to the finding that SAC1 displays polyphosphoinositide phosphatase activity. Moreover, it is still an open question whether SAC1 functions similarly in mammalian cells, since some phenotypes are yeast-specific. Potential protein interaction partners and, connected to that, possible regulatory circuits have not been described. Therefore, we have cloned human SAC1 (hSAC1), show that it behaves similar to ySac1p in terms of substrate specificity, demonstrate that the endogenous protein localizes to the ER and Golgi, and identify for the first time members of the coatomer I (COPI) complex as interaction partners of hSAC1. Mutation of a putative COPI interaction motif (KXKXX) at its C terminus abolishes interaction with COPI and causes accumulation of hSAC1 in the Golgi. In addition, we generated a catalytically inactive mutant, demonstrate that its lipid binding capacity is unaltered, and show that it accumulates in the Golgi, incapable of interacting with the COPI complex despite the presence of the KXKXX motif. These results open the possibility that the enzymatic function of hSAC1 provides a switch for accessibility of the COPI interaction motif.

Pooled analysis of loss of heterozygosity in breast cancer: a genome scan provides comparative evidence for multiple tumor suppressors and identifies novel candidate regions

Somatic loss of heterozygosity (LOH) has been widely reported in breast cancer as a means of identifying putative tumor-suppressor genes. However, individual studies have rarely spanned more than a single chromosome, and the varying criteria used to declare LOH complicate efforts to formally differentiate regions of consistent versus sporadic (random) loss. We report here the compilation of an extensive database from 151 published LOH studies of breast cancer, with summary data from >15,000 tumors and primary allelotypes from >4,300 tumors. Allelic loss was evaluated at 1,168 marker loci, with large variation in the density of informative observations across the genome. Using studies in which primary allelotype information was available, we employed a likelihood-based approach with a formal chromosomal instability and selection model. The approach seeks direct evidence for preferential loss at each locus compared with nearby loci, accounts for heterogeneity across studies, and enables the direct comparison of candidate regions across the genome. Striking preferential loss was observed (in descending order of significance) in specific regions of chromosomes 7q, 16q, 13q, 17p, 8p, 21q, 3p, 18q, 2q, and 19p, as well as other regions, in many cases coinciding with previously identified candidate genes or known fragile sites. Many of these observations were not possible from any single LOH study, and our results suggest that many previously reported LOH results are not systematic or reproducible. Our approach provides a comparative framework for further investigation of regions exhibiting LOH and identifies broad genomic regions for which there exist few data.

Coincidence of synteny breakpoints with malignancy-related deletions on human chromosome 3

We have found previously that during tumor growth intact human chromosome 3 transferred into tumor cells regularly looses certain 3p regions, among them the approximately 1.4-Mb common eliminated region 1 (CER1) at 3p21.3. Fluorescence in situ hybridization analysis of 12 mouse orthologous loci revealed that CER1 splits into two segments in mouse and therefore contains a murine/human conservation breakpoint region (CBR). Several breaks occurred in tumors within the region surrounding the CBR, and this sequence has features that characterize unstable chromosomal regions: deletions in yeast artificial chromosome clones, late replication, gene and segment duplications, and pseudogene insertions. Sequence analysis of the entire 3p12-22 revealed that other cancer-associated deletions (regions eliminated from monochromosomal hybrids carrying an intact chromosome 3 during tumor growth and homozygous deletions found in human tumors) colocalized nonrandomly with murine/human CBRs and were characterized by an increased number of local gene duplications and murine/human conservation mismatches (single genes that do not match into the conserved chromosomal segment). The CBR within CER1 contains a simple tandem TATAGA repeat capable of forming a 40-bp-long secondary hairpin-like structure. This repeat is nonrandomly localized within the other tumor-associated deletions and in the vicinity of 3p12-22 CBRs.

Tumor suppressor genes on chromosome 3p involved in the pathogenesis of lung and other cancers

Loss of heterozygosity (LOH) involving several chromosome 3p regions accompanied by chromosome 3p deletions are detected in almost 100% of small (SCLCs) and more than 90% of non-small (NSCLCs) cell lung cancers. In addition, these changes appear early in the pathogenesis of lung cancer and are found as clonal lesions in the smoking damaged respiratory epithelium including histologically normal epithelium as well as in epithelium showing histologic changes of preneoplasia. These 3p genetic alterations lead to the conclusion that the short arm of human chromosome 3 contains several tumor suppressor gene(s) (TSG(s)). Although the first data suggesting that 3p alterations were involved in lung carcinogenesis were published more than 10 years ago, only recently has significant progress been achieved in identifying the candidate TSGs and beginning to demonstrate their functional role in tumor pathogenesis. Some of the striking results of these findings has been the discovery of multiple 3p TSGs and the importance of tumor acquired promoter DNA methylation as an epigenetic mechanism for inactivating the expression of these genes in lung cancer. This progress, combined with the well known role of smoking as an environmental causative risk factor in lung cancer pathogenesis, is leading to the development of new diagnostic and therapeutic strategies which can be translated into the clinic to combat and prevent the lung cancer epidemic. It is clear now that genetic and epigenetic abnormalities of several genes residing in chromosome region 3p are important for the development of lung cancers but it is still obscure how many of them exist and which of the numerous candidate TSGs are the key players in lung cancer pathogenesis. We review herein our current knowledge and describe the most credible candidate genes.