Long read isoform sequencing reveals hidden transcriptional complexity between cattle subspecies.

The Iso-Seq method of full-length cDNA sequencing is suitable to quantify differentially expressed genes (DEGs), transcripts (DETs) and transcript usage (DTU). However, the higher cost of Iso-Seq relative to RNA-seq has limited the comparison of both methods. Transcript abundance estimated by RNA-seq and deep Iso-Seq data for fetal liver from two cattle subspecies were compared to evaluate concordance. Inter-sample correlation of gene- and transcript-level abundance was higher within technology than between technologies. Identification of DEGs between the cattle subspecies depended on sequencing method with only 44 genes identified by both that included 6 novel genes annotated by Iso-Seq. There was a pronounced difference between Iso-Seq and RNA-seq results at transcript-level wherein Iso-Seq revealed several magnitudes more transcript abundance and usage differences between subspecies. Factors influencing DEG identification included size selection during Iso-Seq library preparation, average transcript abundance, multi-mapping of RNA-seq reads to the reference genome, and overlapping coordinates of genes. Some DEGs called by RNA-seq alone appear to be sequence duplication artifacts. Among the 44 DEGs identified by both technologies some play a role in immune system, thyroid function and cell growth. Iso-Seq revealed hidden transcriptional complexity in DEGs, DETs and DTU genes between cattle subspecies previously missed by RNA-seq.

Genome-wide association study for antibody response to Mycobacterium avium ssp. paratubeculosis in goats.

Mycobacterium avium ssp. paratuberculosis (MAP), causes Johne's disease (JD), or paratuberculosis, a chronic enteritis of ruminants, which in goats is characterized by ileal lesions. The work described here is a case-control association study using the Illumina Caprine SNP50 BeadChip to unravel the genes involved in susceptibility of goats to JD. Goats in herds with a high occurrence of Johne's disease were classified as healthy or infected based on the level of serum antibodies against MAP, and 331 animals were selected for the association study. Goats belonged to the Jonica (157) and Siriana breeds (174). Whole-genome association analysis identified one region suggestive of significance associated with an antibody response to MAP on chromosome 7 (p-value = 1.23 × 10-5 ). These results provide evidence for genetic loci involved in the antibody response to MAP in goats.

Analysis of immune-related microRNAs in cows and newborn calves.

Bovine colostrum contains a high concentration of immune-related microRNAs (miRNAs) that are packaged in exosomes and are very stable. In this study, 5 immune-related miRNAs (miR-142-5p, miR-150, miR-155, miR-181a, and miR-223) were quantified in dam blood, colostrum, and calf blood using reverse transcription quantitative PCR. Their levels in calf blood after colostrum ingestion were investigated to assess whether miRNAs are transferred from the dam to newborn calves. Three groups of Holstein-Friesian bull calves were bottle-fed 2 L of colostrum or milk from different sources twice per day. The group A calves received colostrum from their own dam and the group B calves were fed foster dam colostrum. Each pair of group A and group B calves were fed identical colostrum from the same milking of the corresponding group A dam for 3 d and then bulk tank milk for 7 d after birth. Group C calves were fed only 2L of "pooled colostrum" from multiple dams d 0 to 4 postpartum, and then fed bulk tank milk thereafter for 7 d after birth. The groups were fed colostrum from different sources and different amounts to assess possible miRNA absorption from the colostrum. All miRNAs were at the highest level in colostrum at d 0 and then decreased rapidly after d 1. The level of miR-150 had the largest decrease from 489 × 106 copies/µL (d 0) to 78 × 106 copies/µL (d 1). MicroRNA-223 and miR-155 were the most abundant in both colostrum and milk. Dam colostrum had significantly higher levels of miR-142-5p, miR-155, and miR-181a than the bulk tank milk. However, only the miR-155 concentration was significantly higher in the dam colostrum than in the pooled colostrum. The concentrations of miRNAs in the colostrum were less than in the cow blood (100- to 1,000-fold less). There was no significant correlation between the level of miRNAs in the dam blood and their colostrum, suggesting that miRNA is synthesized locally by the mammary gland rather than being transferred from the blood. MicroRNA-223 had the highest level in both calf and cow blood compared with the other 4 immune-related miRNAs. Calves were born with high levels of immune-related miRNAs in their blood, and there were no significant differences in miRNA levels between the 3 calf groups at birth or after they were fed different colostrum. This suggests that these miRNAs were not transferred from the colostrum to the newborn calves.

Gaur genome reveals expansion of sperm odorant receptors in domesticated cattle.

BACKGROUND: The gaur (Bos gaurus) is the largest extant wild bovine species, native to South and Southeast Asia, with unique traits, and is listed as vulnerable by the International Union for Conservation of Nature (IUCN). RESULTS: We report the first gaur reference genome and identify three biological pathways including lysozyme activity, proton transmembrane transporter activity, and oxygen transport with significant changes in gene copy number in gaur compared to other mammals. These may reflect adaptation to challenges related to climate and nutrition. Comparative analyses with domesticated indicine (Bos indicus) and taurine (Bos taurus) cattle revealed genomic signatures of artificial selection, including the expansion of sperm odorant receptor genes in domesticated cattle, which may have important implications for understanding selection for male fertility. CONCLUSIONS: Apart from aiding dissection of economically important traits, the gaur genome will also provide the foundation to conserve the species.

Adaptive selection signatures in river buffalo with emphasis on immune and major histocompatibility complex genes.

River buffalo is an agriculturally important species with many traits, such as disease tolerance, which promote its use worldwide. Highly contiguous genome assemblies of the river buffalo, goat, pig, human and two cattle subspecies were aligned to study gene gains and losses and signs of positive selection. The gene families that have changed significantly in river buffalo since divergence from cattle play important roles in protein degradation, the olfactory receptor system, detoxification and the immune system. We used the branch site model in PAML to analyse single-copy orthologs to identify positively selected genes that may be involved in skin differentiation, mammary development and bone formation in the river buffalo branch. The high contiguity of the genomes enabled evaluation of differences among species in the major histocompatibility complex. We identified a Babesia-like L1 LINE insertion in the DRB1-like gene in the river buffalo and discuss the implication of this finding.

Distinctive gene expression patterns and imprinting signatures revealed in reciprocal crosses between cattle sub-species.

BACKGROUND: There are two genetically distinct subspecies of cattle, Bos taurus taurus and Bos taurus indicus, which arose from independent domestication events. The two types of cattle show substantial phenotypic differences, some of which emerge during fetal development and are reflected in birth outcomes, including birth weight. We explored gene expression profiles in the placenta and four fetal tissues at mid-gestation from one taurine (Bos taurus taurus; Angus) and one indicine (Bos taurus indicus; Brahman) breed and their reciprocal crosses. RESULTS: In total 120 samples were analysed from a pure taurine breed, an indicine breed and their reciprocal cross fetuses, which identified 6456 differentially expressed genes (DEGs) between the two pure breeds in at least one fetal tissue of which 110 genes were differentially expressed in all five tissues examined. DEGs shared across tissues were enriched for pathways related to immune and stress response functions. Only the liver had a substantial number of DEGs when reciprocal crossed were compared among which 310 DEGs were found to be in common with DEGs identified between purebred livers; these DEGs were significantly enriched for metabolic process GO terms. Analysis of DEGs across purebred and crossbred tissues suggested an additive expression pattern for most genes, where both paternal and maternal alleles contributed to variation in gene expression levels. However, expression of 5% of DEGs in each tissue was consistent with parent of origin effects, with both paternal and maternal dominance effects identified. CONCLUSIONS: These data identify candidate genes potentially driving the tissue-specific differences between these taurine and indicine breeds and provide a biological insight into parental genome effects underlying phenotypic differences in bovine fetal development.

The distribution of runs of homozygosity in the genome of river and swamp buffaloes reveals a history of adaptation, migration and crossbred events.

BACKGROUND: Water buffalo is one of the most important livestock species in the world. Two types of water buffalo exist: river buffalo (Bubalus bubalis bubalis) and swamp buffalo (Bubalus bubalis carabanensis). The buffalo genome has been recently sequenced, and thus a new 90 K single nucleotide polymorphism (SNP) bead chip has been developed. In this study, we investigated the genomic population structure and the level of inbreeding of 185 river and 153 swamp buffaloes using runs of homozygosity (ROH). Analyses were carried out jointly and separately for the two buffalo types. RESULTS: The SNP bead chip detected in swamp about one-third of the SNPs identified in the river type. In total, 18,116 ROH were detected in the combined data set (17,784 SNPs), and 16,251 of these were unique. ROH were present in both buffalo types mostly detected (~ 59%) in swamp buffalo. The number of ROH per animal was larger and genomic inbreeding was higher in swamp than river buffalo. In the separated datasets (46,891 and 17,690 SNPs for river and swamp type, respectively), 19,760 and 10,581 ROH were found in river and swamp, respectively. The genes that map to the ROH islands are associated with the adaptation to the environment, fitness traits and reproduction. CONCLUSIONS: Analysis of ROH features in the genome of the two water buffalo types allowed their genomic characterization and highlighted differences between buffalo types and between breeds. A large ROH island on chromosome 2 was shared between river and swamp buffaloes and contained genes that are involved in environmental adaptation and reproduction.

Colostrum source and passive immunity transfer in dairy bull calves.

Colostrum is essential for good neonate health; however, it is not known whether different calves absorb the nutrients from colostrum equally well. In this study, the absorption of protein, IgG, and γ-glutamyl transferase was compared in newborn dairy bull calves for 1 wk after feeding colostrum from different sources. Thirty-five Holstein-Friesian bull calves were randomly allocated into 3 groups and fed colostrum within 4 h after birth. Group A calves (n = 12) were bottle fed colostrum from their own dam for 3 d. Colostrum from these group A cows was also used as foster cow colostrum for the group B calves (n = 12), such that each group A and B calf pair received identical colostrum from each milking of the respective group A dam (10% of birth weight per day). The group C calves (n = 11) were fed 1 bottle (2 L) of pooled colostrum and transition milk (referred to as pooled colostrum), as was the standard practice on the dairy farm. The pooled colostrum was collected from the other dairy cows on the farm 0 to 4 d postpartum and stored at 4°C for less than 12 h. Blood was sampled from calves before the first feeding and at 1, 2, 3, and 7 d after birth. Levels of total solids, total protein, and IgG were higher in the dam colostrum than in the pooled colostrum. At birth, there were no differences between the calf groups for any measurements, and all calves had very low IgG levels. After receiving colostrum, the glucose, plasma γ-glutamyl transferase, serum total protein, and IgG concentrations increased significantly in all calves. There were no differences in any blood measurements at any time point between the pairs of group A and group B calves that received colostrum from the same cow except for the IgG concentration 2 d after birth. However, the group A calves had a higher total serum protein level and IgG concentration than the group C calves for all the time points after the first feeding. The group B calves had a higher IgG concentration than the group C calves on d 1, 2, and 7 after birth. Compared with groups A and B, there was no difference in the proportion of calves in group C that failed to have passive immunity transferred adequately based on the IgG threshold (<10 g/L). Thus, the calves receiving identical colostrum from the same cow had the same levels of IgG, and even the pooled colostrum provided sufficient transfer of IgG as the calves were fed within 4 h after birth.

Predicting colostrum and calf blood components based on refractometry.

Provision of good quality colostrum is essential for the passive immunity and nutrition of newborn calves. In order to better predict the quality of colostrum and the transfer of passive immunity, the relationships between colostrum components and between calf serum components were examined in this study. Samples of bulk tank milk, colostrum pooled from several cows 0-4 d postpartum, and colostrum collected from individual cows twice daily for 3 d post-partum were compared. With the exception of fat percentage, there were strong correlations between the levels of the components in the pooled colostrum and in the individual cow colostrum collected 0-1 d postpartum. The correlations between total solids as measured by Brix refractometry and total protein, immunoglobulin G (IgG), lactose % and protein % in colostrum within 1 d postpartum and pooled colostrum were 0.92, 0.90, -0.88 and 0.98, respectively. These high correlations enabled these colostrum components to be accurately predicted from Brix % and therefore, the volume of colostrum required to feed neonate calves can be optimised based on Brix refractometry to avoid failure of passive immunity transfer. To assess whether the components obtained from colostrum were correlated in calf blood, newborn calves were separated from their dams before suckling and blood sampled before feeding (day 0), and on days 1 and 7, after receiving colostrum or milk twice a day. The correlations between glucose, total protein, IgG, and gamma-glutamyl transferase (GGT) levels in the calf blood were lower than the correlations observed between the colostrum components. The highest correlation was between serum protein measured by refractometer and serum IgG within one week postpartum. GGT activity was not a good indicator of serum IgG levels. However, serum protein refractometer measurements predicted serum IgG level with high accuracy, providing an on-farm test to determine that calves have received sufficient passive immunity and colostrum components.

Near-Infrared Reflectance Spectroscopy for Predicting the Phospholipid Fraction and the Total Fatty Acid Composition of Freeze-Dried Beef.

Research on fatty acids (FA) is important because their intake is related to human health. NIRS can be a useful tool to estimate the FA of beef but due to the high moisture and the high absorbance of water makes it difficult to calibrate the analyses. This work evaluated near-infrared reflectance spectroscopy as a tool to assess the total fatty acid composition and the phospholipid fraction of fatty acids of beef using freeze-dried meat. An average of 22 unrelated pure breed young bulls from 15 European breeds were reared on a common concentrate-based diet. A total of 332 longissimus thoracis steaks were analysed for fatty acid composition and a freeze-dried sample was subjected to near-infrared spectral analysis. 220 samples (67%) were used as a calibration set with the remaining 110 (33%) being used for validation of the models obtained. There was a large variation in the total FA concentration across the animals giving a good data set for the analysis and whilst the coefficient of variation was nearly 68% for the monounsaturated FA it was only 27% for the polyunsaturated fatty acids (PUFA). PLS method was used to develop the prediction models. The models for the phospholipid fraction had a low R2p and high standard error, while models for neutral lipid had the best performance, in general. It was not possible to obtain a good prediction of many individual PUFA concentrations being present at low concentrations and less variable than other FA. The best models were developed for Total FA, saturated FA, 9c18:1 and 16:1 with R2p greater than 0.76. This study indicates that NIRS is a feasible and useful tool for screening purposes and it has the potential to predict most of the FA of freeze-dried beef.

Maternal peri-conceptional undernourishment perturbs offspring sperm methylome.

The genotype of an organism is stable throughout its life; however, its epigenome is dynamic and can be altered in response to environmental factors, such as diet. Inheritance of acquired epigenetic modifications by the next generation occurs through the germline, although the precise mechanisms remain to be elucidated. Here, we used a sheep model to evaluate if modification of the maternal diet (CTR; control, UND: undernutrition; FA: undernutrition and folic acid supplementation) during the peri-conceptional period affects the genome-wide methylation status of the gametes of male offspring. Sperm DNA methylation, measured by Reduced Representation Bisulfite Sequencing (RRBS), identified Differentially Methylated Regions (DMR) in offspring that experienced in utero undernutrition, both in UND (244) and FA (240), compared with CTR. Gene ontology (GO) analysis identified DMRs in categories related to sperm function, therefore we investigated whether the fertilizing capacity of the semen from the three groups differed in an in vitro fertilization assay. Spermatozoa from the undernourished groups showed lower motility and sperm chromatin structure abnormalities, represented by a higher percentage of DNA fragmentation and an increased number of immature cells, compared with CTR. While good quality blastocysts were obtained from all three groups, the proportion of embryos reaching the blastocyst stage was reduced in the UND vs CTR, an effect partially rescued by the FA treatment. The data reported here show that nutritional stress during early pregnancy leads to epigenetic modifications in the semen of the resulting offspring, the effects of which in next generation remain to be elucidated.

Comparative sequence alignment reveals River Buffalo genomic structural differences compared with cattle.

This study sought to characterize differences in gene content, regulation and structure between taurine cattle and river buffalo (one subspecies of domestic water buffalo) using the extensively annotated UMD3.1 cattle reference genome as a basis for comparisons. We identified 127 deletion CNV regions in river buffalo representing 5 annotated cattle genes. We also characterized 583 merged mobile element insertion (MEI) events within the upstream regions of annotated cattle genes. Transcriptome analysis in various tissue types on river buffalo confirmed the absence of four cattle genes. Four genes which may be related to phenotypic differences in meat quality and color, had upstream MEI predictions and were found to have significantly elevated expression in river buffalo compared with cattle. Our comparative alignment approach and gene expression analyses suggested a functional role for many genomic structural variations, which may contribute to the unique phenotypes of river buffalo.

New insights into mammalian sex chromosome structure and evolution using high-quality sequences from bovine X and Y chromosomes.

BACKGROUND: Mammalian X chromosomes are mainly euchromatic with a similar size and structure among species whereas Y chromosomes are smaller, have undergone substantial evolutionary changes and accumulated male specific genes and genes involved in sex determination. The pseudoautosomal region (PAR) is conserved on the X and Y and pair during meiosis. The structure, evolution and function of mammalian sex chromosomes, particularly the Y chromsome, is still poorly understood because few species have high quality sex chromosome assemblies. RESULTS: Here we report the first bovine sex chromosome assemblies that include the complete PAR spanning 6.84 Mb and three Y chromosome X-degenerate (X-d) regions. The PAR comprises 31 genes, including genes that are missing from the X chromosome in current cattle, sheep and goat reference genomes. Twenty-nine PAR genes are single-copy genes and two are multi-copy gene families, OBP, which has 3 copies and BDA20, which has 4 copies. The Y chromosome X-d1, 2a and 2b regions contain 11, 2 and 2 gametologs, respectively. CONCLUSIONS: The ruminant PAR comprises 31 genes and is similar to the PAR of pig and dog but extends further than those of human and horse. Differences in the pseudoautosomal boundaries are consistent with evolutionary divergence times. A bovidae-specific expansion of members of the lipocalin gene family in the PAR reported here, may affect immune-modulation and anti-inflammatory responses in ruminants. Comparison of the X-d regions of Y chromosomes across species revealed that five of the X-Y gametologs, which are known to be global regulators of gene activity and candidate sexual dimorphism genes, are conserved.

Computational detection and experimental validation of segmental duplications and associated copy number variations in water buffalo ( Bubalus bubalis ).

Duplicated sequences are an important source of gene evolution and structural variation within mammalian genomes. Using a read depth approach based on next-generation sequencing, we performed a genome-wide analysis of segmental duplications (SDs) and associated copy number variations (CNVs) in the water buffalo (Bubalus bubalis). By aligning short reads of Olimpia (the reference water buffalo) to the UMD3.1 cattle genome, we identified 1,038 segmental duplications comprising 44.6 Mb (equivalent to ~1.73% of the cattle genome) of the autosomal and X chromosomal sequence in the buffalo genome. We experimentally validated 70.3% (71/101) of these duplications using fluorescent in situ hybridization. We also detected a total of 1,344 CNV regions across 14 additional water buffaloes, amounting to 59.8 Mb of variable sequence or the equivalent of 2.2% of the cattle genome. The CNV regions overlap 1,245 genes that are significantly enriched for specific biological functions including immune response, oxygen transport, sensory system and signal transduction. Additionally, we performed array Comparative Genomic Hybridization (aCGH) experiments using the 14 water buffaloes as test samples and Olimpia as the reference. Using a linear regression model, a high Pearson correlation (r = 0.781) was observed between the log2 ratios between copy number estimates and the log2 ratios of aCGH probes. We further designed Quantitative PCR assays to confirm CNV regions within or near annotated genes and found 74.2% agreement with our CNV predictions. These results confirm sub-chromosome-scale structural rearrangements present in the cattle and water buffalo. The information on genome variation that will be of value for evolutionary and phenotypic studies, and may be useful for selective breeding of both species.

Variation in the Early Host-Pathogen Interaction of Bovine Macrophages with Divergent Mycobacterium bovis Strains in the United Kingdom.

Bovine tuberculosis has been an escalating animal health issue in the United Kingdom since the 1980s, even though control policies have been in place for over 60 years. The importance of the genetics of the etiological agent, Mycobacterium bovis, in the reemergence of the disease has been largely overlooked. We compared the interaction between bovine monocyte-derived macrophages (bMDM) and two M. bovis strains, AF2122/97 and G18, representing distinct genotypes currently circulating in the United Kingdom. These M. bovis strains exhibited differences in survival and growth in bMDM. Although uptake was similar, the number of viable intracellular AF2122/97 organisms increased rapidly, while G18 growth was constrained for the first 24 h. AF2122/97 infection induced a greater transcriptional response by bMDM than G18 infection with respect to the number of differentially expressed genes and the fold changes measured. AF2122/97 infection induced more bMDM cell death, with characteristics of necrosis and apoptosis, more inflammasome activation, and a greater type I interferon response than G18. In conclusion, the two investigated M. bovis strains interact in significantly different ways with the host macrophage. In contrast to the relatively silent infection by G18, AF2122/97 induces greater signaling to attract other immune cells and induces host cell death, which may promote secondary infections of naive macrophages. These differences may affect early events in the host-pathogen interaction, including granuloma development, which could in turn alter the progression of the disease. Therefore, the potential involvement of M. bovis genotypes in the reemergence of bovine tuberculosis in the United Kingdom warrants further investigation.

A genome-wide scan for signatures of selection in Azeri and Khuzestani buffalo breeds.

BACKGROUND: Identification of genomic regions that have been targets of selection may shed light on the genetic history of livestock populations and help to identify variation controlling commercially important phenotypes. The Azeri and Kuzestani buffalos are the most common indigenous Iranian breeds which have been subjected to divergent selection and are well adapted to completely different regions. Examining the genetic structure of these populations may identify genomic regions associated with adaptation to the different environments and production goals. RESULTS: A set of 385 water buffalo samples from Azeri (N = 262) and Khuzestani (N = 123) breeds were genotyped using the Axiom® Buffalo Genotyping 90 K Array. The unbiased fixation index method (FST) was used to detect signatures of selection. In total, 13 regions with outlier FST values (0.1%) were identified. Annotation of these regions using the UMD3.1 Bos taurus Genome Assembly was performed to find putative candidate genes and QTLs within the selected regions. Putative candidate genes identified include FBXO9, NDFIP1, ACTR3, ARHGAP26, SERPINF2, BOLA-DRB3, BOLA-DQB, CLN8, and MYOM2. CONCLUSIONS: Candidate genes identified in regions potentially under selection were associated with physiological pathways including milk production, cytoskeleton organization, growth, metabolic function, apoptosis and domestication-related changes include immune and nervous system development. The QTL identified are involved in economically important traits in buffalo related to milk composition, udder structure, somatic cell count, meat quality, and carcass and body weight.

Integrative analysis of methylomic and transcriptomic data in fetal sheep muscle tissues in response to maternal diet during pregnancy.

BACKGROUND: Numerous studies have established a link between maternal diet and the physiological and metabolic phenotypes of their offspring. In previous studies in sheep, we demonstrated that different maternal diets altered the transcriptome of fetal tissues. However, the mechanisms underlying transcriptomic changes are poorly understood. DNA methylation is an epigenetic mark regulating transcription and is largely influenced by dietary components of the one-carbon cycle that generate the methyl group donor, SAM. Therefore, in the present study, we tested whether different maternal diets during pregnancy would alter the DNA methylation and gene expression patterns in fetal tissues. RESULTS: Pregnant ewes were randomly divided into two groups which received either hay or corn diet from mid-gestation (day 67 ± 5) until day 131 ± 1 when fetuses were collected by necropsy. A total of 1516 fetal longissimus dorsi (LD) tissues were used for DNA methylation analysis and gene expression profiling. Whole genome DNA methylation using methyl-binding domain enrichment analysis revealed 60 differentially methylated regions (DMRs) between hay and corn fetuses with 39 DMRs more highly methylated in the hay fetuses vs. 21 DMRs more highly methylated in the corn fetuses. Three DMRs (LPAR3, PLIN5-PLIN4, and the differential methylation of a novel lincRNA) were validated using bisulfite sequencing. These DMRs were associated with differential gene expression. Additionally, significant DNA methylation differences were found at the single CpG level. Integrative methylome and transcriptome analysis revealed an association between gene expression and inter-/intragenic methylated regions. Furthermore, intragenic DMRs were found to be associated with expression of neighboring genes. CONCLUSIONS: The findings of this study imply that maternal diet from mid- to late-gestation can shape the epigenome and the transcriptome of fetal tissues, and putatively affect phenotypes of the lambs.

Transcriptional analysis of liver from chickens with fast (meat bird), moderate (F1 layer x meat bird cross) and low (layer bird) growth potential.

BACKGROUND: Divergent selection for meat and egg production in poultry has resulted in strains of birds differing widely in traits related to these products. Modern strains of meat birds can reach live weights of 2 kg in 35 d, while layer strains are now capable of producing more than 300 eggs per annum but grow slowly. In this study, RNA-Seq was used to investigate hepatic gene expression between three groups of birds with large differences in growth potential; meat bird, layer strain as well as an F1 layer x meat bird. The objective was to identify differentially expressed (DE) genes between all three strains to elucidate biological factors underpinning variations in growth performance. RESULTS: RNA-Seq analysis was carried out on total RNA extracted from the liver of meat bird (n = 6), F1 layer x meat bird cross (n = 6) and layer strain (n = 6), males. Differential expression of genes were considered significant at P < 0.05, and a false discovery rate of < 0.05, with any fold change considered. In total, 6278 genes were found to be DE with 5832 DE between meat birds and layers (19%), 2935 DE between meat birds and the cross (9.6%) and 493 DE between the cross and layers (1.6%). Comparisons between the three groups identified 155 significant DE genes. Gene ontology (GO) enrichment and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway analysis of the 155 DE genes showed the FoxO signalling pathway was most enriched (P = 0.001), including genes related to cell cycle regulation and insulin signalling. Significant GO terms included 'positive regulation of glucose import' and 'cellular response to oxidative stress', which is also consistent with FoxOs regulation of glucose metabolism. There were high correlations between FoxO pathway genes and bodyweight, as well as genes related to glycolysis and bodyweight. CONCLUSIONS: This study revealed large transcriptome differences between meat and layer birds. There was significant evidence implicating the FoxO signalling pathway (via cell cycle regulation and altered metabolism) as an active driver of growth variations in chicken. Functional analysis of the FoxO genes is required to understand how they regulate growth and egg production.

Distribution of ncRNAs expression across hypothalamic-pituitary-gonadal axis in Capra hircus.

BACKGROUND: Molecular regulation of the hypothalamic-pituitary-gonadal (HPG) axis plays an essential role in the fine tuning of seasonal estrus in Capra hircus. Noncoding RNAs (ncRNAs) are emerging as key regulators in sexual development and mammalian reproduction. In order to identify ncRNAs and to assess their expression patterns, along the HPG axis, we sequenced ncRNA libraries from hypothalamus, pituitary and ovary of three goats. RESULTS: Among the medium length noncoding RNAs (mncRNAs) identified, small nucleolar RNAs (snoRNAs) and transfer RNAs (tRNAs) were found to be more abundant in ovary and hypothalamus, respectively. The observed GC content was representative for different classes of ncRNAs, allowing the identification of a tRNA-derived RNA fragments (tRFs) subclass, which had a peak distribution around 32-38% GC content in the hypothalamus. Differences observed among organs confirmed the specificity of microRNA (miRNA) profiles for each organ system. CONCLUSIONS: Data on ncRNAs in organs constituting the HPG axis will contribute to understanding their role in the physiological regulation of reproduction in goats.

Genome-wide analysis of DNA methylation in hypothalamus and ovary of Capra hircus.

BACKGROUND: DNA methylation is a frequently studied epigenetic modification due to its role in regulating gene expression and hence in biological processes and in determining phenotypic plasticity in organisms. Rudimentary DNA methylation patterns for some livestock species are publically available: among these, goat methylome deserves to be further explored. RESULTS: Genome-wide DNA methylation maps of the hypothalamus and ovary from Saanen goats were generated using Methyl-CpG binding domain protein sequencing (MBD-seq). Analysis of DNA methylation patterns indicate that the majority of methylation peaks found within genes are located gene body regions, for both organs. Analysis of the distribution of methylated sites per chromosome showed that chromosome X had the lowest number of methylation peaks. The X chromosome has one of the highest percentages of methylated CpG islands in both organs, and approximately 50% of the CpG islands in the goat epigenome are methylated in hypothalamus and ovary. Organ-specific Differentially Methylated Genes (DMGs) were correlated with the expression levels. CONCLUSIONS: The comparison between transcriptome and methylome in hypothalamus and ovary showed that a higher level of methylation is not accompanied by a higher gene suppression. The genome-wide DNA methylation map for two goat organs produced here is a valuable starting point for studying the involvement of epigenetic modifications in regulating goat reproduction performance.