Using online tools at the Bovine Genome Database to manually annotate genes in the new reference genome.

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.

Origins of cattle on Chirikof Island, Alaska, elucidated from genome-wide SNP genotypes.

Feral livestock may harbor genetic variation of commercial, scientific, historical or esthetic value. The origins and uniqueness of feral cattle on Chirikof Island, Alaska, are uncertain. The island is now part of the Alaska Maritime Wildlife Refuge and Federal wildlife managers want grazing to cease, presumably leading to demise of the cattle. Here we characterize the cattle of Chirikof Island relative to extant breeds and discern their origins. Our analyses support the inference that Yakut cattle from Russia arrived first on Chirikof Island, then ~120 years ago the first European taurine cattle were introduced to the island, and finally a large wave of Hereford cattle were introduced on average 40 years ago. In addition, this mixture of European and East-Asian cattle is unique compared with other North American breeds and we find evidence that natural selection in the relatively harsh environment of Chirikof Island has further impacted their genetic architecture. These results provide an objective basis for decisions regarding conservation of the Chirikof Island cattle.

Development and validation of a fast and automated DNA identification line.

The importance of DNA evidence for gaining investigative leads demands a fast workflow for forensic DNA profiling performed in large volumes. Therefore, we developed software solutions for automated DNA profile analysis, contamination check, major donor inference, DNA database (DDB) comparison and reporting of the conclusions. This represents the Fast DNA IDentification Line (FIDL) and this study describes its development, validation and implementation in criminal casework at the authors' institute. This first implementation regards single donor profiles and major contributors to mixtures. The validation included testing of the software components on their own and examination of the performance of different DDB search strategies. Furthermore, end-to-end testing was performed under three conditions: (1) testing of scenarios that can occur in DNA casework practice, (2) tests using three months of previous casework data, and (3) testing in a casework production environment in parallel to standard casework practices. The same DNA database candidates were retrieved by this automated line as by the manual workflow. The data flow was correct, results were reproducible and robust, results requiring manual analysis were correctly flagged, and reported results were as expected. Overall, we found FIDL valid for use in casework practice in our institute. The results from FIDL are automatically reported within three working days from receiving the trace sample. This includes the time needed for registration of the case, DNA extraction, quantification, polymerase chain reaction and capillary electrophoresis. FIDL itself takes less than two hours from intake of the raw CE data to reporting. Reported conclusions are one of five options: (1) candidate retrieved from DDB, (2) no candidate retrieved from DDB, (3) high evidential value with regards to reference within the case, (4) results require examination of expert, or (5) insufficient amount of DNA obtained to generate a DNA profile. In our current process, the automated report is sent within three working days and a complete report, with confirmation of the FIDL results, and signed by a reporting officer is sent at a later time. The signed report may include additional analyses regarding e.g. minor contributors. The automated report with first case results is quickly available to the police enabling them to act upon the DNA results prior to receiving the full DNA report. This line enables a uniform and efficient manner of handling large numbers of traces and cases and provides high value investigative leads in the early stages of the investigation.

Multi-laboratory validation of DNAxs including the statistical library DNAStatistX.

This study describes a multi-laboratory validation of DNAxs, a DNA eXpert System for the data management and probabilistic interpretation of DNA profiles [1], and its statistical library DNAStatistX to which, besides the organising laboratory, four laboratories participated. The software was modified to read multiple data formats and the study was performed prior to the release of the software to the forensic community. The first exercise explored all main functionalities of DNAxs with feedback on user-friendliness, installation and general performance. Next, every laboratory performed likelihood ratio (LR) calculations using their own dataset and a dataset provided by the organising laboratory. The organising laboratory performed LR calculations using all datasets. The datasets were generated with different STR typing kits or analysis systems and consisted of samples varying in DNA amounts, mixture ratios, number of contributors and drop-out level. Hypothesis sets had the correct, under- and over-assigned number of contributors and true and false donors as person of interest. When comparing the results between laboratories, the LRs were foremost within one unit on log10 scale. The few LR results that deviated more had differences for the parameters estimated by the optimizer within DNAStatistX. Some of these were indicated by failed iteration results, others by a failed model validation, since unrealistic hypotheses were included. When these results that do not meet the quality criteria were excluded, as is in accordance with interpretation guidelines, none of the analyses in the different laboratories yielded a different statement in the casework report. Nonetheless, changes in software parameters were sought that minimized differences in outcomes, which made the DNAStatistX module more robust. Overall, the software was found intuitive, user-friendly and valid for use in multiple laboratories.

Quantitative trait loci with effects on feed efficiency traits in Hereford x composite double backcross populations.

Two half-sib families of backcross progeny were produced by mating F(1) Line 1 Hereford (L1) x composite gene combination (CGC) bulls with L1 and CGC cows. Feed intake and periodic weights were measured for 218 backcross progeny. These progenies were genotyped using 232 microsatellite markers that spanned the 29 BTA. Progeny from L1 and CGC females was analysed separately using composite interval mapping to find quantitative trait loci (QTL) affecting daily dry matter intake (DMI), average daily gain (ADG), feed conversion (FCR) and residual feed intake (RFI). Results from both backcrosses were pooled to find additional QTL. In the backcross to L1, QTL were detected for RFI and DMI on BTA11, FCR on BTA16, and ADG on BTA9. In the backcross to CGC, QTL were detected for RFI on BTA10, FCR on BTA12 and 16 and ADG on BTA15 and 17. After pooling, QTL were detected for RFI on BTA 2, 6, 7, 10, 11, 13 and 16; for FCR on BTA 9, 12, 16, 17 and 21; for ADG on BTA 9, 14, 15, 17; and for DMI on BTA 2, 5, 6, 9, 10, 11, 20 and 23.

DNAxs/DNAStatistX: Development and validation of a software suite for the data management and probabilistic interpretation of DNA profiles.

The data management, interpretation and comparison of sets of DNA profiles can be complex, time-consuming and error-prone when performed manually. This, combined with the growing numbers of genetic markers in forensic identification systems calls for expert systems that can automatically compare genotyping results within (large) sets of DNA profiles and assist in profile interpretation. To that aim, we developed a user-friendly software program or DNA eXpert System that is denoted DNAxs. This software includes features to view, infer and match autosomal short tandem repeat profiles with connectivity to up and downstream software programs. Furthermore, DNAxs has imbedded the 'DNAStatistX' module, a statistical library that contains a probabilistic algorithm to calculate likelihood ratios (LRs). This algorithm is largely based on the source code of the quantitative probabilistic genotyping system EuroForMix [1]. The statistical library, DNAStatistX, supports parallel computing which can be delegated to a computer cluster and enables automated queuing of requested LR calculations. DNAStatistX is written in Java and is accessible separately or via DNAxs. Using true and non-contributors to DNA profiles with up to four contributors, the DNAStatistX accuracy and precision were assessed by comparing the DNAStatistX results to those of EuroForMix. Results were the same up to rare differences that could be attributed to the different optimizers used in both software programs. Implementation of dye specific detection thresholds resulted in larger likelihood values and thus a better explanation of the data used in this study. Furthermore, processing time, robustness of DNAStatistX results and the circumstances under which model validations failed were examined. Finally, guidelines for application of the software are shared as an example. The DNAxs software is future-proof as it applies a modular approach by which novel functionalities can be incorporated.

Detection of germline and somatic copy number variations in cattle.

As a complement to the Bovine HapMap Consortium project, we initiated a systematic study of the copy numbervariation (CNV) within the same cattle population using array comparative genomic hybridization (array CGH). Oligonucleotide CGH arrays were designed and fabricated to cover all chromosomes with an average interval of 6 kb using the latest bovine genome assembly. In the initial screening, three Holstein bulls were selected to represent major paternal lineages of the Holstein breed with some maternal linkages between these lines. Dual-label hybridizations were performed using either Hereford L1 Dominette 01449 or L1 Domino 99375 as reference. The CNVs were represented by gains and losses of normalized fluorescence intensities relative to the reference. The data presented here, for the first time, demonstrated that significant amounts of germline and fewer somatic CNVs exist in cattle, that many CNVs are common both across diverse cattle breeds and among individuals within a breed, and that array CGH is an effective tool to systematically detect bovine CNV. Selected CNVs have been confirmed by independent methods using real-time (RT) PCR. The strategy used in this study, based on genome higher-orderarchitecture variation, is a powerful approach to generating resources for the identification of novel genomic variation and candidate genes for economically important traits.

Comparative and physical mapping of 111 previously reported and 105 new porcine microsatellites.

Here we report radiation hybrid mapping of 105 new porcine microsatellite markers on the IMpRH(7000) radiation hybrid panel. In addition, we searched flanking sequences of these markers, as well as 673 previously reported RH-mapped microsatellite markers, for orthology to human sequences. Eighty-seven new and 111 previously mapped sequences exhibited orthology to human sequences. Using a stringent sequence alignment, 25 microsatellite-flanking sequences were found to be highly similar to genic sequences, whereas 173 were similar to non-genic sequences in the human genome. Five markers were located near known breakpoints of synteny between human and swine.

Potential emigration of Siberian cattle germplasm on Chirikof Island, Alaska.

Feral cattle residing in Chirikof Island, Alaska, are relatively distinct from breeds used in commercial production in North America. However, preliminary evidence suggested that they exhibit substantial genetic relationship with cattle from Yakutian region of Siberia. Thus, our objective was to further elucidate quantify the origins, admixture and divergence of the Chirikof Island cattle relative to cattle from Siberia and USA. Subject animals were genotyped at 15 microsatellite loci. Compared with Turano-Mongolian and North American cattle, Chirikof Island cattle had similar variation, with slightly less observed heterozygosity, fewer alleles per locus and a positive fixation index. Analysis of the genetic distances revealed two primary clusters; one that contained the North American breeds and the Kazakh White head, and a second that contained the Yakutian and Kalmyk breeds, and the Chirikof population. Thus, it is suggested that Chirikof Island cattle may be a composite of British breeds emanating from North America and Turano-Mongolian cattle. A potential founder effect, consistent with historical records of the Russian-American period, may contribute to the adaptation of the Chirikof Island cattle to their harsh high-latitude environment. Further study of adaptive mechanisms manifest by these cattle is warranted.

A microsatellite linkage map of the porcine genome.

We report the most extensive genetic linkage map for a livestock species produced to date. We have linked 376 microsatellite (MS) loci with seven restriction fragment length polymorphic loci in a backcross reference population. The 383 markers were placed into 24 linkage groups which span 1997 cM. Seven additional MS did not fall into a linkage group. Linkage groups are assigned to 13 autosomes and the X chromosome (haploid n = 19). This map provides the basis for genetic analysis of quantitative inheritance of phenotypic and physiologic traits in swine.

The sequence of a porcine cDNA encoding alpha-lactalbumin.

A cDNA clone encoding porcine alpha-lactalbumin (alpha LA) was isolated and sequenced. The longest clone was 688 nucleotides (nt) long and encoded a preprotein of 141 amino acids (aa) including a leader peptide of 19 aa. The porcine cDNA exhibited a nt similarity of between 72.2%-83.5% to other alpha LA cDNAs and an aa similarity of between 50.8%-85.2% with other alpha LA aa sequences. The derived aa sequence varied at three positions from a previously reported sequence for porcine alpha LA obtained by direct aa sequencing.

A conceptual database model for genomic research.

We describe a conceptual model for genome databases that facilitates the process of building, maintaining, and disseminating physically anchored genetic linkage maps. The model has been implemented as a relational database at the Roman L. Hruska U.S. Meat Animal Research Center (MARC). Development of consensus maps using disparate data from different reference pedigrees or laboratories is supported. The model is of use to quantitative and population geneticists interested in loci that affect phenotypes and marker-assisted selection, and it is sufficiently flexible for centralized, species genome databases facilitating comparative mapping. The MARC genome database is used to assemble, maintain, and disseminate physically anchored genetic linkage maps for cattle, swine, and sheep currently based on more than 100,000 genotypes from 1,000 markers. Integrated with linkage analysis software, this database permits frequent updates of physically anchored genetic linkage maps.

Evolution of the trappin multigene family in the Suidae.

Trappins are a group of secretory proteins containing a WAP motif with an anchoring domain. Previous studies showed that their genes, especially those of pig, have undergone rapid evolution, which produced trappins with a broad spectrum of actions. To understand the evolution of such a rapidly evolving multigene family, we isolated trappin genes of the Artiodactyla, including pig, wart hog, collared peccary, hippopotamus, and cow, by means of polymerase chain reaction (PCR). Two genes newly isolated from wart hog are orthologs of trappin-1 (SPAI) and trappin-2 (elafin), the others are novel members of the trappin family and named trappins-6 to 11. The divergence of the sequences is greatest in the region that encodes the reactive site, and intron sequences appear to be more highly conserved than the protein-coding sequences, especially among the pig paralogs. Phylogenetic analysis showed that the trappin multigene family members of pig were generated through gene duplication after the divergence of the Suidae (pig and wart hog) and Tayassuidae (collared peccary). Similarities in the gene structure with seminal vesicle clotting proteins (REST) and WAP motif-containing proteins suggest that trappins are naturally occurring fusion proteins created through exon shuffling between ancestral REST and WAP motif-coding genes.

The roles of semantic networks and search efficiency in verbal fluency performance in intractable temporal lobe epilepsy.

Two competing hypotheses (i.e., disruption of semantic networks vs. search inefficiency) concerning the mechanisms underlying impaired semantic verbal fluency in temporal lobe epilepsy (TLE) were tested within a single paradigm. Reports that semantic verbal fluency is more impaired in left than right TLE groups were confirmed by the findings that the left TLE group produced fewer words on a supermarket fluency task than did the normal control (NC) group, and that the performance of the right TLE group was intermediate to that of the left TLE and NC groups. Because both TLE groups generated fewer words per category of supermarket items sampled, and produced a higher ratio of category labels relative to category exemplars than did the NC group, it can be surmised that TLE disrupts semantic memory networks. The findings did not support the competing hypothesis that reduced semantic verbal fluency in TLE is a manifestation of inefficient search/retrieval strategies, possibly associated with distal frontal lobe pathophysiology. Specifically, the TLE and NC groups did not differ significantly in their mean number of perseverations, intrusions, or search efficiency (operationalized as the ratio of the number of shifts between categories to the number of categories sampled).

Pregnancy establishment and maintenance in cattle.

A single ovulation, reciprocal embryo transfer study was used to investigate effects of oocyte competence and maternal environment on pregnancy establishment and maintenance in beef cows. Estrous cycles were synchronized in suckled beef cows and embryo donors were inseminated on d 0 (n = 810). Cows were classified on d 0 as having a small (<12.5 mm) or large (≥12.5 mm) ovulatory follicle and randomly chosen as donors or recipients to remove confounding effects of ovulatory follicle size on fertility. Embryos (n = 393) or oocytes (n = 44) were recovered on d 7, and all viable embryos were transferred into recipients (n = 354). All statistical analyses were conducted using the GLM procedure of SAS. Path analysis (with significance set at P < 0.10) was used to examine potential cause-effect relationships among the measured variables. Greater donor cow BW, circulating estradiol concentration at insemination, postpartum interval, and ovulatory follicle size directly increased (P < 0.10) fertilization success. Greater donor cow age was the only factor that directly decreased (P < 0.10) fertilization success. Viability of d-7 embryos was directly inhibited (P < 0.10) by rapid follicular growth rate from d -2 to 0 and heavier BW. Direct beneficial effects to embryo viability were increased serum progesterone concentration on d -2 and ovulatory follicle size. Pregnancy maintenance from d 7 to 27 was enhanced (P < 0.10) by increased serum estradiol concentration on d 0 and progesterone concentration on d 7 in the recipient cow. Increased follicular diameter in the recipient cow on d 0 was detrimental to pregnancy maintenance from d 7 to 27. This manuscript defines the complex interplay and relative contributions of endocrine and physical factors both prior and subsequent to fertilization that influence both oocyte competence and maternal environment and their roles in establishment and maintenance of pregnancy.

Preovulatory estradiol and the establishment and maintenance of pregnancy in suckled beef cows.

In postpartum beef cows, GnRH-induced ovulation of small dominant follicles decreased pregnancy rates and increased late embryonic/fetal mortality. In Exp. 1, single ovulation reciprocal embryo transfer (ET) was used to examine the relationship between preovulatory serum concentrations of estradiol at GnRH-induced ovulation in donor and recipient cows and establishment and maintenance of pregnancy. Suckled beef cows (n = 1,164) were administered GnRH (GnRH1, 100 µg) on d -9 (GnRH1), PGF(2α) on d -2, and GnRH2 (GnRH2, 100 µg) on d 0 (CO-Synch protocol) either with (donors; n = 810) or without (recipients; n = 354) AI. Single embryos (n = 394) or oocytes (n = 45) were recovered from the donor cows (d 7; ET) and all live embryos were transferred into recipients. Serum concentration of estradiol at GnRH2 was positively correlated with follicle size at GnRH2 (r = 0.45, P < 0.01) and progesterone at ET (r =0.34, P < 0.01). Donor cows with greater estradiol at GnRH2 were more likely to yield an embryo than an unfertilized oocyte (P < 0.01). Donor and recipient cows were retrospectively divided into 4 groups [low estradiol (<8.4 pg/mL) or high estradiol (≥8.4 pg/mL)] based on serum concentration of estradiol at GnRH2. Pregnancy rate at d 27 for low-low (n = 78), low-high (n = 80), high-low (n = 91), and high-high (n = 101) groups (donor-recipient, respectively) was 45, 65, 43, and 61% respectively (P < 0.02). Because recipient cows with greater estradiol concentration at GnRH2 had greater pregnancy rates in Exp. 1, the objective of Exp. 2 was to evaluate the effect of estradiol supplementation on pregnancy rate. Ovulation was synchronized in suckled beef cows (n = 600) using the CO-Synch protocol with the insertion of a controlled internal drug release (CIDR; intravaginal progesterone supplement) from d -9 until d -2. Approximately one-half of the cows (n = 297) received an injection of estradiol cypionate (ECP; 0.5 mg intramuscularly) 24 h before AI. Compared with the no treatment (Control) cows, ECP treatment increased (P < 0.01) pregnancy rates of cows induced to ovulate smaller dominant follicles (<12.2 mm). In conclusion, GnRH-induced ovulation of small dominant follicles was associated with reduced serum estradiol, fertilization rate (donor cows), and pregnancy establishment (recipient cows). Furthermore, ECP supplementation during the preovulatory period increased pregnancy rates in cows induced to ovulate smaller dominant follicles.

Using 50 k single nucleotide polymorphisms to elucidate genomic architecture of line 1 hereford cattle.

Hereford is a major beef breed in the USA, and a sub-population, known as Line 1 (L1), was established in 1934 using two paternal half-sib bulls and 50 unrelated females. L1 has since been maintained as a closed population and selected for growth to 1 year of age. Objectives were to characterize the molecular genetic architecture of L1 (n = 240) by comparing a cross-section of L1 with the general US. Hereford population (AHA, n = 311), estimating effects of imposed selection within L1 based on allele frequencies at 50 K SNP loci, and examining loci-specific effects of heterozygosity on the selection criterion. Animals were genotyped using the Illumina BovineSNP50 Beadchip, and SNP were mapped to UMD3.0 assembly of the bovine genome sequence. Average linkage disequilibrium (LD), measured by square of Pearson correlation, of adjacent SNP was 0.36 and 0.16 in L1 and AHA, respectively. Difference in LD between L1 and AHA decreased as SNP spacing increased. Persistence of phase between L1 and AHA decreased from 0.45 to 0.14 as SNP spacing increased from 50 to 5,000 kb. Extended haplotype homozygosity was greater in L1 than in AHA for 95.6% of the SNP. Knowledge of selection applied to L1 facilitated a novel approach to QTL discovery. Minor allele frequency was (FDR < 0.01) affected by cumulative selection differential at 191 out of 25,901 SNP. With the FDR relaxed to 0.05, 13 regions on BTA2, 5, 6, 9, 11, 14, 15, 18, 23, and 26 are co-located with previously identified QTL for growth. After adjustment of postweaning gain phenotypes for fixed effects and direct additive genetic effects, regression of residuals on genome-wide heterozygosity was -235.3 ± 91.6 kg. However, no SNP-specific loci where heterozygotes were significantly superior to the average of homozygotes were revealed (FDR ≥ 0.17). In conclusion, genome-wide SNP genotypes clarified effects of selection and inbreeding within L1 and differences in genomic architecture between the population segment L1 and the AHA population.

Effects of reduced panel, reference origin, and genetic relationship on imputation of genotypes in Hereford cattle.

The objective of this study was to investigate alternative methods of designing and using reduced SNP panels for imputing SNP genotypes. Two purebred Hereford populations, an experimental population known as Line 1 Hereford (L1, n = 240) and registered Hereford with American Hereford Association (AHA, n = 311), were used. Using different reference samples of 62 to 311 animals with 39,497 SNP on 29 autosomes and study samples of 57 or 62 animals for which genotypes were available for ~2,600 SNP (reduced panels), imputations were performed to predict the other ~36,900 loci that had been masked. An imputation package, including LinkPHASE and DAGPHASE, was used for imputation. Four reduced panels differing in minor allele frequency (MAF) and marker spacing were evaluated. Reduced panels included every 15th SNP across the genome (SNP_space), commercial Illumina Bovine3K Beadchip (SNP_3K), SNP with the highest MAF (SNP_MAF), and SNP with high MAF that were also evenly spaced across the genome (SNP_MS). Imputation accuracy was defined as the correlation of imputed genotypes and real genotypes. Reference samples were either from L1 or AHA. Among animals with genotypes, genetic relationships were estimated based on molecular marker genotypes or pedigree. Reduced panel design, number of animals in the reference sample, reference origin and genetic relationship between animals in the reference, and study samples all affected imputation accuracy (P < 0.001). Across genotyping schemes, imputed genotypes from SNP_MS had the greatest accuracy. A 0.1 increase in average pedigree relationship or average molecular relationship between reference and study samples increased imputation accuracy 10 to 20%. Using reference samples from the L1 population resulted in lower imputation accuracy than using reference samples from the admixed population AHA (P < 0.001). Increasing the number of animals in the reference panel by 100 individuals increased imputation accuracy by 8% when pedigree relationship was used as a covariate and 6% when molecular relationship was used as a covariate. We concluded that imputation accuracy would be increased through optimization of reduced panel design and genotyping strategy.