SEARCH · Search PubMed

Results for “Cell annotation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47Linked to original sources

PubMed2002

Protein domain analysis in the era of complete genomes.

Domains present one of the most useful levels at which to understand protein function, and domain family-based analysis has had a profound impact on the study of individual proteins. Protein domain discovery has been progressing steadily over the past 30 years. What are the realistically achievable goals of sequence-based domain analysis, and how far off are they for the sequences encoded in eukaryotic genomes? Here we address some of the issues involved in better coverage of sequence-based domain annotation, and the integration of these results within the wider context of genomes, structures and function.

Amino Acid Sequence↗

PubMed2003

Reconstruction of the central carbon metabolism of Aspergillus niger.

The topology of central carbon metabolism of Aspergillus niger was identified and the metabolic network reconstructed, by integrating genomic, biochemical and physiological information available for this microorganism and other related fungi. The reconstructed network may serve as a valuable database for annotation of genes identified in future genome sequencing projects on aspergilli. Based on the metabolic reconstruction, a stoichiometric model was set up that includes 284 metabolites and 335 reactions, of which 268 represent biochemical conversions and 67 represent transport processes between the different intracellular compartments and between the cell and the extracellular medium. The stoichiometry of the metabolic reactions was used in combination with biosynthetic requirements for growth and pseudo-steady state mass balances over intracellular metabolites for the quantification of metabolic fluxes using metabolite balancing. This framework was employed to perform an in silico characterisation of the phenotypic behaviour of A. niger grown on different carbon sources. The effects on growth of single reaction deletions were assessed and essential biochemical reactions were identified for different carbon sources. Furthermore, application of the stoichiometric model for assessing the metabolic capabilities of A. niger to produce metabolites was evaluated by using succinate production as a case study.

Aspergillus niger↗

PubMed2001

Comprehensive analysis of a large genomic sequence at the putative B-cell chronic lymphocytic leukaemia (B-CLL) tumour suppresser gene locus.

In many haematological diseases, and more particularly in B-cell chronic lymphocytic leukaemia (B-CLL), the existence of a tumour suppressor gene located within the frequently deleted region 13q14.3, has been put forward. A wide candidate region spanning from marker D13S273 to D13S25 has been proposed and an extensive physical map has been constructed by several teams. In this study, we sequenced a minimal core deleted region that we have previously defined and annotated it with flanking available public sequences. Our analysis shows that this region is gene-poor. Furthermore, our work allowed us to identify new alternative transcripts, spanning core regions, of the previously defined candidate genes DLEU1 and DLEU2. Since their putative involvement in B-CLL was controversial, our present study provide support for reconsidering the DLEU1 and DLEU2 genes as B-CLL candidate genes, with a new definition of their organisation and context.

Alternative Splicing↗

PubMed2003

Structural proteome of human colostral fat globule membrane proteins.

Milk fat globule membrane (MFGM) contains proteins derived from the apical membrane of secreting epithelial cells of the mammary gland. Between 2-4% of total human milk protein content is associated with the fat globule fraction, as MFGM proteins. While MFGM proteins have very low classical nutritional value, they play important roles in various cell processes and defence mechanisms for the newborn. To date, fewer than 30 human MFGM proteins have been identified and characterized, either by immunological methods or by Edman sequencing and mass spectrometry. This study aimed to update the structural proteome of human colostral MFGM proteins and to create an annotated two-dimensional electrophoresis (2-DE) MFGM protein database available on-line. More than one hundred 2-DE spots derived from human colostral MFGM proteins were investigated by matrix-assisted laser desorption/ionization-time of flight mass spectrometry and proteins were identified by three different software packages available on the web (PeptIdent, MS-Fit and ProFound); uncertain identifications were solved by nanoelectrospray ionization-ion trap mass spectrometry using SEQUEST software.

Databases, Protein↗

PubMedDate not supplied

LNCIB human full-length cDNAs collection: towards a better comprehension of the human transcriptome.

LNCIB has been producing a variety of human full-length-enriched, normalized and subtracted cDNA libraries from various cell lines and tissues in different developmental stages by using the CAP-Trapper method. By sequencing 23000 clones of these libraries we identified a pool of about 5800 good quality unique cDNAs. After BLAST analysis on Human RefSeq/Unigene databases, 1717 of these sequences remained with no or poor annotation. We show that cross-species comparative BLAST resulted as a valid tool for the annotation of orthologous genes.

Databases, Factual↗

PubMed2000

Theileria parva genomics reveals an atypical apicomplexan genome.

The discipline of genomics is setting new paradigms in research approaches to resolving problems in human and animal health. We propose to determine the genome sequence of Theileria parva, a pathogen of cattle, using the random shotgun approach pioneered at The Institute for Genomic Research (TIGR). A number of features of the T. parva genome make it particularly suitable for this approach. The G+C content of genomic DNA is about 31%, non-coding repetitive DNA constitutes less than 1% of total DNA and a framework for the 10-12 Mbp genome is available in the form of a physical map for all four chromosomes. Minisatellite sequences are the only dispersed repetitive sequences identified so far, but they are limited in distribution to 13 of 33 SfiI fragments. Telomere and sub-telomeric non-coding sequences occupy less than 10 kbp at each chromosomal end and there are only two units encoding cytoplasmic rRNAs. Three sets of distinct multicopy sequences encoding ORFs have been identified but it is not known if these are associated with expression of parasite antigenic diversity. Protein coding genes exhibit a bias in codon usage and introns when present are unusually short. Like other apicomplexan organisms, T. parva contains two extrachromosomal DNAs, a mitochondrial DNA and a plastid DNA molecule. By annotating the genome sequence, in combination with the use of microarray technology and comparative genomics, we expect to gain significant insights into unique aspects of the biology of T. parva. We believe that the data will underpin future research to aid in the identification of targets of protective CD8+ cell mediated immune responses, and parasite molecules involved in inducing reversible host leukocyte transformation and tumour-like behaviour of transformed parasitised cells.

Animals↗

PubMed1999

Erythropoietin treatment for non-uremic patients: a personal view.

The correction of anemia in patients with chronic renal failure (CRF) has become the most important application of recombinant human erythropoietin (rHuEpo). The merits of rHuEpo therapy in patients with CRF are overt. Firstly, patients with CRF have an absolute deficiency in endogenous erythropoietin production and a relatively low maintenance dose of rHuEpo (often less than 100 IU/kg body weight per week) is effective in avoiding regular transfusions in the majority of the patients with CRF. Secondly, rHuEpo is able to avoid long-term complications of frequent transfusions (hemochromatosis, transfusion-transmissible diseases). Thirdly, patients with uremia notice a considerable improvement in quality of life (QOL) after initiation of rHuEpo. These advantages justify administration of this costly drug in CRF patients. The use of rHuEpo outside the setting of uremia do, however, not cover the complete spectrum of beneficial effects as compared to its use in (pre)dialysis patients. The aim of this overview is to provide some annotations on recently approved (cisplatin-induced anemia, preoperative anemia, zidovudine-related anemia) and possibly future (several types of malignancy and inflammation) indications for rHuEpo in non-uremic patients, leaving out the correction of anemia due to relatively uncommon disorders in the Dutch population (such as sickle cell anemia and thalassemia).

Anemia↗

PubMed2002

Mitogen-activated protein kinase cascades in plants: a new nomenclature.

Mitogen-activated protein kinase (MAPK) cascades are universal signal transduction modules in eukaryotes, including yeasts, animals and plants. These protein phosphorylation cascades link extracellular stimuli to a wide range of cellular responses. In plants, MAPK cascades are involved in responses to various biotic and abiotic stresses, hormones, cell division and developmental processes. Completion of the Arabidopsis genome-sequencing project has revealed the existence of 20 MAPKs, 10 MAPK kinases and 60 MAPK kinase kinases. Here, we propose a simplified nomenclature for Arabidopsis MAPKs and MAPK kinases that might also serve as a basis for standard annotation of these gene families in all plants.

Arabidopsis↗

PubMed2003

Hyper-gravity effects on the Arabidopsis transcriptome.

Callus cultures of Arabidopsis thaliana (cv. Columbia) in Petri dishes were exposed to altered g-forces by centrifugation (1-10 g). Using semi-quantitative RT-PCR transcripts of genes coding for metabolic key enzymes (ADP-glucose pyrophosphorylase, ADPG-PP; beta-amylase, fructose-1,6-bisphosphatase, FBPase; glyceraldehyde-P dehydrogenase, GAPDH; hydroxymethylglutaryl-CoA reductase, HMG; phenylalanine-ammonium-lyase, PAL; PEP carboxylase, PEPC) were used to monitor threshold conditions for g-number (all) and time of exposure (beta-amylase) which led to altered amounts of the gene product. Exposure to approximately 5 g and higher for 1 h resulted in altered transcript levels: transcripts of beta-amylase, PAL, and PEPC were increased, those of ADPG-PP decreased, while those of FBPase, GAPDH, and HMG were not affected. This probably indicates a shift from starch synthesis to starch degradation and increased rates of anaplerosis (PEPC: supply of ketoacids for amino acid synthesis). In order to get more information about g-related effects on gene expression, we used a 1-h exposure to 7 g for a microarray analysis, using a commercial A. thaliana chip with 4105 unique annotated clusters/genes (IncyteGenomics). Transcripts of more than 200 genes were significantly increased in amount (ratio 7 g/1 g control; 2(1.6) and larger). They fall into several categories. Transcripts coding for enzymes of major pathways form the largest group (25%), followed by gene products involved in cellular organization and cell wall formation/rearrangement (17%), signalling, phosphorylation/dephosphorylation (12%), proteolysis and transport (10% each), hormone synthesis plus related events (8%), defense (4%), stress-response (2%), and gravi-sensing (2%). Many of the alterations are part of a general stress response, but some changes related to the synthesis/rearrangement of cell wall components could be more hyper-g-specific. We only found few gene products, which were decreased in relation to 1 g controls, and these were less significant (ratio < 2(1.6)). We thus assume that g-forces above a threshold of about 5 g for 1 h are sensed by plant cells in general, causing distinct metabolic responses, which obviously in part, are regulated by gene expression.

Arabidopsis↗

PubMed2005

Human histopathology of electroanatomic mapping after cooled-tip radiofrequency ablation to treat ventricular tachycardia in remote myocardial infarction.

INTRODUCTION: Catheter ablation of ventricular tachycardia (VT) in remote myocardial infarction (MI) often requires excessive mapping procedures. Documentation of the electrical substrate via electrogram amplitude may help to identify regions of altered myocardium resembling exit areas of reentrant VTs. METHODS AND RESULTS: A patient with multiple symptomatic monomorphic VTs (biventricular ICD, remote MI) underwent electroanatomic substrate mapping (CARTOtrade mark) for VT ablation. Regions of scar (bipolar electrogram amplitudes or=1.5 mV), and "altered" myocardium (0.5-1.5 mV) were identified. Ablation was directed to regions with "altered" myocardium based on pace map correlation. After ablation the clinical VT did not reoccur. The patient died due to worsening of heart failure 7 days afterward. During postmortal evaluation specified sites of electroanatomic mapping were correlated to histopathological findings. Annotated scar areas were documented to consist of areas with massive fibrosis (>or=80% of mural composition). Ablations were found to span through regions with intermediate fibrosis (21-79%) mapped as "altered" myocardium. Ablation produced transmural coagulation necrosis of mesh-like fibrotic tissue with interspersed remnants of myocardial cells up to a maximum depth of 7.0 mm. Subendocardial intramural bleedings were universal findings 7 days after ablation. CONCLUSIONS: Electroanatomic substrate mapping for VT ablation sufficiently identified regions of scar and normal myocardium. Regions with bipolar electrogram amplitudes between 0.5 and 1.5 mV were found to correlate to areas of "intermediate" fibrosis (21-79%) with only remnant strands of myocardial cells and were identified as target region for ablation. Cooled-tip endocardial radiofrequency ablation lead to transmural coagulation necrosis up to a depth of 7.0 mm.

Aged↗

PubMed2003

Detection of a luxS-signaling molecule in Bacillus anthracis.

Quorum-sensing regulation of density-dependent genes has been described for numerous bacterial species. The partially annotated genome sequence of Bacillus anthracis contains an open reading frame (BA5047) predicted to encode an ortholog of luxS, required for synthesis of the quorum-sensing signaling molecule autoinducer-2 (AI-2). To determine whether B. anthracis produces AI-2, the Vibrio harveyi luminescence bioassay was used. Cell-free conditioned media from vaccine (Sterne) strain 34F(2) induced luminescence in V. harveyi reporter strain BB170, indicating its production of AI-2. Cloned BA5047, expressed in Escherichia coli DH5 alpha cells, restored AI-2 activity to these cells. To evaluate whether BA5047 is essential for AI-2 synthesis, it was deleted through allelic exchange with marker rescue; the resulting mutant had no functional luxS activity and had reduced growth in vitro. In the wild-type strain, AI-2 activity was greatest during the exponential phase of growth. In total, these data indicate that BA5047 is a functional luxS ortholog in B. anthracis necessary for growth-phase-specific AI-2 expression. Thus, B. anthracis may utilize extracellular signaling molecules to regulate density-dependent gene expression.

Bacillus anthracis↗

PubMed2005

Comparison of the oxidative phosphorylation (OXPHOS) nuclear genes in the genomes of Drosophila melanogaster, Drosophila pseudoobscura and Anopheles gambiae.

BACKGROUND: In eukaryotic cells, oxidative phosphorylation (OXPHOS) uses the products of both nuclear and mitochondrial genes to generate cellular ATP. Interspecies comparative analysis of these genes, which appear to be under strong functional constraints, may shed light on the evolutionary mechanisms that act on a set of genes correlated by function and subcellular localization of their products. RESULTS: We have identified and annotated the Drosophila melanogaster, D. pseudoobscura and Anopheles gambiae orthologs of 78 nuclear genes encoding mitochondrial proteins involved in oxidative phosphorylation by a comparative analysis of their genomic sequences and organization. We have also identified 47 genes in these three dipteran species each of which shares significant sequence homology with one of the above-mentioned OXPHOS orthologs, and which are likely to have originated by duplication during evolution. Gene structure and intron length are essentially conserved in the three species, although gain or loss of introns is common in A. gambiae. In most tissues of D. melanogaster and A. gambiae the expression level of the duplicate gene is much lower than that of the original gene, and in D. melanogaster at least, its expression is almost always strongly testis-biased, in contrast to the soma-biased expression of the parent gene. CONCLUSIONS: Quickly achieving an expression pattern different from the parent genes may be required for new OXPHOS gene duplicates to be maintained in the genome. This may be a general evolutionary mechanism for originating phenotypic changes that could lead to species differentiation.

Animals↗

PubMed2004

Genomic annotation of the meningioma tumor suppressor locus on chromosome 1p34.

Meningioma is a frequently occurring tumor of the meninges surrounding the central nervous system. Loss of the short arm of chromosome 1 (1p) is the second most frequent chromosomal abnormality observed in these tumors. Previously, we identified a 3.7 megabase (Mb) region of consistent deletion on 1p33-p34 in a panel of 157 tumors. Loss of this region was associated with advanced disease and predictive for tumor relapse. In this report, a high-resolution integrated map of the region was constructed (CompView) to identify all markers in the smallest region of overlapping deletion (SRO). A regional somatic cell hybrid panel was used to more precisely localize those markers identified in CompView as within or overlapping the region. Additional deletion mapping using microsatellites localized to the region narrowed the SRO to approximately 2.8 Mb. The 88 markers remaining in the SRO were used to screen genomic databases to identify large-insert clones. Clones were assembled into a physical map of the region by PCR-based, sequence-tagged site (STS) content mapping. A sequence from clones was used to validate STS content by electronic PCR and to identify transcripts. A minimal tiling path of 43 clones was constructed across the SRO. Sequence data from the most current sequence assembly were used for further validation. A total of 59 genes were ordered within the SRO. In all, 17 of these were selected as likely candidates based on annotation using Gene Ontology Consortium terms, including the MUTYH, PRDX1, FOXD2, FOXE3, PTCH2, and RAD54L genes. This annotation of a putative tumor suppressor locus provides a resource for further analysis of meningioma candidate genes.

Chromosomes, Human, Pair 1↗

PubMed2004

Proteome-wide functional classification and identification of prokaryotic transmembrane proteins by transmembrane topology similarity comparison.

We propose a new method for classifying and identifying transmembrane (TM) protein functions in proteome-scale by applying a single-linkage clustering method based on TM topology similarity, which is calculated simply from comparing the lengths of loop regions. In this study, we focused on 87 prokaryotic TM proteomes consisting of 31 proteobacteria, 22 gram-positive bacteria, 19 other bacteria, and 15 archaea. Prior to performing the clustering, we first categorized individual TM protein sequences as "known," "putative" (similar to "known" sequences), or "unknown" by using the homology search and the sequence similarity comparison against SWISS-PROT to assess the current status of the functional annotation of the TM proteomes based on sequence similarity only. More than three-quarters, that is, 75.7% of the TM protein sequences are functionally "unknown," with only 3.8% and 20.5% of them being classified as "known" and "putative," respectively. Using our clustering approach based on TM topology similarity, we succeeded in increasing the rate of TM protein sequences functionally classified and identified from 24.3% to 60.9%. Obtained clusters correspond well to functional superfamilies or families, and the functional classification and identification are successfully achieved by this approach. For example, in an obtained cluster of TM proteins with six TM segments, 109 sequences out of 119 sequences annotated as "ATP-binding cassette transporter" are properly included and 122 "unknown" sequences are also contained.

Algorithms↗

PubMed2005

Generation, annotation, analysis and database integration of 16,500 white spruce EST clusters.

BACKGROUND: The sequencing and analysis of ESTs is for now the only practical approach for large-scale gene discovery and annotation in conifers because their very large genomes are unlikely to be sequenced in the near future. Our objective was to produce extensive collections of ESTs and cDNA clones to support manufacture of cDNA microarrays and gene discovery in white spruce (Picea glauca [Moench] Voss). RESULTS: We produced 16 cDNA libraries from different tissues and a variety of treatments, and partially sequenced 50,000 cDNA clones. High quality 3' and 5' reads were assembled into 16,578 consensus sequences, 45% of which represented full length inserts. Consensus sequences derived from 5' and 3' reads of the same cDNA clone were linked to define 14,471 transcripts. A large proportion (84%) of the spruce sequences matched a pine sequence, but only 68% of the spruce transcripts had homologs in Arabidopsis or rice. Nearly all the sequences that matched the Populus trichocarpa genome (the only sequenced tree genome) also matched rice or Arabidopsis genomes. We used several sequence similarity search approaches for assignment of putative functions, including blast searches against general and specialized databases (transcription factors, cell wall related proteins), Gene Ontology term assignation and Hidden Markov Model searches against PFAM protein families and domains. In total, 70% of the spruce transcripts displayed matches to proteins of known or unknown function in the Uniref100 database (blastx e-value < 1e-10). We identified multigenic families that appeared larger in spruce than in the Arabidopsis or rice genomes. Detailed analysis of translationally controlled tumour proteins and S-adenosylmethionine synthetase families confirmed a twofold size difference. Sequences and annotations were organized in a dedicated database, SpruceDB. Several search tools were developed to mine the data either based on their occurrence in the cDNA libraries or on functional annotations. CONCLUSION: This report illustrates specific approaches for large-scale gene discovery and annotation in an organism that is very distantly related to any of the fully sequenced genomes. The ArboreaSet sequences and cDNA clones represent a valuable resource for investigations ranging from plant comparative genomics to applied conifer genetics.

Arabidopsis↗

PubMed2006

Analysis of 10,000 ESTs from lymphocytes of the cynomolgus monkey to improve our understanding of its immune system.

BACKGROUND: The cynomolgus monkey (Macaca fascicularis) is one of the most widely used surrogate animal models for an increasing number of human diseases and vaccines, especially immune-system-related ones. Towards a better understanding of the gene expression background upon its immunogenetics, we constructed a cDNA library from Epstein-Barr virus (EBV)-transformed B lymphocytes of a cynomolgus monkey and sequenced 10,000 randomly picked clones. RESULTS: After processing, 8,312 high-quality expressed sequence tags (ESTs) were generated and assembled into 3,728 unigenes. Annotations of these uniquely expressed transcripts demonstrated that out of the 2,524 open reading frame (ORF) positive unigenes (mitochondrial and ribosomal sequences were not included), 98.8% shared significant similarities (E-value less than 1e-10) with the NCBI nucleotide (nt) database, while only 67.7% (E-value less than 1e-5) did so with the NCBI non-redundant protein (nr) database. Further analysis revealed that 90.0% of the unigenes that shared no similarities to the nr database could be assigned to human chromosomes, in which 75 did not match significantly to any cynomolgus monkey and human ESTs. The mapping regions to known human genes on the human genome were described in detail. The protein family and domain analysis revealed that the first, second and fourth of the most abundantly expressed protein families were all assigned to immunoglobulin and major histocompatibility complex (MHC)-related proteins. The expression profiles of these genes were compared with that of homologous genes in human blood, lymph nodes and a RAMOS cell line, which demonstrated expression changes after transformation with EBV. The degree of sequence similarity of the MHC class I and II genes to the human reference sequences was evaluated. The results indicated that class I molecules showed weak amino acid identities (<90%), while class II showed slightly higher ones. CONCLUSION: These results indicated that the genes expressed in the cynomolgus monkey could be used to identify novel protein-coding genes and revise those incomplete or incorrect annotations in the human genome by comparative methods, since the old world monkeys and humans share high similarities at the molecular level, especially within coding regions. The identification of multiple genes involved in the immune response, their sequence variations to the human homologues, and their responses to EBV infection could provide useful information to improve our understanding of the cynomolgus monkey immune system.

5' Untranslated Regions↗

PubMed2026

Bacterial polar metabolites modulate β-amyloid toxicity and cholinergic dysfunction in models of Alzheimer's disease.

Alzheimer's disease is characterized by progressive neurodegeneration driven by β-amyloid (Aβ) toxicity, oxidative stress, and cholinergic dysfunction. In this study, we investigated whether polar metabolites derived from a cultivable bacterial isolate could modulate Aβ-associated neurodegenerative phenotypes in complementary experimental models. A bioactivity-guided approach identified an aqueous fraction with high antioxidant capacity in DPPH, FRAP, and ORAC assays. In a transgenic Drosophila melanogaster model expressing human Aβ, treatment with this fraction significantly reduced amyloid accumulation and attenuated neurodegenerative histopathological alterations. In human SH-SY5Y neuronal cultures, the metabolites improved cell viability under therapeutic, but not preventive, conditions following exposure to aggregated Aβ. The aqueous fraction also exhibited significant inhibitory activity against acetylcholinesterase and butyrylcholinesterase. Whole-genome sequencing assigned the bioactive isolate to the genus Providencia, with comparative genomic analyses suggesting its placement within a distinct taxonomic lineage. Metabolomic profiling by LC-ESI-MS/MS revealed a diverse set of polar metabolites, including metabolites putatively annotated based on spectral matching, previously associated with neuroprotective and cholinesterase-modulating activities. Collectively, these findings demonstrate that bacterial polar metabolites can modulate key pathological features of Alzheimer's disease, supporting their relevance for mechanistic studies of Aβ toxicity and cholinergic dysfunction.

Animals↗

PubMed2007

Immunoaffinity purification of plasma membrane with secondary antibody superparamagnetic beads for proteomic analysis.

Plasma membrane (PM) has very important roles in cell-cell interaction and signal transduction, and it has been extensively targeted for drug design. A major prerequisite for the analysis of PM proteome is the preparation of PM with high purity. Density gradient centrifugation has been commonly employed to isolate PM, but it often occurred with contamination of internal membrane. Here we describe a method for plasma membrane purification using second antibody superparamagnetic beads that combines subcellular fractionation and immunoisolation strategies. Four methods of immunoaffinity were compared, and the variation of crude plasma membrane (CPM), superparamagnetic beads, and antibodies was studied. The optimized method and the number of CPM, beads, and antibodies suitable for proteome analysis were obtained. The PM of mouse liver was enriched 3-fold in comparison with the density gradient centrifugation method, and contamination from mitochondria was reduced 2-fold. The PM protein bands were extracted and trypsin-digested, and the resulting peptides were resolved and characterized by MALDI-TOF-TOF and ESI-Q-TOF, respectively. Mascot software was used to analyze the data against IPI-mouse protein database. Nonredundant proteins (248) were identified, of which 67% are PM or PM-related proteins. No endoplasmic reticulum (ER) or nuclear proteins were identified according to the GO annotation in the optimized method. Our protocol represents a simple, economic, and reproducible tool for the proteomic characterization of liver plasma membrane.

Animals↗