Enteric fermentation by farmed ruminant animals is a major source of methane and constitutes the second largest anthropogenic contributor to global warming. Reducing methane emissions from ruminants is needed to ensure sustainable animal production in the future. Methane yield varies naturally in sheep and is a heritable trait that can be used to select animals that yield less methane per unit of feed eaten. We previously demonstrated elevated expression of hydrogenotrophic methanogenesis pathway genes of methanogenic archaea in the rumens of high methane yield (HMY) sheep compared to their low methane yield (LMY) counterparts. Methane production in the rumen is strongly connected to microbial hydrogen production through fermentation processes. In this study, we investigate the contribution that rumen bacteria make to methane yield phenotypes in sheep.

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.

Seasonal variation in the faecal microbiota of forage-fed horses was investigated over a 12-month period to determine whether the bacterial diversity fluctuated over time. Horses (n = 10) were maintained on pasture for one year, with hay supplemented from June to October. At monthly intervals, data were recorded on pasture availability and climate (collected continuously and averaged on monthly basis), pasture and hay samples were collected for nutrient analysis, and faecal samples were collected from all horses to investigate the diversity of faecal microbiota using next-generation sequencing on the Illumina MiSeq platform. The alpha diversity of bacterial genera was high in all samples (n = 118), with significantly higher Simpson's (p < 0.001) and Shannon-Wiener (p < 0.001) diversity indices observed during the months when horses were kept exclusively on pasture compared to the months when pasture was supplemented with hay. There were significant effects of diet, season, and month (ANOSIM, p < 0.01 for each comparison) on the beta diversity of bacterial genera identified in the faeces. While there was some inter-horse variation, hierarchical clustering of beta diversity indices showed separate clades originating for samples obtained during May, June, and July (late-autumn to winter period), and January, February, and March (a period of drought), with a strong association between bacterial taxa and specific nutrients (dry matter, protein, and structural carbohydrates) and climate variables (rainfall and temperature). Our study supports the hypothesis that the diversity and community structure of the faecal microbiota of horses kept on pasture varied over a 12-month period, and this variation reflects changes in the nutrient composition of the pasture, which in turn is influenced by climatic conditions. The findings of this study may have implications for grazing management and the preparation of conserved forages for those horses susceptible to perturbations of the hindgut microbiota.

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.

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.

Molecular microbial ecology techniques are widely used to study the composition of the rumen microbiota and to increase understanding of the roles they play. Therefore, sampling and DNA extraction methods that result in adequate yields of microbial DNA that also accurately represents the microbial community are crucial. Fifteen different methods were used to extract DNA from cow and sheep rumen samples. The DNA yield and quality, and its suitability for downstream PCR amplifications varied considerably, depending on the DNA extraction method used. DNA extracts from nine extraction methods that passed these first quality criteria were evaluated further by quantitative PCR enumeration of microbial marker loci. Absolute microbial numbers, determined on the same rumen samples, differed by more than 100-fold, depending on the DNA extraction method used. The apparent compositions of the archaeal, bacterial, ciliate protozoal, and fungal communities in identical rumen samples were assessed using 454 Titanium pyrosequencing. Significant differences in microbial community composition were observed between extraction methods, for example in the relative abundances of members of the phyla Bacteroidetes and Firmicutes. Microbial communities in parallel samples collected from cows by oral stomach-tubing or through a rumen fistula, and in liquid and solid rumen digesta fractions, were compared using one of the DNA extraction methods. Community representations were generally similar, regardless of the rumen sampling technique used, but significant differences in the abundances of some microbial taxa such as the Clostridiales and the Methanobrevibacter ruminantium clade were observed. The apparent microbial community composition differed between rumen sample fractions, and Prevotellaceae were most abundant in the liquid fraction. DNA extraction methods that involved phenol-chloroform extraction and mechanical lysis steps tended to be more comparable. However, comparison of data from studies in which different sampling techniques, different rumen sample fractions or different DNA extraction methods were used should be avoided.

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.

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.

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

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.

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.

The management of competition horses in New Zealand often involves rotations of short periods of stall confinement and concentrate feeding, with periods of time at pasture. Under these systems, horses may undergo abrupt dietary changes, with the incorporation of grains or concentrate feeds to the diet to meet performance needs, or sudden changes in the type of forage fed in response to a lack of fresh or conserved forage. Abrupt changes in dietary management are a risk factor for gastrointestinal (GI) disturbances, potentially due to the negative effects observed on the population of GI microbiota. In the present study, the faecal microbiota of horses was investigated to determine how quickly the bacterial communities; (1) responded to dietary change, and (2) stabilised following abrupt dietary transition. Six Thoroughbred mares were stabled for six weeks, consuming freshly cut pasture (weeks 1, 3 and 5), before being abruptly transitioned to conserved forage-based diets, both offered ad libitum. Intestinal markers were administered to measure digesta transit time immediately before each diet change. The conserved forage-based diets were fed according to a 3 × 3 Latin square design (weeks 2, 4 and 6), and comprised a chopped ensiled forage fed exclusively (Diet FE) or with whole oats (Diet FE + O), and perennial ryegrass hay fed with whole oats (Diet H + O). Faecal samples were collected at regular intervals from each horse following the diet changes. High throughput 16S rRNA gene sequencing was used to evaluate the faecal microbiota. There were significant differences in alpha diversity across diets (p < 0.001), and a significant effect of diet on the beta diversity (ANOSIM, p = 0.001), with clustering of samples observed by diet group. There were differences in the bacterial phyla across diets (p < 0.003), with the highest relative abundances observed for Firmicutes (62-64%) in the two diets containing chopped ensiled forage, Bacteroidetes (32-38%) in the pasture diets, and Spirochaetes (17%) in the diet containing hay. Major changes in relative abundances of faecal bacteria appeared to correspond with the cumulative percentage of intestinal markers retrieved in the faeces as the increasing amounts of digesta from each new diet transited the animals. A stable faecal microbiota profile was observed in the samples from 96 h after abrupt transition to the treatment diets containing ensiled chopped forage. The present study confirmed that the diversity and community structure of the faecal bacteria in horses is diet-specific and resilient following dietary transition and emphasised the need to have modern horse feeding management that reflects the ecological niche, particularly by incorporating large proportions of forage into equine diets.

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.

We report the whole genome of a strain of Ligilactobacillus faecis. The complete circular chromosome and plasmid of strain WILCCON 0062 were obtained through a combination of short- and long-read sequencing and may be used to derive unprecedented insights into the genome-level phylogeny and functional capacities of Ligilactobacillus faecis.

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.

The effects of abrupt dietary transition on the faecal microbiota of forage-fed horses over a 3-week period were investigated. Yearling Thoroughbred fillies reared as a cohort were exclusively fed on either an ensiled conserved forage-grain diet ("Group A"; n = 6) or pasture ("Group B"; n = 6) for three weeks prior to the study. After the Day 0 faecal samples were collected, horses of Group A were abruptly transitioned to pasture. Both groups continued to graze similar pasture for three weeks, with faecal samples collected at 4-day intervals. DNA was isolated from the faeces and microbial 16S and 18S rRNA gene amplicons were generated and analysed by pyrosequencing. The faecal bacterial communities of both groups of horses were highly diverse (Simpson's index of diversity > 0.8), with differences between the two groups on Day 0 (P < 0.017 adjusted for multiple comparisons). There were differences between Groups A and B in the relative abundances of four genera, BF311 (family Bacteroidaceae; P = 0.003), CF231 (family Paraprevotellaceae; P = 0.004), and currently unclassified members within the order Clostridiales (P = 0.003) and within the family Lachnospiraceae (P = 0.006). The bacterial community of Group A horses became similar to Group B within four days of feeding on pasture, whereas the structure of the archaeal community remained constant pre- and post-dietary change. The community structure of the faecal microbiota (bacteria, archaea and ciliate protozoa) of pasture-fed horses was also identified. The initial differences observed appeared to be linked to recent dietary history, with the bacterial community of the forage-fed horses responding rapidly to abrupt dietary change.

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

Proteases catalyze the hydrolysis of peptide bonds of proteins, thereby improving dietary protein digestibility, nutrient availability, as well as flavor and texture of fermented food and feed products. The lactobacilli Lactiplantibacillus plantarum (formerly Lactobacillus plantarum) and Pediococcus acidilactici are widely used in food and feed fermentations due to their broad metabolic capabilities and safe use. However, extracellular protease activity in these two species is low. Here, we optimized protease expression and secretion in L. plantarum and P. acidilactici via a genetic engineering strategy.

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.