A Comprehensive Allele Specific Expression Resource for the Equine Transcriptome.

Background: Allele-specific expression (ASE) analysis provides a nuanced view of cis-regulatory mechanisms affecting gene expression. Results: An equine ASE analysis was performed, using integrated Iso-seq and short-read RNA sequencing data from four healthy Thoroughbreds (2 mares and 2 stallions) across 9 tissues from the Functional Annotation of Animal Genomes (FAANG) project. Allele expression was quantified by haplotypes from long-read data, with 42,900 allele expression events compared. Within these events, 635 (1.48%) demonstrated ASE, with liver tissue containing the highest proportion. Genetic variants within ASE events were in histone modified regions 64.2% of the time. Validation of allele-specific variants, using a set of 66 equine liver samples from multiple breeds, confirmed that 97% of variants demonstrated ASE. Conclusions: This valuable publicly accessible resource is poised to facilitate investigations into regulatory variation in equine tissues. Our results highlight the tissue-specific nature of allelic imbalance in the equine genome.

A Comprehensive Allele Specific Expression Resource for the Equine Transcriptome.

Background: Allele-specific expression (ASE) analysis provides a nuanced view of cis-regulatory mechanisms affecting gene expression. Results: In this work, we introduce and highlight the significance of an equine ASE analysis, containing integrated long- and short-read RNA sequencing data, along with insight from histone modification data, from four healthy Thoroughbreds (2 mares and 2 stallions) across 9 tissues. Conclusions: This valuable publicly accessible resource is poised to facilitate investigations into regulatory variation in equine tissues and foster a deeper understanding of the impact of allelic imbalance in equine health and disease at the molecular level.

Trio-binning of a hinny refines the comparative organization of the horse and donkey X chromosomes and reveals novel species-specific features.

We generated single haplotype assemblies from a hinny hybrid which significantly improved the gapless contiguity for horse and donkey autosomal genomes and the X chromosomes. We added over 15 Mb of missing sequence to both X chromosomes, 60 Mb to donkey autosomes and corrected numerous errors in donkey and some in horse reference genomes. We resolved functionally important X-linked repeats: the DXZ4 macrosatellite and ampliconic Equine Testis Specific Transcript Y7 (ETSTY7). We pinpointed the location of the pseudoautosomal boundaries (PAB) and determined the size of the horse (1.8 Mb) and donkey (1.88 Mb) pseudoautosomal regions (PARs). We discovered distinct differences in horse and donkey PABs: a testis-expressed gene, XKR3Y, spans horse PAB with exons1-2 located in Y and exon3 in the X-Y PAR, whereas the donkey XKR3Y is Y-specific. DXZ4 had a similar ~ 8 kb monomer in both species with 10 copies in horse and 20 in donkey. We assigned hundreds of copies of ETSTY7, a sequence horizontally transferred from Parascaris and massively amplified in equids, to horse and donkey X chromosomes and three autosomes. The findings and products contribute to molecular studies of equid biology and advance research on X-linked conditions, sex chromosome regulation and evolution in equids.

First gene-edited calf with reduced susceptibility to a major viral pathogen.

Bovine viral diarrhea virus (BVDV) is one of the most important viruses affecting the health and well-being of bovine species throughout the world. Here, we used CRISPR-mediated homology-directed repair and somatic cell nuclear transfer to produce a live calf with a six amino acid substitution in the BVDV binding domain of bovine CD46. The result was a gene-edited calf with dramatically reduced susceptibility to infection as measured by reduced clinical signs and the lack of viral infection in white blood cells. The edited calf has no off-target edits and appears normal and healthy at 20 months of age without obvious adverse effects from the on-target edit. This precision bred, proof-of-concept animal provides the first evidence that intentional genome alterations in the CD46 gene may reduce the burden of BVDV-associated diseases in cattle and is consistent with our stepwise, in vitro and ex vivo experiments with cell lines and matched fetal clones.

The Assembled Genome of the Stroke-Prone Spontaneously Hypertensive Rat.

BACKGROUND: We report the creation and evaluation of a de novo assembly of the genome of the spontaneously hypertensive rat, the most widely used model of human cardiovascular disease. METHODS: The genome is assembled from long read sequencing (PacBio HiFi and continuous long read data [CLR]) and scaffolded with long-range structural information obtained from Bionano optical maps and proximity ligation sequencing proximity analysis of the genome. The genome assembly was polished with Illumina short reads. Completeness of the assembly was investigated using Benchmarking Universal Single Copy Orthologs analysis. The genome assembly was also evaluated with the rat reference gene set, using NCBI automated protocols. We also generated orthogonal single molecule transcript sequence reads (Iso-Seq) from 8 tissues and used them to validate the coding assembly, to annotate the assembly with RNA transcripts representing unique full length transcript isoforms for each gene and to determine whether divergences between RefSeq sequences and the assembly were attributable to assembly errors or polymorphisms. RESULTS: The assembly analysis indicates that this assembly is comparable in contiguity and completeness to the current rat reference assembly, while the use of HiFi sequencing yields an assembly that is more correct at the single base level. Synteny analysis was performed to uncover the extent of synteny and the presence and distribution of chromosomal rearrangements between the reference and this assembly. CONCLUSION: The resulting genome assembly is reference quality and captures significant structural variation.

Enhancing the compatibility of BioCaRGOS silica sol-gel technology with ctDNA extraction and droplet digital PCR (ddPCR) analysis.

Previously, our group had demonstrated long term stabilization of protein biomarkers using BioCaRGOS, a silica sol-gel technology. Herein, we describe workflow modifications to allow for extraction of cell free DNA (cfDNA) from primary samples containing working concentrations of BioCaRGOS, as well as the compatibility of BioCaRGOS with droplet digital PCR (ddPCR) analysis for pancreatic cancer biomarkers i.e., KRAS circulating tumor DNA (ctDNA). Preliminary attempts to extract ctDNA from BioCaRGOS containing samples demonstrated interference in the extraction of primary samples and the interference with ddPCR analysis when BioCaRGOS was directly introduced to stabilize sample extracts. In our modified technique, we have minimized the interference caused by methanol with ddPCR by complete removal of methanol from the activated BioCaRGOS formulation prior to addition to the biospecimen or ctDNA extract. Interference of the silica matrix present in BioCaRGOS with ctDNA extraction was eliminated through the introduction of invert filtration of the sample prior to extraction. These modifications to the workflow of BioCaRGOS containing samples allow for use of BioCaRGOS for stabilization of trace quantities of nucleic acid biomarkers such as plasma ctDNA, while retaining the capability to extract the biomarker and quantify based on ddPCR.

mRatBN7.2: familiar and unfamiliar features of a new rat genome reference assembly.

Rat genomic tools have been slower to emerge than for those of humans and mice and have remained less thorough and comprehensive. The arrival of a new and improved rat reference genome, mRatBN7.2, in late 2020 is a welcome event. This assembly, like predecessor rat reference assemblies, is derived from an inbred Brown Norway rat. In this "user" survey we hope to provide other users of this assembly some insight into its characteristics and some assessment of its improvements as well as a few caveats that arise from the unique aspects of this assembly. mRatBN7.2 was generated by the Wellcome Sanger Institute as part of the large Vertebrate Genomes Project. This rat assembly has now joined human, mouse, chicken, and zebrafish in the National Center for Biotechnology Information (NCBI)'s Genome Reference Consortium, which provides ongoing curation of the assembly. Here we examine the technical procedures by which the assembly was created and assess how this assembly constitutes an improvement over its predecessor. We also indicate the technical limitations affecting the assembly, providing illustrations of how these limitations arise and the impact that results for this reference assembly.

Delineating the Effects of Passaging and Exposure in a Longitudinal Study of Arsenic-Induced Squamous Cell Carcinoma in a HaCaT Cell Line Model.

Cutaneous squamous cell carcinoma (cSCC) is a major deleterious health effect of chronic arsenic (iAs) exposure. The molecular mechanism of arsenic-induced cSCC remains poorly understood. We recently demonstrated that chronic iAs exposure leads to temporally regulated genome-wide changes in profiles of differentially expressed mRNAs and miRNAs at each stage of carcinogenesis (7, 19, and 28 weeks) employing a well-established passage-matched HaCaT cell line model of arsenic-induced cSCC. Here, we performed longitudinal differential expression analysis (miRNA and mRNA) between the different time points (7 vs 19 weeks and 19 vs 28 weeks) within unexposed and exposed groups, coupled to expression pairing and pathway analyses to differentiate the relative effects of long-term passaging and chronic iAs exposure. Data showed that 66-105 miRNA [p < .05; log2(fold change) > I1I] and 2826-4079 mRNA [p < .001; log2(fold change) > I1I] molecules were differentially expressed depending on the longitudinal comparison. Several mRNA molecules differentially expressed as a function of time, independent of iAs exposure were being targeted by miRNA molecules which were also differentially expressed in a time-dependent manner. Distinct pathways were predicted to be modulated as a function of time or iAs exposure. Some pathways were also modulated both by time and exposure. Thus, the HaCaT model can distinguish between the effects of passaging and chronic iAs exposure individually and corroborate our previously published data on effects of iAs exposure compared with unexposed passage matched HaCaT cells. In addition, this work provides a template for cell line-based longitudinal chronic exposure studies to follow for optimal efficacy.

Genetics of Thoroughbred Racehorse Performance.

Thoroughbred horses have been selected for racing performance for more than 400 years. Despite continued selection, race times have not improved significantly during the past 60 years, raising the question of whether genetic variation for racing performance still exists. Studies using phenotypes such as race time, money earned, and handicapping, however, demonstrate that there is extensive variation within these traits and that they are heritable. Even so, these are poor measures of racing success since Thoroughbreds race at different ages and distances and on different types of tracks, and some may not race at all. With the advent of genomic tools, DNA variants are being identified that contribute to racing success. Aside from strong associations for myostatin variants with best racing distance, weak to modest associations with racing phenotypes are reported for other genomic regions. These data suggest that diverse genetic strategies have contributed to producing a successful racehorse, and genetic variation contributing to athleticism remains important.

Using whole genome sequence to compare variant callers and breed differences of US sheep.

As whole genome sequence (WGS) data sets have become abundant and widely available, so has the need for variant detection and scoring. The aim of this study was to compare the accuracy of commonly used variant calling programs, Freebayes and GATK HaplotypeCaller (GATK-HC), and to use U.S. sheep WGS data sets to identify novel breed-associated SNPs. Sequence data from 145 sheep consisting of 14 U.S. breeds were filtered and biallelic single nucleotide polymorphisms (SNPs) were retained for genotyping analyses. Genotypes from both programs were compared to each other and to genotypes from bead arrays. The SNPs from WGS were compared to the bead array data with breed heterozygosity, principal component analysis and identifying breed associated SNPs to analyze genetic diversity. The average sequence read depth was 2.78 reads greater with 6.11% more SNPs being identified in Freebayes compared to GATK-HC. The genotype concordance of the variant callers to bead array data was 96.0% and 95.5% for Freebayes and GATK-HC, respectively. Genotyping with WGS identified 10.5 million SNPs from all 145 sheep. This resulted in an 8% increase in measured heterozygosity and greater breed separation in the principal component analysis compared to the bead array analysis. There were 1,849 SNPs identified in only the Romanov sheep where all 10 rams were homozygous for one allele and the remaining 135 sheep from 13 breeds were homozygous for the opposite allele. Both variant calling programs had greater than 95% concordance of SNPs with bead array data, and either was suitably accurate for ovine WGS data sets. The use of WGS SNPs improved the resolution of PCA analysis and was critical for identifying Romanov breed-associated SNPs. Subsets of such SNPs could be used to estimate germplasm composition in animals without pedigree information.

Association of ARRDC3 and NFIA variants with bovine congestive heart failure in feedlot cattle.

Background: Bovine congestive heart failure (BCHF) has become increasingly prevalent among feedlot cattle in the Western Great Plains of North America with up to 7% mortality in affected herds. BCHF is an untreatable complex condition involving pulmonary hypertension that culminates in right ventricular failure and death. Genes associated with BCHF in feedlot cattle have not been previously identified. Our aim was to search for genomic regions associated with this disease. Methods: A retrospective, matched case-control design with 102 clinical BCHF cases and their unaffected pen mates was used in a genome-wide association study. Paired nominal data from approximately 560,000 filtered single nucleotide polymorphisms (SNPs) were analyzed with McNemar's test. Results: Two independent genomic regions were identified as having the most significant association with BCHF: the arrestin domain-containing protein 3 gene ( ARRDC3), and the nuclear factor IA gene ( NFIA, mid- p-values, 1x10 -8 and 2x10 -7, respectively). Animals with two copies of risk alleles at either gene were approximately eight-fold more likely to have BCHF than their matched pen mates with either one or zero risk alleles at both genes (CI 95 = 3-17). Further, animals with two copies of risk alleles at both genes were 28-fold more likely to have BCHF than all others ( p-value = 1×10 -7, CI 95 = 4-206). A missense variant in ARRDC3 (C182Y) represents a potential functional variant since the C182 codon is conserved among all other jawed vertebrate species observed. A two-SNP test with markers in both genes showed 29% of 273 BCHF cases had homozygous risk genotypes in both genes, compared to 2.5% in 198 similar unaffected feedlot cattle. This and other DNA tests may be useful for identifying feedlot animals with the highest risk for BCHF in the environments described here. Conclusions: Although pathogenic roles for variants in the ARRDC3 and NFIA genes are unknown, their discovery facilitates classifying animals by genetic risk and allows cattle producers to make informed decisions for selective breeding and animal health management.

Transcriptional and Histochemical Signatures of Bone Marrow Mononuclear Cell-Mediated Resolution of Synovitis.

Osteoarthritis (OA) may result from impaired ability of synovial macrophages to resolve joint inflammation. Increasing macrophage counts in inflamed joints through injection with bone marrow mononuclear cells (BMNC) induces lasting resolution of synovial inflammation. To uncover mechanisms by which BMNC may affect resolution, in this study, differential transcriptional signatures of BMNC in response to normal (SF) and inflamed synovial fluid (ISF) were analyzed. We demonstrate the temporal behavior of co-expressed gene networks associated with traits from related in vivo and in vitro studies. We also identified activated and inhibited signaling pathways and upstream regulators, further determining their protein expression in the synovium of inflamed joints treated with BMNC or DPBS controls. BMNC responded to ISF with an early pro-inflammatory response characterized by a short spike in the expression of a NF-ƙB- and mitogen-related gene network. This response was associated with sustained increased expression of two gene networks comprising known drivers of resolution (IL-10, IGF-1, PPARG, isoprenoid biosynthesis). These networks were common to SF and ISF, but more highly expressed in ISF. Most highly activated pathways in ISF included the mevalonate pathway and PPAR-γ signaling, with pro-resolving functional annotations that improve mitochondrial metabolism and deactivate NF-ƙB signaling. Lower expression of mevalonate kinase and phospho-PPARγ in synovium from inflamed joints treated with BMNC, and equivalent IL-1ß staining between BMNC- and DPBS-treated joints, associates with accomplished resolution in BMNC-treated joints and emphasize the intricate balance of pro- and anti-inflammatory mechanisms required for resolution. Combined, our data suggest that BMNC-mediated resolution is characterized by constitutively expressed homeostatic mechanisms, whose expression are enhanced following inflammatory stimulus. These mechanisms translate into macrophage proliferation optimizing their capacity to counteract inflammatory damage and improving their general and mitochondrial metabolism to endure oxidative stress while driving tissue repair. Such effect is largely achieved through the synthesis of several lipids that mediate recovery of homeostasis. Our study reveals candidate mechanisms by which BMNC provide lasting improvement in patients with OA and suggests further investigation on the effects of PPAR-γ signaling enhancement for the treatment of arthritic conditions.

Evaluating Large Spontaneous Deletions in a Bovine Cell Line Selected for Bovine Viral Diarrhea Virus Resistance.

Bovine viral diarrhea virus's (BVDV) entry into bovine cells involves attachment of virions to cellular receptors, internalization, and pH-dependent fusion with endosomal membranes. The primary host receptor for BVDV is CD46; however, the complete set of host factors required for virus entry is unknown. The Madin-Darby bovine kidney (MDBK) cell line is susceptible to BVDV infection, while a derivative cell line (CRIB) is resistant at the level of virus entry. We performed complete genome sequencing of each to identify genomic variation underlying the resistant phenotype with the aim of identifying host factors essential for BVDV entry. Three large compound deletions in the BVDV-resistant CRIB cell line were identified and predicted to disrupt the function or expression of the genes PTPN12, GRID2, and RABGAP1L. However, CRISPR/Cas9 mediated knockout of these genes, individually or in combination, in the parental MDBK cell line did not impact virus entry or replication. Therefore, resistance to BVDV in the CRIB cell line is not due to the apparent spontaneous loss of PTPN12, GRID2, or RABGAP1L gene function. Identifying the functional cause of BVDV resistance in the CRIB cell line may require more detailed comparisons of the genomes and epigenomes.

Successful ATAC-Seq From Snap-Frozen Equine Tissues.

An assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) has become an increasingly popular method to assess genome-wide chromatin accessibility in isolated nuclei from fresh tissues. However, many biobanks contain only snap-frozen tissue samples. While ATAC-seq has been applied to frozen brain tissues in human, its applicability in a wide variety of tissues in horse remains unclear. The Functional Annotation of Animal Genome (FAANG) project is an international collaboration aimed to provide high quality functional annotation of animal genomes. The equine FAANG initiative has generated a biobank of over 80 tissues from two reference female animals and experiments to begin to characterize tissue specificity of genome function for prioritized tissues have been performed. Due to the logistics of tissue collection and storage, extracting nuclei from a large number of tissues for ATAC-seq at the time of collection is not always practical. To assess the feasibility of using stored frozen tissues for ATAC-seq and to provide a guideline for the equine FAANG project, we compared ATAC-seq results from nuclei isolated from frozen tissue to cryopreserved nuclei (CN) isolated at the time of tissue harvest in liver, a highly cellular homogenous tissue, and lamina, a relatively acellular tissue unique to the horse. We identified 20,000-33,000 accessible chromatin regions in lamina and 22-61,000 in liver, with consistently more peaks identified using CN isolated at time of tissue collection. Our results suggest that frozen tissues are an acceptable substitute when CN are not available. For more challenging tissues such as lamina, nuclei extraction at the time of tissue collection is still preferred for optimal results. Therefore, tissue type and accessibility to intact nuclei should be considered when designing ATAC-seq experiments.

Transcriptomic analysis of equine chorioallantois reveals immune networks and molecular mechanisms involved in nocardioform placentitis.

Nocardioform placentitis (NP) continues to result in episodic outbreaks of abortion and preterm birth in mares and remains a poorly understood disease. The objective of this study was to characterize the transcriptome of the chorioallantois (CA) of mares with NP. The CA were collected from mares with confirmed NP based upon histopathology, microbiological culture and PCR for Amycolatopsis spp. Samples were collected from the margin of the NP lesion (NPL, n = 4) and grossly normal region (NPN, n = 4). Additionally, CA samples were collected from normal postpartum mares (Control; CRL, n = 4). Transcriptome analysis identified 2892 differentially expressed genes (DEGs) in NPL vs. CRL and 2450 DEGs in NPL vs. NPN. Functional genomics analysis elucidated that inflammatory signaling, toll-like receptor signaling, inflammasome activation, chemotaxis, and apoptosis pathways are involved in NP. The increased leukocytic infiltration in NPL was associated with the upregulation of matrix metalloproteinase (MMP1, MMP3, and MMP8) and apoptosis-related genes, such as caspases (CASP3 and CASP7), which could explain placental separation associated with NP. Also, NP was associated with downregulation of several placenta-regulatory genes (ABCG2, GCM1, EPAS1, and NR3C1), angiogenesis-related genes (VEGFA, FLT1, KDR, and ANGPT2), and glucose transporter coding genes (GLUT1, GLUT10, and GLUT12), as well as upregulation of hypoxia-related genes (HIF1A and EGLN3), which could elucidate placental insufficiency accompanying NP. In conclusion, our findings revealed for the first time, the key regulators and mechanisms underlying placental inflammation, separation, and insufficiency during NP, which might lead to the development of efficacious therapies or diagnostic aids by targeting the key molecular pathways.

Dynamic alteration in miRNA and mRNA expression profiles at different stages of chronic arsenic exposure-induced carcinogenesis in a human cell culture model of skin cancer.

Chronic arsenic exposure causes skin cancer, although the underlying molecular mechanisms are not well defined. Altered microRNA and mRNA expression likely play a pivotal role in carcinogenesis. Changes in genome-wide differential expression of miRNA and mRNA at 3 strategic time points upon chronic sodium arsenite (As3+) exposure were investigated in a well-validated HaCaT cell line model of arsenic-induced cutaneous squamous cell carcinoma (cSCC). Quadruplicate independent HaCaT cell cultures were exposed to 0 or 100 nM As3+ for up to 28-weeks (wk). Cell growth was monitored throughout the course of exposure and epithelial-mesenchymal transition (EMT) was examined employing immunoblot. Differentially expressed miRNA and mRNA profiles were generated at 7, 19, and 28-wk by RNA-seq, followed by identification of differentially expressed mRNA targets of differentially expressed miRNAs through expression pairing at each time point. Pathway analyses were performed for total differentially expressed mRNAs and for the miRNA targeted mRNAs at each time point. RNA-seq predictions were validated by immunoblot of selected target proteins. While the As3+-exposed cells grew slower initially, growth was equal to that of unexposed cells by 19-wk (transformation initiation), and exposed cells subsequently grew faster than passage-matched unexposed cells. As3+-exposed cells had undergone EMT at 28-wk. Pathway analyses demonstrate dysregulation of carcinogenesis-related pathways and networks in a complex coordinated manner at each time point. Immunoblot data largely corroborate RNA-seq predictions in the endoplasmic reticulum stress (ER stress) pathway. This study provides a detailed molecular picture of changes occurring during the arsenic-induced transformation of human keratinocytes.

Decreased Tumoral Expression of Colon-Specific Water Channel Aquaporin 8 Is Associated With Reduced Overall Survival in Colon Adenocarcinoma.

BACKGROUND: Colon cancer survival is dependent on metastatic potential and treatment. Large RNA-sequencing data sets may assist in identifying colon cancer-specific biomarkers to improve patient outcomes. OBJECTIVE: This study aimed to identify a highly specific biomarker for overall survival in colon adenocarcinoma by using an RNA-sequencing data set. DESIGN: Raw RNA-sequencing and clinical data for patients with colon adenocarcinoma (n = 271) were downloaded from The Cancer Genome Atlas. A binomial regression model was used to calculate differential RNA expression between paired colon cancer and normal epithelium samples (n = 40). Highly differentially expressed RNAs were examined. SETTINGS: This study was conducted at the University of Louisville using data acquired by The Cancer Genome Atlas. PATIENTS: Patients from US accredited cancer centers between 1998 and 2013 were analyzed. MAIN OUTCOME MEASURES: The primary outcome measures were recurrence-free and overall survival. RESULTS: The median age was 66 years (147/271 men, 180/271 White patients). Thirty RNAs were differentially expressed in colon adenocarcinoma compared with paired normal epithelium, using a log-fold change cutoff of ±6. Using median expression as a cutoff, 4 RNAs were associated with worse overall survival: decreased ZG16 (log-rank = 0.023), aquaporin 8 (log-rank = 0.023), and SLC26A3 (log-rank = 0.098), and increased COL1A1 (log-rank = 0.105). On multivariable analysis, low aquaporin 8 expression (HR, 1.748; 95% CI, 1.016-3.008; p = 0.044) was a risk factor for worse overall survival. Our final aquaporin 8 model had an area under the curve of 0.85 for overall survival. On subgroup analysis, low aquaporin 8 was associated with worse overall survival in patients with high microsatellite instability and in patients with stage II disease. Low aquaporin 8 expression was associated with KRAS and BRAF mutations. Aquaporin 8 immunohistochemistry was optimized for clinical application. LIMITATIONS: This was a retrospective study. CONCLUSION: Aquaporin 8 is a water channel selectively expressed in normal colon tissue. Low aquaporin 8 expression is a risk factor for worse overall survival in patients who have colon cancer. Aquaporin 8 measurement may have a role as a colon-specific prognostic biomarker and help in patient risk stratification for increased surveillance. See Video Abstract at http://links.lww.com/DCR/B603. LA DISMINUCIN DE LA EXPRESIN TUMORAL DE LA ACUAPORINA DEL CANAL DE AGUA ESPECFICO DEL COLON SE ASOCIA CON UNA REDUCCIN DE LA SUPERVIVENCIA GENERAL EN EL ADENOCARCINOMA DE COLON: ANTECEDENTES:La supervivencia del cáncer de colon depende del potencial metastásico y del tratamiento. Grandes conjuntos de datos de secuenciación de ARN pueden ayudar a identificar biomarcadores específicos del cáncer de colon para mejorar los resultados de los pacientes.OBJETIVO:Identificar un biomarcador altamente específico para la supervivencia general en el adenocarcinoma de colon utilizando un conjunto de datos de secuenciación de ARN.DISEÑO:La secuenciación de ARN sin procesar y los datos clínicos para pacientes con adenocarcinoma de colon (n = 271) se descargaron de The Cancer Genome Atlas. Se utilizó un modelo de regresión binomial para calcular la expresión diferencial de ARN entre muestras de cáncer de colon emparejadas y muestras de epitelio normal (n = 40). Se examinaron los ARN expresados de forma altamente diferencial.ENTORNO CLINICO:Este estudio se realizó en la Universidad de Louisville utilizando datos adquiridos por The Cancer Genome Atlas.PACIENTES:Se analizaron pacientes de centros oncológicos acreditados en Estados Unidos entre 1998-2013.PRINCIPALES MEDIDAS DE VALORACION:Las principales medidas de valoración fueron la supervivencia general y libre de recurrencia.RESULTADOS:La mediana de edad fue de 66 años (147/271 hombres, 180/271 caucásicos). Treinta ARN se expresaron diferencialmente en el adenocarcinoma de colon en comparación con el epitelio normal emparejado, utilizando un límite de cambio logarítmico de ± 6. Utilizando la expresión mediana como punto de corte, cuatro ARN se asociaron con una peor supervivencia general: disminución de ZG16 (rango logarítmico = 0,023), acuaporina8 (rango logarítmico = 0,023) y SLC26A3 (rango logarítmico = 0,098) y aumento de COL1A1 (log -rango = 0,105). En el análisis multivariable, la baja expresión de acuaporina8 (HR = 1,748, IC del 95%: 1,016-3,008, p = 0,044) fue un factor de riesgo para una peor supervivencia global. Nuestro modelo de aquaporin8 final tuvo un AUC de 0,85 para la supervivencia global. En el análisis de subgrupos, la acuaporina8 baja se asoció con una peor supervivencia general en pacientes con MSI-H y en pacientes en estadio II. La baja expresión de acuaporina8 se asoció con mutaciones de KRAS y BRAF. La inmunohistoquímica de aquaporina8 se optimizó para su aplicación clínica.LIMITACIONES:Este fue un estudio retrospectivo.CONCLUSIÓN:La acuaporina8 es un canal de agua expresado selectivamente en el tejido normal del colon. La baja expresión de AQP8 es un factor de riesgo de peor supervivencia global en pacientes con cáncer de colon. La medición de aquaporina8 puede tener un papel como un biomarcador de pronóstico específico del colon y ayudar en la estratificación del riesgo del paciente para una mayor vigilancia. Consulte Video Resumen en http://links.lww.com/DCR/B603.

Parental bias in expression and interaction of genes in the equine placenta.

Most autosomal genes in the placenta show a biallelic expression pattern. However, some genes exhibit allele-specific transcription depending on the parental origin of the chromosomes on which the copy of the gene resides. Parentally expressed genes are involved in the reciprocal interaction between maternal and paternal genes, coordinating the allocation of resources between fetus and mother. One of the main challenges of studying parental-specific allelic expression (allele-specific expression [ASE]) in the placenta is the maternal cellular remnant at the fetomaternal interface. Horses (Equus caballus) have an epitheliochorial placenta in which both the endometrial epithelium and the epithelium of the chorionic villi are juxtaposed with minimal extension into the uterine mucosa, yet there is no information available on the allelic gene expression of equine chorioallantois (CA). In the current study, we present a dataset of 1,336 genes showing ASE in the equine CA (https://pouya-dini.github.io/equine-gene-db/) along with a workflow for analyzing ASE genes. We further identified 254 potentially imprinted genes among the parentally expressed genes in the equine CA and evaluated the expression pattern of these genes throughout gestation. Our gene ontology analysis implies that maternally expressed genes tend to decrease the length of gestation, while paternally expressed genes extend the length of gestation. This study provides fundamental information regarding parental gene expression during equine pregnancy, a species with a negligible amount of maternal cellular remnant in its placenta. This information will provide the basis for a better understanding of the role of parental gene expression in the placenta during gestation.

Long non-coding RNA ZFAS1 is a major regulator of epithelial-mesenchymal transition through miR-200/ZEB1/E-cadherin, vimentin signaling in colon adenocarcinoma.

Colon adenocarcinoma is a common cause of cancer-related deaths worldwide. Epithelial-mesenchymal transition is a major regulator of cancer metastasis, and increased understanding of this process is essential to improve patient outcomes. Long non-coding RNA (lncRNA) are important regulators of carcinogenesis. To identify lncRNAs associated with colon carcinogenesis, we performed an exploratory differential gene expression analysis comparing paired colon adenocarcinoma and normal colon epithelium using an RNA-sequencing data set. This analysis identified lncRNA ZFAS1 as significantly increased in colon cancer compared to normal colon epithelium. This finding was validated in an institutional cohort using laser capture microdissection. ZFAS1 was also found to be principally located in the cellular cytoplasm. ZFAS1 knockdown was associated with decreased cellular proliferation, migration, and invasion in two colon cancer cell lines (HT29 and SW480). MicroRNA-200b and microRNA-200c (miR-200b and miR-200c) are experimentally validated targets of ZFAS1, and this interaction was confirmed using reciprocal gene knockdown. ZFAS1 knockdown regulated ZEB1 gene expression and downstream targets E-cadherin and vimentin. Knockdown of miR-200b or miR-200c reversed the effect of ZFAS1 knockdown in the ZEB1/E-cadherin, vimentin signaling cascade, and the effects of cellular migration and invasion, but not cellular proliferation. ZFAS1 knockdown was also associated with decreased tumor growth in an in vivo mouse model. These results demonstrate the critical importance of ZFAS1 as a regulator of the miR-200/ZEB1/E-cadherin, vimentin signaling cascade.