Cross-Species Comparative DNA Methylation Reveals Novel Insights into Complex Trait Genetics among Cattle, Sheep, and Goats.

The cross-species characterization of evolutionary changes in the functional genome can facilitate the translation of genetic findings across species and the interpretation of the evolutionary basis underlying complex phenotypes. Yet, this has not been fully explored between cattle, sheep, goats, and other mammals. Here, we systematically characterized the evolutionary dynamics of DNA methylation and gene expression in 3 somatic tissues (i.e. brain, liver, and skeletal muscle) and sperm across 7 mammalian species, including 3 ruminant livestock species (cattle, sheep, and goats), humans, pigs, mice, and dogs, by generating and integrating 160 DNA methylation and transcriptomic data sets. We demonstrate dynamic changes of DNA hypomethylated regions and hypermethylated regions in tissue-type manner across cattle, sheep, and goats. Specifically, based on the phylo-epigenetic model of DNA methylome, we identified a total of 25,074 hypomethylated region extension events specific to cattle, which participated in rewiring tissue-specific regulatory network. Furthermore, by integrating genome-wide association studies of 50 cattle traits, we provided novel insights into the genetic and evolutionary basis of complex phenotypes in cattle. Overall, our study provides a valuable resource for exploring the evolutionary dynamics of the functional genome and highlights the importance of cross-species characterization of multiomics data sets for the evolutionary interpretation of complex phenotypes in cattle livestock.

Y-Shaped Backbone-Rigidified DNA Tiles for the Construction of Supersized Nondeformable Tetrahedrons for Precise Cancer Therapies.

While engineered DNA nanoframeworks have been extensively exploited for delivery of diagnostic and therapeutic regents, DNA tiling-based DNA frameworks amenable to applications in living systems lag much behind. In this contribution, by developing a Y-shaped backbone-based DNA tiling technique, we assemble Y-shaped backbone-rigidified supersized DNA tetrahedrons (RDT) with 100% efficiency for precisely targeted tumor therapy. RDT displays unparalleled rigidness and unmatched resistance to nuclease degradation so that it almost does not deform under the force exerted by the atomic force microscopy tip, and the residual amount is not less than 90% upon incubating in biological media for 24 h, displaying at least 11.6 times enhanced degradation resistance. Without any targeting ligand, RDT enters the cancer cell in a targeted manner, and internalization specificity is up to 15.8. Moreover, 77% of RDT objects remain intact within living cells for 14 h. The drug loading content of RDT is improved by 4-8 times, and RDT almost 100% eliminates the unintended drug leakage in a stimulated physiological medium. Once systemically administrated into HeLa tumor-bearing mouse models, doxorubicin-loaded RDTs preferentially accumulate in tumor sites and efficiently suppress tumor growth without detectable off-target toxicity. The Y-DNA tiling technique offers invaluable insights into the development of structural DNA nanotechnology for precise medicine.

Quantifying Information of Dynamical Biochemical Reaction Networks.

A large number of complex biochemical reaction networks are included in the gene expression, cell development, and cell differentiation of in vivo cells, among other processes. Biochemical reaction-underlying processes are the ones transmitting information from cellular internal or external signaling. However, how this information is measured remains an open question. In this paper, we apply the method of information length, based on the combination of Fisher information and information geometry, to study linear and nonlinear biochemical reaction chains, respectively. Through a lot of random simulations, we find that the amount of information does not always increase with the length of the linear reaction chain; instead, the amount of information varies significantly when this length is not very large. When the length of the linear reaction chain reaches a certain value, the amount of information hardly changes. For nonlinear reaction chains, the amount of information changes not only with the length of this chain, but also with reaction coefficients and rates, and this amount also increases with the length of the nonlinear reaction chain. Our results will help to understand the role of the biochemical reaction networks in cells.

Differential Membrane Protein Profile in Bovine X- and Y-Sperm.

The aim of this study was to understand the molecular mechanisms behind the biological differences of X- and Y-sperm and to screen the sex-specific candidate antigen proteins for sexed semen production. To this end, we investigated differential expression of total membrane proteins of the two sperm types by using high-purity X- and Y-sperm from 20 Holstein bulls and applying the label-free proteomic technique; 1521 proteins were identified. In the X-sperm group, 8 and 23 proteins were significantly up- and down-regulated, respectively. In the X- and the Y-sperm group, 151 and 88 proteins were specifically expressed, respectively. These were overexpressed in the dynamic changes of the actin cytoskeleton, and cell senescence/apoptosis induced by the immune response, and could result in differences in the state, size, and immune sensitivity of the X-/Y-sperm membranes. The prediction of transmembrane structure, subcellular localization, and Western blotting validation results showed that the CLRN3 and SCAMP1 proteins were cell surface specific antigens of X- and Y-sperm, respectively. Our findings help explain the molecular mechanism behind the biological differences of X-/Y-sperm and lay the foundation for application of immunological methods to produce sex-sorted semen and control livestock sex. Data are available via ProteomeXchange with identifier PXD019435.

A Pilot Feasibility Study of Yttrium-90 Liver Radioembolization Followed by Durvalumab and Tremelimumab in Patients with Microsatellite Stable Colorectal Cancer Liver Metastases.

LESSONS LEARNED: Radioembolization with yttrium-90 resin microspheres can be combined safely with full doses of durvalumab and tremelimumab in patients with metastatic colorectal cancer. Regional radioembolization with yttrium-90 resin microspheres did not result in any hepatic or extrahepatic responses to a combination of durvalumab and tremelimumab. The lack of immunomodulatory responses to yttrium-90 on biopsies before and after treatment rules out a potential role for this strategy in converting a "cold tumor" into an "inflamed," immune responsive tumor. BACKGROUND: PD-1 inhibitors have been ineffective in microsatellite stable (MSS) metastatic colorectal cancer (CRC). Preclinical models suggest that radiation therapy may sensitize MSS CRC to PD-1 blockade. METHODS: Patients with MSS metastatic CRC with liver-predominant disease who progressed following at least one prior line of treatment were treated with yttrium-90 (Y90) radioembolization to the liver (SIR-Spheres; Sirtex, Woburn, MA) followed 2-3 weeks later by the combination of durvalumab and tremelimumab. A Simon two-stage design was implemented, with a planned expansion to 18 patients if at least one response was noted in the first nine patients. RESULTS: Nine patients enrolled in the first stage of the study, all with progressive disease (PD) during or after their first two cycles of treatment. Per preplanned design, the study was closed because of futility. No treatment-related grade 3 or greater toxicities were recorded. Correlative studies with tumor biopsies showed low levels of tumor-infiltrating lymphocyte (TIL) infiltration in tumor cancer islands before and after Y90 radioembolization. CONCLUSION: Y90 radioembolization can be added safely to durvalumab and tremelimumab but did not promote tumor-directed immune responses against liver-metastasized MSS CRC.

Vitamin D as a Primer for Oncolytic Viral Therapy in Colon Cancer Models.

Oncolytic viroimmunotherapy is an exciting modality that can offer lasting anti-tumor immunity for aggressive malignancies like colon cancer. The impact of oncolytic viruses may be extended by combining them with agents to prime a tumor for viral susceptibility. This study investigates vitamin D analogue as an adjunct to oncolytic viral therapy for colon cancer. While vitamin D (VD) has historically been viewed as anti-viral, our in vitro investigations using human colon cancer cell lines showed that VD does not directly inhibit replication of recombinant chimeric poxvirus CF33. VD did restrict growth in HT29 but not HCT116 human colon cancer cells. In vivo investigations using HCT116 and HT29 xenograft models of colon cancer demonstrated that a VD analogue, calcipotriol, was additive with CF33-based viral therapy in VD-responsive HT29 but not in HCT116 tumors. Analyses of RNA-sequencing and gene expression data demonstrated a downregulation in the Jak-STAT signaling pathway with the addition of VD to viral therapy in HT29 models suggesting that the anti-inflammatory properties of VD may enhance the effects of viral therapy in some models. In conclusion, VD may prime oncolytic viral therapy in certain colon cancers.

Association of UDP-galactose-4-epimerase with milk protein concentration in the Chinese Holstein population.

OBJECTIVE: An initial RNA-Sequencing study revealed that UDP-galactose-4-epimerase (GALE) was one of the most promising candidates for milk protein concentration in Chinese Holstein cattle. This enzyme catalyzes the interconversion of UDP-galactose and UDP-glucose, an important step in galactose catabolism. To further validate the genetic effect of GALE on milk protein traits, genetic variations were identified, and genotypes-phenotypes associations were performed. METHODS: The entire coding region and the 5'-regulatory region (5'-UTR) of GALE were re-sequenced using pooled DNA of 17 unrelated sires. Association studies for five milk production traits were performed using a mixed linear animal model with a population encompassing 1,027 Chinese Holstein cows. RESULTS: A total of three variants in GALE were identified, including two novel variants (g.2114 A>G and g.2037 G>A) in the 5'-UTR and one previously reported variant (g.3836 G>C) in an intron. All three single nucleotide polymorphisms (SNPs) were associated with milk yield (p<0.0001), fat yield (p = 0.0006 to <0.0001), protein yield (p = 0.0232 to <0.0001) and protein percentage (p<0.0001), while no significant associations were detected between the SNPs and fat percentage. A strong linkage disequilibrium (D' = 0.96 to 1.00) was observed among all three SNPs, and a 5 Kb haplotype block involving three main haplotypes with GAG, AGC, and AGG was formed. The results of haplotype association analyses were consistent with the results of single locus association analysis (p<0.0001). The phenotypic variance ratio above 3.00% was observed for milk protein yield that was explained by SNP-g.3836G >C. CONCLUSION: Overall, our findings provided new insights into the polymorphic variations in bovine GALE gene and their associations with milk protein concentration. The data indicate their potential uses for marker-assisted breeding or genetic selection schemes.

Phase I clinical trial of the combination of eribulin and everolimus in patients with metastatic triple-negative breast cancer.

BACKGROUND: Alteration of the PI3K/AKT/mTOR pathway is a common genomic abnormality detected in triple-negative breast cancer (TNBC). Everolimus acts synergistically with eribulin in TNBC cell lines and xenograft models. This phase I trial was designed to test the safety and tolerability of combining eribulin and everolimus in patients with metastatic TNBC. METHODS: The primary objective of this study was to evaluate the safety and toxicities of the combination. Patients with metastatic TNBC who had up to four lines of prior chemotherapies were enrolled. The combination of eribulin and everolimus was tested using three dosing levels: A1 (everolimus 5 mg daily; eribulin 1.4 mg/m2 days 1 and 8 every 3 weeks), A2 (everolimus 7.5 mg daily; eribulin 1.4 mg/m2, days 1 and 8 every 3 weeks), and B1 (everolimus 5 mg daily; eribulin 1.1 mg/m2 days 1 and 8 every 3 weeks). RESULTS: Twenty-seven patients with median age 55 years were enrolled. Among 8 evaluable patients who received dose level A1, 4 had dose-limiting toxicities (DLTs). Among 3 evaluable patients treated with dose level A2, 2 had DLTs. Among 12 evaluable patients who received dose level B1, 4 had DLTs. The DLTs were neutropenia, stomatitis, and hyperglycemia. Over the study period, 59% had a ≥ grade 3 toxicity, 44% had ≥ grade 3 hematologic toxicities, and 22% had grade 4 hematologic toxicities. The most common hematological toxicities were neutropenia, leukopenia, and lymphopenia. Thirty-three percent had grade 3 non-hematologic toxicities. The most common non-hematological toxicities were stomatitis, hyperglycemia, and fatigue. The median number of cycles completed was 4 (range 0-8). Among 25 eligible patients, 9 patients (36%) achieved the best response as partial response, 9 (36%) had stable disease, and 7 (28%) had progression. The median time to progression was 2.6 months (95% CI [2.1, 4.0]), and median overall survival (OS) was 8.3 months (95% CI [5.5, undefined]). CONCLUSION: Eribulin 1.1 mg/m2 days 1 and 8 every 3 weeks with everolimus 5 mg daily was defined as the highest dose with acceptable toxicity (RP2D). The combination is safe, and efficacy is modest. A post hoc analysis showed that participants that used dexamethasone mouthwash stayed on treatment for one additional cycle. TRIAL REGISTRATION: ClinicalTrials.gov, NCT02120469. Registered 18 April 2014.

Estimation of the genetic parameters for semen traits in Chinese Holstein bulls.

BACKGROUND: Semen traits are important for the widespread use of superior bulls. Thus, the objective of this study was to estimate the heritability of five semen traits, ejaculate volume (VE), progressive sperm motility (SM), sperm concentration (SC), number of sperm (NSP), and number of progressive motile sperm (NMSP), and their genetic correlations (rg). The dataset being studied consisted of 1450 Chinese Holstein bulls with full pedigree information, born between 1996 and 2012, representing 11 AI centers. Genetic parameters were estimated in a multivariate analysis using the average information restricted maximum likelihood estimation of variance (AI-REML). RESULTS: The estimates of heritability for VE, SM, SC, NSP, and NMSP were 0.15, 0.12, 0.22, 0.16 and 0.12, respectively. The genetic correlations among the five semen traits ranged from 0.02 (VE and SC) to 0.99 (NSP and NMSP). CONCLUSIONS: Our findings provide useful information on the heritability of semen traits in Holstein bulls and the relationships among them, and should assist in selection for improvement of semen traits in Chinese Holstein bulls.

B Cell Lymphoma Immunotherapy Using TLR9-Targeted Oligonucleotide STAT3 Inhibitors.

Growing evidence links the aggressiveness of non-Hodgkin's lymphoma, especially the activated B cell-like type diffuse large B cell lymphomas (ABC-DLBCLs) to Toll-like receptor 9 (TLR9)/MyD88 and STAT3 transcription factor signaling. Here, we describe a dual-function molecule consisting of a clinically relevant TLR9 agonist (CpG7909) and a STAT3 inhibitor in the form of a high-affinity decoy oligodeoxynucleotide (dODN). The CpG-STAT3dODN blocked STAT3 DNA binding and activity, thus reducing expression of downstream target genes, such as MYC and BCL2L1, in human and mouse lymphoma cells. We further demonstrated that injections (i.v.) of CpG-STAT3dODN inhibited growth of human OCI-Ly3 lymphoma in immunodeficient mice. Moreover, systemic CpG-STAT3dODN administration induced complete regression of the syngeneic A20 lymphoma, resulting in long-term survival of immunocompetent mice. Both TLR9 stimulation and concurrent STAT3 inhibition were critical for immune-mediated therapeutic effects, since neither CpG7909 alone nor CpG7909 co-injected with unconjugated STAT3dODN extended mouse survival. The CpG-STAT3dODN induced expression of genes critical to antigen-processing/presentation and Th1 cell activation while suppressing survival signaling. These effects resulted in the generation of lymphoma cell-specific CD8/CD4-dependent T cell immunity protecting mice from tumor rechallenge. Our results suggest that CpG-STAT3dODN as a systemic/local monotherapy or in combination with PD1 blockade can provide an opportunity for treating patients with B cell NHL.

RNA interference screening identifies lenalidomide sensitizers in multiple myeloma, including RSK2.

To identify molecular targets that modify sensitivity to lenalidomide, we measured proliferation in multiple myeloma (MM) cells transfected with 27 968 small interfering RNAs in the presence of increasing concentrations of drug and identified 63 genes that enhance activity of lenalidomide upon silencing. Ribosomal protein S6 kinase (RPS6KA3 or RSK2) was the most potent sensitizer. Other notable gene targets included 5 RAB family members, 3 potassium channel proteins, and 2 peroxisome family members. Single genes of interest included I-κ-B kinase-α (CHUK), and a phosphorylation dependent transcription factor (CREB1), which associate with RSK2 to regulate several signaling pathways. RSK2 knockdown induced cytotoxicity across a panel of MM cell lines and consistently increased sensitivity to lenalidomide. Accordingly, 3 small molecular inhibitors of RSK2 demonstrated synergy with lenalidomide cytotoxicity in MM cells even in the presence of stromal contact. Both RSK2 knockdown and small molecule inhibition downregulate interferon regulatory factor 4 and MYC, and provides an explanation for the synergy between lenalidomide and RSK2 inhibition. Interestingly, RSK2 inhibition also sensitized MM cells to bortezomib, melphalan, and dexamethasone, but did not downregulate Ikaros or influence lenalidomide-mediated downregulation of tumor necrosis factor-α or increase lenalidomide-induced IL-2 upregulation. In summary, inhibition of RSK2 may prove a broadly useful adjunct to MM therapy.

Genetic effects of PDGFRB and MARCH1 identified in GWAS revealing strong associations with semen production traits in Chinese Holstein bulls.

BACKGROUND: Using a genome-wide association study strategy, our previous study discovered 19 significant single-nucleotide polymorphisms (SNPs) related to semen production traits in Chinese Holstein bulls. Among them, three SNPs were within or close to the phosphodiesterase 3A (PDE3A), membrane associated ring-CH-type finger 1 (MARCH1) and platelet derived growth factor receptor beta (PDGFRB) genes. The present study was designed with the objectives of identifying genetic polymorphism of the PDE3A, PDGFRB and MARCH1 genes and their effects on semen production traits in a Holstein bull population. RESULTS: A total of 20 SNPs were detected and genotyped in 730 bulls. Association analyses using de-regressed estimated breeding values of each semen production trait revealed four statistically significant SNPs for one or more semen production traits (P < 0.05): one SNP was located downstream of PDGFRB and three SNPs were located in the promoter of MARCH1. Interestingly, for MARCH1, haplotype-based analysis revealed significant associations of haplotypes with semen volume per ejaculate. Furthermore, high expression of the MARCH1 gene was observed in sperm cells. One SNP (rs43445726) in the regulatory region of MARCH1 had a significant effect on gene expression. CONCLUSION: Our study demonstrated the significant associations of genetic variants of the PDGFRB and MARCH1 genes with semen production traits. The identified SNPs may serve as genetic markers to optimize breeding programs for semen production traits in Holstein bull populations.

SALL4 suppresses PTEN expression to promote glioma cell proliferation via PI3K/AKT signaling pathway.

Spalt-like transcription factor 4 (SALL4), a oncogene, is known to participate in multiple carcinomas, and is up-regulated in glioma. However, its actual role and underlying mechanisms in the development of glioma remain unclear. The present study explored the molecular functions of SALL4 in promoting cell proliferation in glioma. The expression level of SALL4 in 69 human glioma samples and six non-tumor brain tissues was determined using real-time polymerase chain reaction (PCR). Then, we transfected U87 and U251 cell lines with siRNA, and assessed cellular proliferation and cell cycle to understand the function of SALL4, and the relationship between SALL4, PTEN and PI3K/AKT pathway. PCR confirmed that the expression of SALL4 was higher in the glioma samples than non-tumor brain tissues. Cellular growth and proliferation were dramatically reduced following inhibition of SALL4 expression. Western blot showed increase in PTEN expression when SALL4 was silenced, which in turn depressed the activation of PI3K/AKT pathway, suggesting that PTEN was a downstream target of SALL4 in glioma development. Therefore, SALL4 could act as a proto-oncogene by regulating the PTEN/PI3K/AKT signaling pathway, thereby facilitating proliferation of glioma cells.

Genome-wide association study for semen traits of the bulls in Chinese Holstein.

A genome-wide association study (GWAS) was performed to identify markers and candidate genes for five semen traits in the Holstein bull population in China. The analyzed dataset consisted of records from 692 bulls from eight bull stations; each bull was genotyped using the Illumina BovineSNP50 BeadChip. Association tests between each trait and the 41 188 informative high-quality SNPs were achieved with gapit software. In total, 19 suggestive significant SNPs, partly located within the reported QTL regions or within or close to the reported candidate genes, associated with five semen traits were detected. By combining our GWAS results with the biological functions of these genes, eight novel promising candidate genes, including ETNK1, PDE3A, PDGFRB, CSF1R, WT1, DSCAML1, SOD1 and RUNX2, were identified that potentially relate to semen traits. Our findings may provide a basis for further research on the genetic mechanism of semen traits and marker-assisted selection of such traits in Holstein bulls.

Top-hat random fiber Bragg grating.

We examined the possibility of using noise or pseudo-random variations of the refractive index in the design of fiber Bragg gratings (FBGs). We demonstrated theoretically and experimentally that top-hat FBGs may be designed and fabricated using this approach. The reflectivity of the fabricated top-hat FBG matches quite well with that of the designed one.

RNAi screening of the kinome with cytarabine in leukemias.

To identify rational therapeutic combinations with cytarabine (Ara-C), we developed a high-throughput, small-interference RNA (siRNA) platform for myeloid leukemia cells. Of 572 kinases individually silenced in combination with Ara-C, silencing of 10 (1.7%) and 8 (1.4%) kinases strongly increased Ara-C activity in TF-1 and THP-1 cells, respectively. The strongest molecular concepts emerged around kinases involved in cell-cycle checkpoints and DNA-damage repair. In confirmatory siRNA assays, inhibition of WEE1 resulted in more potent and universal sensitization across myeloid cell lines than siRNA inhibition of PKMYT1, CHEK1, or ATR. Treatment of 8 acute myeloid leukemia (AML), acute lymphoblastic leukemia (ALL), and chronic myeloid leukemia (CML) cell lines with commercial and the first-in-class clinical WEE1 kinase inhibitor MK1775 confirmed sensitization to Ara-C up to 97-fold. Ex vivo, adding MK1775 substantially reduced viability in 13 of 14 AML, CML, and myelodysplastic syndrome patient samples compared with Ara-C alone. Maximum sensitization occurred at lower to moderate concentrations of both drugs. Induction of apoptosis was increased using a combination of Ara-C and MK1775 compared with using either drug alone. WEE1 is expressed in primary AML, ALL, and CML specimens. Data from this first siRNA-kinome sensitizer screen suggests that inhibiting WEE1 in combination with Ara-C is a rational combination for the treatment of myeloid and lymphoid leukemias.

Identification of novel variants in BRF1 gene from patient with developmental delay, hearing abnormality, and nervous system anomalies.

Cerebellofaciodental syndrome characterized with dysmorphic features, intellectual disability, and brain anomalies. Now its clinical spectrum expanded more manifestations including bilateral sensorineural hearing impairment and inner ear malformation. Here, we report a 14-month-old boy with global developmental delay and hearing disorder. Whole exome sequencing (WES) revealed the compound heterozygous variants [NM_001519.4: c.652 T > G (p.W218G); c.915 + 1G > T] in the BRF1 gene which inherited from his parents, respectively. The MRI results showed hypoplastic cerebellar vermis, enlarged cisterna magna, and prominent fourth ventricle, the rehabilitation therapy failed to improve the symptoms for our patient. Our finding expands the genetic spectrum of BRF1 variants, which indicates patients with the developmental delay caused by BRF1 variants require other treatments instead of rehabilitation.

Pharmacological profiling of a berbamine derivative for lymphoma treatment.

ABSTRACT: Ca2+/calmodulin-dependent protein kinase II γ (CAMKIIγ) has been identified as a potential target for treating cancer. Based on our previous study of berbamine (BBM) as a CAMKIIγ inhibitor, we have synthesized a new BBM derivative termed PA4. Compared with BBM, PA4 showed improved potency and specificity and was more cytotoxic against lymphoma and leukemia than against other types of cancer. In addition to indirectly targeting c-Myc protein stability, we demonstrated that its cytotoxic effects were also mediated via increased reactive oxygen species production in lymphoma cells. PA4 significantly impeded tumor growth in vivo in a xenograft T-cell lymphoma mouse model. Pharmacokinetics studies demonstrated quick absorption into plasma after oral administration, with a maximum concentration of 1680 ± 479 ng/mL at 5.33 ± 2.31 hours. The calculated oral absolute bioavailability was 34.1%. Toxicity assessment of PA4 showed that the therapeutic window used in our experiments was safe for future development. Given its efficacy, safety, and favorable pharmacokinetic profile, PA4 is a potential lead candidate for treating lymphoma.

Mapping and functional characterization of structural variation in 1060 pig genomes.

BACKGROUND: Structural variations (SVs) have significant impacts on complex phenotypes by rearranging large amounts of DNA sequence. RESULTS: We present a comprehensive SV catalog based on the whole-genome sequence of 1060 pigs (Sus scrofa) representing 101 breeds, covering 9.6% of the pig genome. This catalog includes 42,487 deletions, 37,913 mobile element insertions, 3308 duplications, 1664 inversions, and 45,184 break ends. Estimates of breed ancestry and hybridization using genotyped SVs align well with those from single nucleotide polymorphisms. Geographically stratified deletions are observed, along with known duplications of the KIT gene, responsible for white coat color in European pigs. Additionally, we identify a recent SINE element insertion in MYO5A transcripts of European pigs, potentially influencing alternative splicing patterns and coat color alterations. Furthermore, a Yorkshire-specific copy number gain within ABCG2 is found, impacting chromatin interactions and gene expression across multiple tissues over a stretch of genomic region of ~200 kb. Preliminary investigations into SV's impact on gene expression and traits using the Pig Genotype-Tissue Expression (PigGTEx) data reveal SV associations with regulatory variants and gene-trait pairs. For instance, a 51-bp deletion is linked to the lead eQTL of the lipid metabolism regulating gene FADS3, whose expression in embryo may affect loin muscle area, as revealed by our transcriptome-wide association studies. CONCLUSIONS: This SV catalog serves as a valuable resource for studying diversity, evolutionary history, and functional shaping of the pig genome by processes like domestication, trait-based breeding, and adaptive evolution.

A compendium of genetic regulatory effects across pig tissues.

The Farm Animal Genotype-Tissue Expression (FarmGTEx) project has been established to develop a public resource of genetic regulatory variants in livestock, which is essential for linking genetic polymorphisms to variation in phenotypes, helping fundamental biological discovery and exploitation in animal breeding and human biomedicine. Here we show results from the pilot phase of PigGTEx by processing 5,457 RNA-sequencing and 1,602 whole-genome sequencing samples passing quality control from pigs. We build a pig genotype imputation panel and associate millions of genetic variants with five types of transcriptomic phenotypes in 34 tissues. We evaluate tissue specificity of regulatory effects and elucidate molecular mechanisms of their action using multi-omics data. Leveraging this resource, we decipher regulatory mechanisms underlying 207 pig complex phenotypes and demonstrate the similarity of pigs to humans in gene expression and the genetic regulation behind complex phenotypes, supporting the importance of pigs as a human biomedical model.