Comprehensive characterization of horse genome variation by whole‐genome sequencing of 88 horses

Whole‐genome sequencing studies are vital to gain a thorough understanding of genomic variation. Here, we summarize the results of a whole‐genome sequencing study comprising 88 horses and ponies from diverse breeds at 19.1× average coverage. The paired‐end reads were mapped to the current EquCab3.0 horse reference genome assembly, and we identified approximately 23.5 million single nucleotide variants and 2.3 million short indel variants. Our dataset included at least 7 million variants that were not previously reported. On average, each individual horse genome carried ～5.7 million single nucleotides and 0.8 million small indel variants with respect to the reference genome assembly. The variants were functionally annotated. We provide two examples for potentially deleterious recessive alleles that were identified in a heterozygous state in individual genome sequences. Appropriate management of such deleterious recessive alleles in horse breeding programs should help to improve fertility and reduce the prevalence of heritable diseases. This comprehensive dataset has been made publicly available, will represent a valuable resource for future horse genetic studies and supports the goal of accelerating the rates of genetic gain in domestic horse.     

Targeted next-generation sequencing for identifying genes related to horse temperament

It is a fundamental challenge to discover the association of genetic traits with phenotypic traits. In this study, we aimed to identify possible genetic traits related to horse temperament. Based on previous findings, we selected 71 candidate genes related to temperamental trait and examined them in the human and horse reference genomes (hg38 and equCab2, respectively). We found 16 orthologous genes and, by comparing with the human reference genome, 17 homologous genes in the horse reference genome. We designed probes specific for the 33 horse genes. Using the probes, we built sequencing libraries of the genomic DNA samples from eight aggressive and eight docile horses, and sequenced the constructed libraries using the Illumina Hiseq2500 platform. Through the analysis of the targeted exome sequences, we identified single nucleotide polymorphisms (SNPs) in the genes. SNPs could be served as genetic markers to evaluate aggressive or docile levels of horses. To examine whether any genetic variants are associated with horse temperament, we performed genome-wide association study (GWAS) using the SNP data. GWAS analysis identified ten variants (p-value?<0.05) which could be related to horse temperament. We validated the variants using Sanger sequencing. The most significant variants were found in MAOA (c.1164+41T>C) and AR (c.1047+27G>T) genes with 8.09?×?10^?4 p-value. We suggest that the variants might be used to assess horse temperament and to determine superior horses for riding or racing.     

Genome mining reveals the evolutionary origin and biosynthetic potential of basidiomycete polyketide synthases

Numerous polyketides are known from bacteria, plants, and fungi. However, only a few have been isolated from basidiomycetes. Large scale genome sequencing projects now help anticipate the capacity of basidiomycetes to synthesize polyketides. In this study, we identified and annotated 111 type I and three type III polyketide synthase (PKS) genes from 35 sequenced basidiomycete genomes. Phylogenetic analysis of PKS genes suggests that all main types of fungal iterative PKS had already evolved before the Ascomycota and Basidiomycota diverged. A comparison of genomic and metabolomic data shows that the number of polyketide genes exceeds the number of known polyketide structures by far. Exploiting these results to design degenerate PCR primers, we amplified and cloned the complete sequence of armB, a PKS gene from the melleolide producer Armillaria mellea. We expect this study will serve as a guide for future genomic mining projects to discover structurally diverse mushroom-derived polyketides.     

Comparison of the Equine Reference Sequence with Its Sanger Source Data and New Illumina Reads

The reference assembly for the domestic horse, EquCab2, published in 2009, was built using approximately 30 million Sanger reads from a Thoroughbred mare named Twilight. Contiguity in the assembly was facilitated using nearly 315 thousand BAC end sequences from Twilight’s half brother Bravo. Since then, it has served as the foundation for many genome-wide analyses that include not only the modern horse, but ancient horses and other equid species as well. As data mapped to this reference has accumulated, consistent variation between mapped datasets and the reference, in terms of regions with no read coverage, single nucleotide variants, and small insertions/deletions have become apparent. In many cases, it is not clear whether these differences are the result of true sequence variation between the research subjects’ and Twilight’s genome or due to errors in the reference. EquCab2 is regarded as “The Twilight Assembly.” The objective of this study was to identify inconsistencies between the EquCab2 assembly and the source Twilight Sanger data used to build it. To that end, the original Sanger and BAC end reads have been mapped back to this equine reference and assessed with the addition of approximately 40X coverage of new Illumina Paired-End sequence data. The resulting mapped datasets identify those regions with low Sanger read coverage, as well as variation in genomic content that is not consistent with either the original Twilight Sanger data or the new genomic sequence data generated from Twilight on the Illumina platform. As the haploid EquCab2 reference assembly was created using Sanger reads derived largely from a single individual, the vast majority of variation detected in a mapped dataset comprised of those same Sanger reads should be heterozygous. In contrast, homozygous variations would represent either errors in the reference or contributions from Bravo''s BAC end sequences. Our analysis identifies 720,843 homozygous discrepancies between new, high throughput genomic sequence data generated for Twilight and the EquCab2 reference assembly. Most of these represent errors in the assembly, while approximately 10,000 are demonstrated to be contributions from another horse. Other results are presented that include the binary alignment map file of the mapped Sanger reads, a list of variants identified as discrepancies between the source data and resulting reference, and a BED annotation file that lists the regions of the genome whose consensus was likely derived from low coverage alignments.     

LH pulse frequency and the emergence and growth of ovarian antral follicular waves in the ewe during the luteal phase of the estrous cycle

Background

Oviduct epithelial cells (OEC) co-culture promotes in vitro fertilization (IVF) in human, bovine and porcine species, but no data are available from equine species. Yet, despite numerous attempts, equine IVF rates remain low. Our first aim was to verify a beneficial effect of the OEC on equine IVF. In mammals, oviductal proteins have been shown to interact with gametes and play a role in fertilization. Thus, our second aim was to identify the proteins involved in fertilization in the horse.

Methods & results

In the first experiment, we co-incubated fresh equine spermatozoa treated with calcium ionophore and in vitro matured equine oocytes with or without porcine OEC. We showed that the presence of OEC increases the IVF rates. In the subsequent experiments, we co-incubated equine gametes with OEC and we showed that the IVF rates were not significantly different between 1) gametes co-incubated with equine vs porcine OEC, 2) intact cumulus-oocyte complexes vs denuded oocytes, 3) OEC previously stimulated with human Chorionic Gonadotropin, Luteinizing Hormone and/or oestradiol vs non stimulated OEC, 4) in vivo vs in vitro matured oocytes. In order to identify the proteins responsible for the positive effect of OEC, we first searched for the presence of the genes encoding oviductin, osteopontin and atrial natriuretic peptide A (ANP A) in the equine genome. We showed that the genes coding for osteopontin and ANP A are present. But the one for oviductin either has become a pseudogene during evolution of horse genome or has been not well annotated in horse genome sequence. We then showed that osteopontin and ANP A proteins are present in the equine oviduct using a surface plasmon resonance biosensor, and we analyzed their expression during oestrus cycle by Western blot. Finally, we co-incubated equine gametes with or without purified osteopontin or synthesized ANP A. No significant effect of osteopontin or ANP A was observed, though osteopontin slightly increased the IVF rates.

Conclusion

Our study shows a beneficial effect of homologous and heterologous oviduct cells on equine IVF rates, though the rates remain low. Furthers studies are necessary to identify the proteins involved. We showed that the surface plasmon resonance technique is efficient and powerful to analyze molecular interactions during fertilization.     

Decoding the genome of superior chapatti quality Indian wheat variety ‘C 306’ unravelled novel genomic variants for chapatti and nutrition quality related genes

An Indian wheat variety, ‘C 306’ has good chapatti quality, which is controlled by multiple genes that have not been explored. We report the high quality de novo assembled genome of ‘C 306’ by combining short and long read sequencing data. The hybrid assembly covered 93% of gene space and identified about 142 K coding genes, 34% repetitive DNA and ~ 501 K SSR motifs. The phylogenetic analysis of about 83 K orthologous protein groups suggested the closest relationship with T. turgidum, T. aestivum and Ae. tauschii. Genome wide analysis annotated 69,217,536 genomic variants. Out of them, 1423 missense and 117 deleterious variants identified in processing, nutrition, and chapatti quality related genes such as alpha- and beta-gliadin, SSI, SSIII, SUT1, SBEI, CHS, YSL, DMAS, and NAS encoded proteins. These variants may affect quality genes. The genomic data will be potential genomic resources in wheat breeding programs for quality improvement.