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. 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. 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.

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 , A BCA12 , 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.

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.

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.

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.

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.