Functional annotation of the animal genomes: An integrated annotation resource for the horse.

The genomic sequence of the horse has been available since 2009, providing critical resources for discovering important genomic variants regarding both animal health and population structures. However, to fully understand the functional implications of these variants, detailed annotation of the horse genome is required. Due to the limited availability of functional data for the equine genome, as well as the technical limitations of short-read RNA-seq, existing annotation of the equine genome contains limited information about important aspects of gene regulation, such as alternate isoforms and regulatory elements, which are either not transcribed or transcribed at a very low level. To solve above problems, the Functional Annotation of the Animal Genomes (FAANG) project proposed a systemic approach to tissue collection, phenotyping, and data generation, adopting the blueprint laid out by the Encyclopedia of DNA Elements (ENCODE) project. Here we detail the first comprehensive overview of gene expression and regulation in the horse, presenting 39,625 novel transcripts, 84,613 candidate cis-regulatory elements (CRE) and their target genes, 332,115 open chromatin regions genome wide across a diverse set of tissues. We showed substantial concordance between chromatin accessibility, chromatin states in different genic features and gene expression. This comprehensive and expanded set of genomics resources will provide the equine research community ample opportunities for studies of complex traits in the horse.

Allelic Variation in CXCL16 Determines CD3+ T Lymphocyte Susceptibility to Equine Arteritis Virus Infection and Establishment of Long-Term Carrier State in the Stallion.

Equine arteritis virus (EAV) is the causative agent of equine viral arteritis (EVA), a respiratory, systemic, and reproductive disease of horses and other equid species. Following natural infection, 10-70% of the infected stallions can become persistently infected and continue to shed EAV in their semen for periods ranging from several months to life. Recently, we reported that some stallions possess a subpopulation(s) of CD3+ T lymphocytes that are susceptible to in vitro EAV infection and that this phenotypic trait is associated with long-term carrier status following exposure to the virus. In contrast, stallions not possessing the CD3+ T lymphocyte susceptible phenotype are at less risk of becoming long-term virus carriers. A genome wide association study (GWAS) using the Illumina Equine SNP50 chip revealed that the ability of EAV to infect CD3+ T lymphocytes and establish long-term carrier status in stallions correlated with a region within equine chromosome 11. Here we identified the gene and mutations responsible for these phenotypes. Specifically, the work implicated three allelic variants of the equine orthologue of CXCL16 (EqCXCL16) that differ by four non-synonymous nucleotide substitutions (XM_00154756; c.715 A → T, c.801 G → C, c.804 T → A/G, c.810 G → A) within exon 1. This resulted in four amino acid changes with EqCXCL16S (XP_001504806.1) having Phe, His, Ile and Lys as compared to EqCXL16R having Tyr, Asp, Phe, and Glu at 40, 49, 50, and 52, respectively. Two alleles (EqCXCL16Sa, EqCXCL16Sb) encoded identical protein products that correlated strongly with long-term EAV persistence in stallions (P<0.000001) and are required for in vitro CD3+ T lymphocyte susceptibility to EAV infection. The third (EqCXCL16R) was associated with in vitro CD3+ T lymphocyte resistance to EAV infection and a significantly lower probability for establishment of the long-term carrier state (viral persistence) in the male reproductive tract. EqCXCL16Sa and EqCXCL16Sb exert a dominant mode of inheritance. Most importantly, the protein isoform EqCXCL16S but not EqCXCL16R can function as an EAV cellular receptor. Although both molecules have equal chemoattractant potential, EqCXCL16S has significantly higher scavenger receptor and adhesion properties compared to EqCXCL16R.

Developing a 670k genotyping array to tag ~2M SNPs across 24 horse breeds.

BACKGROUND: To date, genome-scale analyses in the domestic horse have been limited by suboptimal single nucleotide polymorphism (SNP) density and uneven genomic coverage of the current SNP genotyping arrays. The recent availability of whole genome sequences has created the opportunity to develop a next generation, high-density equine SNP array. RESULTS: Using whole genome sequence from 153 individuals representing 24 distinct breeds collated by the equine genomics community, we cataloged over 23 million de novo discovered genetic variants. Leveraging genotype data from individuals with both whole genome sequence, and genotypes from lower-density, legacy SNP arrays, a subset of ~5 million high-quality, high-density array candidate SNPs were selected based on breed representation and uniform spacing across the genome. Considering probe design recommendations from a commercial vendor (Affymetrix, now Thermo Fisher Scientific) a set of ~2 million SNPs were selected for a next-generation high-density SNP chip (MNEc2M). Genotype data were generated using the MNEc2M array from a cohort of 332 horses from 20 breeds and a lower-density array, consisting of ~670 thousand SNPs (MNEc670k), was designed for genotype imputation. CONCLUSIONS: Here, we document the steps taken to design both the MNEc2M and MNEc670k arrays, report genomic and technical properties of these genotyping platforms, and demonstrate the imputation capabilities of these tools for the domestic horse.

Nuclear respiratory factor-1 and bioenergetics in tamoxifen-resistant breast cancer cells.

Acquired tamoxifen (TAM) resistance is a significant clinical problem in treating patients with estrogen receptor α (ERα)+ breast cancer. We reported that ERα increases nuclear respiratory factor-1 (NRF-1), which regulates nuclear-encoded mitochondrial gene transcription, in MCF-7 breast cancer cells and NRF-1 knockdown stimulates apoptosis. Whether NRF-1 and target gene expression is altered in endocrine resistant breast cancer cells is unknown. We measured NRF-1and metabolic features in a cell model of progressive TAM-resistance. NRF-1 and its target mitochondrial transcription factor A (TFAM) were higher in TAM-resistant LCC2 and LCC9 cells than TAM-sensitive MCF-7 cells. Using extracellular flux assays we observed that LCC1, LCC2, and LCC9 cells showed similar oxygen consumption rate (OCR), but lower mitochondrial reserve capacity which was correlated with lower Succinate Dehydrogenase Complex, Subunit B in LCC1 and LCC2 cells. Complex III activity was lower in LCC9 than MCF-7 cells. LCC1, LCC2, and LCC9 cells had higher basal extracellular acidification (ECAR), indicating higher aerobic glycolysis, relative to MCF-7 cells. Mitochondrial bioenergetic responses to estradiol and 4-hydroxytamoxifen were reduced in the endocrine-resistant cells compared to MCF-7 cells. These results suggest the acquisition of altered metabolic phenotypes in response to long term antiestrogen treatment may increase vulnerability to metabolic stress.

Prehistoric genomes reveal the genetic foundation and cost of horse domestication.

The domestication of the horse ∼ 5.5 kya and the emergence of mounted riding, chariotry, and cavalry dramatically transformed human civilization. However, the genetics underlying horse domestication are difficult to reconstruct, given the near extinction of wild horses. We therefore sequenced two ancient horse genomes from Taymyr, Russia (at 7.4- and 24.3-fold coverage), both predating the earliest archeological evidence of domestication. We compared these genomes with genomes of domesticated horses and the wild Przewalski's horse and found genetic structure within Eurasia in the Late Pleistocene, with the ancient population contributing significantly to the genetic variation of domesticated breeds. We furthermore identified a conservative set of 125 potential domestication targets using four complementary scans for genes that have undergone positive selection. One group of genes is involved in muscular and limb development, articular junctions, and the cardiac system, and may represent physiological adaptations to human utilization. A second group consists of genes with cognitive functions, including social behavior, learning capabilities, fear response, and agreeableness, which may have been key for taming horses. We also found that domestication is associated with inbreeding and an excess of deleterious mutations. This genetic load is in line with the "cost of domestication" hypothesis also reported for rice, tomatoes, and dogs, and it is generally attributed to the relaxation of purifying selection resulting from the strong demographic bottlenecks accompanying domestication. Our work demonstrates the power of ancient genomes to reconstruct the complex genetic changes that transformed wild animals into their domesticated forms, and the population context in which this process took place.

Characterization of ovine Nectin-4, a novel peste des petits ruminants virus receptor.

Small ruminants infected with peste des petits ruminants virus exhibit lesions typical of epithelial infection and necrosis. However, the only established host receptor for this virus is the immune cell marker signaling lymphocyte activation molecule (SLAM). We have confirmed that the ovine Nectin-4 protein, when overexpressed in epithelial cells, permits efficient replication of PPRV. Furthermore, this gene was predominantly expressed in epithelial tissues and encoded by multiple haplotypes in sheep breeds from around the world.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats.

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Evidence for origin of lavender foal syndrome among Egyptian Arabian horses in Egypt.

BACKGROUND: Lavender foal syndrome (LFS) is a fatal hereditary condition that is inherited in an autosomal recessive pattern. This detrimental mutation is more common in Arabian foals of Egyptian origin than foals from other bloodlines. Heterozygous horses are carriers of the LFS trait and appear normal, while recessive homozygous foals died shortly after birth due to serious complications. In Egypt, in 2014, an Egyptian foal died after manifestations of neurological signs and abnormal coat colour as LFS signs. Therefore, it is important to identify LFS carriers in the population of Arabian horses in Egypt and to encourage improvement of the Arabian horse industry in Egypt by constructing a breeding system based on genetic background in order to avoid mating between carriers and reduce financial losses from deaths of affected foals. OBJECTIVES: To establish a PCR-based test for detecting the MYO5A gene mutation causing LFS in the registered Arabian horse population in Egypt prior to breeding. STUDY DESIGN: Cross sectional survey (n = 170) plus targeted sampling (n = 30). METHODS: A total of 200 samples were collected from an Arabian farm in Egypt and some of them were traced for LFS based on the farm records. The LFS genotypes were identified using the PCR-RFLP technique, fragment analysis followed by sequence analysis. RESULTS: The overall mutated allele and genotype frequencies (N/L) were 0.08 and 16%, respectively. CONCLUSION: The observed frequency of heterozygotes suggests foals affected with LFS will be produced among Arabian horses in Egypt. Therefore, screening of the entire population for this mutation should be undertaken in the breeding program.

Imputed genomes of historical horses provide insights into modern breeding.

Historical genomes can provide important insights into recent genomic changes in horses, especially the development of modern breeds. In this study, we characterized 8.7 million genomic variants from a panel of 430 horses from 73 breeds, including newly sequenced genomes from 20 Clydesdales and 10 Shire horses. We used this modern genomic variation to impute the genomes of four historically important horses, consisting of publicly available genomes from 2 Przewalski's horses, 1 Thoroughbred, and a newly sequenced Clydesdale. Using these historical genomes, we identified modern horses with higher genetic similarity to those in the past and unveiled increased inbreeding in recent times. We genotyped variants associated with appearance and behavior to uncover previously unknown characteristics of these important historical horses. Overall, we provide insights into the history of Thoroughbred and Clydesdale breeds and highlight genomic changes in the endangered Przewalski's horse following a century of captive breeding.

The localization of centromere protein A is conserved among tissues.

Centromeres are epigenetically specified by the histone H3 variant CENP-A. Although mammalian centromeres are typically associated with satellite DNA, we previously demonstrated that the centromere of horse chromosome 11 (ECA11) is completely devoid of satellite DNA. We also showed that the localization of its CENP-A binding domain is not fixed but slides within an about 500 kb region in different individuals, giving rise to positional alleles. These epialleles are inherited as Mendelian traits but their position can move in one generation. It is still unknown whether centromere sliding occurs during meiosis or during development. Here, we first improve the sequence of the ECA11 centromeric region in the EquCab3.0 assembly. Then, to test whether centromere sliding may occur during development, we map the CENP-A binding domains of ECA11 using ChIP-seq in five tissues of different embryonic origin from the four horses of the equine FAANG (Functional Annotation of ANimal Genomes) consortium. Our results demonstrate that the centromere is localized in the same region in all tissues, suggesting that the position of the centromeric domain is maintained during development.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats.

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared to its predecessor. Gene annotations are now more complete, significantly improving the mapping precision of genomic, transcriptomic, and proteomics data sets. We jointly analyzed 163 short-read whole genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ~20.0 million sequence variations, of which 18.7 thousand are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Decoding the Equine Genome: Lessons from ENCODE.

The horse reference genome assemblies, EquCab2.0 and EquCab3.0, have enabled great advancements in the equine genomics field, from tools to novel discoveries. However, significant gaps of knowledge regarding genome function remain, hindering the study of complex traits in horses. In an effort to address these gaps and with inspiration from the Encyclopedia of DNA Elements (ENCODE) project, the equine Functional Annotation of Animal Genome (FAANG) initiative was proposed to bridge the gap between genome and gene expression, providing further insights into functional regulation within the horse genome. Three years after launching the initiative, the equine FAANG group has generated data from more than 400 experiments using over 50 tissues, targeting a variety of regulatory features of the equine genome. In this review, we examine how valuable lessons learned from the ENCODE project informed our decisions in the equine FAANG project. We report the current state of the equine FAANG project and discuss how FAANG can serve as a template for future expansion of functional annotation in the equine genome and be used as a reference for studies of complex traits in horse. A well-annotated reference functional atlas will also help advance equine genetics in the pan-genome and precision medicine era.

Generation of a Biobank From Two Adult Thoroughbred Stallions for the Functional Annotation of Animal Genomes Initiative.

Following the successful creation of a biobank from two adult Thoroughbred mares, this study aimed to recapitulate sample collection in two adult Thoroughbred stallions as part of the Functional Annotation of the Animal Genome (FAANG) initiative. Both stallions underwent thorough physical, lameness, neurologic, and ophthalmic (including electroretinography) examinations prior to humane euthanasia. Epididymal sperm was recovered from both stallions immediately postmortem and cryopreserved. Aseptically collected full thickness skin biopsies were used to isolate, culture and cryopreserve dermal fibroblasts. Serum, plasma, cerebrospinal fluid, urine, and gastrointestinal content from various locations were collected and cryopreserved. Under guidance of a board-certified veterinary anatomic pathologist, 102 representative tissue samples were collected from both horses. Whole tissue samples were flash-frozen and prioritized tissues had nuclei isolated and cryopreserved. Spatially contemporaneous samples of each tissue were submitted for histologic examination. Antemortem and gross pathologic examination revealed mild abnormalities in both stallions. One stallion (ECA_UCD_AH3) had unilateral thoracic limb lameness and bilateral chorioretinal scars. The second stallion (ECA_UCD_AH4) had subtle symmetrical pelvic limb ataxia, symmetrical prostatomegally, and moderate gastrointestinal nematodiasis. DNA from each was whole-genome sequenced and genotyped using the GGP Equine 70K SNP array. The genomic resources and banked biological samples from these animals augments the existing resource available to the equine genomics community. Importantly we may now improve the resolution of tissue-specific gene regulation as affected by sex, as well as add sex-specific tissues and gametes.

Assessment of genetic variation for the LINE-1 retrotransposon from next generation sequence data.

BACKGROUND: In humans, copies of the Long Interspersed Nuclear Element 1 (LINE-1) retrotransposon comprise 21% of the reference genome, and have been shown to modulate expression and produce novel splice isoforms of transcripts from genes that span or neighbor the LINE-1 insertion site. RESULTS: In this work, newly released pilot data from the 1000 Genomes Project is analyzed to detect previously unreported full length insertions of the retrotransposon LINE-1. By direct analysis of the sequence data, we have identified 22 previously unreported LINE-1 insertion sites within the sequence data reported for a mother/father/daughter trio. CONCLUSIONS: It is demonstrated here that next generation sequencing data, as well as emerging high quality datasets from individual genome projects allow us to assess the amount of heterogeneity with respect to the LINE-1 retrotransposon amongst humans, and provide us with a wealth of testable hypotheses as to the impact that this diversity may have on the health of individuals and populations.

Interaction among variant vascular endothelial growth factor (VEGF) and its receptor in relation to prostate cancer risk.

BACKGROUND: Prostate cancer (PCa) incidence and mortality are disproportionately high among African-American (AA) men. Its detection and perhaps its disparities could be improved through the identification of genetic susceptibility biomarkers within essential biological pathways. Interactions among highly variant genes, central to angiogenesis, may modulate susceptibility for prostate cancer, as previous demonstrated. This study evaluates the interplay among three highly variant genes (i.e., IL-10, TGFbetaR-1, VEGF), their receptors and their influence on PCa within a case-control study consisting of an under-served population. METHODS: This study evaluated single gene and joint modifying effects on PCa risk in a case-control study comprised of 859 AA men (193 cases and 666 controls) using TaqMan qPCR. Interaction among polymorphic IL-10, TGFbetaR-1 and VEGF was analyzed using conventional logistic regression analysis (LR) models, multi-dimensionality reduction (MDR) and interaction entropy graphs. Symbolic modeling allowed validation of gene-gene interaction findings identified by MDR. RESULTS: No significant single gene effects were demonstrated in relation to PCa risk. However, carriers of the VEGF 2482T allele had a threefold increase in the risk of developing aggressive PCa. The presence of VEGF 2482T combined with VEGFR IVS6 + 54 loci were highly significant for the risk of PCa based on MDR and symbolic modeling analyses. These findings were substantiated by 1,000-fold cross validation permutation testing (P = 0.04), respectively. CONCLUSION: These findings suggest the inheritance of VEGF and VEGFR IVS6 + 54 sequence variants may jointly modify PCa susceptibility through their influence on angiogenesis. Larger sub-population studies are needed to validate these findings and evaluate whether the VEGF-VEGR axis may serve as predictors of disease prognosis and ultimately clinical response to available treatment strategies.

Using triallelic SNPs for determining parentage in North American yak ( Bos grunniens) and estimating cattle ( B. taurus) introgression.

Background: Genetic testing for pedigree accuracy is critical for managing genetic diversity in North American (NA) yak ( Bos grunniens), a population expanded mostly from imported zoological park specimens.  DNA testing also enhances species conservation by identifying recent B. taurus F1 hybrid ancestors (within three generations).  Biallelic single nucleotide polymorphisms (SNPs) can accomplish either task, but increases the marker count and costs necessary to achieve both.  Our aim was to identify novel, multifunctional, triallelic yak SNPs (tySNPs), with each having two alleles for yak parentage testing, and a third allele for identifying recent cattle introgression.  Methods:  Genome sequences were aligned to the cattle UMD3.1 assembly and SNPs were screened for 1) heterozygosity in a NA and a Chinese yak, 2) a third allele at high frequency in cattle, and 3) flanking sequences conserved in both species.  Subsequently, tySNPs were filtered for unique alignment to the haplotype-resolved F1 yak assembly.  Allele frequencies were estimated in a subset of 87 tySNPs by genotyping 170 NA yak. Results:  We identified 610 autosomal tySNPs, distributed in 441 clusters with 5 Mb average genome spacing.  The average NA yak minor allele frequency was high (0.296), while average introgressed cattle alleles were low (0.004).  In simulations with tySNPs, 28 were sufficient for globally-unique animal identification (P I=5.81x10 -12), 87 were able to exclude 19 random bulls from parentage at the 99% level without using the dam's genotype (P E=5.3x10 -4), and 87 were able to detect F1 hybridization events after three generations of yak backcrosses (1/16th B. taurus germplasm). Conclusions:  Identifying animals, determining parentage and detecting recent hybridization events was efficient with as few as 87 tySNPs.  A similar triallelic approach could be used with other bottlenecked Bos species that hybridize with cattle, such as NA plains bison ( B. bison).

A MicroRNA gene is hosted in an intron of a schizophrenia-susceptibility gene.

Schizophrenia (SZ) is a neuropsychiatric disorder that affects about 1% of the adult population. Numerous genes have been implicated in SZ susceptibility. MicroRNAs (miRNA) are small RNA molecules that regulate the translation of mRNAs via interactions with their 3' untranslated regions. Identification of known miRNA targets on all human genes indicated that miRNA-346 targets SZ susceptibility genes listed in the SchizophreniaGene database twice as frequently as expected relative to other genes in the genome. The gene encoding this miRNA, miR-346, is located in intron 2 of the glutamate receptor ionotropic delta 1 (GRID1) gene, which has been previously implicated in SZ susceptibility. We used quantitative real-time PCR to determine the expression levels of miR-346 and GRID1 using brain RNA samples from the Stanley Array Collection, Stanley Medical Research Institute. Expression of both miR-346 and GRID1 is lower in SZ patients than that in normal controls (P=0.017 and 0.086, respectively). However, the expression of miR-346 and GRID1 is less correlated in SZ patients than in bipolar patients or in normal controls. This study implicates the importance of a miRNA in SZ.

Epigenetic control of mammalian LINE-1 retrotransposon by retinoblastoma proteins.

Long interspersed nuclear elements (LINEs or L1 elements) are targeted for epigenetic silencing during early embryonic development and remain inactive in most cells and tissues. Here we show that E2F-Rb family complexes participate in L1 elements epigenetic regulation via nucleosomal histone modifications and recruitment of histone deacetylases (HDACs) HDAC1 and HDAC2. Our experiments demonstrated that (i) Rb and E2F interact with human and mouse L1 elements, (ii) L1 elements are deficient in both heterochromatin-associated histone marks H3 tri methyl K9 and H4 tri methyl K20 in Rb family triple knock out (Rb, p107, and p130) fibroblasts (TKO), (iii) L1 promoter exhibits increased histone H3 acetylation in the absence of HDAC1 and HDAC2 recruitment, (iv) L1 expression in TKO fibroblasts is upregulated compared to wild type counterparts, (v) L1 expression increases in the presence of the HDAC inhibitor TSA. On the basis of these findings we propose a model in which L1 sequences throughout the genome serve as centers for heterochromatin formation in an Rb family-dependent manner. As such, Rb proteins and L1 elements may play key roles in heterochromatin formation beyond pericentromeric chromosomal regions. These findings describe a novel mechanism of L1 reactivation in mammalian cells mediated by failure of corepressor protein recruitment by Rb, loss of histone epigenetic marks, heterochromatin formation, and increased histone H3 acetylation.

Coding RNA Sequencing of Equine Endometrium during Maternal Recognition of Pregnancy.

Equine maternal recognition of pregnancy (MRP) is a process whose signal remains unknown. During MRP the conceptus and endometrium communicate to attenuate prostaglandin F2α (PGF) secretion, sparing the corpus luteum and maintaining progesterone production. Recognition of a mobile conceptus by the endometrium is critical by days 14-16 post-ovulation (PO), when endometrium produces PGF, initiating luteolysis. The objective of this study was to evaluate endometrial gene expression changes based upon pregnancy status via RNA sequencing. This experiment utilized a cross-over design with each mare serving as both a pregnant and non-mated control on days nine, 11, and 13 PO (n = 3/status/day). Mares were randomly assigned to collection day and pregnancy confirmed by terminal uterine lavage at the time of endometrial biopsy. Total RNA was isolated and libraries prepared using Illumina TruSeq RNA sample preparation kit. Reads were mapped and annotated using HISAT2 and Stringtie. Expression values were evaluated with DESEQ2 (P ≤ 0.05 indicated significance). On day nine, 11, and 13 there were 1435, 1435 and 916 significant transcripts, respectively. Multiple genes with splice variants had different expression patterns within the same day. These are the first data to evaluate the endometrial transcriptome during MRP on days nine, 11, and 13.