Search

Agricultural data are crucial to many aspects of production, commerce, and research involved in feeding the global community. However, in most agricultural research disciplines standard best practices for data management and publication do not exist. Here we propose a set of best practices in the areas of peer review, minimal dataset development, data repositories, citizen science initiatives, and support for best data management. We illustrate some of these best practices with a case study in dairy agroecosystems research. While many common, and increasingly disparate data management and publication practices are entrenched in agricultural disciplines, opportunities are readily available for promoting and adopting best practices that better enable and enhance data-intensive agricultural research and production. National Institute of Food and AgricultureUnited States Department of Agriculture (USDA)National Institute of Food and Agriculture [2018-67023-27843]; Agricultural Research ServiceUnited States Department of Agriculture (USDA)USDA Agricultural Research Service [8042-31000-001-00-D, 8260-88888-001-00-D] Published version National Institute of Food and Agriculture, Grant/Award Number: 2018-67023-27843; Agricultural Research Service, Grant/Award Numbers: 8042-31000-001-00-D, 8260-88888-001-00-D

We have estimated the average genetic diversity of two Glycine annual and six perennial species based upon 76 orthologous gene sets and performed phylogenetic analysis, divergence analysis and tests for departure from neutrality of the eight species using 52 orthologous gene sets. In addition, 367 orthologous gene sets were used to estimate the relationships of 11 G. canescens accessions. Among the perennials, G. canescens showed the highest nucleotide diversity. The other perennials, except for G. tomentella, had higher nucleotide diversity than the two annuals. Phylogenetic analysis of the Glycine showed a similar genome grouping with the previous report except for G. cyrtoloba and G. stenophita which formed a sister clade in the study. Divergence analysis supported the phylogenetic relationships that G. falcata was the most divergent from G. max, followed by G. cyrtoloba, G. syndetika, G. tomentella D3, G. stenophita and G. canescens. Most genic sequences were homogeneous in the levels of polymorphism and divergence between G. max and other Glycine species based on the HKA test, thus, Glycine perennials may have experienced a very similar evolution as inferred by trans-specific mutation analysis. The greater genetic diversity of most perennial Glycine species and their origins from the warmer and drier climates of Australia suggests the perennials maybe a potential source of heat and drought resistance that will be of value in the face of climate change.

By uniformly analyzing 723 RNA-seq data from 91 tissues and cell types, we built a comprehensive gene atlas and studied tissue specificity of genes in cattle. We demonstrated that tissue-specific genes significantly reflected the tissue-relevant biology, showing distinct promoter methylation and evolution patterns (e.g., brain-specific genes evolve slowest, whereas testis-specific genes evolve fastest). Through integrative analyses of those tissue-specific genes with large-scale genome-wide association studies, we detected relevant tissues/cell types and candidate genes for 45 economically important traits in cattle, including blood/immune system (e.g., CCDC88C) for male fertility, brain (e.g., TRIM46 and RAB6A) for milk production, and multiple growth-related tissues (e.g., FGF6 and CCND2) for body conformation. We validated these findings by using epigenomic data across major somatic tissues and sperm. Collectively, our findings provided novel insights into the genetic and biological mechanisms underlying complex traits in cattle, and our transcriptome atlas can serve as a primary source for biological interpretation, functional validation, studies of adaptive evolution, and genomic improvement in livestock.

A study was conducted to determine the effects of a diet supplemented with fruits and vegetables (FV) on the host whole blood cell (WBC) transcriptome and the composition and function of the intestinal microbiome. Nine six-week-old pigs were fed a pig grower diet alone or supplemented with lyophilized FV equivalent to half the daily recommended amount prescribed for humans by the Dietary Guideline for Americans (DGA) for two weeks. Host transcriptome changes in the WBC were evaluated by RNA sequencing. Isolated DNA from the fecal microbiome was used for 16S rDNA taxonomic analysis and prediction of metabolomic function. Feeding an FV-supplemented diet to pigs induced differential expression of several genes associated with an increase in B-cell development and differentiation and the regulation of cellular movement, inflammatory response, and cell-to-cell signaling. Linear discriminant analysis effect size (LEfSe) in fecal microbiome samples showed differential increases in genera from Lachnospiraceae and Ruminococcaceae families within the order Clostridiales and Erysipelotrichaceae family with a predicted reduction in rgpE-glucosyltransferase protein associated with lipopolysaccharide biosynthesis in pigs fed the FV-supplemented diet. These results suggest that feeding an FV-supplemented diet for two weeks modulated markers of cellular inflammatory and immune function in the WBC transcriptome and the composition of the intestinal microbiome by increasing the abundance of bacterial taxa that have been associated with improved intestinal health.

Background Efforts to improve animal health, and understand genetic bases for production, may benefit from a comprehensive analysis of animal genomes and epigenomes. Although DNA methylation has been well studied in humans and other model species, its distribution patterns and regulatory impacts in cattle are still largely unknown. Here, we present the largest collection of cattle DNA methylation epigenomic data to date. Results Using Holstein cattle, we generated 29 whole genome bisulfite sequencing (WGBS) datasets for 16 tissues, 47 corresponding RNA-seq datasets, and 2 whole genome sequencing datasets. We did read mapping and DNA methylation calling based on two different cattle assemblies, demonstrating the high quality of the long-read-based assembly markedly improved DNA methylation results. We observed large differences across cattle tissues in the methylation patterns of global CpG sites, partially methylated domains (PMDs), hypomethylated regions (HMRs), CG islands (CGIs), and common repeats. We detected that each tissue had a distinct set of PMDs, which showed tissue-specific patterns. Similar to human PMD, cattle PMDs were often linked to a general decrease of gene expression and a decrease in active histone marks and related to long-range chromatin organizations, like topologically associated domains (TADs). We tested a classification of the HMRs based on their distributions relative to transcription start sites (TSSs) and detected tissue-specific TSS-HMRs and genes that showed strong tissue effects. When performing cross-species comparisons of paired genes (two opposite strand genes with their TSS located in the same HMR), we found out they were more consistently co-expressed among human, mouse, sheep, goat, yak, pig, and chicken, but showed lower consistent ratios in more divergent species. We further used these WGBS data to detect 50,023 experimentally supported CGIs across bovine tissues and found that they might function as a guard against C-to-T mutations for TSS-HMRs. Although common repeats were often heavily methylated, some young Bov-A2 repeats were hypomethylated in sperm and could affect the promoter structures by exposing potential transcription factor binding sites. Conclusions This study provides a comprehensive resource for bovine epigenomic research and enables new discoveries about DNA methylation and its role in complex traits.

Summary With the availability of a new highly contiguous Bos taurus reference genome assembly (ARS‐UCD1.2), it is the opportune time to upgrade the bovine gene set by seeking input from researchers. Furthermore, advances in graphical genome annotation tools now make it possible for researchers to leverage sequence data generated with the latest technologies to collaboratively curate genes. For many years the Bovine Genome Database (BGD) has provided tools such as the apollo genome annotation editor to support manual bovine gene curation. The goal of this paper is to explain the reasoning behind the decisions made in the manual gene curation process while providing examples using the existing BGD tools. We will describe the sources of gene annotation evidence provided at the BGD, including RNA‐seq and Iso‐Seq data. We will also explain how to interpret various data visualizations when curating gene models, and will demonstrate the value of manual gene annotation. The process described here can be applied to manual gene curation for other species with similar tools. With a better understanding of manual gene annotation, researchers will be encouraged to edit gene models and contribute to the enhancement of livestock gene sets.

Author(s): Rosen, Benjamin D; Bickhart, Derek M; Schnabel, Robert D; Koren, Sergey; Elsik, Christine G; Tseng, Elizabeth; Rowan, Troy N; Low, Wai Y; Zimin, Aleksey; Couldrey, Christine; Hall, Richard; Li, Wenli; Rhie, Arang; Ghurye, Jay; McKay, Stephanie D; Thibaud-Nissen, Francoise; Hoffman, Jinna; Murdoch, Brenda M; Snelling, Warren M; McDaneld, Tara G; Hammond, John A; Schwartz, John C; Nandolo, Wilson; Hagen, Darren E; Dreischer, Christian; Schultheiss, Sebastian J; Schroeder, Steven G; Phillippy, Adam M; Cole, John B; Van Tassell, Curtis P; Liu, George; Smith, Timothy PL; Medrano, Juan F | Abstract: BackgroundMajor 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.ResultsWe 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 g250× more continuous than the original assembly, with contig N50 g25 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.ConclusionsWe 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.

A targeted metabolomic analysis was performed on tissues derived from pigs fed diets supplemented with white button mushrooms (WBM) to determine the effect on the liver and brain metabolome. Thirty-one pigs were fed a grower diet alone or supplemented with either three or six servings of freeze-dried WBM for six weeks. Tissue metabolomes were analyzed using targeted liquid chromatography-mass spectrometry (LC-MS) combined with chemical similarity enrichment analysis (ChemRICH) and correlated to WBM-induced changes in fecal microbiome composition. Results indicated that WBM can differentially modulate metabolites in liver, brain cortex and hippocampus of healthy pigs. Within the glycero-phospholipids, there was an increase in alkyl-acyl-phosphatidyl-cholines (PC-O 40:3) in the hippocampus of pigs fed six servings of WBM. A broader change in glycerophospholipids and sphingolipids was detected in the liver with a reduction in several lipid species in pigs fed both WBM diets but with an increase in amino acids known as precursors of neurotransmitters in the cortex of pigs fed six servings of WBM. Metabolomic changes were positively correlated with increased abundance of Cryomorphaceae, Lachnospiraceae, Flammeovirgaceae and Ruminococcaceae in the microbiome suggesting that WBM can also positively impact tissue metabolite composition.

Background The domestic pig (Sus scrofa) is important both as a food source and as a biomedical model given its similarity in size, anatomy, physiology, metabolism, pathology, and pharmacology to humans. The draft reference genome (Sscrofa10.2) of a purebred Duroc female pig established using older clone-based sequencing methods was incomplete, and unresolved redundancies, short-range order and orientation errors, and associated misassembled genes limited its utility. Results We present 2 annotated highly contiguous chromosome-level genome assemblies created with more recent long-read technologies and a whole-genome shotgun strategy, 1 for the same Duroc female (Sscrofa11.1) and 1 for an outbred, composite-breed male (USMARCv1.0). Both assemblies are of substantially higher (>90-fold) continuity and accuracy than Sscrofa10.2. Conclusions These highly contiguous assemblies plus annotation of a further 11 short-read assemblies provide an unprecedented view of the genetic make-up of this important agricultural and biomedical model species. We propose that the improved Duroc assembly (Sscrofa11.1) become the reference genome for genomic research in pigs.

Sperm DNA methylation is crucial for fertility and viability of offspring but epigenome evolution in mammals is largely understudied. By comparing sperm DNA methylomes and large-scale genome-wide association study (GWAS) signals between human and cattle, we aimed to examine the DNA methylome evolution and its associations with complex phenotypes in mammals. Our analysis revealed that genes with conserved non-methylated promoters (e.g., ANKS1A and WNT7A) among human and cattle were involved in common system and embryo development, and enriched for GWAS signals of body conformation traits in both species, while genes with conserved hypermethylated promoters (e.g., TCAP and CD80) were engaged in immune responses and highlighted by immune-related traits. On the other hand, genes with human-specific hypomethylated promoters (e.g., FOXP2 and HYDIN) were engaged in neuron system development and enriched for GWAS signals of brain-related traits, while genes with cattle-specific hypomethylated promoters (e.g., LDHB and DGAT2) mainly participated in lipid storage and metabolism. We validated our findings using sperm-retained nucleosome, preimplantation transcriptome, and adult tissue transcriptome data, as well as sequence evolutionary features, including motif binding sites, mutation rates, recombination rates and evolution signatures. In conclusion, our results demonstrate important roles of epigenome evolution in shaping the genetic architecture underlying complex phenotypes, hence enhance signal prioritization in GWAS and provide valuable information for human neurological disorders and livestock genetic improvement.

We estimated average genetic diversity of twoGlycineannual and six perennial species based upon 76 orthologous gene sets and performed phylogenetic analysis, divergence analysis and tests for departure from neutrality of the eight species using 52 orthologous gene sets. In addition, 367 orthologous gene sets were used to estimate the relationships of 11G. canescensaccessions. Among the perennials,G. canescensshowed the highest nucleotide diversity and the other perennials exceptG. tomentellahad higher nucleotide diversity than the two annuals. Phylogenetic analysis of theGlycineshowed a similar genome grouping with the previous report exceptG. cyrtolobaandG. stenophitaformed a sister clade in the study. Divergence analysis supported the phylogenetic relationships,G. falcatawas the most divergent fromG. max, followed byG. cyrtoloba,G. syndetika,G. tomentellaD3,G. stenophitaandG. canescens. Neutrality selection tests within species showed that most genes were subjected to a recent directional selection due to a selective sweep or rapid population expansion. Although most gene sequence had negative and significant Tajima’s D, the sequences were homogeneous in the levels of polymorphism and divergence betweenG. maxand otherGlycinespecies based on the HKA test, thus,Glycineperennials may have experienced very similar evolutionary selection as inferred by trans-specific mutation analysis. The greater genetic diversity of most perennialGlycinespecies and their origins from the warmer and drier climates of Australia suggested the perennials as potential sources of heat and drought resistance that will be of value in the face of climate change.

Aeromonas bestiarum is a Gram-negative mesophilic motile bacterium causing acute hemorrhagic septicemia or chronic skin ulcers in fish. Here, we report the draft genome sequence of A. bestiarum strain GA97-22, which was isolated from rainbow trout in 1997. This genome sequence will improve our understanding of the complex taxonomy of motile aeromonads.

Background Copy number variation (CNV) is an important type of genetic variation contributing to phenotypic differences among mammals and may serve as an alternative molecular marker to single nucleotide polymorphism (SNP) for genome-wide association study (GWAS). Recently, GWAS analysis using CNV has been applied in livestock, although few studies have focused on Holstein cattle. Results We describe 191 CNV detected using intensity data from over 700,000 SNP genotypes generated with the BovineHD Genotyping BeadChip (Illumina, San Diego, CA) in 528 Holstein cows. The CNV were used for GWAS analysis of 10 important production traits of 473 cattle related to feed intake, milk quality, and female fertility, as well as 2 composite traits of net merit and productive life. In total, we detected 57 CNV associated (P

Background Although sperm DNA methylation has been studied in humans and other species, its status in cattle is largely unknown. Results Using whole-genome bisulfite sequencing (WGBS), we profiled the DNA methylome of cattle sperm through comparison with three somatic tissues (mammary gland, brain, and blood). Large differences between cattle sperm and somatic cells were observed in the methylation patterns of global CpGs, pericentromeric satellites, partially methylated domains (PMDs), hypomethylated regions (HMRs), and common repeats. As expected, we observed low methylation in the promoter regions and high methylation in the bodies of active genes. We detected selective hypomethylation of megabase domains of centromeric satellite clusters, which may be related to chromosome segregation during meiosis and their rapid transcriptional activation upon fertilization. We found more PMDs in sperm cells than in somatic cells and identified meiosis-related genes such asKIF2B and REPIN1, which are hypomethylated in sperm but hypermethylated in somatic cells. In addition to the common HMRs around gene promoters, which showed substantial differences between sperm and somatic cells, the sperm-specific HMRs also targeted to distinct spermatogenesis-related genes, including BOLL, MAEL, ASZ1, SYCP3, CTCFL, MND1, SPATA22, PLD6, DDX4, RBBP8, FKBP6, and SYCE1. Although common repeats were heavily methylated in both sperm and somatic cells, some young Bov-A2 repeats, which belong to the SINE family, were hypomethylated in sperm and could affect the promoter structures by introducing new regulatory elements. Conclusions Our study provides a comprehensive resource for bovine sperm epigenomic research and enables new discoveries about DNA methylation and its role in male fertility.

Made available in DSpace on 2018-12-11T17:16:50Z (GMT). No. of bitstreams: 0 Previous issue date: 2017-12-01 The recent evolution of cattle is marked by fluctuations in body size. Height in the Bos taurus lineage was reduced by a factor of ∼1.5 from the Neolithic to the Middle Ages, and increased again only during the Early Modern Ages. Using haplotype analysis, we found evidence that the bovine PLAG1 mutation (Q) with major effects on body size, weight and reproduction is a >1,000 years old derived allele that increased rapidly in frequency in Northwestern European B. taurus between the 16th and 18th centuries. Towards the 19th and 20th centuries, Q was introgressed into non-European B. taurus and Bos indicus breeds. These data implicate a major role of Q in recent changes in body size in modern cattle, and represent one of the first examples of a genomic sweep in livestock that was driven by selection on a complex trait. São Paulo State University (Unesp) School of Agricultural and Veterinarian Sciences Jaboticabal Department of Preventive Veterinary Medicine and Animal Reproduction International Atomic Energy Agency (IAEA) Collaborating Centre on Animal Genomics and Bioinformatics Istituto di Zootecnica BioDNA Centro di Ricerca sulla Biodiversità e sul DNA Antico Università Cattolica del Sacro Cuore São Paulo State University (Unesp) School of Veterinary Medicine Araçatuba Department of Support Production and Animal Health Recombinetics Inc. UFPI - Universidade Federal do Piauí Animal Genomics and Improvement Laboratory Agricultural Research Service USDA Escola de Veterinária e Zootecnia (EVZ) Universidade Federal de Goiás - UFG Campus Samambaia São Paulo State University (Unesp) School of Agricultural and Veterinarian Sciences Jaboticabal Department of Animal Science GenSys Consultores Associados Centro Universitário Católico Salesiano Department of Anthropology University of Maryland Dipartimento di Storia Archeologia e Storia dell'Arte Facoltà di Lettere e Filosofia Università Cattolica del Sacro Cuore Società Friulana di Archeologia Istituto di Genomica Applicata Dipartamento di Scienze Agrarie ed Ambientali Università di Udine BOKU University of Natural Resources and Life Sciences Department of Sustainable Agricultural Systems Division of Livestock Sciences São Paulo State University (Unesp) School of Agricultural and Veterinarian Sciences Jaboticabal Department of Preventive Veterinary Medicine and Animal Reproduction São Paulo State University (Unesp) School of Veterinary Medicine Araçatuba Department of Support Production and Animal Health São Paulo State University (Unesp) School of Agricultural and Veterinarian Sciences Jaboticabal Department of Animal Science

To elucidate the molecular mechanism of host resistance, we characterized the jejunal transcriptome of Angus cattle selected for parasite resistance for over 20 years in response to infection caused by the intestinal worm Cooperia oncophora. The transcript abundance of 56 genes, such as that of mucin 12 (MUC12) and intestinal alkaline phosphatase (ALPI), was significantly higher in resistant cattle. Novel splicing variants, exon skipping events, and gene fusion events, were also detected. An algorithm for the reconstruction of accurate cellular networks (ARACNE) was used to infer de novo regulatory molecular networks in the interactome between the parasite and host. Under a combined cutoff of an error tolerance (ϵ = 0.10) and a stringent P-value threshold of mutual information (1.0 × 10(-5)), a total of 229,100 direct interactions controlled by 20,288 hub genes were identified. Among these hub genes, 7651 genes had ≥ 100 direct neighbors while the top 9778 hub genes controlled more than 50% of total direct interactions. Three lysozyme genes (LYZ1, LYZ2, and LYZ3), which are co-located in bovine chromosome 5 in tandem and are strongly upregulated in resistant cattle, shared a common regulatory network of 55 genes. These ancient antimicrobials were likely involved in regulating host-parasite interactions by affecting host gut microbiome. Notably, ALPI, known as a gut mucosal defense factor, controlled a molecular network consisting 410 genes, including 14 transcription factors (TF) and 10 genes that were significantly regulated in resistant cattle. Several large regulatory networks were controlled by TF, such as STAT6, SREBF1, and ELF4. Gene ontology (GO) processes significantly enriched in the regulatory network controlled by STAT6 included lipid metabolism. Our findings provide insights into the immune regulation of host-parasite interactions and the molecular mechanisms of host resistance in cattle.

Motivation: Identification of structural variants (SVs) in sequence data results in a large number of false positive calls using existing software, which overburdens subsequent validation. Results: Simulations using RAPTR-SV and other, similar algorithms for SV detection revealed that RAPTR-SV had superior sensitivity and precision, as it recovered 66.4% of simulated tandem duplications with a precision of 99.2%. When compared with calls made by Delly and LUMPY on available datasets from the 1000 genomes project, RAPTR-SV showed superior sensitivity for tandem duplications, as it identified 2-fold more duplications than Delly, while making ∼85% fewer duplication predictions. Availability and implementation: RAPTR-SV is written in Java and uses new features in the collections framework in the latest release of the Java version 8 language specifications. A compiled version of the software, instructions for usage and test results files are available on the GitHub repository page: https://github.com/njdbickhart/RAPTR-SV. Contact: derek.bickhart@ars.usda.gov

We investigated diverse genomic selections using high-densit ys ingle nucleotide polymorphism data of fi ve distinct cattle breeds. Based on allele frequency differences, we detected hundreds of candidate regions under positive selection across Holstein, Angus, Charolais, Brahman, and N’Dama. In addition to well-known genes such as KIT, MC1R, ASIP, GHR, LCORL, NCAPG, WIF1 ,a ndABCA12, we found evidence for a variety of novel and less-known genes under selection in cattle, such as LAP3, SAR1B, LRIG3, FGF5 ,a ndNUDCD3. Selective sweeps near LAP3 were then validated by next-generation sequencing. Genome-wide association analysis involving 26,362 Holsteins confirmed that LAP3 and SAR1B were related to milk production traits, suggesting that our candidate regions were likely functional. In addition, haplotype network analyses further revealed distinct selective pressures and evolution patterns across these five cattle breeds. Our results provided a glimpse into diverse genomic selection during cattle domestication, breed formation, and recent genetic improvement. These findings will facilitate genome-assisted breeding to improve animal production and health.

The decrease in sequencing cost and increased sophistication of assembly algorithms for short-read platforms has resulted in a sharp increase in the number of species with genome assemblies. However, these assemblies are highly fragmented, with many gaps, ambiguities, and errors, impeding downstream applications. We demonstrate current state of the art for de novo assembly using the domestic goat (Capra hircus) based on long reads for contig formation, short reads for consensus validation, and scaffolding by optical and chromatin interaction mapping. These combined technologies produced what is, to our knowledge, the most continuous de novo mammalian assembly to date, with chromosome-length scaffolds and only 649 gaps. Our assembly represents a ∼400-fold improvement in continuity due to properly assembled gaps, compared to the previously published C. hircus assembly, and better resolves repetitive structures longer than 1 kb, representing the largest repeat family and immune gene complex yet produced for an individual of a ruminant species.

Aeromonas hydrophila is a Gram-negative bacterium that is particularly adapted to freshwater environments and can cause severe infections in fish and humans. Here, we report the draft genomes of three A. hydrophila catfish and tilapia isolates.

Made available in DSpace on 2018-12-11T16:49:56Z (GMT). No. of bitstreams: 0 Previous issue date: 2017-10-01 Ministero dell’Istruzione, dell’Università e della Ricerca Background: The availability of the bovine genome sequence and SNP panels has improved various genomic analyses, from exploring genetic diversity to aiding genetic selection. However, few of the SNP on the bovine chips are polymorphic in buffalo, therefore a panel of single nucleotide DNA markers exclusive for buffalo was necessary for molecular genetic analyses and to develop genomic selection approaches for water buffalo. The creation of a 90K SNP panel for river buffalo and testing in a genome wide association study for milk production is described here. Methods: The genomes of 73 buffaloes of 4 different breeds were sequenced and aligned against the bovine genome, which facilitated the identification of 22 million of sequence variants among the buffalo genomes. Based on frequencies of variants within and among buffalo breeds, and their distribution across the genome, inferred from the bovine genome sequence, 90,000 putative single nucleotide polymorphisms were selected to create an Axiom® Buffalo Genotyping Array 90K. Results: This 90K “SNP-Chip” was tested in several river buffalo populations and found to have *70% high quality and polymorphic SNPs. Of the 90K SNPs about 24K were also found to be polymorphic in swamp buffalo. The SNP chip was used to investigate the structure of buffalo populations, and could distinguish buffalo from different farms. A Genome Wide Association Study identified genomic regions on 5 chromosomes putatively involved in milk production. Conclusion: The 90K buffalo SNP chip described here is suitable for the analysis of the genomes of river buffalo breeds, and could be used for genetic diversity studies and potentially as a starting point for genome-assisted selection programmes. This SNP Chip could also be used to analyse swamp buffalo, but many loci are not informative and creation of a revised SNP set specific for swamp buffalo would be advised. AIA-LGS Associazione Italiana Allevatori–Laboratorio Genetica e Servizi Fondazione Parco Tecnologico Padano Bioinformatics core Parco Tecnologico Padano Animal Genomics and Improvement Laboratory BARC USDA-ARS Department of Animal Science Iowa State University Institute of Zootechnics Università Cattolica del S. Cuore ANASB-Associazione Nazionale Allevatori Specie Bufalina Universidade Estadual Paulista (UNESP) Câmpus de Araçatuba COMSATS Institute of Information Technology Dipartimento di Scienze Zootecniche ed Ispezione degli Alimenti Facoltà di Agraria Università degli Studi di Napoli Federico II Universidade Federal De Minas Gerais (UFMG) Affymetrix UK Ltd Davies Research Centre School of Animal and Veterinary Sciences University of Adelaide Department of Animal Science University of Arkansas Recombinetics Inc. Council on Dairy Cattle Breeding, 4201 Northview Dr. One Town Centre Universidade Estadual Paulista (UNESP) Câmpus de Araçatuba

Decreased male fertility is a big concern in both human society and the livestock industry. Sperm DNA methylation is commonly believed to be associated with male fertility. However, due to the lack of accurate male fertility records (i.e., limited mating times), few studies have investigated the comprehensive impacts of sperm DNA methylation on male fertility in mammals. In this study, we generated 10 sperm DNA methylomes and performed a preliminary correlation analysis between signals from sperm DNA methylation and signals from large-scale (n = 27,214) genome-wide association studies (GWAS) of 35 complex traits (including 12 male fertility-related traits). We detected genomic regions, which experienced DNA methylation alterations in sperm and were associated with aging and extreme fertility phenotypes (e.g., sire-conception rate or SCR). In dynamic hypomethylated regions (HMRs) and partially methylated domains (PMDs), we found genes (e.g., HOX gene clusters and microRNAs) that were involved in the embryonic development. We demonstrated that genomic regions, which gained rather than lost methylations during aging, and in animals with low SCR were significantly and selectively enriched for GWAS signals of male fertility traits. Our study discovered 16 genes as the potential candidate markers for male fertility, including SAMD5 and PDE5A. Collectively, this initial effort supported a hypothesis that sperm DNA methylation may contribute to male fertility in cattle and revealed the usefulness of functional annotations in enhancing biological interpretation and genomic prediction for complex traits and diseases.

Since 2009, a clonal group of virulent Aeromonas hydrophila strains has been causing severe disease in the catfish aquaculture industry in the southeastern United States. Here, we report draft genomes of four A. hydrophila isolates from catfish aquaculture that represent this clonal group.

The decrease in sequencing cost and increased sophistication of assembly algorithms for short-read platforms has resulted in a sharp increase in the number of species with genome assemblies. However, these assemblies are highly fragmented, with many gaps, ambiguities, and errors, impeding downstream applications. We demonstrate current state of the art for de novo assembly using the domestic goat (Capra hircus), based on long reads for contig formation, short reads for consensus validation, and scaffolding by optical and chromatin interaction mapping. These combined technologies produced the most contiguous de novo mammalian assembly to date, with chromosome-length scaffolds and only 663 gaps. Our assembly represents a >250-fold improvement in contiguity compared to the previously published C. hircus assembly, and better resolves repetitive structures longer than 1 kb, supporting the most complete repeat family and immune gene complex representation ever produced for a ruminant species.

Made available in DSpace on 2018-12-11T17:28:48Z (GMT). No. of bitstreams: 0 Previous issue date: 2016-06-15 The diversity and population genetics of copy number variation (CNV) in domesticated animals are not well understood. In this study, we analysed 75 genomes of major taurine and indicine cattle breeds (including Angus, Brahman, Gir, Holstein, Jersey, Limousin, Nelore, and Romagnola), sequenced to 11-fold coverage to identify 1,853 non-redundant CNV regions. Supported by high validation rates in array comparative genomic hybridization (CGH) and qPCR experiments, these CNV regions accounted for 3.1% (87.5 Mb) of the cattle reference genome, representing a significant increase over previous estimates of the area of the genome that is copy number variable (∼2%). Further population genetics and evolutionary genomics analyses based on these CNVs revealed the population structures of the cattle taurine and indicine breeds and uncovered potential diversely selected CNVs near important functional genes, including AOX1, ASZ1, GAT, GLYAT, and KRTAP9-1. Additionally, 121 CNV gene regions were found to be either breed specific or differentially variable across breeds, such as RICTOR in dairy breeds and PNPLA3 in beef breeds. In contrast, clusters of the PRP and PAG genes were found to be duplicated in all sequenced animals, suggesting that subfunctionalization, neofunctionalization, or overdominance play roles in diversifying those fertility-related genes. These CNV results provide a new glimpse into the diverse selection histories of cattle breeds and a basis for correlating structural variation with complex traits in the future. USDA-ARS Animal Genomics and Improvement Laboratory Department of Animal and Avian Sciences University of Maryland Universidade Estadual Paulista (UNESP) Division of Animal Sciences University of Missouri Informatics Institute University of Missouri Department of Evolution and Ecology University of California Universidade Estadual Paulista (UNESP)

Made available in DSpace on 2018-12-11T17:28:22Z (GMT). No. of bitstreams: 0 Previous issue date: 2016-06-01 Background: Apart from single nucleotide polymorphism (SNP), copy number variation (CNV) is another important type of genetic variation, which may affect growth traits and play key roles for the production of beef cattle. To date, no genome-wide association study (GWAS) for CNV and body traits in beef cattle has been reported, so the present study aimed to investigate this type of association in one of the most important cattle subspecies: Bos indicus (Nellore breed). Results: We have used intensity data from over 700,000 SNP probes across the bovine genome to detect common CNVs in a sample of 2230 Nellore cattle, and performed GWAS between the detected CNVs and nine growth traits. After filtering for frequency and length, a total of 231 CNVs ranging from 894bp to 4,855,088bp were kept and tested as predictors for each growth trait using linear regression analysis with principal components correction. There were 49 significant associations identified among 17 CNVs and seven body traits after false discovery rate correction (P < 0.05). Among the 17 CNVs, three were significant or marginally significant for all the traits. We have compared the locations of associated CNVs with quantitative trait locus and the RefGene database, and found two sets of 9 CNVs overlapping with either known QTLs or genes, respectively. The gene overlapping with CNV100, KCNJ12, is a functional candidate for muscle development and plays critical roles in muscling traits. Conclusion: This study presents the first CNV-based GWAS of growth traits using high density SNP microarray data in cattle. We detected 17 CNVs significantly associated with seven growth traits and one of them (CNV100) may be involved in growth traits through KCNJ12. BARC USDA-ARS Animal Genomics and Improvement Laboratory, BARC-East Northwest A and F University College of Animal Science and Technology Faculdade de Ci�ncias Agr�rias e Veterin�rias UNESP - Univ Estadual Paulista Departamento de Medicina Veterin�ria Preventiva e Reprodu��o Animal Institute of Animal Science Chinese Academy of Agricultural Science College of Animal Science and Food Engineering University of S�o Paulo Department of Veterinary Medicine UNESP - Univ Estadual Paulista Departamento de Zootecnia Faculdade de Ci�ncias Agr�rias e Veterin�rias UNESP - Univ Estadual Paulista Departamento de Apoio Produ��o e Sa�de Animal Faculdade de Medicina Veterin�ria de Ara�atuba Recombinetics Inc. International Atomic Energy Agency (IAEA) Collaborating Centre on Animal Genomics and Bioinformatics Faculdade de Ci�ncias Agr�rias e Veterin�rias UNESP - Univ Estadual Paulista Departamento de Medicina Veterin�ria Preventiva e Reprodu��o Animal UNESP - Univ Estadual Paulista Departamento de Zootecnia Faculdade de Ci�ncias Agr�rias e Veterin�rias UNESP - Univ Estadual Paulista Departamento de Apoio Produ��o e Sa�de Animal Faculdade de Medicina Veterin�ria de Ara�atuba

Aeromonas hydrophila is an opportunistic pathogen residing in freshwater environments that causes infection in fish and mammals. Here, we report the draft genome sequence of A. hydrophila strain TN97-08 isolated from a diseased bluegill ( Lepomis macrochirus ) in 1997.

The feasibility of short-read sequencing for genomic analysis was demonstrated for Fibroporia radiculosa , a copper-tolerant fungus that causes brown rot decay of wood. The effect of read quality on genomic assembly was assessed by filtering Illumina GAIIx reads from a single run of a paired-end library (75-nucleotide read length and 300-bp fragment size) at three different stringency levels and then assembling each data set with Velvet. A simple approach was devised to determine which filter stringency was “best.” Venn diagrams identified the regions containing reads that were used in an assembly but were of a low-enough quality to be removed by a filter. By plotting base quality histograms of reads in this region, we judged whether a filter was too stringent or not stringent enough. Our best assembly had a genome size of 33.6 Mb, an N50 of 65.8 kb for a k -mer of 51, and a maximum contig length of 347 kb. Using GeneMark, 9,262 genes were predicted. TargetP and SignalP analyses showed that among the 1,213 genes with secreted products, 986 had motifs for signal peptides and 227 had motifs for signal anchors. Blast2GO analysis provided functional annotation for 5,407 genes. We identified 29 genes with putative roles in copper tolerance and 73 genes for lignocellulose degradation. A search for homologs of these 102 genes showed that F. radiculosa exhibited more similarity to Postia placenta than Serpula lacrymans . Notable differences were found, however, and their involvements in copper tolerance and wood decay are discussed.

Copy number variations (CNVs) affect a wide range of phenotypic traits; however, CNVs in or near segmental duplication regions are often intractable. Using a read depth approach based on next-generation sequencing, we examined genome-wide copy number differences among five taurine (three Angus, one Holstein, and one Hereford) and one indicine (Nelore) cattle. Within mapped chromosomal sequence, we identified 1265 CNV regions comprising ∼55.6-Mbp sequence—476 of which (∼38%) have not previously been reported. We validated this sequence-based CNV call set with array comparative genomic hybridization (aCGH), quantitative PCR (qPCR), and fluorescent in situ hybridization (FISH), achieving a validation rate of 82% and a false positive rate of 8%. We further estimated absolute copy numbers for genomic segments and annotated genes in each individual. Surveys of the top 25 most variable genes revealed that the Nelore individual had the lowest copy numbers in 13 cases (∼52%, χ2 test; P-value CATHL4 and ULBP17, were highly duplicated in the Nelore individual relative to the taurine cattle, while genes involved in lipid transport and metabolism, including APOL3 and FABP2, were highly duplicated in the beef breeds. These CNV regions also harbor genes like BPIFA2A (BSP30A) and WC1, suggesting that some CNVs may be associated with breed-specific differences in adaptation, health, and production traits. By providing the first individualized cattle CNV and segmental duplication maps and genome-wide gene copy number estimates, we enable future CNV studies into highly duplicated regions in the cattle genome.

Background Next-generation sequencing projects commonly commence by aligning reads to a reference genome assembly. While improvements in alignment algorithms and computational hardware have greatly enhanced the efficiency and accuracy of alignments, a significant percentage of reads often remain unmapped. Results We generated de novo assemblies of unmapped reads from the DNA and RNA sequencing of the Bos taurus reference individual and identified the closest matching sequence to each contig by alignment to the NCBI non-redundant nucleotide database using BLAST. As expected, many of these contigs represent vertebrate sequence that is absent, incomplete, or misassembled in the UMD3.1 reference assembly. However, numerous additional contigs represent invertebrate species. Most prominent were several species of Spirurid nematodes and a blood-borne parasite, Babesia bigemina. These species are either not present in the US or are not known to infect taurine cattle and the reference animal appears to have been host to unsequenced sister species. Conclusions We demonstrate the importance of exploring unmapped reads to ascertain sequences that are either absent or misassembled in the reference assembly and for detecting sequences indicative of parasitic or commensal organisms. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2313-7) contains supplementary material, which is available to authorized users.

Background Milk exosomes are a rich source of microRNAs (miRNAs) that are protected from degradation. Ingestion of milk and subsequent absorption of miRNAs into recipient cells by endocytosis may play a role in the regulation of neonatal innate and adaptive immunity. In contrast, the miRNA content of milk exosomes may also be indicative of a lactating animal’s health; whereby, the presence or absence of specific miRNAs could serve as biomarkers for early detection of bacterial infection that can lead to mastitis. In the present study, we therefore analyzed and compared miRNA expression profiles of milk exosomes from four Holstein cows obtained during mid-lactation prior to and after infection (48 h) of the mammary gland with Staphylococcus aureus. Methods Milk exosomes, purified from control and S. aureus infected cows, were extracted for RNA. Following preparation indexed libraries from both groups the samples were subjected to next generation sequencing. Results Next generation sequencing of eight, unpooled small RNA libraries derived from milk exosomes produced about 60.5 million high-quality, bovine-specific sequence reads for comparison of miRNA expression between treatments. Sequence identity analysis showed the miRNAs make up about 13 % of the average RNA content of these exosomes. Although 417 known bovine miRNAs were identified, miRNAs represented the least diverse class of RNA accounting for only 1 % of all unique sequences. The 20 most prevalent unique sequences within this class accounted for about 90 % of the total miRNA-associated reads across samples. Non-annotated, unique reads provided evidence for another 303 previously unknown bovine miRNAs. Expression analyses found 14 known bovine microRNAs significantly differed in frequency between exosomes from infected and control animals. Conclusions Our survey of miRNA expression from uninfected milk exosomes and those produced in response to infection provides new and comprehensive information supporting a role for delivery into milk of specific miRNAs involved in immune response. In particular, bta-miR-142-5p, and −223 are potential biomarkers for early detection of bacterial infection of the mammary gland. Additionally, 22 mammary-expressed genes involved in regulation of host immune processes and response to inflammation were identified as potential binding targets of the differentially expressed miRNAs. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2044-9) contains supplementary material, which is available to authorized users.

A total of 992,682 single-nucleotide polymorphisms (SNPs) was identified as ideal for Illumina Infinium II BeadChip design after sequencing a diverse set of 17 common bean (Phaseolus vulgaris L) varieties with the aid of next-generation sequencing technology. From these, two BeadChips each with >5000 SNPs were designed. The BARCBean6K_1 BeadChip was selected for the purpose of optimizing polymorphism among market classes and, when possible, SNPs were targeted to sequence scaffolds in the Phaseolus vulgaris 14× genome assembly with sequence lengths >10 kb. The BARCBean6K_2 BeadChip was designed with the objective of anchoring additional scaffolds and to facilitate orientation of large scaffolds. Analysis of 267 F2 plants from a cross of varieties Stampede × Red Hawk with the two BeadChips resulted in linkage maps with a total of 7040 markers including 7015 SNPs. With the linkage map, a total of 432.3 Mb of sequence from 2766 scaffolds was anchored to create the Phaseolus vulgaris v1.0 assembly, which accounted for approximately 89% of the 487 Mb of available sequence scaffolds of the Phaseolus vulgaris v0.9 assembly. A core set of 6000 SNPs (BARCBean6K_3 BeadChip) with high genotyping quality and polymorphism was selected based on the genotyping of 365 dry bean and 134 snap bean accessions with the BARCBean6K_1 and BARCBean6K_2 BeadChips. The BARCBean6K_3 BeadChip is a useful tool for genetics and genomics research and it is widely used by breeders and geneticists in the United States and abroad.

Background Domestication of the now-extinct wild aurochs, Bos primigenius, gave rise to the two major domestic extant cattle taxa, B. taurus and B. indicus. While previous genetic studies have shed some light on the evolutionary relationships between European aurochs and modern cattle, important questions remain unanswered, including the phylogenetic status of aurochs, whether gene flow from aurochs into early domestic populations occurred, and which genomic regions were subject to selection processes during and after domestication. Here, we address these questions using whole-genome sequencing data generated from an approximately 6,750-year-old British aurochs bone and genome sequence data from 81 additional cattle plus genome-wide single nucleotide polymorphism data from a diverse panel of 1,225 modern animals. Results Phylogenomic analyses place the aurochs as a distinct outgroup to the domestic B. taurus lineage, supporting the predominant Near Eastern origin of European cattle. Conversely, traditional British and Irish breeds share more genetic variants with this aurochs specimen than other European populations, supporting localized gene flow from aurochs into the ancestors of modern British and Irish cattle, perhaps through purposeful restocking by early herders in Britain. Finally, the functions of genes showing evidence for positive selection in B. taurus are enriched for neurobiology, growth, metabolism and immunobiology, suggesting that these biological processes have been important in the domestication of cattle. Conclusions This work provides important new information regarding the origins and functional evolution of modern cattle, revealing that the interface between early European domestic populations and wild aurochs was significantly more complex than previously thought. Electronic supplementary material The online version of this article (doi:10.1186/s13059-015-0790-2) contains supplementary material, which is available to authorized users.