Introduction to DOK2 and its potential role in cancer

Cis- and trans-eQTL TWASs of breast and ovarian cancer identify more than 100 susceptibility genes in the BCAC and OCAC consortia

Transcriptome-wide association studies (TWASs) have investigated the role of genetically regulated transcriptional activity in the etiologies of breast and ovarian cancer. However, methods performed to date have focused on the regulatory effects of risk-associated SNPs thought to act in cis on a nearby target gene. With growing evidence for distal (trans) regulatory effects of variants on gene expression, we performed TWASs of breast and ovarian cancer using a Bayesian genome-wide TWAS method (BGW-TWAS) that considers effects of both cis- and trans-expression quantitative trait loci (eQTLs). We applied BGW-TWAS to whole-genome and RNA sequencing data in breast and ovarian tissues from the Genotype-Tissue Expression project to train expression imputation models. We applied these models to large-scale GWAS summary statistic data from the Breast Cancer and Ovarian Cancer Association Consortia to identify genes associated with risk of overall breast cancer, non-mucinous epithelial ovarian cancer, and 10 cancer subtypes. We identified 101 genes significantly associated with risk with breast cancer phenotypes and 8 with ovarian phenotypes. These loci include established risk genes and several novel candidate risk loci, such as ACAP3, whose associations are predominantly driven by trans-eQTLs. We replicated several associations using summary statistics from an independent GWAS of these cancer phenotypes. We further used genotype and expression data in normal and tumor breast tissue from the Cancer Genome Atlas to examine the performance of our trained expression imputation models. This work represents an in-depth look into the role of trans eQTLs in the complex molecular mechanisms underlying these diseases.

Epigenetic disruptions in the offspring hypothalamus in response to maternal infection

DNA methylation is an epigenetic modification that can alter gene expression, and the incidence can vary across developmental stages, inflammatory conditions, and sexes. The effects of viral maternal viral infection and sex on the DNA methylation patterns were studied in the hypothalamus of a pig model of immune activation during development. DNA methylation at single-base resolution in regions of high CpG density was measured on 24 individual hypothalamus samples using reduced representation bisulfite sequencing. Differential over- and under-methylated sites were identified and annotated to proximal genes and corresponding biological processes. A total of 120 sites were differentially methylated (FDR-adjusted p-value < 0.05) between maternal infection or sex groups. Among the 66 sites differentially methylated between groups exposed to inflammatory signals and control, most sites were over-methylated in the challenged group and included sites in the promoter regions of genes SIRT3 and NRBP1. Among the 54 differentially methylated sites between females and males, most sites were over-methylated in females and included sites in the promoter region of genes TNC and EIF4G1. The analysis of the genes proximal to the differentially methylated sites suggested that biological processes potentially impacted include immune response, neuron migration and ensheathment, peptide signaling, adaptive thermogenesis, and tissue development. These results suggest that translational studies should consider that the prolonged effect of maternal infection during gestation may be enacted through epigenetic regulatory mechanisms that may differ between sexes.

Expression and Mutation of SLC45A2 Affects Iris Color in Quail

Iris color is a prominent phenotypic feature of quail. To understand the mechanism of melanin deposition related to quail iris color, iris tissues were selected from Beijing white and Chinese yellow quail for transcriptome analysis. Differentially expressed genes (DEGs) associated with pigmentation were identified using RNA sequencing and validated by quantitative real-time polymerase chain reaction (RT-qPCR). The identified single nucleotide polymorphisms were studied using bioinformatics and iris color correlation analyses. A total of 485 DEGs were obtained, with 223 upregulated and 262 downregulated. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed. Thirty-two genes were annotated using the GO database. Three important pigment synthesis pathways (Notch signaling, melanogenesis, and tyrosine metabolism) were identified in quail iris tissue (P < 0.05). The expression levels of solute carrier family 45 member 2 (SLC45A2), tyrosinase-related protein 1, vitamin D receptor, opsin 5, and docking protein 5 were significantly different between Beijing white and Chinese yellow quail, as verified by RT-qPCR. The c.1061C>T mutation in SLC45A2, which caused a single amino acid change at position 354 (threonine to methionine), was significantly associated with iris color in Beijing white and Chinese yellow quail, and might be the main reason for the different iris colors between these two quail species.

Cis- and trans-eQTL TWAS of breast and ovarian cancer identify more than 100 risk associated genes in the BCAC and OCAC consortia

Transcriptome-wide association studies (TWAS) have investigated the role of genetically regulated transcriptional activity in the etiologies of breast and ovarian cancer. However, methods performed to date have only considered regulatory effects of risk associated SNPs thought to act in cis on a nearby target gene. With growing evidence for distal (trans) regulatory effects of variants on gene expression, we performed TWAS of breast and ovarian cancer using a Bayesian genome-wide TWAS method (BGW-TWAS) that considers effects of both cis- and trans-expression quantitative trait loci (eQTLs). We applied BGW-TWAS to whole genome and RNA sequencing data in breast and ovarian tissues from the Genotype-Tissue Expression project to train expression imputation models. We applied these models to large-scale GWAS summary statistic data from the Breast Cancer and Ovarian Cancer Association Consortia to identify genes associated with risk of overall breast cancer, non-mucinous epithelial ovarian cancer, and 10 cancer subtypes. We identified 101 genes significantly associated with risk with breast cancer phenotypes and 8 with ovarian phenotypes. These loci include established risk genes and several novel candidate risk loci, such as ACAP3, whose associations are predominantly driven by trans-eQTLs. We replicated several associations using summary statistics from an independent GWAS of these cancer phenotypes. We further used genotype and expression data in normal and tumor breast tissue from the Cancer Genome Atlas to examine the performance of our trained expression imputation models. This work represents a first look into the role of trans-eQTLs in the complex molecular mechanisms underlying these diseases.

The transcriptome signature analysis of the epithelial-mesenchymal transition and immune cell infiltration in colon adenocarcinoma

The epithelial-mesenchymal transition (EMT) process is tightly connected to tumors' immune microenvironment. In colon adenocarcinoma (COAD), both the EMT and immune cell infiltration contribute to tumor progression; however, several questions regarding the mechanisms governing the interaction between EMT and the immune response remain unanswered. Our study aims to investigate the cross-talk between these two processes in cases of COAD and identify the key regulators involved. We utilized the EMT and immune signatures of samples from the COAD-TCGA database to identify three subtypes of COAD: high mesenchymal, medium mesenchymal, and low mesenchymal. We observed that EMT was associated with increased tumor immune response and infiltration mediated by pro-inflammatory cytokines. However, EMT was also linked to immunosuppressive activity that involved regulatory T cells, dendritic cells, and the upregulated expression of multiple immune checkpoints, such as PD-1, PDL-1, CTLA-4, and others. Finally, we employed the multivariate random forest feature importance method to identify key genes, such as DOK2 and MSRB3, that may play crucial roles in both EMT and the intratumoral immune response.

Identification of Candidate Genes for Min Pig Villi Hair Traits by Genome-Wide Association of Copy Number Variation

The Min pig is a famous native pig breed in northeast China, which has the special genetic character of villi hair growth in cold seasons. At present, little research has focused on the genetic mechanism of villi hair growth in Min pigs. Copy number variations (CNVs) are a type of variant that may influence many traits. In this study, we first investigated the phenotype of Large White × Min pigs' F2 pig villi hair in detail and then performed a CNV-based genome-wide association study (GWAS) between CNVs and pig villi hair appearance. Finally, a total number of 15 significant CNVRs were found to be associated with Min pig villi hair. The most significant CNVR was located on chromosome 1. Nearest gene annotation analysis indicated that the pig villi hair traits may be associated with the biological process of the G-protein-coupled receptor signaling pathway. QTL overlapping analysis found that among the CNVRs, 14 CNVRs could be co-located with known QTLs. Some genes such as MCHR2, LTBP2, and GFRA2 may be candidate genes for pig villi traits and are worth further study. Our study may provide a basic reference for the selection and breeding of cold-resistant pigs and outdoor breeding.

Genome-Wide Population Structure and Selection Signatures of Yunling Goat Based on RAD-seq

Simple Summary

Goats are important domestic animals that provide meat, milk, fur, and other products for humans. The demand for these products has increased in recent years. Disease resistance among goat breeds is different, but the genetic basis of the differences in resistance to diseases is still unclear and needs to be further studied. In this study, many genes and pathways related to immunity and diseases were identified to be under positive selection between Yunling and Nubian goats using RAD-seq technology. This study on the selection signatures of Yunling goats provides the scientific basis and technical support for the breeding of domestic goats for disease resistance, which has important social and economic significance.

Abstract

Animal diseases impose a huge burden on the countries where diseases are endemic. Conventional control strategies of vaccines and veterinary drugs are to control diseases from a pharmaceutical perspective. Another alternative approach is using pre-existing genetic disease resistance or tolerance. We know that the Yunling goat is an excellent local breed from Yunnan, southwestern China, which has characteristics of strong disease resistance and remarkable adaptability. However, genetic information about the selection signatures of Yunling goats is limited. We reasoned that the genes underlying the observed difference in disease resistance might be identified by investigating selection signatures between two different goat breeds. Herein, we selected the Nubian goat as the reference group to perform the population structure and selection signature analysis by using RAD-seq technology. The results showed that two goat breeds were divided into two clusters, but there also existed gene flow. We used Fst (F-statistics) and π (pi/θπ) methods to carry out selection signature analysis. Eight selected regions and 91 candidate genes were identified, in which some genes such as DOK2, TIMM17A, MAVS, and DOCK8 related to disease and immunity and some genes such as SPEFI, CDC25B, and MIR103 were associated with reproduction. Four GO (Gene Ontology) terms (GO:0010591, GO:001601, GO:0038023, and GO:0017166) were associated with cell migration, signal transduction, and immune responses. The KEGG (Kyoto Encyclopedia of Genes and Genomes) signaling pathways were mainly associated with immune responses, inflammatory responses, and stress reactions. This study preliminarily revealed the genetic basis of strong disease resistance and adaptability of Yunling goats. It provides a theoretical basis for the subsequent genetic breeding of disease resistance of goats.

Identification of genes associated with feather color in Liancheng white duck using FST analysis

Liancheng white duck has two phenotypic traits: white feather and black beak-black foot, but the genes controlling these phenotypic traits are unknown. The objective of this study is to identify various candidate genes related to the plumage of Liancheng white duck. This study used F₂ population construction generated between white Kaiya duck and Liancheng white duck and F_ST analysis between the dominant and recessive loci associated with the Liancheng white duck white feather in order to identify specific gene regions. As per the feather color statistics of the F₂ population, it is estimated that there are about three or four genes controlling the white feather of Liancheng white ducks, and the F_ST results showed that four significant signals were found on chromosomes 4, 12, 13, and 21. Further annotation of these regions led to the identification of five genes involved in the melanin pathway, namely, KIT, CLOCK, MITF, CEBPA, and DOK5. Among them, CEBPA and DOK5 might be affecting the white feather traits of Liancheng white duck by regulating the melanin production and its transfer to the feather. The results provide insightful understanding into the genetic mechanisms of white feather in Liancheng white duck.