Francisco Peñagaricano

Insights into the Maize Pan-Genome and Pan-Transcriptome

Transcriptome sequencing of diverse maize inbreds provided insights into the nature of the maize pan-genome, including identification of 8681 loci absent in the B73 reference sequence. Genome-wide association studies using single nucleotide polymorphisms and transcript abundance variants in the maize pan-genome identified loci associated with traits important for fitness and adaptation. 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.

RNA-Seq analysis uncovers transcriptomic variations between morphologically similar in vivo- and in vitro-derived bovine blastocysts

BackgroundA 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.ResultsA 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.ConclusionsThus, our results support that IVF may influence at the transcriptomic level and that morphology is limited in full characterization of bovine preimplantation embryos.

Maternal Diet during Pregnancy Induces Gene Expression and DNA Methylation Changes in Fetal Tissues in Sheep

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 days 67 ± 3 of gestation until necropsy (days 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 methyltransferase (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.

Effect of Maternal Methionine Supplementation on the Transcriptome of Bovine Preimplantation Embryos

Maternal nutrition exclusively during the periconceptional period can induce remarkable effects on both oocyte maturation and early embryo development, which in turn can have lifelong consequences. The objective of this study was to evaluate the effect of maternal methionine supplementation on the transcriptome of bovine preimplantation embryos. Holstein cows were randomly assigned to one of two treatments differing in level of dietary methionine (1.89 Met vs. 2.43 Met % of metabolizable protein) from calving until embryo flushing. High quality preimplantation embryos from individual cows were pooled and then analyzed by RNA sequencing. Remarkably, a subtle difference in methionine supplementation in maternal diet was sufficient to cause significant changes in the transcriptome of the embryos. A total of 276 genes out of 10,662 showed differential expression between treatments (FDR <0.10). Interestingly, several of the most significant genes are related to embryonic development (e.g., VIM, IFI6, BCL2A1, and TBX15) and immune response (e.g., NKG7, TYROBP, SLAMF7, LCP1, and BLA-DQB). Likewise, gene set enrichment analysis revealed that several Gene Ontology terms, InterPro entries, and KEGG pathways were enriched (FDR <0.05) with differentially expressed genes involved in embryo development and immune system. The expression of most genes was decreased by maternal methionine supplementation, consistent with reduced transcription of genes with increased methylation of specific genes by increased methionine. Overall, our findings provide evidence that supplementing methionine to dams prior to conception and during the preimplantation period can modulate gene expression in bovine blastocysts. The ramifications of the observed gene expression changes for subsequent development of the pregnancy and physiology of the offspring warrant further investigation in future studies.

Genome-wide association study identifies candidate markers for bull fertility in Holstein dairy cattle

The decline in the reproductive efficiency of dairy cattle has become a challenging problem worldwide. Female fertility is now taken into account in breeding goals while generally less attention is given to male fertility. The objective of this study was to perform a genome-wide association study in Holstein bulls to identify genetic variants significantly related to sire conception rate (SCR), a new phenotypic evaluation of bull fertility. The analysis included 1755 sires with SCR data and 38,650 single nucleotide polymorphisms (SNPs) spanning the entire bovine genome. Associations between SNPs and SCR were analyzed using a mixed linear model that included a random polygenic effect and SNP genotype either as a linear covariate or as a categorical variable. A multiple testing correction approach was used to account for the correlation between SNPs because of linkage disequilibrium. After genome-wide correction, eight SNPs showed significant association with SCR. Some of these SNPs are located close to or in the middle of genes with functions related to male fertility, such as the sperm acrosome reaction, chromatin remodeling during the spermatogenesis, and the meiotic process during male germ cell maturation. Some SNPs showed marked dominance effects, which provide more evidence for the relevance of non-additive effects in traits closely related to fitness such as fertility. The results could contribute to the identification of genes and pathways associated with male fertility in dairy cattle.

Inferring quantitative trait pathways associated with bull fertility from a genome-wide association study.

Whole-genome association studies typically focus on genetic markers with the strongest evidence of association. However, single markers often explain only a small component of the genetic variance and hence offer a limited understanding of the trait under study. As such, the objective of this study was to perform a pathway-based association analysis in Holstein dairy cattle in order to identify relevant pathways involved in bull fertility. The results of a single-marker association analysis, using 1,755 bulls with sire conception rate data and genotypes for 38,650 single nucleotide polymorphisms (SNPs), were used in this study. A total of 16,819 annotated genes, including 2,767 significantly associated with bull fertility, were used to interrogate a total of 662 Gene Ontology (GO) terms and 248 InterPro (IP) entries using a test of proportions based on the cumulative hypergeometric distribution. After multiple-testing correction, 20 GO categories and one IP entry showed significant overrepresentation of genes statistically associated with bull fertility. Several of these functional categories such as small GTPases mediated signal transduction, neurogenesis, calcium ion binding, and cytoskeleton are known to be involved in biological processes closely related to male fertility. These results could provide insight into the genetic architecture of this complex trait in dairy cattle. In addition, this study shows that quantitative trait pathways inferred from single-marker analyses could enhance our interpretations of the results of genome-wide association studies.

Association between milk protein gene variants and protein composition traits in dairy cattle

The objective of this study was to identify DNA markers in the 4 casein genes (CSN1S1, CSN1S2, CSN2, and CSN3) and the 2 major whey protein genes (LALBA and LGB) that show associations with milk protein profile measured by reverse-phase HPLC. Fifty-three single nucleotide polymorphisms (SNP) were genotyped for cows in a unique resource population consisting of purebred Holstein and (Holstein × Jersey) × Holstein crossbred animals. Seven traits were analyzed, including concentrations of α(S)-casein (CN), β-CN, κ-CN, α-lactalbumin, β-lactoglobulin, and 2 additional secondary traits, the total concentration of the above 5 milk proteins and the α(S)-CN to β-CN ratio. A substantial fraction of phenotypic variation could be explained by the additive genetic component for the 7 milk protein composition traits studied. Moreover, several SNP were significantly associated with all examined traits at an experiment-wise error rate of 0.05, except for α-lactalbumin. Importantly, the significant SNP explained a large proportion of the phenotypic variation of milk protein composition. Our findings could be used for selecting animals that produce milk with desired composition or desired processing and manufacturing properties.

Maternal nutrition induces gene expression changes in fetal muscle and adipose tissues in sheep

BackgroundMaternal nutrition during different stages of pregnancy can induce significant changes in the structure, physiology, and metabolism of the offspring. These changes could have important implications on food animal production especially if these perturbations impact muscle and adipose tissue development. Here, we evaluated the impact of different maternal isoenergetic diets, alfalfa haylage (HY; fiber), corn (CN; starch), and dried corn distillers grains (DG; fiber plus protein plus fat), on the transcriptome of fetal muscle and adipose tissues in sheep.ResultsPrepartum diets were associated with notable gene expression changes in fetal tissues. In longissimus dorsi muscle, a total of 224 and 823 genes showed differential expression (FDR ≤0.05) in fetuses derived from DG vs. CN and HY vs. CN maternal diets, respectively. Several of these significant genes affected myogenesis and muscle differentiation. In subcutaneous and perirenal adipose tissues, 745 and 208 genes were differentially expressed (FDR ≤0.05), respectively, between CN and DG diets. Many of these genes are involved in adipogenesis, lipogenesis, and adipose tissue development. Pathway analysis revealed that several GO terms and KEGG pathways were enriched (FDR ≤0.05) with differentially expressed genes associated with tissue and organ development, chromatin biology, and different metabolic processes.ConclusionsThese findings provide evidence that maternal nutrition during pregnancy can alter the programming of fetal muscle and fat tissues in sheep. The ramifications of the observed gene expression changes, in terms of postnatal growth, body composition, and meat quality of the offspring, warrant future investigation.

Short communication: A missense mutation in the PROP1 (prophet of Pit 1) gene affects male fertility and milk production traits in the US Holstein population

In previous studies, we reported significant associations of the POU1F1 pathway genes with reproduction and production traits in several dairy cattle populations. Polymorphisms in genes of this pathway were found to be associated with both female and male fertility traits in dairy cattle. The POU1F1 gene is a direct downstream target for the regulation of the prophet of Pit1 (PROP1) gene, also known as PROP paired-like homeobox 1. Interestingly, the position of PROP1 coincides with a quantitative trait locus affecting ovulation rate in cattle. Therefore, the objective of this study was to investigate whether PROP1 affects fertility and milk production traits in Holstein cattle. Using the DNA pooling sequencing approach, a missense single nucleotide polymorphism that replaces a histidine amino acid with an arginine was detected in exon 3 of PROP1. The arginine allele was found to be associated with a decrease in sire conception rate and an increase in productive life, protein yield, and net merit index in a population of 1,951 Holstein bulls. The transcription factors produced from the histidine and arginine isoforms are known to have different transcription, DNA binding, and regulation activities. As such, we propose that the association of the arginine isoform with decreased bull fertility is likely caused by reduced activity of this allele in male functions. The findings of this study suggest PROP1 polymorphisms as candidates in selection programs for fertility, health, and milk production traits in dairy cattle.

An application of MeSH enrichment analysis in livestock

Summary An integral part of functional genomics studies is to assess the enrichment of specific biological terms in lists of genes found to be playing an important role in biological phenomena. Contrasting the observed frequency of annotated terms with those of the background is at the core of overrepresentation analysis (ORA). Gene Ontology (GO) is a means to consistently classify and annotate gene products and has become a mainstay in ORA. Alternatively, Medical Subject Headings (MeSH) offers a comprehensive life science vocabulary including additional categories that are not covered by GO. Although MeSH is applied predominantly in human and model organism research, its full potential in livestock genetics is yet to be explored. In this study, MeSH ORA was evaluated to discern biological properties of identified genes and contrast them with the results obtained from GO enrichment analysis. Three published datasets were employed for this purpose, representing a gene expression study in dairy cattle, the use of SNPs for genome‐wide prediction in swine and the identification of genomic regions targeted by selection in horses. We found that several overrepresented MeSH annotations linked to these gene sets share similar concepts with those of GO terms. Moreover, MeSH yielded unique annotations, which are not directly provided by GO terms, suggesting that MeSH has the potential to refine and enrich the representation of biological knowledge. We demonstrated that MeSH can be regarded as another choice of annotation to draw biological inferences from genes identified via experimental analyses. When used in combination with GO terms, our results indicate that MeSH can enhance our functional interpretations for specific biological conditions or the genetic basis of complex traits in livestock species.