Current challenges and future of agricultural genomes to phenomes in the USA.

Dramatic improvements in measuring genetic variation across agriculturally relevant populations (genomics) must be matched by improvements in identifying and measuring relevant trait variation in such populations across many environments (phenomics). Identifying the most critical opportunities and challenges in genome to phenome (G2P) research is the focus of this paper. Previously (Genome Biol, 23(1):1-11, 2022), we laid out how Agricultural Genome to Phenome Initiative (AG2PI) will coordinate activities with USA federal government agencies expand public-private partnerships, and engage with external stakeholders to achieve a shared vision of future the AG2PI. Acting on this latter step, AG2PI organized the "Thinking Big: Visualizing the Future of AG2PI" two-day workshop held September 9-10, 2022, in Ames, Iowa, co-hosted with the United State Department of Agriculture's National Institute of Food and Agriculture (USDA NIFA). During the meeting, attendees were asked to use their experience and curiosity to review the current status of agricultural genome to phenome (AG2P) work and envision the future of the AG2P field. The topic summaries composing this paper are distilled from two 1.5-h small group discussions. Challenges and solutions identified across multiple topics at the workshop were explored. We end our discussion with a vision for the future of agricultural progress, identifying two areas of innovation needed: (1) innovate in genetic improvement methods development and evaluation and (2) innovate in agricultural research processes to solve societal problems. To address these needs, we then provide six specific goals that we recommend be implemented immediately in support of advancing AG2P research.

Heritability estimates and predictive ability for pig meat quality traits using identity-by-state and identity-by-descent relationships in an F2 population.

Genomic relationships can be computed with dense genome-wide genotypes through different methods, either based on identity-by-state (IBS) or identity-by-descent (IBD). The latter has been shown to increase the accuracy of both estimated relationships and predicted breeding values. However, it is not clear whether an IBD approach would achieve greater heritability ( h 2 ) and predictive ability ( r ̂ y , y ̂ ) than its IBS counterpart for data with low-depth pedigrees. Here, we compare both approaches in terms of the estimated of h 2 and r ̂ y , y ̂ , using data on meat quality and carcass traits recorded in experimental crossbred pigs, with a pedigree constrained to only three generations. Three animal models were fitted which differed on the relationship matrix: an IBS model ( G IBS ), an IBD (defined within the known pedigree) model ( G IBD ), and a pedigree model ( A 22 ). In 9 of 20 traits, the range of increase for the estimates of σ u 2 and h 2 was 1.2-2.9 times greater with G IBS and G IBD models than with A 22 . Whereas for all traits, both parameters were similar between genomic models. The r ̂ y , y ̂ of the genomic models was higher compared to A 22 . A scarce increment in r ̂ y , y ̂ was found with G IBS when compared to G IBD , most likely due to the former recovering sizeable relationships among founder F0 animals.

Altered Larval Yellow Perch Swimming Behavior Due to Methylmercury and PCB126 Detected Using Hidden Markov Chain Models.

Fish swimming behavior is a commonly measured response in aquatic ecotoxicology because behavior is considered a whole organism-level effect that integrates many sensory systems. Recent advancements in animal behavior models, such as hidden Markov chain models (HMM), suggest an improved analytical approach for toxicology. Using both new and traditional approaches, we examined the sublethal effects of PCB126 and methylmercury on yellow perch (YP) larvae (Perca flavescens) using three doses. Both approaches indicate larvae increase activity after exposure to either chemical. The middle methylmercury-dosed larvae showed multiple altered behavior patterns. First, larvae had a general increase in activity, typically performing more behavior states, more time swimming, and more swimming bouts per second. Second, when larvae were in a slow or medium swimming state, these larvae tended to switch between these states more often. Third, larvae swam slower during the swimming bouts. The upper PCB126-dosed larvae exhibited a higher proportion and a fast swimming state, but the total time spent swimming fast decreased. The middle PCB126-dosed larvae transitioned from fast to slow swimming states less often than the control larvae. These results indicate that developmental exposure to very low doses of these neurotoxicants alters YP larvae overall swimming behaviors, suggesting neurodevelopment alteration.

The Social Life of Pigs: Changes in Affiliative and Agonistic Behaviors following Mixing.

This study investigated potentially affiliative behaviors in grow-finish pigs, how these behaviors changed over time and their relationship to agonistic behaviors. A total of 257 Yorkshire barrows were observed for agonistic (reciprocal fights, attacks) and affiliative (nosing, play, non-agonistic contact) behaviors after mixing (at 10 weeks of age), and weeks 3, 6, and 9 after mix. The least square means of affiliative behaviors were compared across time points. Relationships among affiliative and agonistic behaviors were assessed using generalized linear mixed models. Non-agonistic contact with conspecifics increased until week 6 then remained stable between weeks 6 and 9. Nosing was highest at mix, then decreased in the following weeks. Play was lowest at mix and highest at week 3. Affiliative behaviors were negatively related with aggression at mix (p < 0.001). Pigs who engaged in play and nosing behaviors were more likely to be involved in agonistic interactions in the weeks after mixing (p < 0.05), while pigs engaging in non-agonistic contact were less likely to be involved in agonistic interactions (p < 0.001). There appear to be relationships between affiliative and agonistic behaviors in pigs, with contact being the most predictive of less aggression. Future studies could focus on promoting positive non-agonistic contact in unfamiliar pigs as a way to mitigate aggressive interactions.

Safety and efficacy of topically administered netarsudil-latanoprost fixed dose combination (FDC; Rocklatan™) in normal and glaucomatous dogs with ADAMTS10-open-angle glaucoma (ADAMTS10-OAG).

OBJECTIVE: The aim of the study was to evaluate safety and efficacy of topically administered 0.02% netarsudil-0.005% latanoprost fixed-dose combination (FDC) (Rocklatan™; Aerie Pharmaceutical) in normal and glaucomatous dogs with ADAMTS10-open-angle glaucoma (ADAMTS10-OAG). ANIMALS STUDIED: Five normal and five glaucomatous beagle dogs with ADAMTS10-OAG were the study animals. PROCEDURES: In each dog, left (OS) or right eye (OD) was randomly selected for netarsudil-latanoprost FDC treatment. Contralateral eyes served as latanoprost-treated controls. The study was divided into four consecutive study periods: following a 4-day baseline period, two sequential 8-day study periods followed with once daily (q24h) and twice daily (q12h) treatments and ending with a washout period. Efficacy was measured by diurnal intraocular pressure (IOP) and pupil diameter. Safety was assessed by routine ophthalmic examination, gonioscopy, and pachymetry. Differences in least square means of quantitative outcome measures were compared between FDC and latanoprost treatments by using the linear Gaussian model. RESULTS: Baseline IOPs were 13.6 ± 0.7 mmHg (mean ± SEM) in normal and 28.3 ± 1.4 mmHg in OAG-affected dogs. There was a significant decrease in mean diurnal IOP following FDC administration in both normal (q24h: -2.1 mmHg; q12h: -4.1 mmHg) and glaucomatous dogs (q24h: -14.2 mmHg; q12h: -17.7 mmHg; p < .0001). There was no significant difference in the treatment effect when comparing FDC to latanoprost. Both FDC and latanoprost administration resulted in similarly significant pupil constriction (p < .0001). The FDC administration was well-tolerated but resulted in conjunctival hyperemia. CONCLUSIONS: Once or twice daily administration of netarsudil-latanoprost FDC (Rocklatan™) and latanoprost was equally effective in lowering IOP in normal and OAG-affected dogs. There was no netarsudil-related added treatment effect.

Heuristic hyperparameter optimization of deep learning models for genomic prediction.

There is a growing interest among quantitative geneticists and animal breeders in the use of deep learning (DL) for genomic prediction. However, the performance of DL is affected by hyperparameters that are typically manually set by users. These hyperparameters do not simply specify the architecture of the model; they are also critical for the efficacy of the optimization and model-fitting process. To date, most DL approaches used for genomic prediction have concentrated on identifying suitable hyperparameters by exploring discrete options from a subset of the hyperparameter space. Enlarging the hyperparameter optimization search space with continuous hyperparameters is a daunting combinatorial problem. To deal with this problem, we propose using differential evolution (DE) to perform an efficient search of arbitrarily complex hyperparameter spaces in DL models, and we apply this to the specific case of genomic prediction of livestock phenotypes. This approach was evaluated on two pig and cattle datasets with real genotypes and simulated phenotypes (N = 7,539 animals and M = 48,541 markers) and one real dataset (N = 910 individuals and M = 28,916 markers). Hyperparameters were evaluated using cross-validation. We compared the predictive performance of DL models using hyperparameters optimized by DE against DL models with "best practice" hyperparameters selected from published studies and baseline DL models with randomly specified hyperparameters. Optimized models using DE showed a clear improvement in predictive performance across all three datasets. DE optimized hyperparameters also resulted in DL models with less overfitting and less variation in predictive performance over repeated retraining compared to non-optimized DL models.

Integrated Genome-Wide Analysis of MicroRNA Expression Quantitative Trait Loci in Pig Longissimus Dorsi Muscle.

Determining mechanisms regulating complex traits in pigs is essential to improve the production efficiency of this globally important protein source. MicroRNAs (miRNAs) are a class of non-coding RNAs known to post-transcriptionally regulate gene expression affecting numerous phenotypes, including those important to the pig industry. To facilitate a more comprehensive understanding of the regulatory mechanisms controlling growth, carcass composition, and meat quality phenotypes in pigs, we integrated miRNA and gene expression data from longissimus dorsi muscle samples with genotypic and phenotypic data from the same animals. We identified 23 miRNA expression Quantitative Trait Loci (miR-eQTL) at the genome-wide level and examined their potential effects on these important production phenotypes through miRNA target prediction, correlation, and colocalization analyses. One miR-eQTL miRNA, miR-874, has target genes that colocalize with phenotypic QTL for 12 production traits across the genome including backfat thickness, dressing percentage, muscle pH at 24 h post-mortem, and cook yield. The results of our study reveal genomic regions underlying variation in miRNA expression and identify miRNAs and genes for future validation of their regulatory effects on traits of economic importance to the global pig industry.

Time budgets of group-housed pigs in relation to social aggression and production.

Commercial producers house growing pigs by sex and weight to allow for efficient use of resources and provide pigs the welfare benefits of interacting with their conspecifics and more freedom of movement. However, the introduction of unfamiliar pigs can cause increased aggression for 24 to 48 h as pigs establish social relationships. To address this issue, a better understanding of pig behavior is needed. The objectives of this study were to quantify time budgets of pigs following introduction into a new social group and how these changed over time and to investigate how social aggression influences the overall time budgets and production parameters. A total of 257 grow-finish Yorkshire barrows across 20 pens were introduced into new social groups at 10 wk of age (~23 kg) and observed for aggression and time budgets of behavior at four periods: immediately after introduction and 3, 6, and 9 wk later. Pigs were observed for the duration of total aggression and initiated aggression (s) for 9 h after introduction and for 4 h at 3, 6, and 9 wk later. Time budgets were created by scan sampling inactive, movement, ingestion, social, and exploration behaviors every 2 min for 4 h in the afternoon and summarizing the proportion of time each behavior was performed by period. The least square means of each behavior were compared across time points. Pigs spent most of their time inactive. In general, the greatest change in pig behavior was observed between introduction and week 3 (P < 0.003), with gradual changes throughout the study period as pigs became more inactive (week 3 vs. week 6: P = 0.209; week 6 vs. week 9: P = 0.007) and spent less time on other behaviors. Pigs' nonaggressive behavior and production parameters were compared with aggression using generalized linear mixed models. The time pigs spent on nonaggressive behaviors was negatively related to aggression (P < 0.045) with few exceptions. Initiated aggression after introduction was negatively related to loin muscle area (P = 0.003). These results show how finishing pigs spend their time in commercial facilities and indicate that behavior continues to change for up to 9 wk after introduction into a new social group. Efforts to reduce chronic levels of aggression should focus on promoting nonaggressive behaviors, such as exploration and movement, after the initial fighting that occurs immediately after introduction has waned, and should be implemented for up to 9 wk after introduction into new social groups.

Phenomes: the current frontier in animal breeding.

Improvements in genomic technologies have outpaced the most optimistic predictions, allowing industry-scale application of genomic selection. However, only marginal gains in genetic prediction accuracy can now be expected by increasing marker density up to sequence, unless causative mutations are identified. We argue that some of the most scientifically disrupting and industry-relevant challenges relate to 'phenomics' instead of 'genomics'. Thanks to developments in sensor technology and artificial intelligence, there is a wide range of analytical tools that are already available and many more will be developed. We can now address some of the pressing societal demands on the industry, such as animal welfare concerns or efficiency in the use of resources. From the statistical and computational point of view, phenomics raises two important issues that require further work: penalization and dimension reduction. This will be complicated by the inherent heterogeneity and 'missingness' of the data. Overall, we can expect that precision livestock technologies will make it possible to collect hundreds of traits on a continuous basis from large numbers of animals. Perhaps the main revolution will come from redesigning animal breeding schemes to explicitly allow for high-dimensional phenomics. In the meantime, phenomics data will definitely enlighten our knowledge on the biological basis of phenotypes.

Safety and efficacy of topically administered netarsudil (Rhopressa™) in normal and glaucomatous dogs with ADAMTS10-open-angle glaucoma (ADAMTS10-OAG).

OBJECTIVE: To evaluate safety and efficacy of topically administered 0.02% netarsudil ophthalmic solution (Rhopressa™; Aerie Pharmaceutical) in normal and glaucomatous dogs with ADAMTS10-open-angle glaucoma (ADAMTS10-OAG). ANIMALS STUDIED: Five normal and 5 glaucomatous Beagle dogs with ADAMTS10-OAG. PROCEDURES: In each dog, left or right eye was randomly selected for netarsudil treatment. Contralateral eyes were sham-treated with balanced salt solution (BSS). Following a 1-week baseline period, dogs were treated once daily (q24h) during week 2, and twice daily (q12h) during week 3; week 4 served as washout period. Efficacy was measured by diurnal intraocular pressure (IOP) and pupil diameter. Safety was assessed by routine ophthalmic examination, gonioscopy, and pachymetry. Differences in least square means of quantitative outcome measures were compared between netarsudil and BSS sham-treated eyes by linear Gaussian model. RESULTS: Baseline IOPs were 18.5 ± 0.5 mm Hg (mean ± SEM) in normal and 27.8 ± 1.0 mm Hg in OAG dogs. Even though mean IOPs were lower in netarsudil- vs sham-treated eyes, the overall differences were neither significant nor clinically relevant, regardless of treatment frequency (q24h-normal: sham 16.4 ± 1.1 mm Hg vs treatment 15.6 ± 1.0 mm Hg; q24hr-OAG: sham 25.8 ± 2.3 mm Hg vs. treatment 25.7 ± 2.4 mm Hg; q12hr-normal: sham 15.4 ± 0.8 mm Hg vs. treatment 14.4 ± 0.8 mm Hg; q12hr-OAG: sham 26.3 ± 1.7 mm Hg vs. treatment 25.4 ± 1.8 mm Hg). Netarsudil administration was well tolerated but resulted in significant, moderate-to-severe conjunctival hyperemia (P < .001). CONCLUSIONS: Once or twice daily administration of netarsudil resulted in marginal and clinically irrelevant IOP decreases in normal and OAG-affected dogs. Except for conjunctival hyperemia, the drug was well tolerated.

Deciphering Sex-Specific Genetic Architectures Using Local Bayesian Regressions.

Many complex human traits exhibit differences between sexes. While numerous factors likely contribute to this phenomenon, growing evidence from genome-wide studies suggest a partial explanation: that males and females from the same population possess differing genetic architectures. Despite this, mapping gene-by-sex (G×S) interactions remains a challenge likely because the magnitude of such an interaction is typically and exceedingly small; traditional genome-wide association techniques may be underpowered to detect such events, due partly to the burden of multiple test correction. Here, we developed a local Bayesian regression (LBR) method to estimate sex-specific SNP marker effects after fully accounting for local linkage-disequilibrium (LD) patterns. This enabled us to infer sex-specific effects and G×S interactions either at the single SNP level, or by aggregating the effects of multiple SNPs to make inferences at the level of small LD-based regions. Using simulations in which there was imperfect LD between SNPs and causal variants, we showed that aggregating sex-specific marker effects with LBR provides improved power and resolution to detect G×S interactions over traditional single-SNP-based tests. When using LBR to analyze traits from the UK Biobank, we detected a relatively large G×S interaction impacting bone mineral density within ABO, and replicated many previously detected large-magnitude G×S interactions impacting waist-to-hip ratio. We also discovered many new G×S interactions impacting such traits as height and body mass index (BMI) within regions of the genome where both male- and female-specific effects explain a small proportion of phenotypic variance (R2 < 1 × 10-4), but are enriched in known expression quantitative trait loci.

Dysregulated FXR-FGF19 signaling and choline metabolism are associated with gut dysbiosis and hyperplasia in a novel pig model of pediatric NASH.

To investigate the role of bile acids (BAs) in the pathogenesis of diet-induced nonalcoholic steatohepatitis (NASH), we fed a "Western-style diet" [high fructose, high fat (HFF)] enriched with fructose, cholesterol, and saturated fat for 10 wk to juvenile Iberian pigs. We also supplemented probiotics with in vitro BA deconjugating activity to evaluate their potential therapeutic effect in NASH. Liver lipid and function, cytokines, and hormones were analyzed using commercially available kits. Metabolites, BAs, and fatty acids were measured by liquid chromatography-mass spectrometry. Histology and gene and protein expression analyses were performed using standard protocols. HFF-fed pigs developed NASH, cholestasis, and impaired enterohepatic Farnesoid-X receptor (FXR)-fibroblast growth factor 19 (FGF19) signaling in the absence of obesity and insulin resistance. Choline depletion in HFF livers was associated with decreased lipoprotein and cholesterol in serum and an increase of choline-containing phospholipids in colon contents and trimethylamine-N-oxide in the liver. Additionally, gut dysbiosis and hyperplasia increased with the severity of NASH, and were correlated with increased colonic levels of choline metabolites and secondary BAs. Supplementation of probiotics in the HFF diet enhanced NASH, inhibited hepatic autophagy, increased excretion of taurine and choline, and decreased gut microbial diversity. In conclusion, dysregulation of BA homeostasis was associated with injury and choline depletion in the liver, as well as increased biliary secretion, gut metabolism and excretion of choline-based phospholipids. Choline depletion limited lipoprotein synthesis, resulting in hepatic steatosis, whereas secondary BAs and choline-containing phospholipids in colon may have promoted dysbiosis, hyperplasia, and trimethylamine synthesis, causing further damage to the liver.NEW & NOTEWORTHY Impaired Farnesoid-X receptor (FXR)-fibroblast growth factor 19 (FGF19) signaling and cholestasis has been described in nonalcoholic fatty liver disease (NAFLD) patients. However, therapeutic interventions with FXR agonists have produced contradictory results. In a swine model of pediatric nonalcoholic steatohepatitis (NASH), we show that the uncoupling of intestinal FXR-FGF19 signaling and a decrease in FGF19 levels are associated with a choline-deficient phenotype of NASH and increased choline excretion in the gut, with the subsequent dysbiosis, colonic hyperplasia, and accumulation of trimethylamine-N-oxide in the liver.

Estimation of indirect social genetic effects for skin lesion count in group-housed pigs by quantifying behavioral interactions1.

Mixing of pigs into new social groups commonly induces aggressive interactions that result in skin lesions on the body of the animals. The relationship between skin lesions and aggressive behavioral interactions in group-housed pigs can be analyzed within the framework of social genetic effects (SGE). This study incorporates the quantification of aggressive interactions between pairs of animals in the modeling of SGE for skin lesions in different regions of the body in growing pigs. The dataset included 792 pigs housed in 59 pens. Skin lesions in the anterior, central, and caudal regions of the body were counted 24 h after pig mixing. Animals were video-recorded for 9 h postmixing and trained observers recorded the type and duration of aggressive interactions between pairs of animals. The number of seconds that pairs of pigs spent engaged in reciprocal fights and unilateral attack behaviors were used to parametrize the intensity of social interactions (ISI). Three types of models were fitted: direct genetic additive model (DGE), traditional social genetic effect model (TSGE) assuming uniform interactions between dyads, and an intensity-based social genetic effect model (ISGE) that used ISI to parameterize SGE. All models included fixed effects of sex, replicate, lesion scorer, weight at mixing, premixing lesion count, and the total time that the animal spent engaged in aggressive interactions (reciprocal fights and unilateral attack behaviors) as a covariate; a random effect of pen; and a random direct genetic effect. The ISGE models recovered more direct genetic variance than DGE and TSGE, and the estimated heritabilities (h^D2) were highest for all traits (P < 0.01) for the ISGE with ISI parametrized with unilateral attack behavior. The TSGE produced estimates that did not differ significantly from DGE (P > 0.5). Incorporating the ISI into ISGE, even in a small dataset, allowed separate estimation of the genetic parameters for direct and SGE, as well as the genetic correlation between direct and SGE (rs), which was positive for all lesion traits. The estimates from ISGE suggest that if behavioral observations are available, selection incorporating SGE may reduce the consequences of aggressive behaviors after mixing pigs.

Genetic control of longissimus dorsi muscle gene expression variation and joint analysis with phenotypic quantitative trait loci in pigs.

BACKGROUND: Economically important growth and meat quality traits in pigs are controlled by cascading molecular events occurring during development and continuing throughout the conversion of muscle to meat. However, little is known about the genes and molecular mechanisms involved in this process. Evaluating transcriptomic profiles of skeletal muscle during the initial steps leading to the conversion of muscle to meat can identify key regulators of polygenic phenotypes. In addition, mapping transcript abundance through genome-wide association analysis using high-density marker genotypes allows identification of genomic regions that control gene expression, referred to as expression quantitative trait loci (eQTL). In this study, we perform eQTL analyses to identify potential candidate genes and molecular markers regulating growth and meat quality traits in pigs. RESULTS: Messenger RNA transcripts obtained with RNA-seq of longissimus dorsi muscle from 168 F2 animals from a Duroc x Pietrain pig resource population were used to estimate gene expression variation subject to genetic control by mapping eQTL. A total of 339 eQTL were mapped (FDR ≤ 0.01) with 191 exhibiting local-acting regulation. Joint analysis of eQTL with phenotypic QTL (pQTL) segregating in our population revealed 16 genes significantly associated with 21 pQTL for meat quality, carcass composition and growth traits. Ten of these pQTL were for meat quality phenotypes that co-localized with one eQTL on SSC2 (8.8-Mb region) and 11 eQTL on SSC15 (121-Mb region). Biological processes identified for co-localized eQTL genes include calcium signaling (FERM, MRLN, PKP2 and CHRNA9), energy metabolism (SUCLG2 and PFKFB3) and redox hemostasis (NQO1 and CEP128), and results support an important role for activation of the PI3K-Akt-mTOR signaling pathway during the initial conversion of muscle to meat. CONCLUSION: Co-localization of eQTL with pQTL identified molecular markers significantly associated with both economically important phenotypes and gene transcript abundance. This study reveals candidate genes contributing to variation in pig production traits, and provides new knowledge regarding the genetic architecture of meat quality phenotypes.

Induction of maxillary and mandibular squamous epithelial cell proliferation in mink (Neovison vison) by ß-naphthoflavone.

A jaw lesion reported in mink exposed to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and TCDD-like chemicals is considered a potential indicator of exposure to these chemicals. Many of the effects of TCDD-like chemicals are induced through interaction with the aryl hydrocarbon receptor. The present study indicates that mink dosed with ß-naphthoflavone, which is an aryl hydrocarbon receptor ligand but not a TCDD-like chemical, also develop the lesion. Environ Toxicol Chem 2019;38:460-463. © 2018 SETAC.

Whole-Genome Multi-omic Study of Survival in Patients with Glioblastoma Multiforme.

Glioblastoma multiforme (GBM) has been recognized as the most lethal type of malignant brain tumor. Despite efforts of the medical and research community, patients' survival remains extremely low. Multi-omic profiles (including DNA sequence, methylation and gene expression) provide rich information about the tumor. These profiles are likely to reveal processes that may be predictive of patient survival. However, the integration of multi-omic profiles, which are high dimensional and heterogeneous in nature, poses great challenges. The goal of this work was to develop models for prediction of survival of GBM patients that can integrate clinical information and multi-omic profiles, using multi-layered Bayesian regressions. We apply the methodology to data from GBM patients from The Cancer Genome Atlas (TCGA, n = 501) to evaluate whether integrating multi-omic profiles (SNP-genotypes, methylation, copy number variants and gene expression) with clinical information (demographics as well as treatments) leads to an improved ability to predict patient survival. The proposed Bayesian models were used to estimate the proportion of variance explained by clinical covariates and omics and to evaluate prediction accuracy in cross validation (using the area under the Receiver Operating Characteristic curve, AUC). Among clinical and demographic covariates, age (AUC = 0.664) and the use of temozolomide (AUC = 0.606) were the most predictive of survival. Among omics, methylation (AUC = 0.623) and gene expression (AUC = 0.593) were more predictive than either SNP (AUC = 0.539) or CNV (AUC = 0.547). While there was a clear association between age and methylation, the integration of age, the use of temozolomide, and either gene expression or methylation led to a substantial increase in AUC in cross-validaton (AUC = 0.718). Finally, among the genes whose methylation was higher in aging brains, we observed a higher enrichment of these genes being also differentially methylated in cancer.

Sex-Specific Modulation of Fetal Adipogenesis by Gestational Bisphenol A and Bisphenol S Exposure.

The endocrine-disrupting chemical bisphenol A (BPA) increases adipose tissue mass in vivo and promotes adipogenesis in vitro; however, mechanisms explaining BPA's obesogenic effect remain unknown. We investigated the effects of gestational BPA and its analog, bisphenol S (BPS), exposure on the adipogenic differentiation ability of fetal preadipocytes and the role of endoplasmic reticulum stress in regulating this process. Pregnant sheep (n = 7 to 8 per group) mated to the same male were exposed to BPA or BPS from days 30 to 100 of gestation; pregnancies were terminated 20 days later. Adipose tissue was harvested and fetal preadipocytes isolated. Adipose tissue gene expression, adipocyte size, preadipocyte gene expression, adipogenic differentiation, and dynamic expression of genes involved in adipogenesis and endoplasmic reticulum stress were assessed. Gestational BPA enhanced adipogenic differentiation in female, but not male, preadipocytes. The unfolded protein response (UPR) pathway was upregulated in BPA-exposed female preadipocytes supportive of a higher endoplasmic reticulum stress. Increased expression of estradiol receptor 1 and glucocorticoid receptor in female preadipocytes suggests that this may be a potential cause behind the sex-specific effects observed upon BPA exposure. Gestational BPS affected adipogenic terminal differentiation gene expression in male preadipocytes, but not adipogenic differentiation potential. We demonstrate that gestational BPA exposure can modulate the differentiation ability of fetal preadipocytes. UPR upregulation in gestationally BPA-exposed female preadipocytes may contribute to the increased preadipocyte's adipogenic ability. The marked sex-specific effect of BPA highlights higher susceptibility of females to bisphenol A and potentially, a higher risk to develop obesity in adulthood.

Profiling and characterization of a longissimus dorsi muscle microRNA dataset from an F2 Duroc × Pietrain pig resource population.

To elucidate the effects of microRNA (miRNA) regulation in skeletal muscle of adult pigs, miRNA expression profiling was performed with RNA extracted from longissimus dorsi (LD) muscle samples from 174 F2 pigs (~ 5.5 months of age) from a Duroc × Pietrain resource population. Total RNA was extracted from LD samples, and libraries were sequenced on an Illumina HiSeq 2500 platform in 1 × 50 bp format. After processing, 232,826,977 total reads were aligned to the Sus scrofa reference genome (v10.2.79), with 74.8% of total reads mapping successfully. The miRDeep2 software package was utilized to quantify annotated Sus scrofa mature miRNAs from miRBase (Release 21) and to predict candidate novel miRNA precursors. Among the retained 295 normalized mature miRNA expression profiles ssc-miR-1, ssc-miR-133a-3p, ssc-miR-378, ssc-miR-206, and ssc-miR-10b were the most abundant, all of which have previously been shown to be expressed in pig skeletal muscle. Additionally, 27 unique candidate novel miRNA precursors were identified exhibiting homologous sequence to annotated human miRNAs. The composition of classes of small RNA present in this dataset was also characterized; while the majority of unique expressed sequence tags were not annotated in any of the queried databases, the most abundantly expressed class of small RNA in this dataset was miRNAs. This data provides a resource to evaluate miRNA regulation of gene expression and effects on complex trait phenotypes in adult pig skeletal muscle. The raw sequencing data were deposited in the Sequence Read Archive, BioProject PRJNA363073.

Genome-Wide Association Analyses Based on Broadly Different Specifications for Prior Distributions, Genomic Windows, and Estimation Methods.

A currently popular strategy (EMMAX) for genome-wide association (GWA) analysis infers association for the specific marker of interest by treating its effect as fixed while treating all other marker effects as classical Gaussian random effects. It may be more statistically coherent to specify all markers as sharing the same prior distribution, whether that distribution is Gaussian, heavy-tailed (BayesA), or has variable selection specifications based on a mixture of, say, two Gaussian distributions [stochastic search and variable selection (SSVS)]. Furthermore, all such GWA inference should be formally based on posterior probabilities or test statistics as we present here, rather than merely being based on point estimates. We compared these three broad categories of priors within a simulation study to investigate the effects of different degrees of skewness for quantitative trait loci (QTL) effects and numbers of QTL using 43,266 SNP marker genotypes from 922 Duroc-Pietrain F2-cross pigs. Genomic regions were based either on single SNP associations, on nonoverlapping windows of various fixed sizes (0.5-3 Mb), or on adaptively determined windows that cluster the genome into blocks based on linkage disequilibrium. We found that SSVS and BayesA lead to the best receiver operating curve properties in almost all cases. We also evaluated approximate maximum a posteriori (MAP) approaches to BayesA and SSVS as potential computationally feasible alternatives; however, MAP inferences were not promising, particularly due to their sensitivity to starting values. We determined that it is advantageous to use variable selection specifications based on adaptively constructed genomic window lengths for GWA studies.