Concurrently mapping quantitative trait loci associations from multiple subspecies within hybrid populations.

Many of the world's agriculturally important plant and animal populations consist of hybrids of subspecies. Cattle in tropical and sub-tropical regions for example, originate from two subspecies, Bos taurus indicus (Bos indicus) and Bos taurus taurus (Bos taurus). Methods to derive the underlying genetic architecture for these two subspecies are essential to develop accurate genomic predictions in these hybrid populations. We propose a novel method to achieve this. First, we use haplotypes to assign SNP alleles to ancestral subspecies of origin in a multi-breed and multi-subspecies population. Then we use a BayesR framework to allow SNP alleles originating from the different subspecies differing effects. Applying this method in a composite population of B. indicus and B. taurus hybrids, our results show that there are underlying genomic differences between the two subspecies, and these effects are not identified in multi-breed genomic evaluations that do not account for subspecies of origin effects. The method slightly improved the accuracy of genomic prediction. More significantly, by allocating SNP alleles to ancestral subspecies of origin, we were able to identify four SNP with high posterior probabilities of inclusion that have not been previously associated with cattle fertility and were close to genes associated with fertility in other species. These results show that haplotypes can be used to trace subspecies of origin through the genome of this hybrid population and, in conjunction with our novel Bayesian analysis, subspecies SNP allele allocation can be used to increase the accuracy of QTL association mapping in genetically diverse populations.

Using prior information from humans to prioritize genes and gene-associated variants for complex traits in livestock.

Genome-Wide Association Studies (GWAS) in large human cohorts have identified thousands of loci associated with complex traits and diseases. For identifying the genes and gene-associated variants that underlie complex traits in livestock, especially where sample sizes are limiting, it may help to integrate the results of GWAS for equivalent traits in humans as prior information. In this study, we sought to investigate the usefulness of results from a GWAS on human height as prior information for identifying the genes and gene-associated variants that affect stature in cattle, using GWAS summary data on samples sizes of 700,000 and 58,265 for humans and cattle, respectively. Using Fisher's exact test, we observed a significant proportion of cattle stature-associated genes (30/77) that are also associated with human height (odds ratio = 5.1, p = 3.1e-10). Result of randomized sampling tests showed that cattle orthologs of human height-associated genes, hereafter referred to as candidate genes (C-genes), were more enriched for cattle stature GWAS signals than random samples of genes in the cattle genome (p = 0.01). Randomly sampled SNPs within the C-genes also tend to explain more genetic variance for cattle stature (up to 13.2%) than randomly sampled SNPs within random cattle genes (p = 0.09). The most significant SNPs from a cattle GWAS for stature within the C-genes did not explain more genetic variance for cattle stature than the most significant SNPs within random cattle genes (p = 0.87). Altogether, our findings support previous studies that suggest a similarity in the genetic regulation of height across mammalian species. However, with the availability of a powerful GWAS for stature that combined data from 8 cattle breeds, prior information from human-height GWAS does not seem to provide any additional benefit with respect to the identification of genes and gene-associated variants that affect stature in cattle.

Imputation accuracy to whole-genome sequence in Nellore cattle.

BACKGROUND: A cost-effective strategy to explore the complete DNA sequence in animals for genetic evaluation purposes is to sequence key ancestors of a population, followed by imputation mechanisms to infer marker genotypes that were not originally reported in a target population of animals genotyped with single nucleotide polymorphism (SNP) panels. The feasibility of this process relies on the accuracy of the genotype imputation in that population, particularly for potential causal mutations which may be at low frequency and either within genes or regulatory regions. The objective of the present study was to investigate the imputation accuracy to the sequence level in a Nellore beef cattle population, including that for variants in annotation classes which are more likely to be functional. METHODS: Information of 151 key sequenced Nellore sires were used to assess the imputation accuracy from bovine HD BeadChip SNP (~ 777 k) to whole-genome sequence. The choice of the sires aimed at optimizing the imputation accuracy of a genotypic database, comprised of about 10,000 genotyped Nellore animals. Genotype imputation was performed using two computational approaches: FImpute3 and Minimac4 (after using Eagle for phasing). The accuracy of the imputation was evaluated using a fivefold cross-validation scheme and measured by the squared correlation between observed and imputed genotypes, calculated by individual and by SNP. SNPs were classified into a range of annotations, and the accuracy of imputation within each annotation classification was also evaluated. RESULTS: High average imputation accuracies per animal were achieved using both FImpute3 (0.94) and Minimac4 (0.95). On average, common variants (minor allele frequency (MAF) > 0.03) were more accurately imputed by Minimac4 and low-frequency variants (MAF ≤ 0.03) were more accurately imputed by FImpute3. The inherent Minimac4 Rsq imputation quality statistic appears to be a good indicator of the empirical Minimac4 imputation accuracy. Both software provided high average SNP-wise imputation accuracy for all classes of biological annotations. CONCLUSIONS: Our results indicate that imputation to whole-genome sequence is feasible in Nellore beef cattle since high imputation accuracies per individual are expected. SNP-wise imputation accuracy is software-dependent, especially for rare variants. The accuracy of imputation appears to be relatively independent of annotation classification.

Pregnancy rate per cycle is heritable and reduces with cycle in naturally mated tropically adapted beef cows.

Pregnancy in cattle is the outcome of the complex process of initiation of cycling, fertilization, maternal recognition of pregnancy and foeto-placental development. Though much is known about initiation of cycling and associated risk factors, there are virtually no data on pregnancy rate per cycle for naturally mated cattle, especially for extensively managed, tropically adapted genotypes, which this study aimed to determine. Tropical composite (Bos indicus and African Sanga crosses with Bos taurus) and Brahman cattle (n = 2,181) of known pedigree in four-year groups at four sites were mated annually for 84 days. Body condition, ovarian function, pregnancies, calving and lactation were monitored through six full reproductive cycles using 4-8 weekly ultrasound of the reproductive tract outside the calving period and daily monitoring during calving. From this, dates of commencement of cycling and conception in each year were estimated for each animal, enabling calculation of established pregnancy for consecutive 21-day periods while cycling and of pregnancies within four months of calving while lactating (P4M). Pregnancy per 21-day period (cycle) during mating for cycling animals averaged 63%, 71%, 41% and 28% in four consecutive cycles. Pregnant per cycle was 2%-11% higher in tropical composites than in Brahmans. The only other consistently significant risk to becoming pregnant was if cycling commenced later than three weeks before mating commenced. P4M averaged 62% and was lower for cows in sub-optimal body condition and in first-parity and later-calving cows. Pregnant per cycle was moderately heritable (~20%), while heritability was moderate to high (33%) for P4M. Selection for pregnant per cycle could be achieved indirectly by selection for P4M, a trait that is readily measured.

Genetic control of temperament traits across species: association of autism spectrum disorder risk genes with cattle temperament.

BACKGROUND: Temperament traits are of high importance across species. In humans, temperament or personality traits correlate with psychological traits and psychiatric disorders. In cattle, they impact animal welfare, product quality and human safety, and are therefore of direct commercial importance. We hypothesized that genetic factors that contribute to variation in temperament among individuals within a species will be shared between humans and cattle. Using imputed whole-genome sequence data from 9223 beef cattle from three cohorts, a series of genome-wide association studies was undertaken on cattle flight time, a temperament phenotype measured as the time taken for an animal to cover a short-fixed distance after release from an enclosure. We also investigated the association of cattle temperament with polymorphisms in bovine orthologs of risk genes for neuroticism, schizophrenia, autism spectrum disorders (ASD), and developmental delay disorders in humans. RESULTS: Variants with the strongest associations were located in the bovine orthologous region that is involved in several behavioural and cognitive disorders in humans. These variants were also partially validated in independent cattle cohorts. Genes in these regions (BARHL2, NDN, SNRPN, MAGEL2, ABCA12, KIFAP3, TOPAZ1, FZD3, UBE3A, and GABRA5) were enriched for the GO term neuron migration and were differentially expressed in brain and pituitary tissues in humans. Moreover, variants within 100 kb of ASD susceptibility genes were associated with cattle temperament and explained 6.5% of the total additive genetic variance in the largest cattle cohort. The ASD genes with the most significant associations were GABRB3 and CUL3. Using the same 100 kb window, a weak association was found with polymorphisms in schizophrenia risk genes and no association with polymorphisms in neuroticism and developmental delay disorders risk genes. CONCLUSIONS: Our analysis showed that genes identified in a meta-analysis of cattle temperament contribute to neuron development functions and are differentially expressed in human brain tissues. Furthermore, some ASD susceptibility genes are associated with cattle temperament. These findings provide evidence that genetic control of temperament might be shared between humans and cattle and highlight the potential for future analyses to leverage results between species.

Use of whole-genome sequence data and novel genomic selection strategies to improve selection for age at puberty in tropically-adapted beef heifers.

BACKGROUND: In tropically-adapted beef heifers, application of genomic prediction for age at puberty has been limited due to low prediction accuracies. Our aim was to investigate novel methods of pre-selecting whole-genome sequence (WGS) variants and alternative analysis methodologies; including genomic best linear unbiased prediction (GBLUP) with multiple genomic relationship matrices (MGRM) and Bayesian (BayesR) analyses, to determine if prediction accuracy for age at puberty can be improved. METHODS: Genotypes and phenotypes were obtained from two research herds. In total, 868 Brahman and 960 Tropical Composite heifers were recorded in the first population and 3695 Brahman, Santa Gertrudis and Droughtmaster heifers were recorded in the second population. Genotypes were imputed to 23 million whole-genome sequence variants. Eight strategies were used to pre-select variants from genome-wide association study (GWAS) results using conditional or joint (COJO) analyses. Pre-selected variants were included in three models, GBLUP with a single genomic relationship matrix (SGRM), GBLUP MGRM and BayesR. Five-way cross-validation was used to test the effect of marker panel density (6 K, 50 K and 800 K), analysis model, and inclusion of pre-selected WGS variants on prediction accuracy. RESULTS: In all tested scenarios, prediction accuracies for age at puberty were highest in BayesR analyses. The addition of pre-selected WGS variants had little effect on the accuracy of prediction when BayesR was used. The inclusion of WGS variants that were pre-selected using a meta-analysis with COJO analyses by chromosome, fitted in a MGRM model, had the highest prediction accuracies in the GBLUP analyses, regardless of marker density. When the low-density (6 K) panel was used, the prediction accuracy of GBLUP was equal (0.42) to that with the high-density panel when only six additional sequence variants (identified using meta-analysis COJO by chromosome) were included. CONCLUSIONS: While BayesR consistently outperforms other methods in terms of prediction accuracies, reasonable improvements in accuracy can be achieved when using GBLUP and low-density panels with the inclusion of a relatively small number of highly relevant WGS variants.

Unraveling genetic sensitivity of beef cattle to environmental variation under tropical conditions.

BACKGROUND: Selection of cattle that are less sensitive to environmental variation in unfavorable environments and more adapted to harsh conditions is of primary importance for tropical beef cattle production systems. Understanding the genetic background of sensitivity to environmental variation is necessary for developing strategies and tools to increase efficiency and sustainability of beef production. We evaluated the degree of sensitivity of beef cattle performance to environmental variation, at the animal and molecular marker levels (412 K single nucleotide polymorphisms), by fitting and comparing the results of different reaction norm models (RNM), using a comprehensive dataset of Nellore cattle raised under diverse environmental conditions. RESULTS: Heteroscedastic RNM (with different residual variances for environmental level) provided better fit than homoscedastic RNM. In addition, spline and quadratic RNM outperformed linear RNM, which suggests the existence of a nonlinear genetic component affecting the performance of Nellore cattle. This nonlinearity indicates that within-animal sensitivity depends on the environmental gradient (EG) level and that animals may present different patterns of sensitivity according to the range of environmental variations. The spline RNM showed that sensitivity to environmental variation from harsh to average EG is lowly correlated with sensitivity from average to good EG, at both the animal and molecular marker levels. Although the genomic regions that affect sensitivity in harsher environments were not the same as those associated with less challenging environments, the candidate genes within those regions participate in common biological processes such as those related to inflammatory and immune response. Some plausible candidate genes were identified. CONCLUSIONS: Sensitivity of tropical beef cattle to environmental variation is not continuous along the environmental gradient, which implies that animals that are less sensitive to harsher conditions are not necessarily less responsive to variations in better environmental conditions, and vice versa. The same pattern was observed at the molecular marker level, i.e. genomic regions and, consequently, candidate genes associated with sensitivity to harsh conditions were not the same as those associated with sensitivity to less challenging conditions.

Bias in relative survival methods when using incorrect life-tables: lung and bladder cancer by smoking status and ethnicity in New Zealand.

Relative survival and excess mortality approaches are commonly used to estimate and compare net survival from cancer. These approaches are based on the assumption that the underlying (non-cancer) mortality rate of cancer patients is the same as that of the general population. This assumption is likely to be violated particularly in the context of smoking-related cancers. The magnitude of this bias has not been estimated. The objective of this article is to estimate the bias in relative survival ratios (RSRs) and excess mortality rate ratios (EMRRs) from using total population compared to correct subpopulation specific life-tables. Analyses were conducted on 1996-2001 linked census-cancer data (including smoking status) for people with lung and bladder cancer, using sex-specific (standard practice), sex- and ethnic-specific, sex- and smoking-specific and sex-, ethnic- and smoking-specific life-tables. Five-year RSRs using sex-specific life-tables, compared to fully stratified life-tables, were underestimated by 10-25% for current smoking and Maori populations. For example, the current smoker male bladder cancer RSR was 0.700 for sex-specific life-tables, compared to 0.838 for fully stratified life-tables. Similarly, EMRRs comparing current to never smokers and Maori to non-Maori were overestimated using sex-specific life-tables only: modestly only for lung cancer, but markedly for bladder cancer. For example, the EMRR comparing current to never smokers with bladder cancer in a fully adjusted regression model was 1.475 when using sex-specific life-tables only, but reduced to 1.098 when using fully stratified life-tables. Substantial bias can occur when estimating relative cancer survival across subpopulations if non-matching life-tables are used.

Genome-wide identification of copy number variation and association with fat deposition in thin and fat-tailed sheep breeds.

Copy number variants (CNVs) are a type of genetic polymorphism which contribute to phenotypic variation in several species, including livestock. In this study, we used genomic data of 192 animals from 3 Iranian sheep breeds including 96 Baluchi sheep and 47 Lori-Bakhtiari sheep as fat-tailed breeds and 47 Zel sheep as thin-tailed sheep breed genotyped with Illumina OvineSNP50K Beadchip arrays. Also, for association test, 70 samples of Valle del Belice sheep were added to the association test as thin-tailed sheep breed. PennCNV and CNVRuler software were, respectively, used to study the copy number variation and genomic association analyses. We detected 573 and 242 CNVs in the fat and thin tailed breeds, respectively. In terms of CNV regions (CNVRs), these represented 328 and 187 CNVRs that were within or overlapping with 790 known Ovine genes. The CNVRs covered approximately 73.85 Mb of the sheep genome with average length 146.88 kb, and corresponded to 2.6% of the autosomal genome sequence. Five CNVRs were randomly chosen for validation, of which 4 were experimentally confirmed using Real time qPCR. Functional enrichment analysis showed that genes harbouring CNVs in thin-tailed sheep were involved in the adaptive immune response, regulation of reactive oxygen species biosynthetic process and response to starvation. In fat-tailed breeds these genes were involved in cellular protein modification process, regulation of heart rate, intestinal absorption, olfactory receptor activity and ATP binding. Association test identified one copy gained CNVR on chromosomes 6 harbouring two protein-coding genes HGFAC and LRPAP1. Our findings provide information about genomic structural changes and their association to the interested traits including fat deposition and environmental compatibility in sheep.

Assigning scores for ordered categorical responses.

Deciding on the best statistical method to apply when the response variable is ordinal is essential because the way the categories are ordered in the data is relevant as it could change the results of the analysis. Although the models for continuous variables have similarities to those for ordinal variables, this paper presents the advantages of the use of the ordering information on the outcomes with methods developed for modeling ordinal data such as the ordered stereotype model. The novelty of this article lies in showing the dangers of assigning equally spaced scores to ordered response categories in statistical analysis, which are illustrated with a simulation study and a case study. We propose a new way to use the score parameters, which incorporates the fitted spacing dictated by the data. Additionally, this article uses score parameter estimates in the ordered stereotype model to propose a new measure to calculate continuous medians in the raw data: the adjusted c-median. It benefits the general audience who can easily understand the median as a summary statistic. Supplementary materials for this article are available online.

A method for ordinal outcomes: The ordered stereotype model.

OBJECTIVE: The collection and use of ordinal variables are common in many psychological and psychiatric studies. Although the models for continuous variables have similarities to those for ordinal variables, there are advantages when a model developed for modeling ordinal data is used such as avoiding "floor" and "ceiling" effects and avoiding to assign scores, as it happens in continuous models, which can produce results sensitive to the score assigned. This paper introduces and focuses on the application of the ordered stereotype model, which was developed for modeling ordinal outcomes and is not so popular as other models such as linear regression and proportional odds models. This paper aims to compare the performance of the ordered stereotype model with other more commonly used models among researchers and practitioners. METHODS: This article compares the performance of the stereotype model against the proportional odd and linear regression models, with three, four, and five levels of ordinal categories and sample sizes 100, 500, and 1000. This paper also discusses the problem of treating ordinal responses as continuous using a simulation study. The trend odds model is also presented in the application. RESULTS: Three types of models were fitted in one real-life example, including ordered stereotype, proportional odds, and trend odds models. They reached similar conclusions in terms of the significance of covariates. The simulation study evaluated the performance of the ordered stereotype model under four cases. The performance varies depending on the scenarios. CONCLUSIONS: The method presented can be applied to several areas of psychiatry dealing with ordinal outcomes. One of the main advantages of this model is that it breaks with the assumption of levels of the ordinal response are equally spaced, which might be not true.

The burden of cancer in New Zealand: a comparison of incidence and DALY metrics and its relevance for ethnic disparities.

AIM: Cancer burden measured in disability adjusted life years (DALYs) captures survival and disability impacts of incident cancers. In this paper, we estimate the prospective burden of disease arising from 27 cancer sites diagnosed in 2006, by sex and ethnicity; and determine how its distribution differs from that for incidence rates alone. METHODS: Using a prospective approach, Markov and cancer disease models were used to estimate DALYs with inputs of population counts, incidence and excess mortality rates, disability weights, and background mortality. DALYs were discounted at 3.5% per year. RESULTS: The age standardised Maori:non-Maori incidence rate ratios were 1.00 for males and 1.19 for females, whereas for DALYs they were greater at 1.42 for males and 1.68 for females. The total burden of cancer for 2006 incident cases (i.e. not age standardised) was estimated to be approximately 127,000 DALYs. Breast (27%), lung (14%) and colorectal (13%) cancers for females and lung (16%), colorectal (14%), and prostate (16%) cancers for males were the top contributors. By ethnicity, Maori experienced a substantially higher burden from lung cancer (around 25% for both sexes). CONCLUSIONS: Due to Maori both having higher rates of cancers with a worse survival (e.g. lung cancer), and tending to have worse survival for each cancer site, ethnic disparities in the age-standardised DALY burden were greater than those for incidence (rate ratios of 1.52 and 1.07 respectively, sexes pooled).