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.

MicroRNA breed and parent-of-origin effects provide insights into biological pathways differentiating cattle subspecies in fetal liver.

Introduction: MicroRNAs (miRNAs) play a crucial role in regulating gene expression during key developmental processes, including fetal development. Brahman (Bos taurus indicus) and Angus (Bos taurus taurus) cattle breeds represent two major cattle subspecies with strikingly different phenotypes. Methods: We analyzed miRNA expression in liver samples of purebred and reciprocal crosses of Angus and Brahman to investigate breed and parent-of-origin effects at the onset of accelerated fetal growth. Results: We identified eight novel miRNAs in fetal liver samples and 14 differentially expressed miRNAs (DEMs) between purebred samples. Correlation of gene expression modules and miRNAs by breed and parent-of-origin effects revealed an enrichment of genes associated with breed-specific differences in traits such as heat tolerance (Brahman) and fat deposition (Angus). We demonstrate that genes predicted to be targets of DEMs were more likely to be differentially expressed than non-targets (p-value < 0.05). We identified several miRNAs (bta-miR-187, bta-miR-216b, bta-miR-2284c, bta-miR-2285c, bta-miR-2285cp, bta-miR-2419-3p, bta-miR-2419-5p, and bta-miR-11984) that showed similar correlation patterns as bta-miR-2355-3p, which has been associated with the glutamatergic synapse pathway, a key facilitator of heat tolerance. Furthermore, we report Angus-breed-specific miRNAs (bta-miR-2313-5p, btamiR-490, bta-miR-2316, and bta-miR-11990) that may be involved in fat deposition. Finally, we showed that the DEMs identified in fetal liver are involved in Rap1, MAPK, and Ras signalling pathways, which are important for fetal development, muscle development and metabolic traits such as fat metabolism. Conclusion: Our work sheds light on the miRNA expression patterns that contribute to gene expression differences driving phenotypic differences in indicine and taurine cattle.

Asymmetric growth-limiting development of the female conceptus.

Introduction: Sex differences in prenatal growth may contribute to sex-dependent programming effects on postnatal phenotype. Methods: We integrated for the first time phenotypic, histomorphological, clinico-chemical, endocrine and gene expression analyses in a single species, the bovine conceptus at mid-gestation. Results: We demonstrate that by mid-gestation, before the onset of accelerated growth, the female conceptus displays asymmetric lower growth compared to males. Female fetuses were smaller with lower ponderal index and organ weights than males. However, their brain:body weight, brain:liver weight and heart:body weight ratios were higher than in males, indicating brain and heart 'sparing'. The female placenta weighed less and had lower volumes of trophoblast and fetal connective tissue than the male placenta. Female umbilical cord vessel diameters were smaller, and female-specific relationships of body weight and brain:liver weight ratios with cord vessel diameters indicated that the umbilico-placental vascular system creates a growth-limiting environment where blood flow is redistributed to protect brain and heart growth. Clinico-chemical indicators of liver perfusion support this female-specific growth-limiting phenotype, while lower insulin-like growth factor 2 (IGF2) gene expression in brain and heart, and lower circulating IGF2, implicate female-specific modulation of key endocrine mediators by nutrient supply. Conclusion: This mode of female development may increase resilience to environmental perturbations in utero and contribute to sex-bias in programming outcomes including susceptibility to non-communicable diseases.

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.

Evaluation of miRNA as Biomarkers of Emotional Valence in Pigs.

The ability to assess the welfare of animals is dependent on our ability to accurately determine their emotional (affective) state, with particular emphasis being placed on the identification of positive emotions. The challenge remains that current physiological and behavioral indices are either unable to distinguish between positive and negative emotional states, or they are simply not suitable for a production environment. Therefore, the development of novel measures of animal emotion is a necessity. Here we investigated the efficacy of microRNA (miRNA) in the brain and blood as biomarkers of emotional state in the pig. Female Large White × Landrace pigs (n = 24) were selected at weaning and trained to perform a judgment bias test (JBT), before being exposed for 5 weeks to either enriched (n = 12) or barren housing (n = 12) conditions. Pigs were tested on the JBT once prior to treatment, and immediately following treatment. MiRNA and neurotransmitters were analyzed in blood and brain tissue after euthanasia. Treatment had no effect on the outcomes of the JBT. There was also no effect of treatment on miRNA expression in blood or the brain (FDR p > 0.05). However, pigs exposed to enriched housing had elevated dopamine within the striatum compared to pigs in barren housing (p = 0.02). The results imply that either (a) miRNAs are not likely to be valid biomarkers of a positive affective state, at least under the type of conditions employed in this study, or (b) that the study design used to modify affective state was not able to create differential affective states, and therefore establish the validity of miRNA as biomarkers.

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.

Mitochondrial DNA Depletion in Granulosa Cell Derived Nuclear Transfer Tissues.

Somatic cell nuclear transfer (SCNT) is a key technology with broad applications that range from production of cloned farm animals to derivation of patient-matched stem cells or production of humanized animal organs for xenotransplantation. However, effects of aberrant epigenetic reprogramming on gene expression compromise cell and organ phenotype, resulting in low success rate of SCNT. Standard SCNT procedures include enucleation of recipient oocytes before the nuclear donor cell is introduced. Enucleation removes not only the spindle apparatus and chromosomes of the oocyte but also the perinuclear, mitochondria rich, ooplasm. Here, we use a Bos taurus SCNT model with in vitro fertilized (IVF) and in vivo conceived controls to demonstrate a ∼50% reduction in mitochondrial DNA (mtDNA) in the liver and skeletal muscle, but not the brain, of SCNT fetuses at day 80 of gestation. In the muscle, we also observed significantly reduced transcript abundances of mtDNA-encoded subunits of the respiratory chain. Importantly, mtDNA content and mtDNA transcript abundances correlate with hepatomegaly and muscle hypertrophy of SCNT fetuses. Expression of selected nuclear-encoded genes pivotal for mtDNA replication was similar to controls, arguing against an indirect epigenetic nuclear reprogramming effect on mtDNA amount. We conclude that mtDNA depletion is a major signature of perturbations after SCNT. We further propose that mitochondrial perturbation in interaction with incomplete nuclear reprogramming drives abnormal epigenetic features and correlated phenotypes, a concept supported by previously reported effects of mtDNA depletion on the epigenome and the pleiotropic phenotypic effects of mtDNA depletion in humans. This provides a novel perspective on the reprogramming process and opens new avenues to improve SCNT protocols for healthy embryo and tissue development.

Extended haplotype-phasing of long-read de novo genome assemblies using Hi-C.

Haplotype-resolved genome assemblies are important for understanding how combinations of variants impact phenotypes. To date, these assemblies have been best created with complex protocols, such as cultured cells that contain a single-haplotype (haploid) genome, single cells where haplotypes are separated, or co-sequencing of parental genomes in a trio-based approach. These approaches are impractical in most situations. To address this issue, we present FALCON-Phase, a phasing tool that uses ultra-long-range Hi-C chromatin interaction data to extend phase blocks of partially-phased diploid assembles to chromosome or scaffold scale. FALCON-Phase uses the inherent phasing information in Hi-C reads, skipping variant calling, and reduces the computational complexity of phasing. Our method is validated on three benchmark datasets generated as part of the Vertebrate Genomes Project (VGP), including human, cow, and zebra finch, for which high-quality, fully haplotype-resolved assemblies are available using the trio-based approach. FALCON-Phase is accurate without having parental data and performance is better in samples with higher heterozygosity. For cow and zebra finch the accuracy is 97% compared to 80-91% for human. FALCON-Phase is applicable to any draft assembly that contains long primary contigs and phased associate contigs.

Haplotype-resolved genomes provide insights into structural variation and gene content in Angus and Brahman cattle.

Inbred animals were historically chosen for genome analysis to circumvent assembly issues caused by haplotype variation but this resulted in a composite of the two genomes. Here we report a haplotype-aware scaffolding and polishing pipeline which was used to create haplotype-resolved, chromosome-level genome assemblies of Angus (taurine) and Brahman (indicine) cattle subspecies from contigs generated by the trio binning method. These assemblies reveal structural and copy number variants that differentiate the subspecies and that variant detection is sensitive to the specific reference genome chosen. Six genes with immune related functions have additional copies in the indicine compared with taurine lineage and an indicus-specific extra copy of fatty acid desaturase is under positive selection. The haplotyped genomes also enable transcripts to be phased to detect allele-specific expression. This work exemplifies the value of haplotype-resolved genomes to better explore evolutionary and functional variations.

Dual-locus DNA metabarcoding reveals southern hairy-nosed wombats (Lasiorhinus latifrons Owen) have a summer diet dominated by toxic invasive plants.

Habitat degradation and summer droughts severely restrict feeding options for the endangered southern hairy-nosed wombat (SHNW; Lasiorhinus latifrons). We reconstructed SHNW summer diets by DNA metabarcoding from feces. We initially validated rbcL and ndhJ diet reconstructions using autopsied and captive animals. Subsequent diet reconstructions of wild wombats broadly reflected vegetative ground cover, implying local rather than long-range foraging. Diets were all dominated by alien invasives. Chemical analysis of alien food revealed Carrichtera annua contains high levels of glucosinolates. Clinical examination (7 animals) and autopsy (12 animals) revealed that the most degraded site also contained most individuals showing signs of glucosinolate poisoning. We infer that dietary poisoning through the ingestion of alien invasives may have contributed to the recent population crashes in the region. In floristically diverse sites, individuals appear to be able to manage glucosinolate intake by avoidance or episodic feeding but this strategy is less tractable in the most degraded sites. We conclude that recovery of the most affected populations may require effective Carrichtera management and interim supplementary feeding. More generally, we argue that protection against population decline by poisoning in territorial herbivores requires knowledge of their diet and of those food plants containing toxic principles.

Pronounced Inter-Individual Variation in Plasma Cortisol Response to Fluoxetine Hydrochloride in the Pig.

Animal welfare assessment requires measures for positive affective state. Pharmacological agents that manipulate affective state can be used to evaluate novel biomarkers for affective state assessment. However, to validate that an agent has modified brain function, a reliable indicator is required. Circulating cortisol has been used as a reporter for effective delivery of the antidepressant selective serotonin reuptake inhibitor (SSRI) fluoxetine hydrochloride to the brain in humans and sheep. Here, we tested cortisol as a reporter for effective delivery of fluoxetine hydrochloride to the pig brain. We hypothesized that following administration of SSRI, innervation of the serotonergic reward pathway would result in activation of the hypothalamic-pituitary-adrenal (HPA) axis, leading to increased circulating cortisol levels. Large White-Landrace gilts received either a single intravenous dose of 100 mg fluoxetine hydrochloride suspended in 10 mL saline (n = 4), or 10 mL saline alone (n = 4). Blood samples were collected every 15 min for one hour prior to, and six hours post-treatment. The interaction of treatment x time on mean plasma cortisol levels between 15-165 min post-treatment was significant (p = 0.048) with highest cortisol concentrations of SSRI treated pigs versus controls (+ 98%) at 135 min post-treatment. However, individual cortisol profiles after SSRI treatment revealed high inter-individual variation in response. We conclude that, while combined data imply that plasma cortisol may be a readout for SSRI efficacy, inter-individual variation in SSRI response may preclude application of this approach in the pig. Given the current limited sample size, further research to confirm these findings is needed.

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.

De novo assembly of haplotype-resolved genomes with trio binning.

Complex allelic variation hampers the assembly of haplotype-resolved sequences from diploid genomes. We developed trio binning, an approach that simplifies haplotype assembly by resolving allelic variation before assembly. In contrast with prior approaches, the effectiveness of our method improved with increasing heterozygosity. Trio binning uses short reads from two parental genomes to first partition long reads from an offspring into haplotype-specific sets. Each haplotype is then assembled independently, resulting in a complete diploid reconstruction. We used trio binning to recover both haplotypes of a diploid human genome and identified complex structural variants missed by alternative approaches. We sequenced an F1 cross between the cattle subspecies Bos taurus taurus and Bos taurus indicus and completely assembled both parental haplotypes with NG50 haplotig sizes of >20 Mb and 99.998% accuracy, surpassing the quality of current cattle reference genomes. We suggest that trio binning improves diploid genome assembly and will facilitate new studies of haplotype variation and inheritance.

Ultrasonographic measurements in first trimester concepti identify predictors of birth weight and postnatal development in cattle.

The placenta is a major driver of prenatal growth and involved in programming of postnatal performance. We therefore determined placental and embryo-fetal ultrasonographic parameters in early pregnancy and their relationships with birth weight and postnatal weights in a Bos indicus-Bos taurus composite beef cattle population. Pregnancies were generated in 2-yr-old Droughtmaster heifers by artificial insemination after estrus synchronization in 2 consecutive years (2009, n = 36 and 2010, n = 57), with a subset of 2010 heifers used again as lactating 3-yr-old cows in 2011 (n = 24). Each cohort was managed as 1 contemporary group for measurements of Corpus luteum diameter, amnion length and width, placentome width and thickness, and embryo-fetal crown-rump length, at 7 and 8 wk of gestation. This was followed by recordings of birth weight, branding weight at 5 to 6 mo of age and weaning weight 2 mo later. At a significance threshold of P < 0.05, placentome thickness at week 7 was negatively correlated with weights at birth (r = -0.23), branding (r = -0.25), and weaning (r = -0.35), whereas placentome width at week 7 (r = 0.24) and thickness at week 8 (r = 0.29) were positively correlated with birth weight. Thicker placentomes in males at week 7 (7%) difference mirrored sex differences in weights at birth (7%), branding (10%), and weaning (6%). The sex difference trend for birth weight was not consistent across sire-year combinations, ranging from -3.2 to +4.7 kg (birth weight of males - females per sire). These results support the hypothesis that placental parameters at the transition from embryo to fetal stage are major predictors of fetal and postnatal growth, albeit with significant environmentally induced plasticity, in stabilized B. indicus-B. taurus composite populations, and suggest that elements of B. indicus-B. taurus reciprocal differences in birth weight persist in composite populations.

Atlas of tissue- and developmental stage specific gene expression for the bovine insulin-like growth factor (IGF) system.

The insulin-like growth factor (IGF) axis is fundamental for mammalian growth and development. However, no comprehensive reference data on gene expression across tissues and pre- and postnatal developmental stages are available for any given species. Here we provide systematic promoter- and splice variant specific information on expression of IGF system components in embryonic (Day 48), fetal (Day 153), term (Day 277, placenta) and juvenile (Day 365-396) tissues of domestic cow, a major agricultural species and biomedical model. Analysis of spatiotemporal changes in expression of IGF1, IGF2, IGF1R, IGF2R, IGFBP1-8 and IR genes, as well as lncRNAs H19 and AIRN, by qPCR, indicated an overall increase in expression from embryo to fetal stage, and decrease in expression from fetal to juvenile stage. The stronger decrease in expression of lncRNAs (average -16-fold) and ligands (average -12.1-fold) compared to receptors (average -5.7-fold) and binding proteins (average -4.3-fold) is consistent with known functions of IGF peptides and supports important roles of lncRNAs in prenatal development. Pronounced overall reduction in postnatal expression of IGF system components in lung (-12.9-fold) and kidney (-13.2-fold) are signatures of major changes in organ function while more similar hepatic expression levels (-2.2-fold) are evidence of the endocrine rather than autocrine/paracrine role of IGFs in postnatal growth regulation. Despite its rapid growth, placenta displayed a more stable expression pattern than other organs during prenatal development. Quantitative analyses of contributions of promoters P0-P4 to global IGF2 transcript in fetal tissues revealed that P4 accounted for the bulk of transcript in all tissues but skeletal muscle. Demonstration of IGF2 expression in fetal muscle and postnatal liver from a promoter orthologous to mouse and human promoter P0 provides further evidence for an evolutionary and developmental shift from placenta-specific P0-expression in rodents and suggests that some aspects of bovine IGF expression may be closer to human than mouse.

Remodelling of the bovine placenta: Comprehensive morphological and histomorphological characterization at the late embryonic and early accelerated fetal growth stages.

INTRODUCTION: Placental function impacts growth and development with lifelong consequences for performance and health. We provide novel insights into placental development in bovine, an important agricultural species and biomedical model. METHODS: Concepti with defined genetics and sex were recovered from nulliparous dams managed under standardized conditions to study placental gross morphological and histomorphological parameters at the late embryo (Day48) and early accelerated fetal growth (Day153) stages. RESULTS: Placentome number increased 3-fold between Day48 and Day153. Placental barrier thickness was thinner, and volume of placental components, and surface areas and densities were higher at Day153 than Day48. We confirmed two placentome types, flat and convex. At Day48, there were more convex than flat placentomes, and convex placentomes had a lower proportion of maternal connective tissue (P < 0.01). However, this was reversed at Day153, where convex placentomes were lower in number and had greater volume of placental components (P < 0.01- P < 0.001) and greater surface area (P < 0.001) than flat placentomes. Importantly, embryo (r = 0.50) and fetal (r = 0.30) weight correlated with total number of convex but not flat placentomes. DISCUSSION: Extensive remodelling of the placenta increases capacity for nutrient exchange to support rapidly increasing embryo-fetal weight from Day48 to Day153. The cellular composition of convex placentomes, and exclusive relationships between convex placentome number and embryo-fetal weight, provide strong evidence for these placentomes as drivers of prenatal growth. The difference in proportion of maternal connective tissue between placentome types at Day48 suggests that this tissue plays a role in determining placentome shape, further highlighting the importance of early placental development.

Widespread differential maternal and paternal genome effects on fetal bone phenotype at mid-gestation.

Parent-of-origin-dependent (epi)genetic factors are important determinants of prenatal development that program adult phenotype. However, data on magnitude and specificity of maternal and paternal genome effects on fetal bone are lacking. We used an outbred bovine model to dissect and quantify effects of parental genomes, fetal sex, and nongenetic maternal effects on the fetal skeleton and analyzed phenotypic and molecular relationships between fetal muscle and bone. Analysis of 51 bone morphometric and weight parameters from 72 fetuses recovered at day 153 gestation (54% term) identified six principal components (PC1-6) that explained 80% of the variation in skeletal parameters. Parental genomes accounted for most of the variation in bone wet weight (PC1, 72.1%), limb ossification (PC2, 99.8%), flat bone size (PC4, 99.7%), and axial skeletal growth (PC5, 96.9%). Limb length showed lesser effects of parental genomes (PC3, 40.8%) and a significant nongenetic maternal effect (gestational weight gain, 29%). Fetal sex affected bone wet weight (PC1, p < 0.0001) and limb length (PC3, p < 0.05). Partitioning of variation explained by parental genomes revealed strong maternal genome effects on bone wet weight (74.1%, p < 0.0001) and axial skeletal growth (93.5%, p < 0.001), whereas paternal genome controlled limb ossification (95.1%, p < 0.0001). Histomorphometric data revealed strong maternal genome effects on growth plate height (98.6%, p < 0.0001) and trabecular thickness (85.5%, p < 0.0001) in distal femur. Parental genome effects on fetal bone were mirrored by maternal genome effects on fetal serum 25-hydroxyvitamin D (96.9%, p < 0.001) and paternal genome effects on alkaline phosphatase (90.0%, p < 0.001) and their correlations with maternally controlled bone wet weight and paternally controlled limb ossification, respectively. Bone wet weight and flat bone size correlated positively with muscle weight (r = 0.84 and 0.77, p < 0.0001) and negatively with muscle H19 expression (r = -0.34 and -0.31, p < 0.01). Because imprinted maternally expressed H19 regulates growth factors by miRNA interference, this suggests muscle-bone interaction via epigenetic factors.

Tissue-specific and minor inter-individual variation in imprinting of IGF2R is a common feature of Bos taurus Concepti and not correlated with fetal weight.

The insulin-like growth factor 2 receptor (IGF2R) is essential for prenatal growth regulation and shows gene dosage effects on fetal weight that can be affected by in-vitro embryo culture. Imprinted maternal expression of murine Igf2r is well documented for all fetal tissues excluding brain, but polymorphic imprinting and biallelic expression were reported for IGF2R in human. These differences have been attributed to evolutionary changes correlated with specific reproductive strategies. However, data from species suitable for testing this hypothesis are lacking. The domestic cow (Bos taurus) carries a single conceptus with a similar gestation length as human. We identified 12 heterozygous concepti informative for imprinting studies among 68 Bos taurus fetuses at Day 80 of gestation (28% term) and found predominantly maternal IGF2R expression in all fetal tissues but brain, which escapes imprinting. Inter-individual variation in allelic expression bias, i.e. expression of the repressed paternal allele relative to the maternal allele, ranged from 4.6-8.9% in heart, 4.3-10.2% in kidney, 6.1-11.2% in liver, 4.6-15.8% in lung and 3.2-12.2% in skeletal muscle. Allelic bias for mesodermal tissues (heart, skeletal muscle) differed significantly (P<0.05) from endodermal tissues (liver, lung). The placenta showed partial imprinting with allelic bias of 22.9-34.7% and differed significantly (P<0.001) from all other tissues. Four informative fetuses were generated by in-vitro fertilization (IVF) with embryo culture and two individuals displayed fetal overgrowth. However, there was no evidence for changes in imprinting or DNA methylation after IVF, or correlations between allelic bias and fetal weight. In conclusion, imprinting of Bos taurus IGF2R is similar to mouse except in placenta, which could indicate an effect of reproductive strategy. Common minor inter-individual variation in allelic bias and absence of imprinting abnormalities in IVF fetuses suggest changes in IGF2R expression in overgrown fetuses could be modulated through other mechanisms than changes in imprinting.

Maternal and paternal genomes differentially affect myofibre characteristics and muscle weights of bovine fetuses at midgestation.

Postnatal myofibre characteristics and muscle mass are largely determined during fetal development and may be significantly affected by epigenetic parent-of-origin effects. However, data on such effects in prenatal muscle development that could help understand unexplained variation in postnatal muscle traits are lacking. In a bovine model we studied effects of distinct maternal and paternal genomes, fetal sex, and non-genetic maternal effects on fetal myofibre characteristics and muscle mass. Data from 73 fetuses (Day153, 54% term) of four genetic groups with purebred and reciprocal cross Angus and Brahman genetics were analyzed using general linear models. Parental genomes explained the greatest proportion of variation in myofibre size of Musculus semitendinosus (80-96%) and in absolute and relative weights of M. supraspinatus, M. longissimus dorsi, M. quadriceps femoris and M. semimembranosus (82-89% and 56-93%, respectively). Paternal genome in interaction with maternal genome (P<0.05) explained most genetic variation in cross sectional area (CSA) of fast myotubes (68%), while maternal genome alone explained most genetic variation in CSA of fast myofibres (93%, P<0.01). Furthermore, maternal genome independently (M. semimembranosus, 88%, P<0.0001) or in combination (M. supraspinatus, 82%; M. longissimus dorsi, 93%; M. quadriceps femoris, 86%) with nested maternal weight effect (5-6%, P<0.05), was the predominant source of variation for absolute muscle weights. Effects of paternal genome on muscle mass decreased from thoracic to pelvic limb and accounted for all (M. supraspinatus, 97%, P<0.0001) or most (M. longissimus dorsi, 69%, P<0.0001; M. quadriceps femoris, 54%, P<0.001) genetic variation in relative weights. An interaction between maternal and paternal genomes (P<0.01) and effects of maternal weight (P<0.05) on expression of H19, a master regulator of an imprinted gene network, and negative correlations between H19 expression and fetal muscle mass (P<0.001), suggested imprinted genes and miRNA interference as mechanisms for differential effects of maternal and paternal genomes on fetal muscle.