Extensively grazed cattle are often mustered only once a year. Therefore, birthdates are typically unknown or inaccurate. Birthdates would be useful for deriving important traits (growth rate; calving interval), breed registrations, and making management decisions. Epigenetic clocks use methylation of DNA to predict an individual’s age. An epigenetic clock for cattle could provide a solution to the challenges of industry birthdate recording. Here we derived the first epigenetic clock for tropically adapted cattle using portable sequencing devices from tail hair, a tissue which is widely used in industry for genotyping. Cattle (n = 66) with ages ranging from 0.35 to 15.7 years were sequenced using Oxford Nanopore Technologies MinION and methylation was called at CpG sites across the genome. Sites were then filtered and used to calculate a covariance relationship matrix based on methylation state. Best linear unbiased prediction was used with 10-fold cross validation to predict age. A second methylation relationship matrix was also calculated that contained sites associated with genes used in the dog and human epigenetic clocks. The correlation between predicted age and actual age was 0.71 for all sites and 0.60 for dog and human gene epigenetic clock sites. The mean absolute deviation was 1.4 years for animals aged less than 3 years of age, and 1.5 years for animals aged 3–10 years. This is the first reported epigenetic clock using industry relevant samples in cattle.

Background Genomic prediction describes the use of SNP genotypes to predict complex phenotypes and has been widely applied in humans and agriculture species. Genotyping-by-sequencing, a method which uses low-coverage sequence data paired with genotype imputation, is becoming increasingly popular for SNP genotyping. The development of Oxford Nanopore Technologies’ (ONT) MinION sequencer has now made genotyping-by-sequencing portable and rapid. Here we evaluate the speed and accuracy of genomic predictions using low-coverage ONT sequence data in a population of cattle using four imputation approaches. We also investigate the effect of SNP reference panel size on their performance. Results SNP array genotypes and ONT sequence data for 64 beef heifers were used to calculate genomic estimated breeding values (GEBVs) from 641k SNP for four traits. Accuracy of the GEBVs was much higher when flanking SNP from sequence data was used to help impute the 641k panel used for genomic predictions. Using the imputation package QUILT, correlations between ONT and low-density SNP array genomic breeding values were greater than 0.91 and up to 0.97 for sequencing coverages as low as 0.1x using a panel of 48 million SNP that flanked the 641k in the prediction equation. Imputation time was significantly reduced by decreasing the number of flanking sequence SNP used in imputation for all methods. Genomic breeding values calculated using QUILT also had higher correlations to high density SNP arrays than genomic breeding values from imputed-low density arrays for coverages as low as 0.5x. Conclusions Here we demonstrated accurate genomic prediction is possible with ONT sequence data from sequencing coverages as low as 0.1x, and imputation time can be as short as 10 minutes per sample.

The northern Australia beef cattle industry operates in harsh environmental conditions which consistently suppress female fertility. To better understand the environmental effect on cattle raised extensively in northern Australia, new environmental descriptors were defined for 54 commercial herds located across the region. Three fertility traits, based on the presence of a corpus luteum at 600 d of age, indicating puberty, (CL Presence, n = 25,176), heifer pregnancy (n = 20,989) and first lactation pregnancy (n = 10,072) were recorded. Temperature, humidity, and rainfall were obtained from publicly available data based on herd location. Being pubertal at 600 d (i.e. CL Presence) increased the likelihood of success at heifer pregnancy and first lactation pregnancy (P < 0.05), underscoring the importance of early puberty in reproductive success. A temperature humidity index (THI) of 65-70 had a significant (P < 0.05) negative effect on first lactation pregnancy rate, heifer pregnancy and puberty at 600 d of age. Area under the curve of daily THI was significant (P < 0.05) and reduced the likelihood of pregnancy at first lactation and puberty at 600 days. Deviation from long-term average rainfall was not significant (P < 0.05) for any trait. Average daily weight gain had a significant and positive relationship (P < 0.05) for heifer and first lactation pregnancy. The results indicate that chronic or cumulative heat load is more determinantal to reproductive performance than acute heat stress. The reason for the lack of a clear relationship between acute heat stress and reproductive performance is unclear but may be partially explained by peak THI and peak nutrition coinciding at the same time. Sufficient evidence was found to justify the use of average daily weight gain and chronic heat load as descriptors to define an environmental gradient.