Suzanne Rowe

Utilization of high throughput genome sequencing technology for large scale single nucleotide polymorphism discovery in red deer and Canadian elk

Deer farming is a significant international industry. For genetic improvement, using genomic tools, an ordered array of DNA variants and associated flanking sequence across the genome is required. This work reports a comparative assembly of the deer genome and subsequent DNA variant identification. Next generation sequencing combined with an existing bovine reference genome enabled the deer genome to be assembled sufficiently for large-scale SNP discovery. In total, 28 Gbp of sequence data were generated from seven Cervus elaphus (European red deer and Canadian elk) individuals. After aligning sequence to the bovine reference genome build UMD 3.0 and binning reads into one Mbp groups; reads were assembled and analyzed for SNPs. Greater than 99% of the non-repetitive fraction of the bovine genome was covered by deer chromosomal scaffolds. We identified 1.8 million SNPs meeting Illumina InfiniumII SNP chip technical threshold. Markers on the published Red x Pere David deer linkage map were aligned to both UMD3.0 and the new deer chromosomal scaffolds. This enabled deer linkage groups to be assigned to deer chromosomal scaffolds, although the mapping locations remain based on bovine order. Genotyping of 270 SNPs on a Sequenom MS system showed that 88% of SNPs identified could be amplified. Also, inheritance patterns showed no evidence of departure from Hardy-Weinberg equilibrium. A comparative assembly of the deer genome, alignment with existing deer genetic linkage groups and SNP discovery has been successfully completed and validated facilitating application of genomic technologies for subsequent deer genetic improvement.

Gene network analysis identifies rumen epithelial cell proliferation, differentiation and metabolic pathways perturbed by diet and correlated with methane production

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.

Linkage Disequilibrium Estimation in Low Coverage High-Throughput Sequencing Data

High-throughput sequencing methods provide a cost-effective approach for genotyping and are commonly used in population genetics studies. A drawback of these methods, however, is that sequencing and genotyping errors can arise... High-throughput sequencing methods that multiplex a large number of individuals have provided a cost-effective approach for discovering genome-wide genetic variation in large populations. These sequencing methods are increasingly being utilized in population genetic studies across a diverse range of species. Two side-effects of these methods, however, are (1) sequencing errors and (2) heterozygous genotypes called as homozygous due to only one allele at a particular locus being sequenced, which occurs when the sequencing depth is insufficient. Both of these errors have a profound effect on the estimation of linkage disequilibrium (LD) and, if not taken into account, lead to inaccurate estimates. We developed a new likelihood method, GUS-LD, to estimate pairwise linkage disequilibrium using low coverage sequencing data that accounts for undercalled heterozygous genotypes and sequencing errors. Our findings show that accurate estimates were obtained using GUS-LD, whereas underestimation of LD results if no adjustment is made for the errors.

Genome-wide association study of lung lesions and pleurisy in New Zealand lambs1

Abstract Pneumonia is an important issue for sheep production, leading to reduced growth rate and a predisposition to pleurisy. The objective of this study was to identify loci associated with pneumonic lesions and pleurisy in New Zealand progeny test lambs. The lungs from 3,572 progeny-test lambs were scored for presence and severity of pneumonic lesions and pleurisy at slaughter. Animals were genotyped using the Illumina Ovine Infinium HD SNP BeadChip (606,006 markers). The heritability of lung lesion score and pleurisy were calculated using the genomic relationship matrix, and genome-wide association analyses were conducted using EMMAX and haplotype trend regression. At slaughter, 35% of lambs had pneumonic lesions, with 9% showing lesions on more than half of any individual lobe. The number of lambs recorded as having pleurisy by the processing plants was 9%. Heritability estimates for pneumonic lesions and pleurisy scores adjusted for heteroscedasticity (CPSa and PLEURa) were 0.16 (± 0.03) and 0.05 (± 0.02), respectively. Five single-nucleotide polymorphisms (SNPs) were significantly associated with pneumonic lesions at the genome-wide level, and additional 37 SNPs were suggestively significant. Four SNPs were significantly associated with pleurisy, with an additional 11 SNPs reaching the suggestive level of significance. There were no regions that overlapped between the 2 traits. Multiple SNPs were in regions that contained genes involved in either the DNA damage response or the innate immune response, including several that had previously been reported to have associations with respiratory disease. Both EMMAX and HTR analyses of pleurisy data showed a significant peak on chromosome 2, located downstream from the transcription factor SP3. SP3 activates or suppresses the expression of numerous genes, including several genes with known functions in the immune system. This study identified several SNPs associated with genes involved in both the innate immune response and the response to DNA damage that are associated with pneumonic lesions and pleurisy in lambs at slaughter. Additionally, the identification in sheep of several SNPs within genes that have previously been associated with the respiratory system in cattle, pigs, rats, and mice indicates that there may be common pathways that underlie the response to invasion by respiratory pathogens in multiple species.

Genetic parameters of methane emissions determined using portable accumulation chambers in lambs and ewes grazing pasture and genetic correlations with emissions determined in respiration chambers1

Abstract Methane (CH4) emission traits were previously found to be heritable and repeatable in sheep fed alfalfa pellets in respiration chambers (RC). More rapid screening methods are, however, required to increase genetic progress and to provide a cost-effective method to the farming industry for maintaining the generation of breeding values in the future. The objective of the current study was to determine CH4 and carbon dioxide (CO2) emissions using several 1-h portable accumulation chamber (PAC) measurements from lambs and again as ewes while grazing ryegrass-based pasture. Many animals with PAC measurements were also measured in RC while fed alfalfa pellets at 2.0 × maintenance metabolizable energy requirements (MEm). Heritability estimates from mixed models for CH4 and CO2 production (g/d) were 0.19 and 0.16, respectively, when measured using PAC with lambs; 0.20 and 0.27, respectively, when measured using PAC with ewes; and 0.23 and 0.34, respectively, when measured using RC with lambs. For measured gas traits, repeatabilities of measurements collected 14 d apart ranged from 0.33 to 0.55 for PAC (combined lambs and ewes) and were greater at 0.65 to 0.76 for the same traits measured using RC. Genetic correlations (rg) between PAC in lambs and ewes were 0.99 for CH4, 0.93 for CH4 + CO2, and 0.85 for CH4/(CH4 + CO2), suggesting that CH4 emissions in lambs and ewes are the same trait. Genetic correlations between PAC and RC measurements were lower, at 0.62 to 0.67 for CH4 and 0.41 to 0.42 for CH4 + CO2, likely reflecting different environmental conditions associated with the protocols used with the 2 measurement methods. The CH4/(CH4 + CO2) ratio was the most similar genetic trait measured using PAC (both lambs and ewes, 63% and 66% selection efficiency, respectively) compared with CH4 yield (g/kg DMI) measured using RC. These results suggest that PAC measurements have considerable value as a rapid low-cost method to estimate breeding values for CH4 emissions in sheep.

Across-Experiment Transcriptomics of Sheep Rumen Identifies Expression of Lipid/Oxo-Acid Metabolism and Muscle Cell Junction Genes Associated With Variation in Methane-Related Phenotypes

Ruminants are significant contributors to the livestock generated component of the greenhouse gas, methane (CH4). The CH4 is primarily produced by the rumen microbes. Although the composition of the diet and animal intake amount have the largest effect on CH4 production and yield (CH4 production/dry matter intake, DMI), the host also influences CH4 yield. Shorter rumen feed mean retention time (MRT) is associated with higher dry matter intake and lower CH4 yield, but the molecular mechanism(s) by which the host affects CH4 production remain unclear. We integrated rumen wall transcriptome data and CH4 phenotypes from two independent experiments conducted with sheep in Australia (AUS, n = 62) and New Zealand (NZ, n = 24). The inclusion of the AUS data validated the previously identified clusters and gene sets representing rumen epithelial, metabolic and muscular functions. In addition, the expression of the cell cycle genes as a group was consistently positively correlated with acetate and butyrate concentrations (p < 0.05, based on AUS and NZ data together). The expression of a group of metabolic genes showed positive correlations in both AUS and NZ datasets with CH4 production (p < 0.05) and yield (p < 0.01). These genes encode key enzymes in the ketone body synthesis pathway and included members of the poorly characterized aldo-keto reductase 1C (AKR1C) family. Several AKR1C family genes appear to have ruminant specific evolution patterns, supporting their specialized roles in the ruminants. Combining differential gene expression in the rumen wall muscle of the shortest and longest MRT AUS animals (no data available for the NZ animals) with correlation and network analysis, we identified a set of rumen muscle genes involved in cell junctions as potential regulators of MRT, presumably by influencing contraction rates of the smooth muscle component of the rumen wall. Higher rumen expression of these genes, including SYNPO (synaptopodin, p < 0.01) and NEXN (nexilin, p < 0.05), was associated with lower CH4 yield in both AUS and NZ datasets. Unlike the metabolic genes, the variations in the expression of which may reflect the availability of rumen metabolites, the muscle genes are currently our best candidates for causal genes that influence CH4 yield.