Paul Maclean

Annotation of sheep keratin intermediate filament genes and their patterns of expression

Abstract:  Keratin IF (KRT) and keratin‐associated protein genes encode the majority of wool and hair proteins. We have identified cDNA sequences representing nine novel sheep KRT genes, increasing the known active genes from eight to 17, a number comparable to that in the human. However, the absence of KRT37 in the type I family and the discovery of type II KRT87 in sheep exemplify species‐specific compositional differences in hair KRT genes. Phylogenetic analysis of hair KRT genes within type I and type II families in the sheep, cattle and human genomes revealed a high degree of consistency in their sequence conservation and grouping. However, there were differences in the fibre compartmentalisation and keratinisation zones for the expression of six ovine KRT genes compared with their human orthologs. Transcripts of three genes (KRT40, KRT82 and KRT84) were only present in the fibre cuticle. KRT32, KRT35 and KRT85 were expressed in both the cuticle and the fibre cortex. The remaining 11 genes (KRT31, KRT33A, KRT33B, KRT34, KRT36, KRT38‐39, KRT81, KRT83 and KRT86‐87) were expressed only in the cortex. Species‐specific differences in the expressed keratin gene sets, their relative expression levels and compartmentalisation are discussed in the context of their underlying roles in wool and hair developmental programmes and the distinctive characteristics of the fibres produced.

Efficient introgression of allelic variants by embryo-mediated editing of the bovine genome

The recent development of designer nucleases allows for the efficient and precise introduction of genetic change into livestock genomes. Most studies so far have focused on the introduction of random mutations in cultured cells and the use of nuclear transfer to generate animals with edited genotypes. To circumvent the intrinsic uncertainties of random mutations and the inefficiencies of nuclear transfer we directed our efforts to the introduction of specific genetic changes by homology-driven repair directly in in vitro produced embryos. Initially, we injected zinc finger nuclease (ZFN)-encoding mRNA or DNA into bovine zygotes to verify cleavage activity at their target site within the gene for beta-lactoglobulin (LGB) and detected ZFN-induced random mutations in 30% to 80% of embryos. Next, to precisely change the LGB sequence, we co-injected ZFNs or transcription activator-like effector nucleases (TALENs) with DNA oligonucleotides (ODNs). Analysis of co-injected embryos showed targeted changes in up to 33% (ZFNs) and 46% (TALENs) of blastocysts. Deep sequence analysis of selected embryos revealed contributions of the targeted LGB allele can reach 100% which implies that genome editing by zygote injections can facilitate the one-step generation of non-mosaic livestock animals with pre-designed biallelic modifications.

Molecular Characterization of the Onset and Progression of Colitis in Inoculated Interleukin-10 Gene-Deficient Mice: A Role for PPARα

The interleukin-10 gene-deficient (Il10 −/−) mouse is a model of human inflammatory bowel disease and Ppara has been identified as one of the key genes involved in regulation of colitis in the bacterially inoculated Il10 −/− model. The aims were to (1) characterize colitis onset and progression using a histopathological, transcriptomic, and proteomic approach and (2) investigate links between PPARα and IL10 using gene network analysis. Bacterial inoculation resulted in severe colitis in Il10 −/− mice from 10 to 12 weeks of age. Innate and adaptive immune responses showed differences in gene expression relating to colitis severity. Actin cytoskeleton dynamics, innate immunity, and apoptosis-linked gene and protein expression data suggested a delayed remodeling process in 12-week-old Il10 −/− mice. Gene expression changes in 12-week-old Il10 −/− mice were related to PPARα signaling likely to control colitis, but how PPARα activation might regulate intestinal IL10 production remains to be determined.

Genome-wide DNA methylation patterns and transcription analysis in sheep muscle.

DNA methylation plays a central role in regulating many aspects of growth and development in mammals through regulating gene expression. The development of next generation sequencing technologies have paved the way for genome-wide, high resolution analysis of DNA methylation landscapes using methodology known as reduced representation bisulfite sequencing (RRBS). While RRBS has proven to be effective in understanding DNA methylation landscapes in humans, mice, and rats, to date, few studies have utilised this powerful method for investigating DNA methylation in agricultural animals. Here we describe the utilisation of RRBS to investigate DNA methylation in sheep Longissimus dorsi muscles. RRBS analysis of ∼1% of the genome from Longissimus dorsi muscles provided data of suitably high precision and accuracy for DNA methylation analysis, at all levels of resolution from genome-wide to individual nucleotides. Combining RRBS data with mRNAseq data allowed the sheep Longissimus dorsi muscle methylome to be compared with methylomes from other species. While some species differences were identified, many similarities were observed between DNA methylation patterns in sheep and other more commonly studied species. The RRBS data presented here highlights the complexity of epigenetic regulation of genes. However, the similarities observed across species are promising, in that knowledge gained from epigenetic studies in human and mice may be applied, with caution, to agricultural species. The ability to accurately measure DNA methylation in agricultural animals will contribute an additional layer of information to the genetic analyses currently being used to maximise production gains in these species.

Analyses of Hypomethylated Oil Palm Gene Space

Demand for palm oil has been increasing by an average of ∼8% the past decade and currently accounts for about 59% of the worlds vegetable oil market. This drives the need to increase palm oil production. Nevertheless, due to the increasing need for sustainable production, it is imperative to increase productivity rather than the area cultivated. Studies on the oil palm genome are essential to help identify genes or markers that are associated with important processes or traits, such as flowering, yield and disease resistance. To achieve this, 294,115 and 150,744 sequences from the hypomethylated or gene-rich regions of Elaeis guineensis and E. oleifera genome were sequenced and assembled into contigs. An additional 16,427 shot-gun sequences and 176 bacterial artificial chromosomes (BAC) were also generated to check the quality of libraries constructed. Comparison of these sequences revealed that although the methylation-filtered libraries were sequenced at low coverage, they still tagged at least 66% of the RefSeq supported genes in the BAC and had a filtration power of at least 2.0. A total 33,752 microsatellites and 40,820 high-quality single nucleotide polymorphism (SNP) markers were identified. These represent the most comprehensive collection of microsatellites and SNPs to date and would be an important resource for genetic mapping and association studies. The gene models predicted from the assembled contigs were mined for genes of interest, and 242, 65 and 14 oil palm transcription factors, resistance genes and miRNAs were identified respectively. Examples of the transcriptional factors tagged include those associated with floral development and tissue culture, such as homeodomain proteins, MADS, Squamosa and Apetala2. The E. guineensis and E. oleifera hypomethylated sequences provide an important resource to understand the molecular mechanisms associated with important agronomic traits in oil palm.

Links between sulphur oxidation and sulphur-oxidising bacteria abundance and diversity in soil microcosms based on soxB functional gene analysis

Sulphur-oxidising bacteria (SOB) play a key role in the biogeochemical cycling of sulphur in soil ecosystems. However, the ecology of SOB is poorly understood, and there is little knowledge about the taxa capable of sulphur oxidation, their distribution, habitat preferences and ecophysiology. Furthermore, as yet there are no conclusive links between SOB community size or structure and rates of sulphur oxidation. We have developed a molecular approach based on primer design targeting the soxB functional gene of nonfilamentous chemolithotrophic SOB that allows assessment of both abundance and diversity. Cloning and sequencing revealed considerable diversity of known soxB genotypes from agricultural soils and also evidence for previously undescribed taxa. In a microcosm experiment, abundance of soxB genes increased with sulphur oxidation rate in soils amended with elemental sulphur. Addition of elemental sulphur to soil had a significant effect in the soxB gene diversity, with the chemolithotrophic Thiobacillus-like Betaproteobacteria sequences dominating clone libraries 6 days after sulphur application. Using culture-independent methodology, the study provides evidence for links between abundance and diversity of SOB and sulphur oxidation. The methodology provides a new tool for investigation of the ecology and role of SOB in soil sulphur biogeochemistry.

Key bacterial families (Clostridiaceae, Erysipelotrichaceae and Bacteroidaceae) are related to the digestion of protein and energy in dogs

Background Much of the recent research in companion animal nutrition has focussed on understanding the role of diet on faecal microbiota composition. To date, diet-induced changes in faecal microbiota observed in humans and rodents have been extrapolated to pets in spite of their very different dietary and metabolic requirements. This lack of direct evidence means that the mechanisms by which microbiota influences health in dogs are poorly understood. We hypothesised that changes in faecal microbiota correlate with physiological parameters including apparent macronutrient digestibility. Methods Fifteen adult dogs were assigned to two diet groups, exclusively fed either a premium kibbled diet (kibble; K; n = 8) or a raw red meat diet (meat; M; n = 7) for nine weeks. Apparent digestibility of macronutrients (protein, fat, gross energy and dry matter), faecal weight, faecal health scores, faecal VFA concentrations and faecal microbial composition were determined. Datasets were integrated using mixOmics in R. Results Faecal weight and VFA levels were lower and the apparent digestibility of protein and energy were higher in dogs on the meat diet. Diet significantly affected 27 microbial families and 53 genera in the faeces. In particular, the abundances of Bacteriodes, Prevotella, Peptostreptococcus and Faecalibacterium were lower in dogs fed the meat diet, whereas Fusobacterium, Lactobacillus and Clostridium were all more abundant. Discussion Our results show clear associations of specific microbial taxa with diet composition. For example, Clostridiaceae, Erysipelotrichaceae and Bacteroidaceae were highly correlated to parameters such as protein and fat digestibility in the dog. By understanding the relationship between faecal microbiota and physiological parameters we will gain better insights into the effects of diet on the nutrition of our pets.

Gene Expression Changes in the Colon Epithelium Are Similar to Those of Intact Colon during Late Inflammation in Interleukin-10 Gene Deficient Mice

In addition to their role in absorption and secretion, epithelial cells play an important role in the protection of the colon mucosa from the resident microbiota and are important for the maintenance of homeostasis. Microarray analysis of intact colon samples is widely used to gain an overview of the cellular pathways and processes that are active in the colon during inflammation. Laser microdissection of colon epithelial cells allows a more targeted analysis of molecular pathways in the mucosa, preceding and during inflammation, with potentially increased sensitivity to changes in specific cell populations. The aim of this study was to investigate the molecular changes that occur in early and late inflammation stages in colon epithelium of a mouse model of inflammatory bowel diseases. Microarray analysis of intact colon samples and microdissected colon epithelial cell samples from interleukin-10 gene deficient and control mice at 6 and 12 weeks of age was undertaken. Results of gene set enrichment analysis showed that more immune-related pathways were identified between interleukin-10 gene deficient and control mice at 6 weeks of age in epithelial cells than intact colon. This suggests that targeting epithelial cells could increase sensitivity for detecting immune changes that occur early in the inflammatory process. However, in the later stages of inflammation, microarray analyses of intact colon and epithelium both provide a similar overview of gene expression changes in the colon mucosa at the pathway level.

Improving the genetic representation of rare taxa within complex microbial communities using DNA normalization methods

Complex microbial communities typically contain a large number of low abundance species, which collectively, comprise a considerable proportion of the community. This ‘rare biosphere’ has been speculated to contain keystone species and act as a repository of genomic diversity to facilitate community adaptation. Many environmental microbes are currently resistant to cultivation, and can only be accessed via culture‐independent approaches. To enhance our understanding of the role of the rare biosphere, we aimed to improve their metagenomic representation using DNA normalization methods, and assess normalization success via shotgun DNA sequencing. A synthetic metagenome was constructed from the genomic DNA of five bacterial species, pooled in a defined ratio spanning three orders of magnitude. The synthetic metagenome was fractionated and thermally renatured, allowing the most abundant sequences to hybridize. Double‐stranded DNA was removed either by hydroxyapatite chromatography, or by a duplex‐specific nuclease (DSN). The chromatographic method failed to enrich for the genomes present in low starting abundance, whereas the DSN method resulted in all genomes reaching near equimolar abundance. The representation of the rarest member was increased by approximately 450‐fold. De novo assembly of the normalized metagenome enabled up to 18.0% of genes from the rarest organism to be assembled, in contrast to the un‐normalized sample, where genes were not able to be assembled at the same sequencing depth. This study has demonstrated that the application of normalization methods to metagenomic samples is a powerful tool to enrich for sequences from rare taxa, which will shed further light on their ecological niches.

KDM4B-mediated reduction of H3K9me3 and H3K36me3 levels improves somatic cell reprogramming into pluripotency

Correct reprogramming of epigenetic marks is essential for somatic cells to regain pluripotency. Repressive histone (H) lysine (K) methylation marks are known to be stable and difficult to reprogram. In this study, we generated transgenic mice and mouse embryonic fibroblasts (MEFs) for the inducible expression of KDM4B, a demethylase that removes H3 K9 and H3K36 trimethylation (me3) marks (H3K9/36me3). Upon inducing Kdm4b, H3K9/36me3 levels significantly decreased compared to non-induced controls. Concurrently, H3K9me1 levels significantly increased, while H3K9me2 and H3K27me3 remained unchanged. The global transcriptional impact of Kdm4b-mediated reduction in H3K9/36me3 levels was examined by comparative microarray analysis and mRNA-sequencing of three independent transgenic MEF lines. We identified several commonly up-regulated targets, including the heterochromatin-associated zinc finger protein 37 and full-length endogenous retrovirus repeat elements. Following optimized zona-free somatic nuclear transfer, reduced H3K9/36me3 levels were restored within hours. Nevertheless, hypo-methylated Kdm4b MEF donors reprogrammed six-fold better into cloned blastocysts than non-induced donors. They also reprogrammed nine-fold better into induced pluripotent stem cells that gave rise to teratomas and chimeras. In summary, we firmly established H3K9/36me3 as a major roadblock to somatic cell reprogramming and identified transcriptional targets of derestricted chromatin that could contribute towards improving this process in mouse.