Divergent gene expression patterns in alcohol and opioid use disorders lead to consistent alterations in functional networks within the Dorsolateral Prefrontal Cortex.

Substance Use Disorders (SUDs) manifest as persistent drug-seeking behavior despite adverse consequences, with Alcohol Use Disorder (AUD) and Opioid Use Disorder (OUD) representing prevalent forms associated with significant mortality rates and economic burdens. The co-occurrence of AUD and OUD is common, necessitating a deeper comprehension of their intricate interactions. While the causal link between these disorders remains elusive, shared genetic factors are hypothesized. Leveraging public datasets, we employed genomic and transcriptomic analyses to explore conserved and distinct molecular pathways within the dorsolateral prefrontal cortex associated with AUD and OUD. Our findings unveil modest transcriptomic overlap at the gene level between the two disorders but substantial convergence on shared biological pathways. Notably, these pathways predominantly involve inflammatory processes, synaptic plasticity, and key intracellular signaling regulators. Integration of transcriptomic data with the latest genome-wide association studies (GWAS) for problematic alcohol use (PAU) and OUD not only corroborated our transcriptomic findings but also confirmed the limited shared heritability between the disorders. Overall, our study indicates that while alcohol and opioids induce diverse transcriptional alterations at the gene level, they converge on select biological pathways, offering promising avenues for novel therapeutic targets aimed at addressing both disorders simultaneously.

Integration of selective sweeps across the sheep genome: understanding the relationship between production and adaptation traits.

BACKGROUND: Livestock populations are under constant selective pressure for higher productivity levels for different selective purposes. This pressure results in the selection of animals with unique adaptive and production traits. The study of genomic regions associated with these unique characteristics has the potential to improve biological knowledge regarding the adaptive process and how it is connected to production levels and resilience, which is the ability of an animal to adapt to stress or an imbalance in homeostasis. Sheep is a species that has been subjected to several natural and artificial selective pressures during its history, resulting in a highly specialized species for production and adaptation to challenging environments. Here, the data from multiple studies that aim at mapping selective sweeps across the sheep genome associated with production and adaptation traits were integrated to identify confirmed selective sweeps (CSS). RESULTS: In total, 37 studies were used to identify 518 CSS across the sheep genome, which were classified as production (147 prodCSS) and adaptation (219 adapCSS) CSS based on the frequency of each type of associated study. The genes within the CSS were associated with relevant biological processes for adaptation and production. For example, for adapCSS, the associated genes were related to the control of seasonality, circadian rhythm, and thermoregulation. On the other hand, genes associated with prodCSS were related to the control of feeding behaviour, reproduction, and cellular differentiation. In addition, genes harbouring both prodCSS and adapCSS showed an interesting association with lipid metabolism, suggesting a potential role of this process in the regulation of pleiotropic effects between these classes of traits. CONCLUSIONS: The findings of this study contribute to a deeper understanding of the genetic link between productivity and adaptability in sheep breeds. This information may provide insights into the genetic mechanisms that underlie undesirable genetic correlations between these two groups of traits and pave the way for a better understanding of resilience as a positive ability to respond to environmental stressors, where the negative effects on production level are minimized.

Genomic screening of allelic and genotypic transmission ratio distortion in horse.

The phenomenon in which the expected Mendelian inheritance is altered is known as transmission ratio distortion (TRD). The TRD analysis relies on the study of the transmission of one of the two alleles from a heterozygous parent to the offspring. These distortions are due to biological mechanisms affecting gametogenesis, embryo development and/or postnatal viability, among others. In this study, TRD phenomenon was characterized in horses using SNP-by-SNP model by TRDscan v.2.0 software. A total of 1,041 Pura Raza Español breed horses were genotyped with 554,634 SNPs. Among them, 277 horses genotyped in trios (stallion-mare-offspring) were used to perform the TRD analysis. Our results revealed 140 and 42 SNPs with allelic and genotypic patterns, respectively. Among them, 63 displayed stallion-TRD and 41 exhibited mare-TRD, while 36 SNPs showed overall TRD. In addition, 42 SNPs exhibited heterosis pattern. Functional analyses revealed that the annotated genes located within the TRD regions identified were associated with biological processes and molecular functions related to spermatogenesis, oocyte division, embryonic development, and hormonal activity. A total of 10 functional candidate genes related to fertility were found. To our knowledge, this is the most extensive study performed to evaluate the presence of alleles and functional candidate genes with transmission ratio distortion affecting reproductive performance in the domestic horse.

Unravelling transmission ratio distortion across the bovine genome: identification of candidate regions for reproduction defects.

BACKGROUND: Biological mechanisms affecting gametogenesis, embryo development and postnatal viability have the potential to alter Mendelian inheritance expectations resulting in observable transmission ratio distortion (TRD). Although the discovery of TRD cases have been around for a long time, the current widespread and growing use of DNA technologies in the livestock industry provides a valuable resource of large genomic data with parent-offspring genotyped trios, enabling the implementation of TRD approach. In this research, the objective is to investigate TRD using SNP-by-SNP and sliding windows approaches on 441,802 genotyped Holstein cattle and 132,991 (or 47,910 phased) autosomal SNPs. RESULTS: The TRD was characterized using allelic and genotypic parameterizations. Across the whole genome a total of 604 chromosomal regions showed strong significant TRD. Most (85%) of the regions presented an allelic TRD pattern with an under-representation (reduced viability) of carrier (heterozygous) offspring or with the complete or quasi-complete absence (lethality) for homozygous individuals. On the other hand, the remaining regions with genotypic TRD patterns exhibited the classical recessive inheritance or either an excess or deficiency of heterozygote offspring. Among them, the number of most relevant novel regions with strong allelic and recessive TRD patterns were 10 and 5, respectively. In addition, functional analyses revealed candidate genes regulating key biological processes associated with embryonic development and survival, DNA repair and meiotic processes, among others, providing additional biological evidence of TRD findings. CONCLUSIONS: Our results revealed the importance of implementing different TRD parameterizations to capture all types of distortions and to determine the corresponding inheritance pattern. Novel candidate genomic regions containing lethal alleles and genes with functional and biological consequences on fertility and pre- and post-natal viability were also identified, providing opportunities for improving breeding success in cattle.

Transcriptome analysis of perirenal fat from Spanish Assaf suckling lamb carcasses showing different levels of kidney knob and channel fat.

Introduction: Suckling lamb meat is highly appreciated in European Mediterranean countries because of its mild flavor and soft texture. In suckling lamb carcasses, perirenal and pelvic fat depots account for a large fraction of carcass fat accumulation, and their proportions are used as an indicator of carcass quality. Material and Methods: This study aimed to characterize the genetic mechanisms that regulate fat deposition in suckling lambs by evaluating the transcriptomic differences between Spanish Assaf lambs with significantly different proportions of kidney knob and channel fat (KKCF) depots in their carcasses (4 High-KKCF lambs vs. 4 Low-KKCF lambs). Results: The analyzed fat tissue showed overall dominant expression of white adipose tissue gene markers, although due to the young age of the animals (17-36 days), the expression of some brown adipose tissue gene markers (e.g., UCP1, CIDEA) was still identified. The transcriptomic comparison between the High-KKCF and Low-KKCF groups revealed a total of 80 differentially expressed genes (DEGs). The enrichment analysis of the 49 DEGs with increased expression levels in the Low-KKCF lambs identified significant terms linked to the biosynthesis of lipids and thermogenesis, which may be related to the higher expression of the UCP1 gene in this group. In contrast, the enrichment analysis of the 31 DEGs with increased expression in the High-KKCF lambs highlighted angiogenesis as a key biological process supported by the higher expression of some genes, such as VEGF-A and THBS1, which encode a major angiogenic factor and a large adhesive extracellular matrix glycoprotein, respectively. Discussion: The increased expression of sestrins, which are negative regulators of the mTOR complex, suggests that the preadipocyte differentiation stage is being inhibited in the High-KKCF group in favor of adipose tissue expansion, in which vasculogenesis is an essential process. All of these results suggest that the fat depots of the High-KKCF animals are in a later stage of development than those of the Low-KKCF lambs. Further genomic studies based on larger sample sizes and complementary analyses, such as the identification of polymorphisms in the DEGs, should be designed to confirm these results and achieve a deeper understanding of the genetic mechanisms underlying fat deposition in suckling lambs.

Intergenerational impact of dietary protein restriction in dairy ewes on epigenetic marks in the perirenal fat of their suckling lambs.

In sheep, nutrition during the prepubertal stage is essential for growth performance and mammary gland development. However, the potential effects of nutrient restriction in a prepuberal stage over the progeny still need to be better understood. Here, the intergenerational effect of maternal protein restriction at prepubertal age (2 months of age) on methylation patterns was evaluated in the perirenal fat of Assaf suckling lambs. In total, 17 lambs from ewes subjected to dietary protein restriction (NPR group, 44% less protein) and 17 lambs from control ewes (C group) were analyzed. These lambs were ranked based on their carcass proportion of perirenal and cavitary fat and classified into HighPCF and LowPCF groups. The perirenal tissue from 4 NPR-LowPCF, 4 NPR-HighPCF, 4 C-LowPCF, and 4 C-HighPCF lambs was subjected to whole-genome bisulfite sequencing and differentially methylated regions (DMRs) were identified. Among other relevant processes, these DMRs were mapped in genes responsible for regulating the transition of brown to white adipose tissue and nonshivering thermoregulation, which might be associated with better adaptation/survival of lambs in the perinatal stage. The current study provides important biological insights about the intergenerational effect on the methylation pattern of an NPR in replacement ewes.

Underlying genetic architecture of resistance to mastitis in dairy cattle: A systematic review and gene prioritization analysis of genome-wide association studies.

Mastitis, the most frequent disease in dairy cattle. Resistance to mastitis is a complex, polygenic trait controlled by several genes, each with small effects. Genome-wide association studies have been widely used to identify genomic variants associated with complex traits, including resistance to mastitis, to elucidate the underlying genetic architecture of the trait. However, no systematic review and gene prioritization analysis have been conducted to date on GWAS results for resistance to mastitis in dairy cattle. Hence, the objective was to perform a systematic review and gene prioritization analysis of GWAS studies to identify potential functional candidate genes associated with resistance to mastitis-related traits in dairy cattle. Four electronic databases were searched from inception to December 2020, supplemented with multiple sources of gray literature, to identify eligible articles. Annotation for genes and quantitative trait loci (QTL), and QTL enrichment analysis were conducted using GALLO. Gene prioritization analysis was performed by a guilty-by-association approach using GUILDify and ToppGene. From 52 articles included within this systematic review, 30 articles were used for further functional analyses. Gene and QTL annotation resulted in 9,125 and 43,646 unique genes and QTL, respectively, from 39 studies. In general, overlapping of genes across studies was very low (mean ± SD = 0.02% ± 0.07%). Most annotated genes were associated with somatic cell count-related traits and the Holstein breed. Within all annotated genes, 74 genes were shared among Holstein, Jersey, and Ayrshire breeds. Approximately 7.5% of annotated QTL were related to QTL class "health." Within the health QTL class, 2.6 and 2.2% of QTL were associated with clinical mastitis and somatic cell count-related traits. Enrichment analysis of QTL demonstrated that many enriched QTL were associated with somatic cell score located in Bos taurus autosomes 5, 6, 16, and 20. The prioritization analysis resulted in 427 significant genes after multiple test correction (false discovery rate of 5%) from 26 studies. Most prioritized genes were located in Bos taurus autosomes 19 and 7, and most top-ranked genes were from the cytokine superfamily (e.g., chemokines, interleukins, transforming growth factors, and tumor necrosis factor genes). Although most prioritized genes (397) were associated with somatic cell count-related traits, only 54 genes were associated with clinical mastitis-related traits. Twenty-four genes (ABCC9, ACHE, ADCYAP1, ARC, BCL2L1, CDKN1A, EPO, GABBR2, GDNF, GNRHR, IKBKE, JAG1, KCNJ8, KCNQ1, LIFR, MC3R, MYOZ3, NFKB1, OSMR, PPP3CA, PRLR, SHARPIN, SLC1A3, and TNFRSF25) were reported for both somatic cell count and clinical mastitis-related traits. Prioritized genes were mainly associated with immune response, regulation of secretion, locomotion, cell proliferation, and development. In conclusion, this study provided a fine-mapping of previously identified genomic regions associated with resistance to mastitis and identified key functional candidate genes for resistance to mastitis, which can be used to develop enhanced genomic strategies to combat mastitis by increasing mastitis resistance through genetic selection.

Unravelling the genetics of non-random fertilization associated with gametic incompatibility.

In the dairy industry, mate allocation is dependent on the producer's breeding goals and the parents' breeding values. The probability of pregnancy differs among sire-dam combinations, and the compatibility of a pair may vary due to the combination of gametic haplotypes. Under the hypothesis that incomplete incompatibility would reduce the odds of fertilization, and complete incompatibility would lead to a non-fertilizing or lethal combination, deviation from Mendelian inheritance expectations would be observed for incompatible pairs. By adding an interaction to a transmission ratio distortion (TRD) model, which detects departure from the Mendelian expectations, genomic regions linked to gametic incompatibility can be identified. This study aimed to determine the genetic background of gametic incompatibility in Holstein cattle. A total of 283,817 genotyped Holstein trios were used in a TRD analysis, resulting in 422 significant regions, which contained 2075 positional genes further investigated for network, overrepresentation, and guilt-by-association analyses. The identified biological pathways were associated with immunology and cellular communication and a total of 16 functional candidate genes were identified. Further investigation of gametic incompatibility will provide opportunities to improve mate allocation for the dairy cattle industry.

Deviations from Mendelian Inheritance on Bovine X-Chromosome Revealing Recombination, Sex-of-Offspring Effects and Fertility-Related Candidate Genes.

Transmission ratio distortion (TRD), or significant deviations from Mendelian inheritance, is a well-studied phenomenon on autosomal chromosomes, but has not yet received attention on sex chromosomes. TRD was analyzed on 3832 heterosomal single nucleotide polymorphisms (SNPs) and 400 pseudoautosomal SNPs spanning the length of the X-chromosome using 436,651 genotyped Holstein cattle. On the pseudoautosomal region, an opposite sire-TRD pattern between male and female offspring was identified for 149 SNPs. This finding revealed unique SNPs linked to a specific-sex (Y- or X-) chromosome and describes the accumulation of recombination events across the pseudoautosomal region. On the heterosomal region, 13 SNPs and 69 haplotype windows were identified with dam-TRD. Functional analyses for TRD regions highlighted relevant biological functions responsible to regulate spermatogenesis, development of Sertoli cells, homeostasis of endometrium tissue and embryonic development. This study uncovered the prevalence of different TRD patterns across both heterosomal and pseudoautosomal regions of the X-chromosome and revealed functional candidate genes for bovine reproduction.

Integrated analyses of the methylome and transcriptome to unravel sex differences in the perirenal fat from suckling lambs.

In sheep, differences were observed regarding fat accumulation and fatty acid (FA) composition between males and females, which may impact the quality and organoleptic characteristics of the meat. The integration of different omics technologies is a relevant approach for investigating biological and genetic mechanisms associated with complex traits. Here, the perirenal tissue of six male and six female Assaf suckling lambs was evaluated using RNA sequencing and whole-genome bisulfite sequencing (WGBS). A multiomic discriminant analysis using multiblock (s)PLS-DA allowed the identification of 314 genes and 627 differentially methylated regions (within these genes), which perfectly discriminate between males and females. These candidate genes overlapped with previously reported QTLs for carcass fat volume and percentage of different FAs in milk and meat from sheep. Additionally, differentially coexpressed (DcoExp) modules of genes between males (nine) and females (three) were identified that harbour 22 of these selected genes. Interestingly, these DcoExp were significantly correlated with fat percentage in different deposits (renal, pelvic, subcutaneous and intramuscular) and were associated with relevant biological processes for adipogenesis, adipocyte differentiation, fat volume and FA composition. Consequently, these genes may potentially impact adiposity and meat quality traits in a sex-specific manner, such as juiciness, tenderness and flavour.

Detection and Visualization of Heterozygosity-Rich Regions and Runs of Homozygosity in Worldwide Sheep Populations.

In this study, we chose 17 worldwide sheep populations of eight breeds, which were intensively selected for different purposes (meat, milk, or wool), or locally-adapted breeds, in order to identify and characterize factors impacting the detection of runs of homozygosity (ROH) and heterozygosity-rich regions (HRRs) in sheep. We also applied a business intelligence (BI) tool to integrate and visualize outputs from complementary analyses. We observed a prevalence of short ROH, and a clear distinction between the ROH profiles across populations. The visualizations showed a fragmentation of medium and long ROH segments. Furthermore, we tested different scenarios for the detection of HRR and evaluated the impact of the detection parameters used. Our findings suggest that HRRs are small and frequent in the sheep genome; however, further studies with higher density SNP chips and different detection methods are suggested for future research. We also defined ROH and HRR islands and identified common regions across the populations, where genes related to a variety of traits were reported, such as body size, muscle development, and brain functions. These results indicate that such regions are associated with many traits, and thus were under selective pressure in sheep breeds raised for different purposes. Interestingly, many candidate genes detected within the HRR islands were associated with brain integrity. We also observed a strong association of high linkage disequilibrium pattern with ROH compared with HRR, despite the fact that many regions in linkage disequilibrium were not located in ROH regions.

GALLO: An R package for genomic annotation and integration of multiple data sources in livestock for positional candidate loci.

BACKGROUND: The development of high-throughput sequencing and genotyping methodologies has enabled the identification of thousands of genomic regions associated with several complex traits. The integration of multiple sources of biological information is a crucial step required to better understand patterns regulating the development of these traits. FINDINGS: Genomic Annotation in Livestock for positional candidate LOci (GALLO) is an R package developed for the accurate annotation of genes and quantitative trait loci (QTLs) located in regions identified in common genomic analyses performed in livestock, such as genome-wide association studies and transcriptomics using RNA sequencing. Moreover, GALLO allows the graphical visualization of gene and QTL annotation results, data comparison among different grouping factors (e.g., methods, breeds, tissues, statistical models, studies), and QTL enrichment in different livestock species such as cattle, pigs, sheep, and chickens. CONCLUSIONS: Consequently, GALLO is a useful package for annotation, identification of hidden patterns across datasets, and data mining previously reported associations, as well as the efficient examination of the genetic architecture of complex traits in livestock.

Identifying functionally relevant candidate genes for inflexible ethanol intake in mice and humans using a guilt-by-association approach.

Gene prioritization approaches are useful tools to explore and select candidate genes in transcriptome studies. Knowing the importance of processes such as neuronal activity, intracellular signal transduction, and synapse plasticity to the development and maintenance of compulsive ethanol drinking, the aim of the present study was to explore and identify functional candidate genes associated with these processes in an animal model of inflexible pattern of ethanol intake. To do this, we applied a guilt-by-association approach, using the GUILDify and ToppGene software, in our previously published microarray data from the prefrontal cortex (PFC) and striatum of inflexible drinker mice. We then tested some of the prioritized genes that showed a tissue-specific pattern in postmortem brain tissue (PFC and nucleus accumbens (NAc)) from humans with alcohol use disorder (AUD). In the mouse brain, we prioritized 44 genes in PFC and 26 in striatum, which showed opposite regulation patterns in PFC and striatum. The most prioritized of them (i.e., Plcb1 and Prkcb in PFC, and Dnm2 and Lrrk2 in striatum) were associated with synaptic neuroplasticity, a neuroadaptation associated with excessive ethanol drinking. The identification of transcription factors among the prioritized genes suggests a crucial role for Irf4 in the pattern of regulation observed between PFC and striatum. Lastly, the differential transcription of IRF4 and LRRK2 in PFC and nucleus accumbens in postmortem brains from AUD compared to control highlights their involvement in compulsive ethanol drinking in humans and mice.

Weighted Gene Correlation Network Meta-Analysis Reveals Functional Candidate Genes Associated with High- and Sub-Fertile Reproductive Performance in Beef Cattle.

Improved reproductive efficiency could lead to economic benefits for the beef industry, once the intensive selection pressure has led to a decreased fertility. However, several factors limit our understanding of fertility traits, including genetic differences between populations and statistical limitations. In the present study, the RNA-sequencing data from uterine samples of high-fertile (HF) and sub-fertile (SF) animals was integrated using co-expression network meta-analysis, weighted gene correlation network analysis, identification of upstream regulators, variant calling, and network topology approaches. Using this pipeline, top hub-genes harboring fixed variants (HF × SF) were identified in differentially co-expressed gene modules (DcoExp). The functional prioritization analysis identified the genes with highest potential to be key-regulators of the DcoExp modules between HF and SF animals. Consequently, 32 functional candidate genes (10 upstream regulators and 22 top hub-genes of DcoExp modules) were identified. These genes were associated with the regulation of relevant biological processes for fertility, such as embryonic development, germ cell proliferation, and ovarian hormone regulation. Additionally, 100 candidate variants (single nucleotide polymorphisms (SNPs) and insertions and deletions (INDELs)) were identified within those genes. In the long-term, the results obtained here may help to reduce the frequency of subfertility in beef herds, reducing the associated economic losses caused by this condition.

Transcriptomic Profiles of Monocyte-Derived Macrophages in Response to Escherichia coli is Associated with the Host Genetics.

Reactive Nitrogen Species (RNS) are a group of bactericidal molecules produced by macrophages in response to pathogens in a process called oxidative burst. Nitric oxide (NO-) is a member of RNS produced from arginine by inducible Nitric Oxide Synthase (iNOS) enzyme. The activity of iNOS and production of NO- by macrophages following stimulation is one of the indicators of macrophage polarization towards M1/proinflammatory. Production of NO- by bovine monocyte-derived macrophage (MDM) and mouse peritoneal macrophages has been shown to be strongly associated with host genetic with the heritability of 0.776 in bovine MDM and 0.8 in mouse peritoneal macrophages. However, the mechanism of genetic regulation of macrophage response has remained less explored. In the current study, the transcriptome of bovine MDMs was compared between two extreme phenotypes that had been classified as high and low responder based on NO- production. The results showed that 179 and 392 genes were differentially expressed (DE) between high and low responder groups at 3 and 18 hours after exposure to Escherichia coli, respectively. A set of 11 Transcription Factors (TFs) (STAT1, IRF7, SPI1, STAT4, IRF1, HIF1A, FOXO3, REL, NFAT5, HIC1, and IRF4) at 3 hours and a set of 13 TFs (STAT1, IRF1, HIF1A, STAT4, ATF4, TP63, EGR1, CDKN2A, RBL1, E2F1, PRDM1, GATA3, and IRF4) at 18 hours after exposure to E. coli were identified to be differentially regulated between the high and low responder phenotypes. These TFs were found to be divided into two clusters of inflammatory- and hypoxia-related TFs. Functional analysis revealed that some key canonical pathways such as phagocytosis, chemotaxis, antigen presentation, and cell-to-cell signalling are enriched among the over-expressed genes by high responder phenotype. Based on the results of this study, it was inferred that the functional characteristics of bovine MDMs are associated with NO-based classification. Since NO- production is strongly associated with host genetics, this study for the first time shows the distinct proinflammatory profiles of macrophages are controlled by the natural genetic polymorphism in an outbred population. In addition, the results suggest that genetics can be considered as a new dimension in the current model of macrophage polarization which is currently described by the combination of stimulants, only.

Reducing cryptic relatedness in genomic data sets via a central node exclusion algorithm.

Cryptic relatedness is a confounding factor in genetic diversity and genetic association studies. Development of strategies to reduce cryptic relatedness in a sample is a crucial step for downstream genetic analyses. This study uses a node selection algorithm, based on network degrees of centrality, to evaluate its applicability and impact on evaluation of genetic diversity and population stratification. 1,036 Guzerá (Bos indicus) females were genotyped using Illumina Bovine SNP50 v2 BeadChip. Four strategies were compared. The first and second strategies consist on a iterative exclusion of most related individuals based on PLINK kinship coefficient (φij) and VanRaden's φij, respectively. The third and fourth strategies were based on a node selection algorithm. The fourth strategy, Network G matrix, preserved the larger number of individuals with a better diversity and representation from the initial sample. Determining the most probable number of populations was directly affected by the kinship metric. Network G matrix was the better strategy for reducing relatedness due to producing a larger sample, with more distant individuals, a more similar distribution when compared with the full data set in the MDS plots and keeping a better representation of the population structure. Resampling strategies using VanRaden's φij as a relationship metric was better to infer the relationships among individuals. Moreover, the resampling strategies directly impact the genomic inflation values in genomewide association studies. The use of the node selection algorithm also implies better selection of the most central individuals to be removed, providing a more representative sample.

Correction: Biotinidase deficiency: Genotype-biochemical phenotype association in Brazilian patients.

[This corrects the article DOI: 10.1371/journal.pone.0177503.].