Genomes at the species level are dynamic, with genes present in every individual (core) and genes in a subset of individuals (dispensable) that collectively constitute the pan-genome. Using transcriptome sequencing of seedling RNA from 503 maize (Zea mays) inbred lines to characterize the maize pan-genome, we identified 8681 representative transcript assemblies (RTAs) with 16.4% expressed in all lines and 82.7% expressed in subsets of the lines. Interestingly, with linkage disequilibrium mapping, 76.7% of the RTAs with at least one single nucleotide polymorphism (SNP) could be mapped to a single genetic position, distributed primarily throughout the nonpericentromeric portion of the genome. Stepwise iterative clustering of RTAs suggests, within the context of the genotypes used in this study, that the maize genome is restricted and further sampling of seedling RNA within this germplasm base will result in minimal discovery. Genome-wide association studies based on SNPs and transcript abundance in the pan-genome revealed loci associated with the timing of the juvenile-to-adult vegetative and vegetative-to-reproductive developmental transitions, two traits important for fitness and adaptation. This study revealed the dynamic nature of the maize pan-genome and demonstrated that a substantial portion of variation may lie outside the single reference genome for a species.

Background: Transforming large amounts of genomic data into valuable knowledge for predicting complex traits has been an important challenge for animal and plant breeders. Prediction of complex traits has not escaped the current excitement on machine-learning, including interest in deep learning algorithms such as multilayer perceptrons (MLP) and convolutional neural networks (CNN). The aim of this study was to compare the predictive performance of two deep learning methods (MLP and CNN), two ensemble learning methods [random forests (RF) and gradient boosting (GB)], and two parametric methods [genomic best linear unbiased prediction (GBLUP) and Bayes B] using real and simulated datasets. Methods: The real dataset consisted of 11,790 Holstein bulls with sire conception rate (SCR) records and genotyped for 58k single nucleotide polymorphisms (SNPs). To support the evaluation of deep learning methods, various simulation studies were conducted using the observed genotype data as template, assuming a heritability of 0.30 with either additive or non-additive gene effects, and two different numbers of quantitative trait nucleotides (100 and 1000). Results: In the bull dataset, the best predictive correlation was obtained with GB (0.36), followed by Bayes B (0.34), GBLUP (0.33), RF (0.32), CNN (0.29) and MLP (0.26). The same trend was observed when using mean squared error of prediction. The simulation indicated that when gene action was purely additive, parametric methods outperformed other methods. When the gene action was a combination of additive, dominance and of two-locus epistasis, the best predictive ability was obtained with gradient boosting, and the superiority of deep learning over the parametric methods depended on the number of loci controlling the trait and on sample size. In fact, with a large dataset including 80k individuals, the predictive performance of deep learning methods was similar or slightly better than that of parametric methods for traits with non-additive gene action. Conclusions: For prediction of traits with non-additive gene action, gradient boosting was a robust method. Deep learning approaches were not better for genomic prediction unless non-additive variance was sizable.

Two retrospective studies examining data of 7,500 lactating cows from a single herd were performed with the objective of evaluating the long-term effects of clinical disease during the early postpartum period on milk production, reproduction, and culling of dairy cows through 305 days in milk (DIM). In the first study, data regarding health, milk production, reproduction, and culling of 5,085 cows were summarized. Cows were classified according to incidence of clinical problem (metritis, mastitis, lameness, digestive problem, or respiratory problem) during the first 21 DIM (ClinD21). During 305 d of lactation, cows that had ClinD21 produced, on average, 410 kg less milk, 17 kg less fat, and 12 kg less protein compared with cows that did not have ClinD21 (NoClinD21). Although the interval to first breeding was not different between groups of interest, pregnancy rate through 305 DIM was lower in cows that had ClinD21 [adjusted hazard ratio (AHR) = 0.81]. When individual breedings were analyzed, cows that had ClinD21 presented lower rates of pregnancy per breeding for breedings performed before 150 DIM, reduced numbers of calving per breeding for breedings performed before 200 DIM, and greater number of pregnancy losses for all breedings performed through 305 DIM. The rate of culling from calving through 305 DIM was higher in cows that had a single ClinD21 (AHR = 1.79) and in cows that had multiple ClinD21 (AHR = 3.06), which resulted in a greater proportion of cows leaving the herd by 305 DIM (NoClinD21 = 22.6%; single ClinD21 = 35.7%; multiple ClinD21 = 53.8%). In the second study, data regarding postpartum health and 305-d yields of milk, fat, and protein were collected from 2,415 primiparous cows that had genomic testing information. Genomic estimated breeding values (EBV) were used to predict 305-d yields of milk, fat, and protein. Genomic EBV and predicted yields of milk, fat, and protein did not differ between cows that had ClinD21 and those that did not have ClinD21. In contrast, the observed 305-d yields of milk, fat, and protein were reduced by 345, 10, and 10 kg, respectively, in cows that had ClinD21 compared with cows that did not have ClinD21. We conclude that clinical disease diagnosed and treated during the first 21 DIM has long-term effects on lactation performance, reproduction, and culling of dairy cows, which contribute to detrimental consequences of health problems on sustainability of dairy herds. Replication of our studies in multiple herds will be important to confirm our findings in a larger population.

Background: A valuable tool for both research and industry, in vitro fertilization (IVF) has applications range from gamete selection and preservation of traits to cloning. Although IVF has achieved worldwide use, with approximately 339,685 bovine embryos transferred in 2010 alone, there are still continuing difficulties with efficiency. It is rare to have more than 40% of fertilized in vitro cattle oocytes reach blastocyst stage by day 8 of culture, and pregnancy rates are reported as less than 45% for in vitro produced embryos. To investigate potential influences in-vitro fertilization (IVF) has on embryonic development, this study compares in vivo- and in vitro-derived bovine blastocysts at a similar stage and quality grade (expanded, excellent quality) to determine the degree of transcriptomic variation beyond morphology using RNA-Seq.Results: A total of 26,906,451 and 38,184,547 fragments were sequenced for in vitro and in vivo embryo pools, respectively. We detected expression for a total of 17,634 genes, with 793 genes showing differential expression between the two embryo populations with false discovery rate (FDR) < 0.05. There were also 395 novel transcribed units found, of which 45 were differentially expressed (FDR < 0.05). In addition, 4,800 genes showed evidence of alternative splicing, with 873 genes displaying differential alternative splicing between the two pools (FDR < 0.05). Using GO enrichment analysis, multiple biological pathways were found to be significantly enriched for differentially expressed genes (FDR < 0.01), including cholesterol and sterol synthesis, system development, and cell differentiation.Conclusions: Thus, our results support that IVF may influence at the transcriptomic level and that morphology is limited in full characterization of bovine preimplantation embryos.

Exposure to intrauterine heat stress during late gestation affects offspring performance into adulthood. However, underlying mechanistic links between thermal insult in fetal life and postnatal outcomes are not completely understood. We examined morphology, DNA methylation, and gene expression of liver and mammary gland for bull calves and heifers that were gestated under maternal conditions of heat stress or cooling (i.e. in utero heat stressedvs. in utero cooled calves). Mammary tissue was harvested from dairy heifers during their first lactation and liver from bull calves at birth. The liver of in utero heat stressed bull calves contained more cells and the mammary glands of in utero heat stressed heifers were comprised of smaller alveoli. We identified more than 1,500 CpG sites differently methylated between maternal treatment groups. These CpGs were associated with approximately 400 genes, which play a role in processes, such asdevelopment,innate immune defense, cell signaling, andtranscription and translation. We also identified over 100 differentially expressed genes in the mammary gland with similar functions. Interestingly, fifty differentially methylated genes were shared by both bull calf liver and heifer mammary gland. Intrauterine heat stress alters the methylation profile of liver and mammary DNA and programs their morphology in postnatal life, which may contribute to the poorer performance of in utero heat stressed calves.

Studies in rats and mice have established that maternal nutrition induces epigenetic modifications, sometimes permanently, that alter gene expression in the fetus, which in turn leads to phenotypic changes. However, limited data is available on the influence of maternal diet on epigenetic modifications and gene expression in sheep. Therefore, the objectives of this study were to investigate the impact of different maternal dietary energy sources on the expression of imprinted genes in fetuses in sheep. Ewes were naturally bred to a single sire and from d 67 ± 3 of gestation until necropsy (d 130 ± 1), they were fed one of three diets of alfalfa haylage (HY; fiber), corn (CN; starch), or dried corn distiller’s grains (DG; fiber plus protein plus fat). A total of 26 fetuses were removed from the dams and longissimus dorsi, semitendinosus, perirenal adipose depot, and subcutaneous adipose depot tissues were collected for expression and DNA methylation analyses. Expression analysis of nine imprinted genes and three DNA methylatransferase (DNMTs) genes showed significant effects of the different maternal diets on the expression of these genes. The methylation levels of CpG islands of both IGF2R and H19 were higher in HY and DG than CN fetuses in both males and females. This result is consistent with the low amino acid content of the CN diet, a source of methyl group donors, compared to HY and DG diets. Thus, results of this study provide evidence of association between maternal nutrition during pregnancy and transcriptomic and epigenomic alterations of imprinted genes and DNMTs in the fetal tissues.