ITGB2 is a central hub-gene associated with inflammation and early fibro-atheroma development in a swine model of atherosclerosis.

Background and aim: The complex dynamic interplay between different biological pathways involved in atherosclerosis development has rendered the identification of specific therapeutic targets a challenging quest. We aimed to identify specific genes and mechanistic pathways associated with the early development of fibro-atheromas in a swine model of atherosclerosis. Methods: The Wisconsin Miniature Swine™ model of Familial Hypercholesterolemia (WMS-FH, n = 11) and genetically related WMS controls (WMS-N, n = 11) were used. The infrarenal aorta was harvested from both groups for histopathologic and transcriptomic profiling at 12 months. Bioinformatic analysis was performed to identify hub genes and pathways central to disease pathophysiology. The expression of ITGB2, the top ranked hub gene, was manipulated in cell culture and the expression of interconnected genes was tested. Results: Fibro-atheromatous lesions were documented in all WMS-FH aortic tissues and displayed internal elastic lamina (IEL) disruption, significant reduction of myofibroblast presence and disorganized collagen deposition. No fibro-atheromas were observed in the control group. A total of 266 differentially expressed genes (DEGs) were upregulated in WMS-FH aortic tissues, while 29 genes were downregulated. Top identified hub genes included ITGB2, C1QA, LCP2, SPI1, CSF1R, C5AR1, CTSS, MPEG1, C1QC, and CSF2RB. Overexpression of ITGB2 resulted in elevated expression of other interconnected genes expressed in porcine endothelial cells. Conclusion: In a swine translational model of atherosclerosis, transcriptomic analysis identified ITGB2 as a central hub gene associated inflammation and early fibroatheroma development making it a potential therapeutic target at this stage of disease.

Disruption of the nuclear localization signal in RBM20 is causative in dilated cardiomyopathy.

Human patients carrying genetic mutations in RNA binding motif 20 (RBM20) develop a clinically aggressive dilated cardiomyopathy (DCM). Genetic mutation knockin (KI) animal models imply that altered function of the arginine-serine-rich (RS) domain is crucial for severe DCM. To test this hypothesis, we generated an RS domain deletion mouse model (Rbm20ΔRS). We showed that Rbm20ΔRS mice manifested DCM with mis-splicing of RBM20 target transcripts. We found that RBM20 was mis-localized to the sarcoplasm in Rbm20ΔRS mouse hearts and formed RBM20 granules similar to those detected in mutation KI animals. In contrast, mice lacking the RNA recognition motif showed similar mis-splicing of major RBM20 target genes but did not develop DCM or exhibit RBM20 granule formation. Using in vitro studies with immunocytochemical staining, we demonstrated that only DCM-associated mutations in the RS domain facilitated RBM20 nucleocytoplasmic transport and promoted granule assembly. Further, we defined the core nuclear localization signal (NLS) within the RS domain of RBM20. Mutation analysis of phosphorylation sites in the RS domain suggested that this modification may be dispensable for RBM20 nucleocytoplasmic transport. Collectively, our findings revealed that disruption of RS domain-mediated nuclear localization is crucial for severe DCM caused by NLS mutations.

The Activation of Protamine 1 Using Epigenome Editing Decreases the Proliferation of Tumorigenic Cells.

DNA methyltransferases (DNMT) and histone deacetylases (HDAC) inhibitors are used as cancer epigenome drugs. However, these epigenetic drugs lack targeting specificity and could risk inducing genome instability and the expression of oncogenes. Therefore, there is a need to develop new therapeutic strategies where specific cancer genes can be targeted for silencing or activation. The CRISPR/dCas9 system represents a promising, powerful therapeutic tool because of its simplicity and specificity. Protamine 1 (PRM1) is exclusively expressed in sperm and has a vital role in the tight packaging of DNA, thus inducing transcriptional silencing in sperm cells. We hypothesized that the activation of the PRM1 gene in tumorigenic cells would lead to DNA condensation and reduce the proliferation of these cells. To test our hypothesis, we transfected human embryonic kidney cells 293T with a dCas9-P300 plasmid that adds acetyl groups to the promoter region of PRM1 via specific gRNAs plasmids. RNA-Seq analysis of transfected cells revealed high specificity of targeted gene activation. PRM1 expression resulted in a significant decrease in cell proliferation as measured by the BrdU ELISA assay. To confirm that the activation of PRM1 was due to acetyl groups deposited to H3K27, a ChIP-qPCR was performed. The acetylation of the PRM1 promoter region targeted by dCas9-p300 in transfected cells was higher than that of the control cells. Interestingly, the targeted promoter region for acetylation showed reduced DNA methylation. These findings demonstrate the efficacy of epigenome editing in activating PRM1 in non-expressing tumorigenic cells, which could be used as a promising therapeutic strategy in cancer treatment.

RBM20 phosphorylation and its role in nucleocytoplasmic transport and cardiac pathogenesis.

Arginine-serine (RS) domain(s) in splicing factors are critical for protein-protein interaction in pre-mRNA splicing. Phosphorylation of RS domain is important for splicing control and nucleocytoplasmic transport in the cell. RNA-binding motif 20 (RBM20) is a splicing factor primarily expressed in the heart. A previous study using phospho-antibody against RS domain showed that RS domain can be phosphorylated. However, its actual phosphorylation sites and function have not been characterized. Using middle-down mass spectrometry, we identified 16 phosphorylation sites, two of which (S638 and S640 in rats, or S637 and S639 in mice) were located in the RSRSP stretch in the RS domain. Mutations on S638 and S640 regulated splicing, promoted nucleocytoplasmic transport and protein-RNA condensates. Phosphomimetic mutations on S638 and S640 indicated that phosphorylation was not the major cause for RBM20 nucleocytoplasmic transport and condensation in vitro. We generated a S637A knock-in (KI) mouse model (Rbm20S637A ) and observed the reduced RBM20 phosphorylation. The KI mice exhibited aberrant gene splicing, protein condensates, and a dilated cardiomyopathy (DCM)-like phenotype. Transcriptomic profiling demonstrated that KI mice had altered expression and splicing of genes involving cardiac dysfunction, protein localization, and condensation. Our in vitro data showed that phosphorylation was not a direct cause for nucleocytoplasmic transport and protein condensation. Subsequently, the in vivo results reveal that RBM20 mutations led to cardiac pathogenesis. However, the role of phosphorylation in vivo needs further investigation.

Investigating potential causal relationships among carcass and meat quality traits using structural equation model in Nellore cattle.

The objective of this study was to investigate potential causal relationships among hot carcass weight (HCW), longissimus muscle area (LMA), backfat thickness (BF), Warner-Bratzler shear force (WBSF), and marbling score (MB) traits in Nellore cattle using structural equation models (SEM). The SEM fitted comprises the following links between traits: WBSF → LMA, WBSF → HCW, HCW → LMA, BF → HCW, and BF → MB, where the arrows indicate the causal direction between traits, with structural coefficients posterior means (posterior standard deviation) equal to -0.29 cm2/kg (0.09), 0.43 kg/kg (0.29), 0.10 cm2/kg (0.006), 1.92 kg/mm (0.28), and 0.03 score-grade/mm (0.006), respectively. The final SEM revealed some important putative causal relationships among the traits studied here. The implied causal effects suggest that interventions on meat tenderness and fat content would affect overall growth and muscle deposition. Knowledge regarding potential causal relationships inferred among the traits studied here can have important implications for the genetic selection and management of Nellore cattle for improvement of carcass and meat quality.

RBM20S639G mutation is a high genetic risk factor for premature death through RNA-protein condensates.

Dilated cardiomyopathy (DCM) is a heritable and genetically heterogenous disease often idiopathic and a leading cause of heart failure with high morbidity and mortality. DCM caused by RNA binding motif protein 20 (RBM20) mutations is diverse and needs a more complete mechanistic understanding. RBM20 mutation S637G (S639G in mice) is linked to severe DCM and early death in human patients. In this study, we generated a RBM20 S639G mutation knock-in (KI) mouse model to validate the function of S639G mutation and examine the underlying mechanisms. KI mice exhibited severe DCM and premature death with a ~ 50% mortality in two months old homozygous (HM) mice. KI mice had enlarged atria and increased ANP and BNP biomarkers. The S639G mutation promoted RBM20 trafficking and ribonucleoprotein (RNP) granules in the sarcoplasm. RNA Seq data revealed differentially expressed and spliced genes were associated with arrhythmia, cardiomyopathy, and sudden death. KI mice also showed a reduction of diastolic stiffness and impaired contractility at both the left ventricular (LV) chamber and cardiomyocyte levels. Our results indicate that the RBM20 S639G mutation leads to RNP granules causing severe heart failure and early death and this finding strengthens the novel concept that RBM20 cardiomyopathy is a RNP granule disease.

Differentially expressed genes identified through RNA-seq with extreme values of principal components for beef fatty acid in Nelore cattle.

The aim of this study was to identify differentially expressed genes (DEG) in the Longissimus thoracis muscle of Nelore cattle related to fatty acid (FA) profile through RNA sequencing and principal component analysis (PCA). Two groups of 10 animals each were selected containing PC1 and PC2 extreme DEG values (HIGH × LOW) for each FA group. The intramuscular fat (IMF) was compared between cluster groups by ANOVA, and only the sum of monounsaturated FA (MUFA) and ω3 showed significant differences (p < .05). Interestingly, the highest percentage (95%) of phenotypic variation explained by the sum of the first two PC was observed for ω3, which also displayed the lowest number of DEG (n = 1). The lowest percentage (59%) was observed for MUFA, which also revealed the largest number of DEG (n = 66). Since only MUFA and ω3 exhibited significant differences between cluster groups, we can conclude that the differences observed for the remaining groups are not due to the percentage of IMF. Several genes that have been previously associated with meat quality and FA traits were identified as DEG in this study. The functional analysis revealed one KEGG pathway and eight GO terms as significant (p < .05), in which we highlighted the purine metabolism, glycolytic process, adenosine triphosphate binding and bone development. These results strongly contribute to the knowledge of the biological mechanisms involved in meat FA profile of Nelore cattle.

Genomic Study of Babesia bovis Infection Level and Its Association With Tick Count in Hereford and Braford Cattle.

Bovine babesiosis is a tick-borne disease caused by intraerythrocytic protozoa and leads to substantial economic losses for the livestock industry throughout the world. Babesia bovis is considered the most pathogenic species, which causes bovine babesiosis in Brazil. Genomic data could be used to evaluate the viability of improving resistance against B. bovis infection level (IB) through genomic selection, and, for that, knowledge of genetic parameters is needed. Furthermore, genome-wide association studies (GWAS) could be conducted to provide a better understanding of the genetic basis of the host response to B. bovis infection. No previous work in quantitative genetics of B. bovis infection was found. Thus, the objective of this study was to estimate the genetic correlation between IB and tick count (TC), evaluate predictive ability and applicability of genomic selection, and perform GWAS in Hereford and Braford cattle. The single-step genomic best linear unbiased prediction method was used, which allows the estimation of both breeding values and marker effects. Standard phenotyping was conducted for both traits. IB quantifications from the blood of 1,858 animals were carried using quantitative PCR assays. For TC, one to three subsequent tick counts were performed by manually counting adult female ticks on one side of each animal's body that was naturally exposed to ticks. Animals were genotyped using the Illumina BovineSNP50 panel. The posterior mean of IB heritability, estimated by the Bayesian animal model in a bivariate analysis, was low (0.10), and the estimations of genetic correlation between IB and TC were also low (0.15). The cross-validation genomic prediction accuracy for IB ranged from 0.18 to 0.35 and from 0.29 to 0.32 using k-means and random clustering, respectively, suggesting that genomic predictions could be used as a tool to improve genetics for IB, especially if a larger training population is developed. The top 10 single nucleotide polymorphisms from the GWAS explained 5.04% of total genetic variance for IB, which were located on chromosomes 1, 2, 5, 6, 12, 17, 18, 16, 24, and 26. Some candidate genes participate in immunity system pathways indicating that those genes are involved in resistance to B. bovis in cattle. Although the genetic correlation between IB and TC was weak, some candidate genes for IB were also reported in tick infestation studies, and they were also involved in biological resistance processes. This study contributes to improving genetic knowledge regarding infection by B. bovis in cattle.

Comparison between haplotype-based and individual snp-based genomic predictions for beef fatty acid profile in Nelore cattle.

The aim of this study was to evaluate the genomic predictions using the single-step genomic best linear unbiased predictor (ssGBLUP) method based on SNPs and haplotype markers associated with beef fatty acids (FAs) profile in Nelore cattle. The data set contained records from 963 Nelore bulls finished in feedlot (±90 days) and slaughtered with approximately 24 months of age. Meat samples from the Longissimus dorsi muscle were taken for FAs profile measurement. FAs were quantified by gas chromatography using a SP-2560 capillary column. Animals were genotyped with the high-density SNP panel (BovineHD BeadChip assay) containing 777,962 markers. SNPs with a minor allele frequency and a call rate lower than 0.05 and 0.90, respectively, monomorphic, located on sex chromosomes, and with unknown position were removed from the data set. After genomic quality control, a total of 469,981 SNPs and 892 samples were available for subsequent analyses. Missing genotypes were imputed and phased using the FImpute software. Haplotype blocks were defined based on linkage disequilibrium using the Haploview software. The model to estimate variance components and genetic parameters and to predict the genomic values included the random genetic additive effects, fixed effects of the contemporary group and the age at slaughter as a linear covariate. Accuracies using the haplotype-based approach ranged from 0.07 to 0.31, and those SNP-based ranged from 0.06 to 0.33. Regression coefficients ranged from 0.07 to 0.74 and from 0.08 to 1.45 using the haplotype- and SNP-based approaches, respectively. Despite the low to moderate accuracies for the genomic values, it is possible to obtain genetic progress trough selection using genomic information based either on SNPs or haplotype markers. The SNP-based approach allows less biased genomic evaluations, and it is more feasible when taking into account the computational and operational cost underlying the haplotypes inference.

Genomic selection for meat quality traits in Nelore cattle.

The objective of this study was to present heritability estimates and accuracy of genomic prediction using different methods for meat quality traits in Nelore cattle. Approximately 5000 animals with phenotypes and genotypes of 412,000 SNPs, were divided into two groups: (1) training population: animals born from 2008 to 2013 and (2) validation population: animals born in 2014. A single-trait animal model was used to estimate heritability and to adjust the phenotype. The methods of GBLUP, Improved Bayesian Lasso and Bayes Cπ were performed to estimate the SNP effects. Accuracy of genomic prediction was calculated using Pearson's correlations between direct genomic values and adjusted phenotypes, divided by the square root of heritability of each trait (0.03-0.19). The accuracies varied from 0.23 to 0.73, with the lowest accuracies estimated for traits associated with fat content and the greatest accuracies observed for traits of meat color and tenderness. There were small differences in genomic prediction accuracy between methods.

Polymorphism analysis in genes of the somatotropic axis in Nellore cattle selected for growth.

The somatotropic axis consists of genes that are involved in muscular development. These genes are potential regions of study to identify possible QTL for economically important traits in beef cattle. The aim of this study was to verify the existence of GH1, POU1F1, and GHR polymorphisms in Nellore cattle to verify the influence of selection in these mutations and to analyse the association between molecular markers and body weight at different ages, yearling hip height, carcass fat thickness and loin eye area. Six hundred forty-five animals from the Centro APTA Bovinos de Corte, were genotyped by PCR-RFLP techniques. The association analyses were performed with general mixed models taking into consideration the effect of one marker, and other model taking into consideration interactions between two molecular markers. Only the molecular markers rs81109601 on GH1 and rs109136815 on GHR were polymorphic; however, they were not found to be under selection. The association of the GHR rs109136815 marker and loin eye area was observed (p<0.05), as well as the effect of interaction between the markers and the female body weight at 550 days of age (p<0.04). The interaction effect should be considered in situations where the interactivity between two genes is known.

First polymorphisms in JY-1 gene in cattle (Bos taurus indicus) and their association with sexual precocity and growth traits.

Protein JY-1 is a bovine oocyte-specific protein that regulates granulosa cell function and is involved in early embryonic development, influencing the chance of pregnancy. This study investigated molecular markers for the JY-1 gene. Seven SNPs were identified in exon 3 of the gene. The positions of the SNPs in the exon and the respective substitutions are: 163 (T/C), 281 (T/C), 321 (T/C), 532 (T/C), 652 (A/G), 679 (T/C), and 722 (G/C) (GenBank: JN592587 and JF262042.2). SNP 163 is located in a coding region and causes a proline-to-leucine substitution. The other SNPs are located in the 3'UTR region. SNPs 163, 281, 321, and 679 were genotyped in 297 Nellore heifers and the haplotypes were constructed. The haplotypes of JY-1 were not correlated with the traits studied at 5 %.