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.

Influence of Maternal BLV Infection on miRNA and tRF Expression in Calves.

Small non-coding RNAs, such as microRNAs (miRNA) and tRNA-derived fragments (tRF), are known to be involved in post-transcriptional gene regulation. Research has provided evidence that small RNAs may influence immune development in calves. Bovine leukosis is a disease in cattle caused by Bovine Leukemia Virus (BLV) that leads to increased susceptibility to opportunistic pathogens. No research has addressed the potential influence that a maternal BLV infection may have on gene regulation through the differential expression of miRNAs or tRFs in progeny. Blood samples from 14-day old Holstein calves born to BLV-infected dams were collected. Antibodies for BLV were assessed using ELISA and levels of BLV provirus were assessed using qPCR. Total RNA was extracted from whole blood samples for small RNA sequencing. Five miRNAs (bta-miR-1, bta-miR-206, bta-miR-133a, bta-miR-133b, and bta-miR-2450d) and five tRFs (tRF-36-8JZ8RN58X2NF79E, tRF-20-0PF05B2I, tRF-27-W4R951KHZKK, tRF-22-S3M8309NF, and tRF-26-M87SFR2W9J0) were dysregulated in calves born to BLV-infected dams. The miRNAs appear to be involved in the gene regulation of immunological responses and muscle development. The tRF subtypes and parental tRNA profiles in calves born to infected dams appear to be consistent with previous publications in adult cattle with BLV infection. These findings offer insight into how maternal BLV infection status may impact the development of offspring.

AgAnimalGenomes: browsers for viewing and manually annotating farm animal genomes.

Current genome sequencing technologies have made it possible to generate highly contiguous genome assemblies for non-model animal species. Despite advances in genome assembly methods, there is still room for improvement in the delineation of specific gene features in the genomes. Here we present genome visualization and annotation tools to support seven livestock species (bovine, chicken, goat, horse, pig, sheep, and water buffalo), available in a new resource called AgAnimalGenomes. In addition to supporting the manual refinement of gene models, these browsers provide visualization tracks for hundreds of RNAseq experiments, as well as data generated by the Functional Annotation of Animal Genomes (FAANG) Consortium. For species with predicted gene sets from both Ensembl and RefSeq, the browsers provide special tracks showing the thousands of protein-coding genes that disagree across the two gene sources, serving as a valuable resource to alert researchers to gene model issues that may affect data interpretation. We describe the data and search methods available in the new genome browsers and how to use the provided tools to edit and create new gene models.

Characteristics of circulating small noncoding RNAs in plasma and serum during human aging.

Objective: Aging is a complicated process that triggers age-related disease susceptibility through intercellular communication in the microenvironment. While the classic secretome of senescence-associated secretory phenotype (SASP) including soluble factors, growth factors, and extracellular matrix remodeling enzymes are known to impact tissue homeostasis during the aging process, the effects of novel SASP components, extracellular small noncoding RNAs (sncRNAs), on human aging are not well established. Methods: Here, by utilizing 446 small RNA-seq samples from plasma and serum of healthy donors found in the Extracellular RNA (exRNA) Atlas data repository, we correlated linear and nonlinear features between circulating sncRNAs expression and age by the maximal information coefficient (MIC) relationship determination. Age predictors were generated by ensemble machine learning methods (Adaptive Boosting, Gradient Boosting, and Random Forest) and core age-related sncRNAs were determined through weighted coefficients in machine learning models. Functional investigation was performed via target prediction of age-related miRNAs. Results: We observed the number of highly expressed transfer RNAs (tRNAs) and microRNAs (miRNAs) showed positive and negative associations with age respectively. Two-variable (sncRNA expression and individual age) relationships were detected by MIC and sncRNAs-based age predictors were established, resulting in a forecast performance where all R 2 values were greater than 0.96 and root-mean-square errors (RMSE) were less than 3.7 years in three ensemble machine learning methods. Furthermore, important age-related sncRNAs were identified based on modeling and the biological pathways of age-related miRNAs were characterized by their predicted targets, including multiple pathways in intercellular communication, cancer and immune regulation. Conclusion: In summary, this study provides valuable insights into circulating sncRNAs expression dynamics during human aging and may lead to advanced understanding of age-related sncRNAs functions with further elucidation.

Implications of tRNA abundance on translation elongation across bovine tissues.

Introduction: Translation is a crucial stage of gene expression. It may also act as an additional layer of regulation that plays an important role in gene expression and function. Highly expressed genes are believed to be codon-biased to support increased protein production, in which quickly translated codons correspond to highly abundant tRNAs. Synonymous SNPs, considered to be silent due to the degeneracy of the genetic code, may shift protein abundance and function through alterations in translational efficiency and suboptimal pairing to lowly abundant tRNAs. Methods: Here, we applied Quantitative Mature tRNA sequencing (QuantM-tRNAseq) and ribosome profiling across bovine tissues in order to investigate the relationship between tRNA expression and slowed translation. Results: Moreover, we have identified genes modulated at transcriptional and/or translational levels underlying tissue-specific biological processes. We have also successfully defined pausing sites that depict the regulatory information encoded within the open reading frame of transcripts, which could be related to translation rate and facilitate proper protein folding. This work offers an atlas of distinctive pausing sites across three bovine tissues, which provides an opportunity to predict codon optimality and understand tissue-specific mechanisms of regulating protein synthesis.

The counterpart congenital overgrowth syndromes Beckwith-Wiedemann Syndrome in human and large offspring syndrome in bovine involve alterations in DNA methylation, transcription, and chromatin configuration.

Beckwith-Wiedemann Syndrome (BWS, OMIM #130650) is a congenital epigenetic disorder in humans which affects approximately 1 in 10,340 children. The incidence is likely an underestimation as the condition is usually recognized based on observable phenotypes at birth. BWS children have up to a 28% risk of developing tumors and currently, only 80% of patients can be corroborated molecularly (epimutations/variants). It is unknown how the subtypes of this condition are molecularly similar/dissimilar globally, therefore there is a need to deeply characterize the syndrome at the molecular level. Here we characterize the methylome, transcriptome and chromatin configuration of 18 BWS individuals together with the animal model of the condition, the bovine large offspring syndrome (LOS). Sex specific comparisons are performed for a subset of the BWS patients and LOS. Given that this epigenetic overgrowth syndrome has been characterized as a loss-of-imprinting condition, parental allele-specific comparisons were performed using the bovine animal model. In general, the differentially methylated regions (DMRs) detected in BWS and LOS showed significant enrichment for CTCF binding sites. Altered chromosome compartments in BWS and LOS were positively correlated with gene expression changes, and the promoters of differentially expressed genes showed significant enrichment for DMRs, differential topologically associating domains, and differential A/B compartments in some comparisons of BWS subtypes and LOS. We show shared regions of dysregulation between BWS and LOS, including several HOX gene clusters, and also demonstrate that altered DNA methylation differs between the clinically epigenetically identified BWS patients and those identified as having DNA variants (i.e. CDKN1C microdeletion). Lastly, we highlight additional genes and genomic regions that have the potential to serve as targets for biomarker development to improve current molecular methodologies. In summary, our results suggest that genome-wide alternation of chromosome architecture, which is partially caused by DNA methylation changes, also contribute to the development of BWS and LOS.

Differentially expressed tRNA-derived fragments in bovine fetuses with assisted reproduction induced congenital overgrowth syndrome.

Background: As couples struggle with infertility and livestock producers wish to rapidly improve genetic merit in their herd, assisted reproductive technologies (ART) have become increasingly popular in human medicine as well as the livestock industry. Utilizing ART can cause an increased risk of congenital overgrowth syndromes, such as Large Offspring Syndrome (LOS) in ruminants. A dysregulation of transcripts has been observed in bovine fetuses with LOS, which is suggested to be a cause of the phenotype. Our recent study identified variations in tRNA expression in LOS individuals, leading us to hypothesize that variations in tRNA expression can influence the availability of their processed regulatory products, tRNA-derived fragments (tRFs). Due to their resemblance in size to microRNAs, studies suggest that tRFs target mRNA transcripts and regulate gene expression. Thus, we have sequenced small RNA isolated from skeletal muscle and liver of day 105 bovine fetuses to elucidate the mechanisms contributing to LOS. Moreover, we have utilized our previously generated tRNA sequencing data to analyze the contribution of tRNA availability to tRF abundance. Results: 22,289 and 7,737 unique tRFs were predicted in the liver and muscle tissue respectively. The greatest number of reads originated from 5' tRFs in muscle and 5' halves in liver. In addition, mitochondrial (MT) and nuclear derived tRF expression was tissue-specific with most MT-tRFs and nuclear tRFs derived from LysUUU and iMetCAU in muscle, and AsnGUU and GlyGCC in liver. Despite variation in tRF abundance within treatment groups, we identified differentially expressed (DE) tRFs across Control-AI, ART-Normal, and ART-LOS groups with the most DE tRFs between ART-Normal and ART-LOS groups. Many DE tRFs target transcripts enriched in pathways related to growth and development in the muscle and tumor development in the liver. Finally, we found positive correlation coefficients between tRNA availability and tRF expression in muscle (R = 0.47) and liver (0.6). Conclusion: Our results highlight the dysregulation of tRF expression and its regulatory roles in LOS. These tRFs were found to target both imprinted and non-imprinted genes in muscle as well as genes linked to tumor development in the liver. Furthermore, we found that tRNA transcription is a highly modulated event that plays a part in the biogenesis of tRFs. This study is the first to investigate the relationship between tRNA and tRF expression in combination with ART-induced LOS.

SARS CoV-2 (Delta Variant) Infection Kinetics and Immunopathogenesis in Domestic Cats.

Continued emergence of SARS-CoV-2 variants highlights the critical need for adaptable and translational animal models for acute COVID-19. Limitations to current animal models for SARS CoV-2 (e.g., transgenic mice, non-human primates, ferrets) include subclinical to mild lower respiratory disease, divergence from clinical COVID-19 disease course, and/or the need for host genetic modifications to permit infection. We therefore established a feline model to study COVID-19 disease progression and utilized this model to evaluate infection kinetics and immunopathology of the rapidly circulating Delta variant (B.1.617.2) of SARS-CoV-2. In this study, specific-pathogen-free domestic cats (n = 24) were inoculated intranasally and/or intratracheally with SARS CoV-2 (B.1.617.2). Infected cats developed severe clinical respiratory disease and pulmonary lesions at 4- and 12-days post-infection (dpi), even at 1/10 the dose of previously studied wild-type SARS-CoV-2. Infectious virus was isolated from nasal secretions of delta-variant infected cats in high amounts at multiple timepoints, and viral antigen was co-localized in ACE2-expressing cells of the lungs (pneumocytes, vascular endothelium, peribronchial glandular epithelium) and strongly associated with severe pulmonary inflammation and vasculitis that were more pronounced than in wild-type SARS-CoV-2 infection. RNA sequencing of infected feline lung tissues identified upregulation of multiple gene pathways associated with cytokine receptor interactions, chemokine signaling, and viral protein-cytokine interactions during acute infection with SARS-CoV-2. Weighted correlation network analysis (WGCNA) of differentially expressed genes identified several distinct clusters of dysregulated hub genes that are significantly correlated with both clinical signs and lesions during acute infection. Collectively, the results of these studies help to delineate the role of domestic cats in disease transmission and response to variant emergence, establish a flexible translational model to develop strategies to prevent the spread of SARS-CoV-2, and identify potential targets for downstream therapeutic development.

RNA-seq reveals insights into molecular mechanisms of metabolic restoration via tryptophan supplementation in low birth weight piglet model.

Low birth weight (LBW) is associated with metabolic disorders in early life. While dietary l-tryptophan (Trp) can ameliorate postprandial plasma triglycerides (TG) disposal in LBW piglets, the genetic and biological basis underlying Trp-caused alterations in lipid metabolism is poorly understood. In this study, we collected 24 liver samples from 1-mo-old LBW and normal birth weight (NBW) piglets supplemented with different concentrations of dietary Trp (NBW with 0% Trp, N0; LBW with 0% Trp, L0; LBW with 0.4% Trp, L4; LBW with 0.8% Trp, L8; N = 6 in each group.) and conducted systematic, transcriptome-wide analysis using RNA sequencing (RNA-seq). We identified 39 differentially expressed genes (DEG) between N0 and L0, and genes within "increased dose effect" clusters based on dose-series expression profile analysis, enriched in fatty acid response of gene ontology (GO) biological process (BP). We then identified RNA-binding proteins including SRSF1, DAZAP1, PUM2, PCBP3, IGF2BP2, and IGF2BP3 significantly (P < 0.05) enriched in alternative splicing events (ASE) in comparison with L0 as control. There were significant positive and negative relationships between candidate genes from co-expression networks (including PID1, ANKRD44, RUSC1, and CYP2J34) and postprandial plasma TG concentration. Further, we determined whether these candidate hub genes were also significantly associated with metabolic and cardiovascular traits in humans via human phenome-wide association study (Phe-WAS), and analysis of mammalian orthologs suggests a functional conservation between human and pig. Our work demonstrates that transcriptomic changes during dietary Trp supplementation in LBW piglets. We detected candidate genes and related BP that may play roles on lipid metabolism restoration. These findings will help to better understand the amino acid support in LBW metabolic complications.

Low birth weight (LBW) has been associated with higher rate of mortality and morbidity and the development of metabolic complications, leaving burdens on livestock production and human health care. The feasibility of LBW metabolic restoration via postnatal nutrition compensation has been verified and the role of one of essential amino acids, l-tryptophan (Trp), on rescuing lipid metabolism in LBW was determined, while the underlying molecular mechanism and key gene regulation is little known. Our study was conducted to identify the unique molecular mechanisms between LBW and normal birth weight (NBW), and to identify the metabolic restoration related genes and biological processes after dietary Trp supplementation in LBW piglet model. We found that differentially expressed genes (DEG) between LBW and NBW were related to fatty acid response based on gene ontology enrichment analysis, and LBW piglets supplemented with Trp showed lower postprandial plasma triglycerides (TG) level as NBW, with similar expression feature of lipid metabolism related genes.

Allele-specific aberration of imprinted domain chromosome architecture associates with large offspring syndrome.

Large offspring syndrome (LOS) and Beckwith-Wiedemann syndrome are similar epigenetic congenital overgrowth conditions in ruminants and humans, respectively. We have reported global loss-of-imprinting, methylome epimutations, and gene misregulation in LOS. However, less than 4% of gene misregulation can be explained with short range (<20kb) alterations in DNA methylation. Therefore, we hypothesized that methylome epimutations in LOS affect chromosome architecture which results in misregulation of genes located at distances >20kb in cis and in trans (other chromosomes). Our analyses focused on two imprinted domains that frequently reveal misregulation in these syndromes, namely KvDMR1 and IGF2R. Using bovine fetal fibroblasts, we identified CTCF binding at IGF2R imprinting control region but not KvDMR1, and allele-specific chromosome architecture of these domains in controls. In LOS, analyses identified erroneous long-range contacts and clustering tendency in the direction of expression of misregulated genes. In conclusion, altered chromosome architecture is associated with LOS.

Characterization of tRNA expression profiles in large offspring syndrome.

BACKGROUND: Assisted Reproductive Technologies (ART) use can increase the risk of congenital overgrowth syndromes, such as large offspring syndrome (LOS) in ruminants. Epigenetic variations are known to influence gene expression and differentially methylated regions (DMRs) were previously determined to be associated with LOS in cattle. We observed DMRs overlapping tRNA clusters which could affect tRNA abundance and be associated with tissue specificity or overgrowth. Variations in tRNA expression have been identified in several disease pathways suggesting an important role in the regulation of biological processes. Understanding the role of tRNA expression in cattle offers an opportunity to reveal mechanisms of regulation at the translational level. We analyzed tRNA expression in the skeletal muscle and liver tissues of day 105 artificial insemination-conceived, ART-conceived with a normal body weight, and ART-conceived bovine fetuses with a body weight above the 97th percentile compared to Control-AI. RESULTS: Despite the centrality of tRNAs to translation, in silico predictions have revealed dramatic differences in the number of tRNA genes between humans and cattle (597 vs 1,659). Consistent with reports in human, only a fraction of predicted tRNA genes are expressed. We detected the expression of 474 and 487 bovine tRNA genes in the muscle and liver with the remainder being unexpressed. 193 and 198 unique tRNA sequences were expressed in all treatment groups within muscle and liver respectively. In addition, an average of 193 tRNA sequences were expressed within the same treatment group in different tissues. Some tRNA isodecoders were differentially expressed between treatment groups. In the skeletal muscle and liver, we categorized 11 tRNA isoacceptors with undetected expression as well as an isodecoder that was unexpressed in the liver (SerGGA). Our results identified variation in the proportion of tRNA gene copies expressed between tissues and differences in the highest contributing tRNA anticodon within an amino acid family due to treatment and tissue type. Out of all amino acid families, roughly half of the most highly expressed tRNA isoacceptors correlated to their most frequent codon in the bovine genome. CONCLUSION: Although the number of bovine tRNA genes is nearly triple of that of the tRNA genes in human, there is a shared occurrence of transcriptionally inactive tRNA genes in both species. We detected differential expression of tRNA genes as well as tissue- and treatment- specific tRNA transcripts with unique sequence variations that could modulate translation during protein homeostasis or cellular stress, and give rise to regulatory products targeting genes related to overgrowth in the skeletal muscle and/or tumor development in the liver of LOS individuals. While the absence of certain isodecoders may be relieved by wobble base pairing, missing tRNA species could increase the likelihood of mistranslation or mRNA degradation.

Effect of food source availability in the salivary gland transcriptome of the unique burying beetle Nicrophorus pustulatus (Coleoptera: Silphidae).

Nicrophorus is a genus of beetles that bury and transform small vertebrate carcasses into a brood ball coated with their oral and anal secretions to prevent decay and that will serve as a food source for their young. Nicrophorus pustulatus is an unusual species with the ability to overtake brood of other burying beetles and whose secretions, unlike other Nicrophorus species, has been reported not to exhibit antimicrobial properties. This work aims to better understand how the presence or absence of a food source influences the expression of genes involved in the feeding process of N. pustulatus. To achieve that, total RNA was extracted from pooled samples of salivary gland tissue from N. pustulatus and sequenced using an Illumina platform. The resulting reads were used to assemble a de novo transcriptome using Trinity. Duplicates with more than 95% similarity were removed to obtain a "unigene" set. Annotation of the unigene set was done using the Trinotate pipeline. Transcript abundance was determined using Kallisto and differential gene expression analysis was performed using edgeR. A total of 651 genes were found to be differentially expressed, including 390 upregulated and 261 downregulated genes in fed insects compared to starved. Several genes upregulated in fed beetles are associated with the insect immune response and detoxification processes with only one transcript encoding for the antimicrobial peptide (AMP) defensin. These results confirm that N. pustulatus does not upregulate the production of genes encoding AMPs during feeding. This study provides a snapshot of the changes in gene expression in the salivary glands of N. pustulatus following feeding while providing a well described transcriptome for the further analysis of this unique burying beetle.

The genomic basis of evolutionary differentiation among honey bees.

In contrast to the western honey bee, Apis mellifera, other honey bee species have been largely neglected despite their importance and diversity. The genetic basis of the evolutionary diversification of honey bees remains largely unknown. Here, we provide a genome-wide comparison of three honey bee species, each representing one of the three subgenera of honey bees, namely the dwarf (Apis florea), giant (A. dorsata), and cavity-nesting (A. mellifera) honey bees with bumblebees as an outgroup. Our analyses resolve the phylogeny of honey bees with the dwarf honey bees diverging first. We find that evolution of increased eusocial complexity in Apis proceeds via increases in the complexity of gene regulation, which is in agreement with previous studies. However, this process seems to be related to pathways other than transcriptional control. Positive selection patterns across Apis reveal a trade-off between maintaining genome stability and generating genetic diversity, with a rapidly evolving piRNA pathway leading to genomes depleted of transposable elements, and a rapidly evolving DNA repair pathway associated with high recombination rates in all Apis species. Diversification within Apis is accompanied by positive selection in several genes whose putative functions present candidate mechanisms for lineage-specific adaptations, such as migration, immunity, and nesting behavior.

De novo assembly of the cattle reference genome with single-molecule sequencing.

BACKGROUND: Major advances in selection progress for cattle have been made following the introduction of genomic tools over the past 10-12 years. These tools depend upon the Bos taurus reference genome (UMD3.1.1), which was created using now-outdated technologies and is hindered by a variety of deficiencies and inaccuracies. RESULTS: We present the new reference genome for cattle, ARS-UCD1.2, based on the same animal as the original to facilitate transfer and interpretation of results obtained from the earlier version, but applying a combination of modern technologies in a de novo assembly to increase continuity, accuracy, and completeness. The assembly includes 2.7 Gb and is >250× more continuous than the original assembly, with contig N50 >25 Mb and L50 of 32. We also greatly expanded supporting RNA-based data for annotation that identifies 30,396 total genes (21,039 protein coding). The new reference assembly is accessible in annotated form for public use. CONCLUSIONS: We demonstrate that improved continuity of assembled sequence warrants the adoption of ARS-UCD1.2 as the new cattle reference genome and that increased assembly accuracy will benefit future research on this species.

Bovine Genome Database: new annotation tools for a new reference genome.

The Bovine Genome Database (BGD) (http://bovinegenome.org) has been the key community bovine genomics database for more than a decade. To accommodate the increasing amount and complexity of bovine genomics data, BGD continues to advance its practices in data acquisition, curation, integration and efficient data retrieval. BGD provides tools for genome browsing (JBrowse), genome annotation (Apollo), data mining (BovineMine) and sequence database searching (BLAST). To augment the BGD genome annotation capabilities, we have developed a new Apollo plug-in, called the Locus-Specific Alternate Assembly (LSAA) tool, which enables users to identify and report potential genome assembly errors and structural variants. BGD now hosts both the newest bovine reference genome assembly, ARS-UCD1.2, as well as the previous reference genome, UMD3.1.1, with cross-genome navigation and queries supported in JBrowse and BovineMine, respectively. Other notable enhancements to BovineMine include the incorporation of genomes and gene annotation datasets for non-bovine ruminant species (goat and sheep), support for multiple assemblies per organism in the Regions Search tool, integration of additional ontologies and development of many new template queries. To better serve the research community, we continue to focus on improving existing tools, developing new tools, adding new datasets and encouraging researchers to use these resources.

Bovine Genome Database: Tools for Mining the Bos taurus Genome.

The Bovine Genome Database (BGD; http://bovinegenome.org ) is a web-accessible resource that supports bovine genomics research by providing genome annotation and data mining tools. BovineMine is a tool within BGD that integrates BGD data, including the genome, genes, precomputed gene expression levels and variant consequences, with external data sources that include quantitative trait loci (QTL), orthologues, Gene Ontology, gene interactions, and pathways. BovineMine enables researchers without programming skills to create custom integrated datasets for use in downstream analyses. This chapter describes how to enhance a bovine genomics project using the Bovine Genome Database, with data mining examples demonstrating BovineMine.

Hymenoptera Genome Database: Using HymenopteraMine to Enhance Genomic Studies of Hymenopteran Insects.

The Hymenoptera Genome Database (HGD; http://hymenopteragenome.org ) is a genome informatics resource for insects of the order Hymenoptera, which includes bees, ants and wasps. HGD provides genome browsers with manual annotation tools (JBrowse/Apollo), BLAST, bulk data download, and a data mining warehouse (HymenopteraMine). This chapter focuses on the use of HymenopteraMine to create annotation data sets that can be exported for use in downstream analyses. HymenopteraMine leverages the InterMine platform to combine genome assemblies and official gene sets with data from OrthoDB, RefSeq, FlyBase, Gene Ontology, UniProt, InterPro, KEGG, Reactome, dbSNP, PubMed, and BioGrid, as well as precomputed gene expression information based on publicly available RNAseq. Built-in template queries provide starting points for data exploration, while the QueryBuilder tool supports construction of complex custom queries. The List Analysis and Genomic Regions search tools execute queries based on uploaded lists of identifiers and genome coordinates, respectively. HymenopteraMine facilitates cross-species data mining based on orthology and supports meta-analyses by tracking identifiers across gene sets and genome assemblies.

Progesterone effects on extracellular vesicles in the sheep uterus.

Progesterone (P4) acts via the endometrium to promote conceptus growth and implantation for pregnancy establishment. Many cells release extracellular vesicles (EVs) that are membrane-bound vesicles of endosomal and plasma membrane origin. In sheep, endometrial-derived EVs were found to traffic to the conceptus trophectoderm. Thus, EVs are hypothesized to be an important mode of intercellular communication by transferring select RNAs, proteins, and lipids between the endometrium and conceptus. Electron microscopy analysis found that the endometrial luminal and glandular epithelia were the primary source of EVs in the uterus of cyclic sheep. Size exclusion chromatography and nanoparticle tracking analysis (NTA) found that total EV number in the uterine lumen increased from day 10 to 14 in cyclic sheep. Next, ewes were ovariectomized and hormone replaced to determine effects of P4 on the endometrium and EVs in the uterine lumen. Transcriptome analyses found that P4 regulated 1611 genes and nine miRNAs in the endometrium. Total EV number in the uterine lumen was increased by P4 treatment. Small RNA sequencing of EVs detected expression of 768 miRNAs and determined that P4 regulated seven of those miRNAs. These studies provide fundamental new information on how P4 influences endometrial function to regulate conceptus growth for pregnancy establishment in sheep.

Global misregulation of genes largely uncoupled to DNA methylome epimutations characterizes a congenital overgrowth syndrome.

Assisted reproductive therapies (ART) have become increasingly common worldwide and numerous retrospective studies have indicated that ART-conceived children are more likely to develop the overgrowth syndrome Beckwith-Wiedemann (BWS). In bovine, the use of ART can induce a similar overgrowth condition, which is referred to as large offspring syndrome (LOS). Both BWS and LOS involve misregulation of imprinted genes. However, it remains unknown whether molecular alterations at non-imprinted loci contribute to these syndromes. Here we examined the transcriptome of skeletal muscle, liver, kidney, and brain of control and LOS bovine fetuses and found that different tissues within LOS fetuses have perturbations of distinct gene pathways. Notably, in skeletal muscle, multiple pathways involved in myoblast proliferation and fusion into myotubes are misregulated in LOS fetuses. Further, characterization of the DNA methylome of skeletal muscle demonstrates numerous local methylation differences between LOS and controls; however, only a small percent of differentially expressed genes (DEGs), including the imprinted gene IGF2R, could be associated with the neighboring differentially methylated regions. In summary, we not only show that misregulation of non-imprinted genes and loss-of-imprinting characterize the ART-induced overgrowth syndrome but also demonstrate that most of the DEGs is not directly associated with DNA methylome epimutations.

Global assessment of imprinted gene expression in the bovine conceptus by next generation sequencing.

Genomic imprinting is an epigenetic mechanism that leads to parental-allele-specific gene expression. Approximately 150 imprinted genes have been identified in humans and mice but less than 30 have been described as imprinted in cattle. For the purpose of de novo identification of imprinted genes in bovine, we determined global monoallelic gene expression in brain, skeletal muscle, liver, kidney and placenta of day ∼105 Bos taurus indicus × Bos taurus taurus F1 conceptuses using RNA sequencing. To accomplish this, we developed a bioinformatics pipeline to identify parent-specific single nucleotide polymorphism alleles after filtering adenosine to inosine (A-to-I) RNA editing sites. We identified 53 genes subject to monoallelic expression. Twenty three are genes known to be imprinted in the cow and an additional 7 have previously been characterized as imprinted in human and/or mouse that have not been reported as imprinted in cattle. Of the remaining 23 genes, we found that 10 are uncharacterized or unannotated transcripts located in known imprinted clusters, whereas the other 13 genes are distributed throughout the bovine genome and are not close to any known imprinted clusters. To exclude potential cis-eQTL effects on allele expression, we corroborated the parental specificity of monoallelic expression in day 86 Bos taurus taurus × Bos taurus taurus conceptuses and identified 8 novel bovine imprinted genes. Further, we identified 671 candidate A-to-I RNA editing sites and describe random X-inactivation in day 15 bovine extraembryonic membranes. Our results expand the imprinted gene list in bovine and demonstrate that monoallelic gene expression can be the result of cis-eQTL effects.