Cysteine and tyrosine-rich 1 (CYYR1), a novel unpredicted gene on human chromosome 21 (21q21.2), encodes a cysteine and tyrosine-rich protein and defines a new family of highly conserved vertebrate-specific genes

CYYR1 promotes the degradation of the E3 ubiquitin ligase WWP1 and is associated with favorable prognosis in breast cancer

Ubiquitination plays a crucial role in cellular homeostasis by regulating the degradation, localization, and activity of proteins, ensuring proper cell function and balance. Among E3 ubiquitin ligases, WW domain-containing protein 1 (WWP1) is implicated in cell proliferation, survival, and apoptosis. Notably WWP1 is frequently amplified in breast cancer and associated with poor prognosis. Here, we identify the protein cysteine and tyrosine-rich protein 1 (CYYR1) that had previously no assigned function, as a regulator of WWP1 activity and stability. We show that CYYR1 binds to the WW domains of the E3 ubiquitin ligase WWP1 through its PPxY motifs. This interaction triggers K63-linked autoubiquitination and subsequent degradation of WWP1. We furthermore demonstrate that CYYR1 localizes to late endosomal vesicles and directs polyubiquitinated WWP1 toward lysosomal degradation through binding to ANKyrin repeat domain-containing protein 13 A (ANKRD13A). Moreover, we found that CYYR1 expression attenuates breast cancer cell growth in anchorage-dependent and independent colony formation assays in a PPxY-dependent manner. Finally, we highlight that CYYR1 expression is significantly decreased in breast cancer and is associated with beneficial clinical outcome. Taken together our study suggests tumor suppressor properties for CYYR1 through regulation of WWP1 autoubiquitination and lysosomal degradation.

Difficulty in obtaining the complete mRNA coding sequence at 5' region (5' end mRNA artifact): Causes, consequences in biology and medicine and possible solutions for obtaining the actual amino acid sequence of proteins (Review)

The known difficulty in obtaining the actual full length, complete sequence of a messenger RNA (mRNA) may lead to the erroneous determination of its coding sequence at the 5' region (5' end mRNA artifact), and consequently to the wrong assignment of the translation start codon, leading to the inaccurate prediction of the encoded polypeptide at its amino terminus. Among the known human genes whose study was affected by this artifact, we can include disco interacting protein 2 homolog A (DIP2A; KIAA0184), Down syndrome critical region 1 (DSCR1), SON DNA binding protein (SON), trefoil factor 3 (TFF3) and URB1 ribosome biogenesis 1 homolog (URB1; KIAA0539) on chromosome 21, as well as receptor for activated C kinase 1 (RACK1, also known as GNB2L1), glutaminyl‑tRNA synthetase (QARS) and tyrosyl-DNA phosphodiesterase 2 (TDP2) along with another 474 loci, including interleukin 16 (IL16). In this review, we discuss the causes of this issue, its quantitative incidence in biomedical research, the consequences in biology and medicine, and the possible solutions for obtaining the actual amino acid sequence of proteins in the post-genomics era.

Characterization of human gene locus CYYR1: a complex multi-transcript system

Cysteine/tyrosine-rich 1 (CYYR1) is a gene we previously identified on human chromosome 21 starting from an in-depth bioinformatics analysis of chromosome 21 segment 40/105 (21q21.3), where no coding region had previously been predicted. CYYR1 was initially characterized as a four-exon gene, whose brain-derived cDNA sequencing predicts a 154-amino acid product. In this study we provide, with in silico and in vitro analyses, the first detailed description of the human CYYR1 locus. The analysis of this locus revealed that it is composed of a multi-transcript system, which includes at least seven CYYR1 alternative spliced isoforms and a new CYYR1 antisense gene (named CYYR1-AS1). In particular, we cloned, for the first time, the following isoforms: CYYR1-1,2,3,4b and CYYR1-1,2,3b, which present a different 3' transcribed region, with a consequent different carboxy-terminus of the predicted proteins; CYYR1-1,2,4 lacks exon 3; CYYR1-1,2,2bis,3,4 presents an additional exon between exon 2 and exon 3; CYYR1-1b,2,3,4 presents a different 5' untranslated region when compared to CYYR1. The complexity of the locus is enriched by the presence of an antisense transcript. We have cloned a long transcript overlapping with CYYR1 as an antisense RNA, probably a non-coding RNA. Expression analysis performed in different normal tissues, tumour cell lines as well as in trisomy 21 and euploid fibroblasts has confirmed a quantitative and qualitative variability in the expression pattern of the multi-transcript locus, suggesting a possible role in complex diseases that should be further investigated.

High throughput gene expression analysis identifies reliable expression markers of human corneal endothelial cells

Considerable interest has been generated for the development of suitable corneal endothelial graft alternatives through cell-tissue engineering, which can potentially alleviate the shortage of corneal transplant material. The advent of less invasive suture-less key-hole surgery options such as Descemet’s Stripping Endothelial Keratoplasty (DSEK) and Descemet’s Membrane Endothelial Keratoplasty (DMEK), which involve transplantation of solely the endothelial layer instead of full thickness cornea, provide further impetus for the development of alternative endothelial grafts for clinical applications. A major challenge for this endeavor is the lack of specific markers for this cell type. To identify genes that reliably mark corneal endothelial cells (CECs) in vivo and in vitro, we performed RNA-sequencing on freshly isolated human CECs (from both young and old donors), CEC cultures, and corneal stroma. Gene expression of these corneal cell types was also compared to that of other human tissue types. Based on high throughput comparative gene expression analysis, we identified a panel of markers that are: i) highly expressed in CECs from both young donors and old donors; ii) expressed in CECs in vivo and in vitro; and iii) not expressed in corneal stroma keratocytes and the activated corneal stroma fibroblasts. These were SLC4A11, COL8A2 and CYYR1. The use of this panel of genes in combination reliably ascertains the identity of the CEC cell type.

Unexpected diversity in Shisa-like proteins suggests the importance of their roles as transmembrane adaptors

The Shisa family of single-transmembrane proteins is characterized by an N-terminal cysteine-rich domain and a proline-rich C-terminal region. Its founding member, Xenopus Shisa, promotes head development by antagonizing Wnt and FGF signaling. Recently, a mouse brain-specific Shisa protein CKAMP44 (Shisa9) was shown to play an important role in AMPA receptor desensitization. We used sequence similarity searches against protein, genome and EST databases to study the evolutionary origin and phylogenetic distribution of Shisa homologs. In addition to nine Shisa subfamilies in vertebrates, we detected distantly related Shisa homologs that possess an N-terminal domain with six conserved cysteines. These Shisa-like proteins include FAM159 and KIAA1644 mainly from vertebrates, and members from various bilaterian invertebrates and Porifera, suggesting their presence in the last common ancestor of Metazoa. Shisa-like genes have undergone large expansions in Branchiostoma floridae and Saccoglossus kowalevskii, and appear to have been lost in certain insects. Pattern-based searches against eukaryotic proteomes also uncovered several other families of predicted single-transmembrane proteins with a similar cysteine-rich domain. We refer to these proteins (Shisa/Shisa-like, WBP1/VOPP1, CX, DUF2650, TMEM92, and CYYR1) as STMC6 proteins (single-transmembrane proteins with conserved 6 cysteines). STMC6 genes are widespread in Metazoa, with the human genome containing 17 members. Frequent occurrences of PY motifs in STMC6 proteins suggest that most of them could interact with WW-domain-containing proteins, such as the NEDD4 family E3 ubiquitin ligases, and could play critical roles in protein degradation and sorting. STMC6 proteins are likely transmembrane adaptors that regulate membrane proteins such as cell surface receptors.

Mapping transcription mechanisms from multimodal genomic data

BACKGROUND

Identification of expression quantitative trait loci (eQTLs) is an emerging area in genomic study. The task requires an integrated analysis of genome-wide single nucleotide polymorphism (SNP) data and gene expression data, raising a new computational challenge due to the tremendous size of data.

RESULTS

We develop a method to identify eQTLs. The method represents eQTLs as information flux between genetic variants and transcripts. We use information theory to simultaneously interrogate SNP and gene expression data, resulting in a Transcriptional Information Map (TIM) which captures the network of transcriptional information that links genetic variations, gene expression and regulatory mechanisms. These maps are able to identify both cis- and trans- regulating eQTLs. The application on a dataset of leukemia patients identifies eQTLs in the regions of the GART, PCP4, DSCAM, and RIPK4 genes that regulate ADAMTS1, a known leukemia correlate.

CONCLUSIONS

The information theory approach presented in this paper is able to infer the dependence networks between SNPs and transcripts, which in turn can identify cis- and trans-eQTLs. The application of our method to the leukemia study explains how genetic variants and gene expression are linked to leukemia.

Identification of human gene products containing Pro-Pro-x-Tyr (PY) motifs that enhance glutathione and endocytotic marker uptake in yeast

In an attempt to identify genes involved in glutathione (GSH) transport, a human mammary gland cDNA library was screened for clones capable of complementing a defect in GSH uptake in yeast cells that lack Hgt1p, the primary yeast GSH uptake transporter. Five genes capable of rescuing growth on sulfur-deficient GSH-containing medium were identified: prostate transmembrane protein, androgen induced 1 (PMEPA1); lysosomal-associated protein transmembrane 4 alpha (LAPTM4alpha); solute carrier family 25, member 1 (SLC25A1); lipopolysaccharide-induced TNF factor (LITAF); and cysteine/tyrosine-rich-1 (CYYR1). All of these genes encode small integral membrane proteins of unknown function, although none appear to encode prototypical GSH transporters. Nevertheless, they all increased both intracellular glutathione levels and [(3)H]GSH uptake rates. [(3)H]GSH uptake was uniformly inhibited by high concentrations of unlabeled GSH, GSSG, and ophthalmic acid. Interestingly, each protein is predicted to contain Pro-Pro-x-Tyr (PY) motifs, which are thought to be important for regulating protein cell surface expression. Uptake of the endocytotic markers lucifer yellow and FM4-64 was also enhanced by each of the five genes. Mutations of the PY motifs in LITAF largely abolished all of its effects. In summary, although the results do not reveal novel GSH transporters, they identify five PY-containing human gene products that may influence plasma membrane transport activity.

Detection of growth hormone doping by gene expression profiling of peripheral blood

CONTEXT

GH abuse is a significant problem in many sports, and there is currently no robust test that allows detection of doping beyond a short window after administration.

OBJECTIVE

Our objective was to evaluate gene expression profiling in peripheral blood leukocytes in-vivo as a test for GH doping in humans.

DESIGN

Seven men and thirteen women were administered GH, 2 mg/d sc for 8 wk. Blood was collected at baseline and at 8 wk. RNA was extracted from the white cell fraction. Microarray analysis was undertaken using Agilent 44K G4112F arrays using a two-color design. Quantitative RT-PCR using TaqMan gene expression assays was performed for validation of selected differentially expressed genes.

RESULTS

GH induced an approximately 2-fold increase in circulating IGF-I that was maintained throughout the 8 wk of the study. GH induced significant changes in gene expression with 353 in women and 41 in men detected with a false discovery rate of less than 5%. None of the differentially expressed genes were common between men and women. The maximal changes were a doubling for up-regulated or halving for down-regulated genes, similar in magnitude to the variation between individuals. Quantitative RT-PCR for seven target genes showed good concordance between microarray and quantitative PCR data in women but not in men.

CONCLUSION

Gene expression analysis of peripheral blood leukocytes is unlikely to be a viable approach for the detection of GH doping.

Sequence, "subtle" alternative splicing and expression of the CYYR1 (cysteine/tyrosine-rich 1) mRNA in human neuroendocrine tumors

Background

CYYR1 is a recently identified gene located on human chromosome 21 whose product has no similarity to any known protein and is of unknown function. Analysis of expressed sequence tags (ESTs) have revealed high human CYYR1 expression in cells belonging to the diffuse neuroendocrine system (DNES). These cells may be the origin of neuroendocrine (NE) tumors. The aim of this study was to conduct an initial analysis of sequence, splicing and expression of the CYYR1 mRNA in human NE tumors.

Methods

The CYYR1 mRNA coding sequence (CDS) was studied in 32 NE tumors by RT-PCR and sequence analysis. A subtle alternative splicing was identified generating two isoforms of CYYR1 mRNA differing in terms of the absence (CAG^- isoform, the first described mRNA for CYYR1 locus) or the presence (CAG⁺ isoform) of a CAG codon. When present, this specific codon determines the presence of an alanine residue, at the exon 3/exon 4 junction of the CYYR1 mRNA. The two mRNA isoform amounts were determined by quantitative relative RT-PCR in 29 NE tumors, 2 non-neuroendocrine tumors and 10 normal tissues. A bioinformatic analysis was performed to search for the existence of the two CYYR1 isoforms in other species.

Results

The CYYR1 CDS did not show differences compared to the reference sequence in any of the samples, with the exception of an NE tumor arising in the neck region. Sequence analysis of this tumor identified a change in the CDS 333 position (T instead of C), leading to the amino acid mutation P111S. NE tumor samples showed no significant difference in either CYYR1 CAG^- or CAG⁺ isoform expression compared to control tissues. CYYR1 CAG^- isoform was significantly more expressed than CAG⁺ isoform in NE tumors as well as in control samples investigated. Bioinformatic analysis revealed that only the genomic sequence of Pan troglodytes CYYR1 is consistent with the possible existence of the two described mRNA isoforms.

Conclusion

A new "subtle" splicing isoform (CAG⁺) of CYYR1 mRNA, the sequence and the expression of this gene were defined in a large series of NE tumors.

mRNA 5′ region sequence incompleteness: a potential source of systematic errors in translation initiation codon assignment in human mRNAs

The amino acid sequence of gene products is routinely deduced from the nucleotide sequence of the relative cloned cDNA, according to the rules for recognition of start codon (first-AUG rule, optimal sequence context) and the genetic code. From this prediction stem most subsequent types of product analysis, although all standard methods for cDNA cloning are affected by a potential inability to effectively clone the 5' region of mRNA. Revision by bioinformatics and cloning methods of 109 known genes located on human chromosome 21 (HC 21) shows that 60 mRNAs lack any in-frame stop upstream of the first-AUG, and that in five cases (DSCR1, KIAA0184, KIAA0539, SON, and TFF3) the coding region at the 5' end was incompletely characterized in the original descriptions. We describe the respective consequences for genomic annotation, domain and ortholog identification, and functional experiments design. We have also analyzed the sequences of 13,124 human mRNAs (RefSeq databank), discovering that in 6448 cases (49%), an in-frame stop codon is present upstream of the initiation codon, while in the other 6676 mRNAs (51%), identification of additional bases at the mRNA 5' region could well reveal some new upstream in-frame AUG codons in the optimal context. Proportionally to the HC 21 data, about 550 known human genes might thus be affected by this 5' end mRNA artifact.