Transcriptome analysis of cynomolgus macaques throughout their lifespan reveals age-related immune patterns.

Despite the different perspectives by diverse research sectors spanning several decades, aging research remains uncharted territory for human beings. Therefore, we investigated the transcriptomic characteristics of eight male healthy cynomolgus macaques, and the annual sampling was designed with two individuals in four age groups. As a laboratory animal, the macaques were meticulously shielded from all environmental factors except aging. The results showed recent findings of certain immune response and the age-associated network of primate immunity. Three important aging patterns were identified and each gene clusters represented a different immune response. The increased expression pattern was predominantly associated with innate immune cells, such as Neutrophils and NK cells, causing chronic inflammation with aging whereas the other two decreased patterns were associated with adaptive immunity, especially "B cell activation" affecting antibody diversity of aging. Furthermore, the hub gene network of the patterns reflected transcriptomic age and correlated with human illness status, aiding in future human disease prediction. Our macaque transcriptome profiling results offer systematic insights into the age-related immunological features of primates.

Sense-oriented AluYRa1 elements provide a lineage-specific transcription environment for polyadenylation.

Transposable elements cause alternative splicing (AS) in different ways, contributing to transcript diversification. Alternative polyadenylation (APA), one of the AS events, is related to the generation of mRNA isoforms in 70% of human genes. In this study, we tried to investigate AluYRa1s located at the terminal region of cynomolgus monkey genes, utilizing both computational analysis and molecular experimentation. We found that ten genes had AluYRa1 at their 3' end, and nine of these AluYRa1s were sense-oriented. Furthermore, in seven genes, AluYRa1s were expected to have a similar consensus sequence for polyadenylation cleavage. Additional computational analysis using the annotation files from the UCSC database showed that AluYRa1 was more involved in polyadenylation than in open reading frame exon splicing. To examine the extent of AluYRa1 involvement in polyadenylation, RNA-seq data from 30 normal cynomolgus monkeys were analyzed using TAPAS, a recently devised software that detects all the promising polyadenylation sites including APA sites. We observed that approximately 74% of possible polyadenylation sites in the analyzed genes were provided by sense-oriented AluYRa1. In conclusion, AluYRa1 is an Old-World monkey-specific TE, and its sense-oriented insertion at the 3'UTR region tends to provide a favorable environment for polyadenylation, diversifying gene transcripts.

Particulate Matter Exposure During Oocyte Maturation: Cell Cycle Arrest, ROS Generation, and Early Apoptosis in Mice.

Particulate matter (PM) is a general atmospheric pollutant released into the air by an anthropogenic and naturally derived mixture of substances. Current studies indicate that fine dust can result in different health defects, including endothelial dysfunction, asthma, lung cancer, cardiovascular diseases, uterine leiomyoma, deterioration in sperm quality, and overall birth impairment. However, the most prominent effects of PM10 (diameter < 10 µM) exposure on the female reproductive system, especially with respect to oocyte maturation, remain unclear. In the present study, maturing mouse oocytes were treated with PM10 and the phenotypes of the resulting toxic effects were investigated. Exposure to PM10 led to impairment of maturation capacity by inducing cell cycle arrest and blocking normal polar body extrusion during in vitro maturation and activation of fertilization of mouse oocytes. Additionally, defects in tubulin formation and DNA alignment were observed in PM10-treated oocytes during metaphase I to anaphase/telophase I transition. Moreover, PM10 induced reactive oxygen species generation, mitochondrial dysfunction, DNA damage, and early apoptosis. Taken together, these results indicate that PM10 exposure leads to a decline in oocyte quality and affects the subsequent embryonic development potential of mammalian oocytes.

Longitudinal profiling of the blood transcriptome in an African green monkey aging model.

African green monkeys (AGMs, Chlorocebus aethiops) are Old World monkeys which are used as experimental models in biomedical research. Recent technological advances in next generation sequencing are useful for unraveling the genetic mechanisms underlying senescence, aging, and age-related disease. To elucidate the normal aging mechanisms in older age, the blood transcriptomes of nine healthy, aged AGMs (15‒23 years old), were analyzed over two years. We identified 910‒1399 accumulated differentially expressed genes (DEGs) in each individual, which increased with age. Aging-related DEGs were sorted across the three time points. A major proportion of the aging-related DEGs belonged to gene ontology (GO) categories involved in translation and rRNA metabolic processes. Next, we sorted common aging-related DEGs across three time points over two years. Common aging-related DEGs belonged to GO categories involved in translation, cellular component biogenesis, rRNA metabolic processes, cellular component organization, biogenesis, and RNA metabolic processes. Furthermore, we identified 29 candidate aging genes that were upregulated across the time series analysis. These candidate aging genes were linked to protein synthesis. This study describes a changing gene expression pattern in AGMs during aging using longitudinal transcriptome sequencing. The candidate aging genes identified here may be potential targets for the treatment of aging.

A single mutation in the ACTR8 gene associated with lineage-specific expression in primates.

BACKGROUND: Alternative splicing (AS) generates various transcripts from a single gene and thus plays a significant role in transcriptomic diversity and proteomic complexity. Alu elements are primate-specific transposable elements (TEs) and can provide a donor or acceptor site for AS. In a study on TE-mediated AS, we recently identified a novel AluSz6-exonized ACTR8 transcript of the crab-eating monkey (Macaca fascicularis). In the present study, we sought to determine the molecular mechanism of AluSz6 exonization of the ACTR8 gene and investigate its evolutionary and functional consequences in the crab-eating monkey. RESULTS: We performed RT-PCR and genomic PCR to analyze AluSz6 exonization in the ACTR8 gene and the expression of the AluSz6-exonized transcript in nine primate samples, including prosimians, New world monkeys, Old world monkeys, and hominoids. AluSz6 integration was estimated to have occurred before the divergence of simians and prosimians. The Alu-exonized transcript obtained by AS was lineage-specific and expressed only in Old world monkeys and apes, and humans. This lineage-specific expression was caused by a single G duplication in AluSz6, which provides a new canonical 5' splicing site. We further identified other alternative transcripts that were unaffected by the AluSz6 insertion. Finally, we observed that the alternative transcripts were transcribed into new isoforms with C-terminus deletion, and in silico analysis showed that these isoforms do not have a destructive function. CONCLUSIONS: The single G duplication in the TE sequence is the source of TE exonization and AS, and this mutation may suffer a different fate of ACTR8 gene expression during primate evolution.

Cooperative evolution of two different TEs results in lineage-specific novel transcripts in the BLOC1S2 gene.

BACKGROUND: The BLOC1S2 gene encodes the multifunctional protein BLOS2, a shared subunit of two lysosomal trafficking complexes: i) biogenesis of lysosome-related organelles complex-1 and i) BLOC-1-related complex. In our previous study, we identified an intriguing unreported transcript of the BLOC1S2 gene that has a novel exon derived from two transposable elements (TEs), MIR and AluSp. To investigate the evolutionary footprint and molecular mechanism of action of this transcript, we performed PCR and RT-PCR experiments and sequencing analyses using genomic DNA and RNA samples from humans and various non-human primates. RESULTS: The results showed that the MIR element had integrated into the genome of our common ancestor, specifically in the BLOC1S2 gene region, before the radiation of all primate lineages and that the AluSp element had integrated into the genome of our common ancestor, fortunately in the middle of the MIR sequences, after the divergence of Old World monkeys and New World monkeys. The combined MIR and AluSp sequences provide a 3' splice site (AG) and 5' splice site (GT), respectively, and generate the Old World monkey-specific transcripts. Moreover, branch point sequences for the intron removal process are provided by the MIR and AluSp combination. CONCLUSIONS: We show for the first time that sequential integration into the same location and sequence divergence events of two different TEs generated lineage-specific transcripts through sequence collaboration during primate evolution.

Potential Applications of Non-thermal Plasma in Animal Husbandry to Improve Infrastructure.

Infrastructure in animal husbandry refers to fundamental facilities and services necessary for better living conditions of animals and its economy to function through better productivity. Mainly, infrastructure can be divided into two categories: hard infrastructure and soft infrastructure. Physical infrastructure, such as buildings, roads, and water supplying systems, belongs to hard infrastructure. Soft infrastructure includes services which are required to maintain economic, health, cultural and social standards of animal husbandry. Therefore, the proper management of infrastructure in animal husbandry is necessary for animal welfare and its economy. Among various technologies to improve the quality of infrastructure, non-thermal plasma (NTP) technology is an effectively applicable technology in different stages of animal husbandry. NTP is mainly helpful in maintaining better health conditions of animals in several ways via decontamination from microorganisms present in air, water, food, instruments and surfaces of animal farming systems. Furthermore, NTP is used in the treatment of waste water, vaccine production, wound healing in animals, odor-free ventilation, and packaging of animal food or animal products. This review summarizes the recent studies of NTP which can be related to the infrastructure in animal husbandry.

Successful application of human-based methyl capture sequencing for methylome analysis in non-human primate models.

BACKGROUND: The characterization of genomic or epigenomic variation in human and animal models could provide important insight into pathophysiological mechanisms of various diseases, and lead to new developments in disease diagnosis and clinical intervention. The African green monkey (AGM; Chlorocebus aethiops) and cynomolgus monkey (CM; Macaca fascicularis) have long been considered important animal models in biomedical research. However, non-human primate-specific methods applicable to epigenomic analyses in AGM and CM are lacking. The recent development of methyl-capture sequencing (MC-seq) has an unprecedented advantage of cost-effectiveness, and further allows for extending the methylome coverage compared to conventional sequencing approaches. RESULTS: Here, we used a human probe-designed MC-seq method to assay DNA methylation in DNA obtained from 13 CM and three AGM blood samples. To effectively adapt the human probe-designed target region for methylome analysis in non-human primates, we redefined the target regions, focusing on regulatory regions and intragenic regions with consideration of interspecific sequence homology and promoter region variation. Methyl-capture efficiency was controlled by the sequence identity between the captured probes based on the human reference genome and the AGM and CM genome sequences, respectively. Using reasonable guidelines, 56 and 62% of the human-based capture probes could be effectively mapped for DNA methylome profiling in the AGM and CM genome, respectively, according to numeric global statistics. In particular, our method could cover up to 89 and 87% of the regulatory regions of the AGM and CM genome, respectively. CONCLUSIONS: Use of human-based MC-seq methods provides an attractive, cost-effective approach for the methylome profiling of non-human primates at the single-base resolution level.

Identification and characterization of the tyrosinase gene (TYR) and its transcript variants (TYR_1 and TYR_2) in the crab-eating macaque (Macaca fascicularis).

Tyrosinase is a copper-containing enzyme that regulates melanin biosynthesis and is encoded by the tyrosinase (TYR) gene. Previous studies demonstrated that mutations in TYR could lead to oculocutaneous albinism type 1 (OCA1) owing to the failure of melanin formation. Although a previous study found that albinism in the rhesus monkey was derived from a mutation in TYR, the identification and characterization of this gene in non-human primates has not been achieved thus far. Thus, using the rapid amplification of cDNA ends (RACE) and internal reverse transcription PCR (RT-PCR) we identified the full-length sequence of TYR in the crab-eating macaque, and two different transcript variants (TYR_1 and TYR_2). While TYR_1 comprised five exons and its coding sequence was highly similar to that of humans, TYR_2 comprised four exons and was generated by a third-exon-skipping event. Interestingly, these two transcripts were also present in the African green monkey (Old World monkey) and the common marmoset (New World monkey). Deduced amino acid sequence analyses revealed that TYR_2 had a shorter C-terminal region than TYR_1 owing to the exon-skipping event. Thus, the present study is the first to identify and characterize a full-length TYR gene in a non-human primate, while the further validation of the third-exon-skipping in TYR indicates that this event is well conserved in the primate lineage. Therefore, this study provides useful and important information for the study of albinism using non-human primate models.

Alu-Derived Alternative Splicing Events Specific to Macaca Lineages in CTSF Gene.

Cathepsin F, which is encoded by CTSF, is a cysteine proteinase ubiquitously expressed in several tissues. In a previous study, novel transcripts of the CTSF gene were identified in the crab-eating monkey deriving from the integration of an Alu element-AluYRa1. The occurrence of AluYRa1-derived alternative transcripts and the mechanism of exonization events in the CTSF gene of human, rhesus monkey, and crab-eating monkey were investigated using PCR and reverse transcription PCR on the genomic DNA and cDNA isolated from several tissues. Results demonstrated that AluYRa1 was only integrated into the genome of Macaca species and this lineage-specific integration led to exonization events by producing a conserved 3' splice site. Six transcript variants (V1-V6) were generated by alternative splicing (AS) events, including intron retention and alternative 5' splice sites in the 5' and 3' flanking regions of CTSF_AluYRa1. Among them, V3-V5 transcripts were ubiquitously expressed in all tissues of rhesus monkey and crab-eating monkey, whereas AluYRa1-exonized V1 was dominantly expressed in the testis of the crab-eating monkey, and V2 was only expressed in the testis of the two monkeys. These five transcript variants also had different amino acid sequences in the C-terminal region of CTSF, as compared to reference sequences. Thus, species-specific Alu-derived exonization by lineage-specific integration of Alu elements and AS events seems to have played an important role during primate evolution by producing transcript variants and gene diversification.

Transcriptome analysis of paired primary colorectal carcinoma and liver metastases reveals fusion transcripts and similar gene expression profiles in primary carcinoma and liver metastases.

BACKGROUND: Despite the clinical significance of liver metastases, the difference between molecular and cellular changes in primary colorectal cancers (CRC) and matched liver metastases is poorly understood. METHODS: In order to compare gene expression patterns and identify fusion genes in these two types of tumors, we performed high-throughput transcriptome sequencing of five sets of quadruple-matched tissues (primary CRC, liver metastases, normal colon, and liver). RESULTS: The gene expression patterns in normal colon and liver were successfully distinguished from those in CRCs; however, RNA sequencing revealed that the gene expression between primary CRCs and their matched liver metastases is highly similar. We identified 1895 genes that were differentially expressed in the primary carcinoma and liver metastases, than that in the normal colon tissues. A major proportion of the transcripts, identified by gene expression profiling as significantly enriched in the primary carcinoma and metastases, belonged to gene ontology categories involved in the cell cycle, mitosis, and cell division. Furthermore, we identified gene fusion events in primary carcinoma and metastases, and the fusion transcripts were experimentally confirmed. Among these, a chimeric transcript resulting from the fusion of RNF43 and SUPT4H1 was found to occur frequently in primary colorectal carcinoma. In addition, knockdown of the expression of this RNF43-SUPT4H1 chimeric transcript was found to have a growth-inhibitory effect in colorectal cancer cells. CONCLUSIONS: The present study reports a high concordance of gene expression in the primary carcinoma and liver metastases, and reveals potential new targets, such as fusion genes, against primary and metastatic colorectal carcinoma.

Identification of Alternative Variants and Insertion of the Novel Polymorphic AluYl17 in TSEN54 Gene during Primate Evolution.

TSEN54 encodes a subunit of the tRNA-splicing endonuclease complex, which catalyzes the identification and cleavage of introns from precursor tRNAs. Previously, we identified an AluSx-derived alternative transcript in TSEN54 of cynomolgus monkey. Reverse transcription-polymerase chain reaction (RT-PCR) amplification and TSEN54 sequence analysis of primate and human samples identified five novel alternative transcripts, including the AluSx exonized transcript. Additionally, we performed comparative expression analysis via RT-qPCR in various cynomolgus, rhesus monkey, and human tissues. RT-qPCR amplification revealed differential expression patterns. Furthermore, genomic PCR amplification and sequencing of primate and human DNA samples revealed that AluSx elements were integrated in human and all of the primate samples tested. Intriguingly, in langur genomic DNA, an additional AluY element was inserted into AluSx of intron eight of TSEN54. The new AluY element showed polymorphic insertion. Using standardized nomenclature for Alu repeats, the polymorphic AluY of the langur TSEN54 was designated as being of the AluYl17 subfamily. Our results suggest that integration of the AluSx element in TSEN54 contributed to diversity in transcripts and induced lineage- or species-specific evolutionary events such as alternative splicing and polymorphic insertion during primate evolution.

Molecular characterization and expression analysis of mouse epidermal growth factor-like domain 8.

Epidermal growth factor (EGF)-like (EGFL) domain, a common structural module in numerous secreted or transmembrane proteins, is generally involved in protein-protein interactions. To date, several EGFL proteins have been identified and characterized, but little is known about EGFL domain 8 (EGFL8). The present study reported the molecular characterization and expression analysis of EGFL8 in mice. Mouse EGFL8 amplified using a reverse transcription-polymerase chain reaction approach was sequenced and characterized. Mouse EGFL8 encodes a protein of 293 amino acids with two EGFL domains, an Emilin-like domain and a Ca(2+)-binding EGFL domain, which has a molecular mass of 32 kDa. The coding sequence has a high degree of amino acid sequence identity across species, and the EGFL domain has been highly conserved in various species during evolutionary radiation. A phylogenetic tree calculated using the neighbor-joining method revealed that EGFL8 and EGFL7 are more closely associated with each other than either is to EGFL3, and they cluster with EGFL6. It was found that mouse EGFL8 protein was highly expressed in diverse mouse tissue types, including the thymus, lymph nodes, testis, ovaries, epididymis, ductus deferens, ileum, colon, stomach, esophagus, lung, uterus, urinary bladder, skin, spleen, adrenal glands and penis. These results are of great use in understanding the biological roles of mouse EGFL8 for further study.

Snail and serpinA1 promote tumor progression and predict prognosis in colorectal cancer.

The role of Snail and serpin peptidase inhibitor clade A member 1 (serpinA1) in tumorigenesis has been previously identified. However, the exact role and mechanism of these proteins in progression of colorectal cancer (CRC) are controversial. In this study, we investigated the role of Snail and serpinA1 in colorectal cancer (CRC) and examined the mechanisms through which these proteins mediate CRC progression. Immunohistochemical analysis of 528 samples from patients with CRC showed that elevated expression of Snail or serpinA1 was correlated with advanced stage, lymph node metastasis, and poor prognosis. Moreover, we detected a correlation between Snail and serpinA1 expression. Functional studies performed using the CRC cell lines DLD-1 and SW-480 showed that overexpression of Snail or serpinA1 significantly increased CRC cell invasion and migration. Conversely, knockdown of Snail or serpinA1 expression suppressed CRC cell invasion and migration. ChIP analysis revealed that Snail regulated serpinA1 by binding to its promoter. In addition, fibronectin mediated Snail and serpinA1 signaling was involved in CRC cell invasion and migration. Taken together, our data showed that Snail and serpinA1 promoted CRC progression through fibronectin. These findings suggested that Snail and serpinA1 were novel prognostic biomarkers and candidate therapeutic targets in CRC.

Genome-wide analysis of DNA methylation before-and after exercise in the thoroughbred horse with MeDIP-Seq.

Athletic performance is an important criteria used for the selection of superior horses. However, little is known about exercise-related epigenetic processes in the horse. DNA methylation is a key mechanism for regulating gene expression in response to environmental changes. We carried out comparative genomic analysis of genome-wide DNA methylation profiles in the blood samples of two different thoroughbred horses before and after exercise by methylated-DNA immunoprecipitation sequencing (MeDIP-Seq). Differentially methylated regions (DMRs) in the pre-and post-exercise blood samples of superior and inferior horses were identified. Exercise altered the methylation patterns. After 30 min of exercise, 596 genes were hypomethylated and 715 genes were hypermethylated in the superior horse, whereas in the inferior horse, 868 genes were hypomethylated and 794 genes were hypermethylated. These genes were analyzed based on gene ontology (GO) annotations and the exercise-related pathway patterns in the two horses were compared. After exercise, gene regions related to cell division and adhesion were hypermethylated in the superior horse, whereas regions related to cell signaling and transport were hypermethylated in the inferior horse. Analysis of the distribution of methylated CpG islands confirmed the hypomethylation in the gene-body methylation regions after exercise. The methylation patterns of transposable elements also changed after exercise. Long interspersed nuclear elements (LINEs) showed abundance of DMRs. Collectively, our results serve as a basis to study exercise-based reprogramming of epigenetic traits.

Genome-wide analysis of DNA methylation patterns in horse.

BACKGROUND: DNA methylation is an epigenetic regulatory mechanism that plays an essential role in mediating biological processes and determining phenotypic plasticity in organisms. Although the horse reference genome and whole transcriptome data are publically available the global DNA methylation data are yet to be known. RESULTS: We report the first genome-wide DNA methylation characteristics data from skeletal muscle, heart, lung, and cerebrum tissues of thoroughbred (TH) and Jeju (JH) horses, an indigenous Korea breed, respectively by methyl-DNA immunoprecipitation sequencing. The analysis of the DNA methylation patterns indicated that the average methylation density was the lowest in the promoter region, while the density in the coding DNA sequence region was the highest. Among repeat elements, a relatively high density of methylation was observed in long interspersed nuclear elements compared to short interspersed nuclear elements or long terminal repeat elements. We also successfully identified differential methylated regions through a comparative analysis of corresponding tissues from TH and JH, indicating that the gene body regions showed a high methylation density. CONCLUSIONS: We provide report the first DNA methylation landscape and differentially methylated genomic regions (DMRs) of thoroughbred and Jeju horses, providing comprehensive DMRs maps of the DNA methylome. These data are invaluable resource to better understanding of epigenetics in the horse providing information for the further biological function analyses.

Structure and Expression Analyses of SVA Elements in Relation to Functional Genes.

SINE-VNTR-Alu (SVA) elements are present in hominoid primates and are divided into 6 subfamilies (SVA-A to SVA-F) and active in the human population. Using a bioinformatic tool, 22 SVA element-associated genes are identified in the human genome. In an analysis of genomic structure, SVA elements are detected in the 5' untranslated region (UTR) of HGSNAT (SVA-B), MRGPRX3 (SVA-D), HYAL1 (SVA-F), TCHH (SVA-F), and ATXN2L (SVA-F) genes, while some elements are observed in the 3'UTR of SPICE1 (SVA-B), TDRKH (SVA-C), GOSR1 (SVA-D), BBS5 (SVA-D), NEK5 (SVA-D), ABHD2 (SVA-F), C1QTNF7 (SVA-F), ORC6L (SVA-F), TMEM69 (SVA-F), and CCDC137 (SVA-F) genes. They could contribute to exon extension or supplying poly A signals. LEPR (SVA-C), ALOX5 (SVA-D), PDS5B (SVA-D), and ABCA10 (SVA-F) genes also showed alternative transcripts by SVA exonization events. Dominant expression of HYAL1_SVA appeared in lung tissues, while HYAL1_noSVA showed ubiquitous expression in various human tissues. Expression of both transcripts (TDRKH_SVA and TDRKH_noSVA) of the TDRKH gene appeared to be ubiquitous. Taken together, these data suggest that SVA elements cause transcript isoforms that contribute to modulation of gene regulation in various human tissues.

Promoter activity analysis and methylation characterization of LTR elements of PERVs in NIH miniature pig.

The potential risk of porcine endogenous retrovirus (PERV) transmission is an important issue in xenotransplantation (pig-to-human transplantation). Long terminal repeats (LTRs) in PERV elements show promoter activity that could affect neighboring functional genes. The methylation status and promoter activities of 3 LTR structures (PERV-LTR1, LTR2, and LTR3 elements) belonging to the PERV-A family were examined using luciferase reporter genes in human liver cell lines (HepG2 and Hep3B). The PERV LTR3 element exhibited hypomethylation and stronger promoter activity than the other LTR elements in human liver cells. We also performed comparative sequences analysis of the PERV LTR elements by using bioinformatics tools. Our findings showed that several transcription factors such as Nkx2-2 and Elk-1 positively influenced the high transcriptional activity of the PERV LTR3 element.

A New Exon Derived from a Mammalian Apparent LTR Retrotransposon of the SUPT16H Gene.

The SUPT16H gene known as FACTP140 is required for the transcription of other genes. For transcription, genes need to be complexed with accessory factors, including transcription factors and RNA polymerase II. One such factor, FACT, interacts with histones H2A/H2B for nucleosome disassembly and transcription elongation. The SUPT16H gene has a transcript and many expressed sequence tags (ESTs). We were especially interested in an MaLR-derived transcript (EST, BX333035) that included a new exon introduced by a transposable element, a mammalian apparent LTR retrotransposon (MaLR). The MaLR was detected ranging from humans to galagos, indicating the MaLR in the SUPT16H gene is integrated into the primate ancestor genome. A new exon was created by alternative donor site provided by the MaLR. The original transcript and the MaLR-derived transcript were expressed in various human, rhesus monkey, and other primate tissues. Additionally, we identified a new alternative transcript that included the MaLR, but there was no significant difference in the expression of the original transcript and the MaLR-derived transcript. Interestingly, the new alternative transcript and the MaLR-derived transcript had the MaLR sequence in the new exon, but they had different structures by adopting different 3' splice sites. From this study, we verified transposable elements that contributed to transcriptome diversity.