Christopher M. Seabury

Resolving the evolution of extant and extinct ruminants with high-throughput phylogenomics

The Pecorans (higher ruminants) are believed to have rapidly speciated in the Mid-Eocene, resulting in five distinct extant families: Antilocapridae, Giraffidae, Moschidae, Cervidae, and Bovidae. Due to the rapid radiation, the Pecoran phylogeny has proven difficult to resolve, and 11 of the 15 possible rooted phylogenies describing ancestral relationships among the Antilocapridae, Giraffidae, Cervidae, and Bovidae have each been argued as representations of the true phylogeny. Here we demonstrate that a genome-wide single nucleotide polymorphism (SNP) genotyping platform designed for one species can be used to genotype ancient DNA from an extinct species and DNA from species diverged up to 29 million years ago and that the produced genotypes can be used to resolve the phylogeny for this rapidly radiated infraorder. We used a high-throughput assay with 54,693 SNP loci developed for Bos taurus taurus to rapidly genotype 678 individuals representing 61 Pecoran species. We produced a highly resolved phylogeny for this diverse group based upon 40,843 genome-wide SNP, which is five times as many informative characters as have previously been analyzed. We also establish a method to amplify and screen genomic information from extinct species, and place Bison priscus within the Bovidae. The quality of genotype calls and the placement of samples within a well-supported phylogeny may provide an important test for validating the fidelity and integrity of ancient samples. Finally, we constructed a phylogenomic network to accurately describe the relationships between 48 cattle breeds and facilitate inferences concerning the history of domestication and breed formation.

Worldwide Patterns of Ancestry, Divergence, and Admixture in Domesticated Cattle

The domestication and development of cattle has considerably impacted human societies, but the histories of cattle breeds and populations have been poorly understood especially for African, Asian, and American breeds. Using genotypes from 43,043 autosomal single nucleotide polymorphism markers scored in 1,543 animals, we evaluate the population structure of 134 domesticated bovid breeds. Regardless of the analytical method or sample subset, the three major groups of Asian indicine, Eurasian taurine, and African taurine were consistently observed. Patterns of geographic dispersal resulting from co-migration with humans and exportation are recognizable in phylogenetic networks. All analytical methods reveal patterns of hybridization which occurred after divergence. Using 19 breeds, we map the cline of indicine introgression into Africa. We infer that African taurine possess a large portion of wild African auroch ancestry, causing their divergence from Eurasian taurine. We detect exportation patterns in Asia and identify a cline of Eurasian taurine/indicine hybridization in Asia. We also identify the influence of species other than Bos taurus taurus and B. t. indicus in the formation of Asian breeds. We detect the pronounced influence of Shorthorn cattle in the formation of European breeds. Iberian and Italian cattle possess introgression from African taurine. American Criollo cattle originate from Iberia, and not directly from Africa with African ancestry inherited via Iberian ancestors. Indicine introgression into American cattle occurred in the Americas, and not Europe. We argue that cattle migration, movement and trading followed by admixture have been important forces in shaping modern bovine genomic variation.

Production of hornless dairy cattle from genome-edited cell lines

1Department of Biochemistry, University of Utah School of Medicine, Salt Lake City, Utah, USA. 2Department of Animal Science, University of California, Davis, California, USA. 3Division of Animal Sciences, University of Missouri, Columbia, Missouri, USA. 4Department of Biology, University of Utah, Salt Lake City, Utah, USA. 5Department of Genetics, Cell Biology, and Development, University of Minnesota, Minneapolis, Minnesota, USA. 6Center for Genome Engineering, University of Minnesota, Minneapolis, Minnesota, USA. e-mail: dana@biochem.utah.edu

Analysis of sequence variability and protein domain architectures for bovine peptidoglycan recognition protein 1 and Toll-like receptors 2 and 6

The mammalian Toll-like receptors (TLRs) recognize invading pathogens, thereafter provoking innate immune responses, whereas peptidoglycan recognition protein 1 (PGLYRP1) is directly microbicidal. The primary objective of this study was to characterize single-nucleotide polymorphisms (SNPs) and insertion-deletion polymorphisms (indels) within bovine TLR2, TLR6, and PGLYRP1, thereby facilitating future TLR signaling, association, and PGLYRP1 microbicidal assays relevant to bovine innate immunity. Comparative sequence analysis for 10 bovine breeds revealed 83 polymorphisms (82 SNPs, 1 indel), with 15 nonsynonymous SNPs located within predicted functional domains. Of the 83 polymorphisms detected, 72 (87%) are reported here for the first time. Several predicted amino acid replacements encoded by bovine TLR2 and TLR6, but not PGLYRP1, resulted in the confident prediction of protein domain alterations. Prediction and comparison of protein domain architectures for TLR2 and TLR6 revealed six regions of leucine-rich-repeat patterning that was conserved among multiple species. Collectively, differences in the patterns and frequencies of polymorphism were noted between bovine TLRs that predominantly recognize viral ligands (TLRs 3, 7, 8, 9) and those that recognize microbial and/or unknown ligands (TLRs 1, 2, 5, 6, 10).

Diversity and evolution of 11 innate immune genes in Bos taurus taurus and Bos taurus indicus cattle

The Toll-like receptor (TLR) and peptidoglycan recognition protein 1 (PGLYRP1) genes play key roles in the innate immune systems of mammals. While the TLRs recognize a variety of invading pathogens and induce innate immune responses, PGLYRP1 is directly microbicidal. We used custom allele-specific assays to genotype and validate 220 diallelic variants, including 54 nonsynonymous SNPs in 11 bovine innate immune genes (TLR1-TLR10, PGLYRP1) for 37 cattle breeds. Bayesian haplotype reconstructions and median joining networks revealed haplotype sharing between Bos taurus taurus and Bos taurus indicus breeds at every locus, and we were unable to differentiate between the specialized B. t. taurus beef and dairy breeds, despite an average polymorphism density of one locus per 219 bp. Ninety-nine tagSNPs and one tag insertion-deletion polymorphism were sufficient to predict 100% of the variation at all 11 innate immune loci in both subspecies and their hybrids, whereas 58 tagSNPs captured 100% of the variation at 172 loci in B. t. taurus. PolyPhen and SIFT analyses of nonsynonymous SNPs encoding amino acid replacements indicated that the majority of these substitutions were benign, but up to 31% were expected to potentially impact protein function. Several diversity-based tests provided support for strong purifying selection acting on TLR10 in B. t. taurus cattle. These results will broadly impact efforts related to bovine translational genomics.

Comparative PRNP genotyping of U.S. cattle sires for potential association with BSE

The recent discovery of significant associations between bovine spongiform encephalopathy (BSE) susceptibility in German cattle and the frequency distributions of insertion/deletion (indel) polymorphisms within the bovine PRNP gene prompted an evaluation of 132 commercial U.S. artificial insemination (AI) sires from 39 breeds. Forward primer sequences from published primer sets targeting indels within the putative bovine PRNP promoter, intron 1, and the 3′ UTR (untranslated region) were synthesized with unique 5′ fluorescent labels and utilized to develop a rapid multiplexed PCR assay for identifying BSE-associated indels as well as facilitating polymorphism analyses and/or marker-assisted selection. Significant differences (p < 0.05 all tests) were detected between the frequencies of bovine PRNP promoter alleles for 48 healthy German cattle previously described and 132 commercial U.S. cattle sires. The frequency of the 23-bp promoter allele observed for commercial U.S. cattle sires strongly resembled that recently described for 43 BSE-affected German cattle. No significant difference (p = 0.051) was detected between the distributions of promoter genotypes for healthy German cattle and our panel of commercial U.S. cattle sires. Interestingly, significant differences (p < 0.01; p < 0.02) were also noted between the frequencies and distributions of intron 1 alleles and genotypes, respectively, for BSE-affected German cattle and our panel of U.S. cattle sires. No significant allelic or genotypic differences were detected for the 14-bp 3′ UTR indel for any given comparison between German cattle and commercial U.S. cattle sires.

A 4,103 marker integrated physical and comparative map of the horse genome

A comprehensive second-generation whole genome radiation hybrid (RH II), cytogenetic and comparative map of the horse genome (2n = 64) has been developed using the 5000rad horse × hamster radiation hybrid panel and fluorescence in situ hybridization (FISH). The map contains 4,103 markers (3,816 RH; 1,144 FISH) assigned to all 31 pairs of autosomes and the X chromosome. The RH maps of individual chromosomes are anchored and oriented using 857 cytogenetic markers. The overall resolution of the map is one marker per 775 kilobase pairs (kb), which represents a more than five-fold improvement over the first-generation map. The RH II incorporates 920 markers shared jointly with the two recently reported meiotic maps. Consequently the two maps were aligned with the RH II maps of individual autosomes and the X chromosome. Additionally, a comparative map of the horse genome was generated by connecting 1,904 loci on the horse map with genome sequences available for eight diverse vertebrates to highlight regions of evolutionarily conserved syntenies, linkages, and chromosomal breakpoints. The integrated map thus obtained presents the most comprehensive information on the physical and comparative organization of the equine genome and will assist future assemblies of whole genome BAC fingerprint maps and the genome sequence. It will also serve as a tool to identify genes governing health, disease and performance traits in horses and assist us in understanding the evolution of the equine genome in relation to other species.

Susceptibility loci revealed for bovine respiratory disease complex in pre-weaned holstein calves.

BackgroundBovine respiratory disease complex (BRDC) is an infectious disease of cattle that is caused by a combination of viral and/or bacterial pathogens. Selection for cattle with reduced susceptibility to respiratory disease would provide a permanent tool for reducing the prevalence of BRDC. The objective of this study was to identify BRDC susceptibility loci in pre-weaned Holstein calves as a prerequisite to using genetic improvement as a tool for decreasing the prevalence of BRDC. High density SNP genotyping with the Illumina BovineHD BeadChip was conducted on 1257 male and 757 female Holstein calves from California (CA), and 767 calves identified as female from New Mexico (NM). Of these, 1382 were classified as BRDC cases, and 1396 were classified as controls, with all phenotypes assigned using the McGuirk health scoring system. During the acquisition of blood for DNA isolation, two deep pharyngeal and one mid-nasal diagnostic swab were obtained from each calf for the identification of bacterial and viral pathogens. Genome-wide association analyses were conducted using four analytical approaches (EIGENSTRAT, EMMAX-GRM, GBLUP and FvR). The most strongly associated SNPs from each individual analysis were ranked and evaluated for concordance. The heritability of susceptibility to BRDC in pre-weaned Holstein calves was estimated.ResultsThe four statistical approaches produced highly concordant results for 373 top ranked SNPs that defined 126 chromosomal regions for the CA population. Similarly, in NM, 370 SNPs defined 138 genomic regions that were identified by all four approaches. When the two populations were combined (i.e., CA + NM) and analyzed, 324 SNPs defined 116 genomic regions that were associated with BRDC across all analytical methods. Heritability estimates for BRDC were 21% for both CA and NM as individual populations, but declined to 13% when the populations were combined.ConclusionsFour analytical approaches utilizing both single and multi-marker association methods revealed common genomic regions associated with BRDC susceptibility that can be further characterized and used for genomic selection. Moderate heritability estimates were observed for BRDC susceptibility in pre-weaned Holstein calves, thereby supporting the application of genomic selection to reduce the prevalence of BRDC in U.S. Holsteins.

Quantum dots trigger immunomodulation of the NFκB pathway in human skin cells.

The immunological effects of quantum dots are dependent on a variety of factors including, but not limited to, exposure time and dosing concentrations. In this study, we investigated the influence of 15 nm CdSe/ZnS-COOH quantum dot nanocrystals (QDs) on cell density, viability, and morphology in human epidermal keratinocytes (HEK) and human dermal fibroblasts (HDF). Furthermore, inflammatory and non-inflammatory immune responses were measured using protein and real time PCR array analysis from HDF cells exposed to predetermined sub-lethal concentrations of QDs. CdSe/ZnS-COOH QDs caused concentration-dependent (1-120 nM exposure concentrations) and time-dependent (8 h or 48 h) cell death, as evidenced by metabolic activity and morphological changes. QD exposure induced upregulation of apoptotic, inflammatory and immunoregulatory proteins such as TNF-α, IL-1B and IL-10. HMOX1, an indicator of stress due to reactive oxygen intermediates (ROIs) and/or metals, was upregulated at the later time point as well. QDs also caused modulation of genes known to be associated with inflammatory (IL1-β, CCL2, IRAK-2), immune (IL-1, IL-6, PGLYRP1, SERPINA1, IL-10), stress due to ROIs and/or heavy metals (HMOX1), and apoptotic (CASP1, ADORA2A) responses. Cellular effects from QD exposure were found to primarily follow the NFκB pathway. In addition, QDs induced a differential cytotoxicity in keratinocytes and fibroblasts at different exposure concentrations and time points, even at physiologically relevant dosing concentrations, thus emphasizing the need to investigate potential mechanisms of action among different cell types within the same target organ.

A Multi-Platform Draft de novo Genome Assembly and Comparative Analysis for the Scarlet Macaw (Ara macao)

Data deposition to NCBI Genomes This Whole Genome Shotgun project has been deposited at DDBJ/EMBL/GenBank under the accession AMXX00000000 (SMACv1.0, unscaffolded genome assembly). The version described in this paper is the first version (AMXX01000000). The scaffolded assembly (SMACv1.1) has been deposited at DDBJ/EMBL/GenBank under the accession AOUJ00000000, and is also the first version (AOUJ01000000). Strong biological interest in traits such as the acquisition and utilization of speech, cognitive abilities, and longevity catalyzed the utilization of two next-generation sequencing platforms to provide the first-draft de novo genome assembly for the large, new world parrot Ara macao (Scarlet Macaw). Despite the challenges associated with genome assembly for an outbred avian species, including 951,507 high-quality putative single nucleotide polymorphisms, the final genome assembly (>1.035 Gb) includes more than 997 Mb of unambiguous sequence data (excluding N’s). Cytogenetic analyses including ZooFISH revealed complex rearrangements associated with two scarlet macaw macrochromosomes (AMA6, AMA7), which supports the hypothesis that translocations, fusions, and intragenomic rearrangements are key factors associated with karyotype evolution among parrots. In silico annotation of the scarlet macaw genome provided robust evidence for 14,405 nuclear gene annotation models, their predicted transcripts and proteins, and a complete mitochondrial genome. Comparative analyses involving the scarlet macaw, chicken, and zebra finch genomes revealed high levels of nucleotide-based conservation as well as evidence for overall genome stability among the three highly divergent species. Application of a new whole-genome analysis of divergence involving all three species yielded prioritized candidate genes and noncoding regions for parrot traits of interest (i.e., speech, intelligence, longevity) which were independently supported by the results of previous human GWAS studies. We also observed evidence for genes and noncoding loci that displayed extreme conservation across the three avian lineages, thereby reflecting their likely biological and developmental importance among birds.