Facial eczema (FE) is a significant metabolic disease that affects New Zealand ruminants. Ingestion of the mycotoxin sporidesmin leads to liver and bile duct damage, which can result in photosensitisation, reduced productivity and death. Strategies used to manage the incidence and severity of the disease include breeding. In sheep, there is considerable genetic variation in the response to FE. A commercial testing program is available for ram breeders who aim to increase tolerance, determined by the concentration of the serum enzyme, gamma-glutamyltransferase 21 days after a measured sporidesmin challenge (GGT21). Genome-wide association studies were carried out to determine regions of the genome associated with GGT21. Two regions on chromosomes 15 and 24 are reported, which explain 5% and 1% of the phenotypic variance in the response to FE, respectively. The region on chromosome 15 contains the β-globin locus. Of the significant SNPs in the region, one is a missense variant within the haemoglobin subunit β (HBB) gene. Mass spectrometry of haemoglobin from animals with differing genotypes at this locus indicated that genotypes are associated with different forms of adult β-globin. Haemoglobin haplotypes have previously been associated with variation in several health-related traits in sheep and warrant further investigation regarding their role in tolerance to FE in sheep. We show a strategic approach to the identification of regions of importance for commercial breeding programs with a combination of discovery, statistical and biological validation. This study highlights the power of using increased density genotyping for the identification of influential genomic regions, combined with subsequent inclusion on lower density genotyping platforms.

Genotyping-by-sequencing (GBS) is becoming an attractive alternative to array-based methods for genotyping individuals for a large number of single nucleotide polymorphisms (SNPs). Costs can be lowered by reducing the mean sequencing depth, but this results in genotype calls of lower quality. A common analysis strategy is to filter SNPs to just those with sufficient depth, thereby greatly reducing the number of SNPs available. We investigate methods for estimating relatedness using GBS data, including results of low depth, using theoretical calculation, simulation and application to a real data set.

To assist cattle producers transition from microsatellite (MS) to single nucleotide polymorphism (SNP) genotyping for parental verification we previously devised an effective and inexpensive method to impute MS alleles from SNP haplotypes. While the reported method was verified with only a limited data set (N = 479) from Brown Swiss, Guernsey, Holstein, and Jersey cattle, some of the MS-SNP haplotype associations were concordant across these phylogenetically diverse breeds. This implied that some haplotypes predate modern breed formation and remain in strong linkage disequilibrium. To expand the utility of MS allele imputation across breeds, MS and SNP data from more than 8000 animals representing 39 breeds (Bos taurus and B. indicus) were used to predict 9410 SNP haplotypes, incorporating an average of 73 SNPs per haplotype, for which alleles from 12 MS markers could be accurately be imputed. Approximately 25% of the MS-SNP haplotypes were present in multiple breeds (N = 2 to 36 breeds). These shared haplotypes allowed for MS imputation in breeds that were not represented in the reference population with only a small increase in Mendelian inheritance inconsistancies. Our reported reference haplotypes can be used for any cattle breed and the reported methods can be applied to any species to aid the transition from MS to SNP genetic markers. While ~91% of the animals with imputed alleles for 12 MS markers had ≤1 Mendelian inheritance conflicts with their parents' reported MS genotypes, this figure was 96% for our reference animals, indicating potential errors in the reported MS genotypes. The workflow we suggest autocorrects for genotyping errors and rare haplotypes, by MS genotyping animals whose imputed MS alleles fail parentage verification, and then incorporating those animals into the reference dataset.

Identification of genomic regions that have been targets of selection for phenotypic traits is one of the most important and challenging areas of research in animal genetics. However, currently there are relatively few genomic regions identified that have been subject to positive selection. In this study, a genome-wide scan using ~50,000 Single Nucleotide Polymorphisms (SNPs) was performed in an attempt to identify genomic regions associated with fat deposition in fat-tail breeds. This trait and its modification are very important in those countries grazing these breeds.

The potent greenhouse gas methane (CH4) is produced in the rumens of ruminant animals from hydrogen produced during microbial degradation of ingested feed. The natural animal-to-animal variation in the amount of CH4 emitted and the heritability of this trait offer a means for reducing CH4 emissions by selecting low-CH4 emitting animals for breeding. We demonstrate that differences in rumen microbial community structure are linked to high and low CH4 emissions in sheep. Bacterial community structures in 236 rumen samples from 118 high- and low-CH4 emitting sheep formed gradual transitions between three ruminotypes. Two of these (Q and S) were linked to significantly lower CH4 yields (14.4 and 13.6 g CH4/kg dry matter intake [DMI], respectively) than the third type (H; 15.9 g CH4/kg DMI; p<0.001). Low-CH4 ruminotype Q was associated with a significantly lower ruminal acetate to propionate ratio (3.7±0.4) than S (4.4±0.7; p<0.001) and H (4.3±0.5; p<0.001), and harbored high relative abundances of the propionate-producing Quinella ovalis. Low-CH4 ruminotype S was characterized by lactate- and succinate-producing Fibrobacter spp., Kandleria vitulina, Olsenella spp., Prevotella bryantii, and Sharpea azabuensis. High-CH4 ruminotype H had higher relative abundances of species belonging to Ruminococcus, other Ruminococcaceae, Lachnospiraceae, Catabacteriaceae, Coprococcus, other Clostridiales, Prevotella, other Bacteroidales, and Alphaproteobacteria, many of which are known to form significant amounts of hydrogen. We hypothesize that lower CH4 yields are the result of bacterial communities that ferment ingested feed to relatively less hydrogen, which results in less CH4 being formed.

Cattle and other ruminants have evolved the ability to derive most of their metabolic energy requirement from otherwise indigestible plant matter through a symbiotic relationship with plant fibre degrading microbes within a specialised fermentation chamber, the rumen. The genetic changes underlying the evolution of the ruminant lifestyle are poorly understood. The BPI-like locus encodes several putative innate immune proteins, expressed predominantly in the oral cavity and airways, which are structurally related to Bactericidal/Permeability Increasing protein (BPI). We have previously reported the expression of variant BPI-like proteins in cattle (Biochim Biophys Acta 2002, 1579, 92-100). Characterisation of the BPI-like locus in cattle would lead to a better understanding of the role of the BPI-like proteins in cattle physiology

Ruminants obtain nutrients from microbial fermentation of plant material, primarily in their rumen, a multilayered forestomach. How the different layers of the rumen wall respond to diet and influence microbial fermentation, and how these process are regulated, is not well understood. Gene expression correlation networks were constructed from full thickness rumen wall transcriptomes of 24 sheep fed two different amounts and qualities of a forage and measured for methane production. The network contained two major negatively correlated gene sub-networks predominantly representing the epithelial and muscle layers of the rumen wall. Within the epithelium sub-network gene clusters representing lipid/oxo-acid metabolism, general metabolism and proliferating and differentiating cells were identified. The expression of cell cycle and metabolic genes was positively correlated with dry matter intake, ruminal short chain fatty acid concentrations and methane production. A weak correlation between lipid/oxo-acid metabolism genes and methane yield was observed. Feed consumption level explained the majority of gene expression variation, particularly for the cell cycle genes. Many known stratified epithelium transcription factors had significantly enriched targets in the epithelial gene clusters. The expression patterns of the transcription factors and their targets in proliferating and differentiating skin is mirrored in the rumen, suggesting conservation of regulatory systems.

Black-billed gulls (Larus bulleri) are endemic to New Zealand and are suspected to be undergoing substantial population declines. They primarily breed on open gravel beds in braided rivers of the South Island-a habitat that is diminishing and becoming increasingly modified. Although management of this species is increasing, little has been published on their movements and demographics. In this study, both mitochondrial DNA (mtDNA) control region domain I and nuclear single nucleotide polymorphisms (SNPs) were examined to help understand the connectivity and population structure of black-billed gulls across the country and to help inform management decisions. Mitochondrial DNA showed no population structure, with high haplotype and low nucleotide diversity, and analyses highlighted mitochondrial introgression with the closely related red-billed gulls (Larus novaehollandiae scopulinus). Nuclear DNA analyses, however, identified two groups, with Rotorua birds in the North Island being distinct from the rest of New Zealand, and isolation-by-distance evident across the South Island populations. Gene flow primarily occurs between nearby colonies with a stepwise movement across the landscape. The importance from a genetic perspective of the more isolated North Island birds (1.6% of total population) needs to be further evaluated. From our results, we infer that the South Island black-billed gull management should focus on maintaining several populations within each region rather than focusing on single specific colonies or river catchments. Future study is needed to investigate the genetic structure of populations at the northern limit of the species' range, and identify the mechanisms behind, and extent of, the hybridisation between red-billed and black-billed gulls.

Microbial community profiles have been associated with a variety of traits, including methane emissions in livestock. These profiles can be difficult and expensive to obtain for thousands of samples (e.g. for accurate association of microbial profiles with traits), therefore the objective of this work was to develop a low-cost, high-throughput approach to capture the diversity of the rumen microbiome. Restriction enzyme reduced representation sequencing (RE-RRS) using ApeKI or PstI, and two bioinformatic pipelines (reference-based and reference-free) were compared to bacterial 16S rRNA gene sequencing using repeated samples collected two weeks apart from 118 sheep that were phenotypically extreme (60 high and 58 low) for methane emitted per kg dry matter intake (n = 236). DNA was extracted from freeze-dried rumen samples using a phenol chloroform and bead-beating protocol prior to RE-RRS. The resulting sequences were used to investigate the repeatability of the rumen microbial community profiles, the effect of laboratory and analytical method, and the relationship with methane production. The results suggested that the best method was PstI RE-RRS analyzed with the reference-free approach, which accounted for 53.3±5.9% of reads, and had repeatabilities of 0.49±0.07 and 0.50±0.07 for the first two principal components (PC1 and PC2), phenotypic correlations with methane yield of 0.43±0.06 and 0.46±0.06 for PC1 and PC2, and explained 41±8% of the variation in methane yield. These results were significantly better than for bacterial 16S rRNA gene sequencing of the same samples (p<0.05) except for the correlation between PC2 and methane yield. A Sensitivity study suggested approximately 2000 samples could be sequenced in a single lane on an Illumina HiSeq 2500, meaning the current work using 118 samples/lane and future proposed 384 samples/lane are well within that threshold. With minor adaptations, our approach could be used to obtain microbial profiles from other metagenomic samples.

Investments in genetic selection have played a major role in the New Zealand sheep industry competitiveness. Selection may erode genetic diversity, which is a crucial factor for the success of breeding programs. Better understanding of linkage disequilibrium (LD) and ancestral effective population size (Ne) through quantifying this diversity and comparison between populations allows for more informed decisions with regards to selective breeding taking population genetic diversity into account. The estimation of N e can be determined via genetic markers and requires knowledge of genetic distances between these markers. Single nucleotide polymorphisms (SNP) data from a sample of 12,597 New Zealand crossbred and purebred sheep genotyped with the Illumina Ovine SNP50 BeadChip was used to perform a genome-wide scan of LD and N e . Three methods to estimate genetic distances were investigated: 1) M1: a ratio fixed across the whole genome of one Megabase per centiMorgan; 2) M2: the ratios of genetic distance (using M3, below) over physical distance fixed for each chromosome; and, 3) M3: a genetic map of inter-SNP distances estimated using CRIMAP software (v2.503).

Ruminant livestock represent the single largest anthropogenic source of the potent greenhouse gas methane, which is generated by methanogenic archaea residing in ruminant digestive tracts. While differences between individual animals of the same breed in the amount of methane produced have been observed, the basis for this variation remains to be elucidated. To explore the mechanistic basis of this methane production, we measured methane yields from 22 sheep, which revealed that methane yields are a reproducible, quantitative trait. Deep metagenomic and metatranscriptomic sequencing demonstrated a similar abundance of methanogens and methanogenesis pathway genes in high and low methane emitters. However, transcription of methanogenesis pathway genes was substantially increased in sheep with high methane yields. These results identify a discrete set of rumen methanogens whose methanogenesis pathway transcription profiles correlate with methane yields and provide new targets for CH4 mitigation at the levels of microbiota composition and transcriptional regulation.

Genotype imputation is a key element of the implementation of genomic selection within the New Zealand sheep industry, but many factors can influence imputation accuracy. Our objective was to provide practical directions on the implementation of imputation strategies in a multi-breed sheep population genotyped with three single nucleotide polymorphism (SNP) panels: 5K, 50K and HD (600K SNPs).

Events of inbreeding are inevitable in critically endangered species. Reduced population sizes and unique life-history traits can increase the severity of inbreeding, leading to declines in fitness and increased risk of extinction. Here, we investigate levels of inbreeding in a critically endangered flightless parrot, the kākāpō (Strigops habroptilus), wherein a highly inbred island population and one individual from the mainland of New Zealand founded the entire extant population. Genotyping-by-sequencing (GBS), and a genotype calling approach using a chromosome-level genome assembly, identified a filtered set of 12,241 single-nucleotide polymorphisms (SNPs) among 161 kākāpō, which together encompass the total genetic potential of the extant population. Multiple molecular-based estimates of inbreeding were compared, including genome-wide estimates of heterozygosity (FH), the diagonal elements of a genomic-relatedness matrix (FGRM), and runs of homozygosity (RoH, FRoH). In addition, we compared levels of inbreeding in chicks from a recent breeding season to examine if inbreeding is associated with offspring survival. The density of SNPs generated with GBS was sufficient to identify chromosomes that were largely homozygous with RoH distributed in similar patterns to other inbred species. Measures of inbreeding were largely correlated and differed significantly between descendants of the two founding populations. However, neither inbreeding nor ancestry was found to be associated with reduced survivorship in chicks, owing to unexpected mortality in chicks exhibiting low levels of inbreeding. Our study highlights important considerations for estimating inbreeding in critically endangered species, such as the impacts of small population sizes and admixture between diverse lineages.

Knowledge about the genetic diversity of a population is a crucial parameter for the implementation of successful genomic selection and conservation of genetic resources. The aim of this research was to establish the scientific basis for the implementation of genomic selection in a composite Terminal sheep breeding scheme by providing consolidated linkage disequilibrium (LD) measures across SNP markers, estimating consistency of gametic phase between breed-groups, and assessing genetic diversity measures, such as effective population size (Ne), and population structure parameters, using a large number of animals (n = 14,845) genotyped with a high density SNP chip (606,006 markers). Information generated in this research will be useful for optimizing molecular breeding values predictions and managing the available genetic resources.

Meat is an important source of high-value protein providing sustainable nutrition for human health. The discolouration of meat results in significant waste, which threatens the sustainability of meat production in terms of availability, affordability, and utilisation. Advancing the knowledge of factors and underlying mechanisms for meat discolouration supports the sustainability transformation of meat production practices. Previous studies found that colour stability may be associated with signature changes in certain metabolites, including NADH, glutamate, methionine, and testosterone. This study aimed to confirm the effect of these metabolites and sex, sire, and muscle type on lamb meat colour. NADH and glutamate improved colour stability as evidenced by the increased metmyoglobin reductase activity, while methionine and testosterone had detrimental effects. Overall, lamb meat was discoloured with retail display for up to 10 days at 4 °C. The semitendinosus muscle had higher L*, b*, and hue angle and lower a* (p < 0.05) than other muscles, especially in ewes. Lamb meat from rams had a higher L* and hue angle and lower a* than the ewes (p < 0.05), especially in the colour-labile group, suggesting an interaction between sex and sire. The outcomes of this study will help make the production of meat more sustainable by assisting the meat industry in improving the selection of animals for meat production and processing practices to reduce meat waste due to discolouration.

Gastrointestinal nematodes are one of the most serious causes of disease in domestic ruminants worldwide. There is considerable variation in resistance to gastrointestinal nematodes within and between sheep breeds, which appears to be due to underlying genetic diversity. Through selection of resistant animals, rapid genetic progress has been demonstrated in both research and commercial flocks. Recent advances in genome sequencing and genomic technologies provide new opportunities to understand the ovine host response to gastrointestinal nematodes at the molecular level, and to identify polymorphisms conferring nematode resistance.

Genotype imputation from low-density (LD) to high-density single nucleotide polymorphism (SNP) chips is an important step before applying genomic selection, since denser chips tend to provide more reliable genomic predictions. Imputation methods rely partially on linkage disequilibrium between markers to infer unobserved genotypes. Bos indicus cattle (e.g. Nelore breed) are characterized, in general, by lower levels of linkage disequilibrium between genetic markers at short distances, compared to taurine breeds. Thus, it is important to evaluate the accuracy of imputation to better define which imputation method and chip are most appropriate for genomic applications in indicine breeds.

Single nucleotide polymorphisms (SNPs) are an abundant form of genetic variation in the genome of every species and are useful for gene mapping and association studies. Of particular interest are non-synonymous SNPs, which may alter protein function and phenotype. We therefore examined bovine expressed sequences for non-synonymous SNPs and validated and tested selected SNPs for their association with measured traits.

Cattle are important agriculturally and relevant as a model organism. Previously described genetic and radiation hybrid (RH) maps of the bovine genome have been used to identify genomic regions and genes affecting specific traits. Application of these maps to identify influential genetic polymorphisms will be enhanced by integration with each other and with bacterial artificial chromosome (BAC) libraries. The BAC libraries and clone maps are essential for the hybrid clone-by-clone/whole-genome shotgun sequencing approach taken by the bovine genome sequencing project.

Arctic charr (Salvelinus alpinus) is a species of high economic value for the aquaculture industry, and of high ecological value due to its Holarctic distribution in both marine and freshwater environments. Novel genome sequencing approaches enable the study of population and quantitative genetic parameters even on species with limited or no prior genomic resources. Low coverage genotyping by sequencing (GBS) was applied in a selected strain of Arctic charr in Sweden originating from a landlocked freshwater population. For the needs of the current study, animals from year classes 2013 (171 animals, parental population) and 2017 (759 animals; 13 full sib families) were used as a template for identifying genome wide single nucleotide polymorphisms (SNPs). GBS libraries were constructed using the PstI and MspI restriction enzymes. Approximately 14.5K SNPs passed quality control and were used for estimating a genomic relationship matrix. Thereafter a wide range of analyses were conducted in order to gain insights regarding genetic diversity and investigate the efficiency of the genomic information for parentage assignment and breeding value estimation. Heterozygosity estimates for both year classes suggested a slight excess of heterozygotes. Furthermore, FST estimates among the families of year class 2017 ranged between 0.009 - 0.066. Principal components analysis (PCA) and discriminant analysis of principal components (DAPC) were applied aiming to identify the existence of genetic clusters among the studied population. Results obtained were in accordance with pedigree records allowing the identification of individual families. Additionally, DNA parentage verification was performed, with results in accordance with the pedigree records with the exception of a putative dam where full sib genotypes suggested a potential recording error. Breeding value estimation for juvenile growth through the usage of the estimated genomic relationship matrix clearly outperformed the pedigree equivalent in terms of prediction accuracy (0.51 opposed to 0.31). Overall, low coverage GBS has proven to be a cost-effective genotyping platform that is expected to boost the selection efficiency of the Arctic charr breeding program.