Single-Step Genome Editing of Small Ruminant Embryos by Electroporation.

We investigated the possibility of single-step genome editing in small ruminants by CRISPR-Cas9 zygote electroporation. We targeted SOCS2 and PDX1 in sheep embryos and OTX2 in goat embryos, utilizing a dual sgRNA approach. Gene editing efficiency was compared between microinjection and three different electroporation settings performed at four different times of embryo development. Electroporation of sheep zygotes 6 h after fertilization with settings that included short high-voltage (poring) and long low-voltage (transfer) pulses was efficient at producing SOCS2 knock-out blastocysts. The mutation rate after CRISPR/Cas9 electroporation was 95.6% ± 8%, including 95.4% ± 9% biallelic mutations; which compared favorably to 82.3% ± 8% and 25% ± 10%, respectively, when using microinjection. We also successfully disrupted the PDX1 gene in sheep and the OTX2 gene in goat embryos. The biallelic mutation rate was 81 ± 5% for PDX1 and 85% ± 6% for OTX2. In conclusion, using single-step CRISPR-Cas9 zygote electroporation, we successfully introduced biallelic deletions in the genome of small ruminant embryos.

Effect of calving season on the parameters and components of the lactation curve in Holstein dairy cows managed in a semi-desert climate.

The lactation curve in dairy cows is influenced by the calving season, which is highly dependent on the warm climate in semi-arid regions. Objective herein was to evaluate effects of calving season on the parameters and components of the lactation curve in Holstein cows. The study included 278,317 milk records collected from 1086 cows from the 1st to 5th calving and good body condition score. The cows were grouped according to the season in which they calved: winter-calved (CS1), summer-calved (CS2), and autumn-calved cows (CS3). Ambient temperature and humidity data were used to calculate the temperature-humidity index (THI). The NLIN procedure was used to estimate the parameters of the lactation curve that served to calculate the components. The mixed procedure was executed to analyze the fixed effect of calving season. Associations between lactation curve traits were tested using correlation and regression analyses. A univariate model was utilized to calculate heritability. Average THI values during the lactation period were 73.5, 68.5, and 69.5 units for CS1, CS2, and CS3 groups, respectively. Initial milk production and increasing rate to the maximum milk yield in CS1 and CS3 groups were higher (P < 0.05) than CS2 cows. However, persistency and total milk yield during the entire lactation period were superior (P < 0.05) for CS2 and CS3 cows compared to CS1 cows, probably due to the moderate heat stress during the lactation period in the CS1 group. In cows from CS2 and CS3 groups, total milk production at 305 days was moderately correlated with initial milk production (r = 0.47; P < 0.05), and highly correlated with milk yield at peak day (r = 0.91; P < 0.05) which resulted as reliable predictor for total milk yield during the entire lactation (R2 = 0.83). In conclusion, the THI prevailing during the different calving seasons appeared to be an important factor influencing the performance of the lactation curve.

A SNP within the PMCH gene as a molecular marker associated with fertility traits in Angus and Brangus beef heifers raised under a desert environment.

Components of the GH/IGF1 endocrine axis regulate growth and reproductive traits in cattle. The pro-melanin-concentrating hormone (PMCH) gene located within chromosome 5 belongs to this axis. Objective herein was to evaluate PMCH single-nucleotide polymorphisms (SNPs) as molecular markers associated with age at first calving, calving interval, and age at second calving in Angus and Brangus beef heifers raised in desert conditions. Five SNPs within the PMCH gene were included in the study. Three of these SNPs had minor allele frequency > 10% and only one SNP did not deviate from Hardy-Weinberg equilibrium. A genotype to phenotype association analyses was performed using a mixed-effects model which included phenotype as the response variable, SNP genotype, breed, year of birth and age of dam as fixed terms, and sire as a random effect. Genotypes from the SNP rs135033882 were found to be associated (P < 0.05) with all evaluated fertility traits, and the term breed resulted as a significant source of variation only for age at second calving. The allele A was the favorable allele because it decreased the age at first calving 98.6 days, the calving interval 85.3 days, and the age at second calving 183.1 days, in Angus and Brangus heifers. In conclusion, we proposed a SNP within the PMCH gene as a potential candidate marker associated with reproductive performance in Angus and Brangus beef heifers raised in a desert climate.

A comparative analysis of chromatin accessibility in cattle, pig, and mouse tissues.

BACKGROUND: Although considerable progress has been made towards annotating the noncoding portion of the human and mouse genomes, regulatory elements in other species, such as livestock, remain poorly characterized. This lack of functional annotation poses a substantial roadblock to agricultural research and diminishes the value of these species as model organisms. As active regulatory elements are typically characterized by chromatin accessibility, we implemented the Assay for Transposase Accessible Chromatin (ATAC-seq) to annotate and characterize regulatory elements in pigs and cattle, given a set of eight adult tissues. RESULTS: Overall, 306,304 and 273,594 active regulatory elements were identified in pig and cattle, respectively. 71,478 porcine and 47,454 bovine regulatory elements were highly tissue-specific and were correspondingly enriched for binding motifs of known tissue-specific transcription factors. However, in every tissue the most prevalent accessible motif corresponded to the insulator CTCF, suggesting pervasive involvement in 3-D chromatin organization. Taking advantage of a similar dataset in mouse, open chromatin in pig, cattle, and mice were compared, revealing that the conservation of regulatory elements, in terms of sequence identity and accessibility, was consistent with evolutionary distance; whereas pig and cattle shared about 20% of accessible sites, mice and ungulates only had about 10% of accessible sites in common. Furthermore, conservation of accessibility was more prevalent at promoters than at intergenic regions. CONCLUSIONS: The lack of conserved accessibility at distal elements is consistent with rapid evolution of enhancers, and further emphasizes the need to annotate regulatory elements in individual species, rather than inferring elements based on homology. This atlas of chromatin accessibility in cattle and pig constitutes a substantial step towards annotating livestock genomes and dissecting the regulatory link between genome and phenome.

In vitro breeding: application of embryonic stem cells to animal production†.

Embryonic stem cells (ESCs) are derived from the inner cell mass of preimplantation blastocysts. For decades, attempts to efficiently derive ESCs in animal livestock species have been unsuccessful, but this goal has recently been achieved in cattle. Together with the recent reconstitution of the germ cell differentiation processes from ESCs in mice, these achievements open new avenues for the development of promising technologies oriented toward improving health, animal production, and the environment. In this article, we present a strategy that will notably accelerate genetic improvement in livestock populations by reducing the generational interval, namely in vitro breeding (IVB). IVB combines genomic selection, a widely used strategy for genetically improving livestock, with ESC derivation and in vitro differentiation of germ cells from pluripotent stem cells. We also review the most recent findings in the fields on which IVB is based. Evidence suggests this strategy will be soon within reach.

Neuropeptidome of the Hypothalamus and Pituitary Gland of Indicine × Taurine Heifers: Evidence of Differential Neuropeptide Processing in the Pituitary Gland before and after Puberty.

Puberty in cattle is regulated by an endocrine axis, which includes a complex milieu of neuropeptides in the hypothalamus and pituitary gland. The neuropeptidome of hypothalamic-pituitary gland tissue of pre- (PRE) and postpubertal (POST) Bos indicus-influenced heifers was characterized, followed by quantitative analysis of 51 fertility-related neuropeptides in these tissues. Comparison of peptide abundances with gene expression levels allowed assessment of post-transcriptional peptide processing. On the basis of classical cleavage, 124 mature neuropeptides from 35 precursor proteins were detected in hypothalamus and pituitary gland tissues of three PRE and three POST Brangus heifers. An additional 19 peptides (cerebellins, PEN peptides) previously reported as neuropeptides that did not follow classical cleavage were also identified. In the pre-pubertal hypothalamus, a greater diversity of neuropeptides (25.8%) was identified relative to post-pubertal heifers, while in the pituitary gland, 38.6% more neuropeptides were detected in the post-pubertal heifers. Neuro-tissues of PRE and POST heifers revealed abundance differences ( p < 0.05) in peptides from protein precursors involved in packaging and processing (e.g., the granin family and ProSAAS) or neuron stimulation (PENK, CART, POMC, cerebellins). On their own, the transcriptome data of the precursor genes could not predict the neuropeptide profile in the exact same tissues in several cases. This provides further evidence of the importance of differential processing of the neuropeptide precursors in the pituitary before and after puberty.

Genome-wide identification of tissue-specific long non-coding RNA in three farm animal species.

BACKGROUND: Numerous long non-coding RNAs (lncRNAs) have been identified and their roles in gene regulation in humans, mice, and other model organisms studied; however, far less research has been focused on lncRNAs in farm animal species. While previous studies in chickens, cattle, and pigs identified lncRNAs in specific developmental stages or differentially expressed under specific conditions in a limited number of tissues, more comprehensive identification of lncRNAs in these species is needed. The goal of the FAANG Consortium (Functional Annotation of Animal Genomes) is to functionally annotate animal genomes, including the annotation of lncRNAs. As one of the FAANG pilot projects, lncRNAs were identified across eight tissues in two adult male biological replicates from chickens, cattle, and pigs. RESULTS: Comprehensive lncRNA annotations for the chicken, cattle, and pig genomes were generated by utilizing RNA-seq from eight tissue types from two biological replicates per species at the adult developmental stage. A total of 9393 lncRNAs in chickens, 7235 lncRNAs in cattle, and 14,429 lncRNAs in pigs were identified. Including novel isoforms and lncRNAs from novel loci, 5288 novel lncRNAs were identified in chickens, 3732 in cattle, and 4870 in pigs. These transcripts match previously known patterns of lncRNAs, such as generally lower expression levels than mRNAs and higher tissue specificity. An analysis of lncRNA conservation across species identified a set of conserved lncRNAs with potential functions associated with chromatin structure and gene regulation. Tissue-specific lncRNAs were identified. Genes proximal to tissue-specific lncRNAs were enriched for GO terms associated with the tissue of origin, such as leukocyte activation in spleen. CONCLUSIONS: LncRNAs were identified in three important farm animal species using eight tissues from adult individuals. About half of the identified lncRNAs were not previously reported in the NCBI annotations for these species. While lncRNAs are less conserved than protein-coding genes, a set of positionally conserved lncRNAs were identified among chickens, cattle, and pigs with potential functions related to chromatin structure and gene regulation. Tissue-specific lncRNAs have potential regulatory functions on genes enriched for tissue-specific GO terms. Future work will include epigenetic data from ChIP-seq experiments to further refine these annotations.

Polymorphisms within the prolactin and growth hormone/insulin-like growth factor-1 functional pathways associated with fertility traits in Holstein cows raised in a hot-humid climate.

Prolactin (PRL), growth hormone (GH), and insulin-like growth factor-1 (IGF-1) are in hormone-response pathways involved in energy metabolism during thermoregulation processes in cattle. Objective herein was to study the association between single nucleotide polymorphisms (SNP) within genes of the PRL and GH/IGF-1 pathways with fertility traits such as services per conception (SPC) and days open (DO) in Holstein cattle lactating under a hot-humid climate. Ambient temperature and relative humidity were used to calculate the temperature-humidity index (THI) which revealed that the cows were exposed to heat stress conditions from June to November of 2012 in southern Sonora, Mexico. Individual blood samples from all cows were collected, spotted on FTA cards, and used to genotype a 179 tag SNP panel within 44 genes from the PRL and GH/IGF-1 pathways. The associative analyses among SNP genotypes and fertility traits were performed using mixed-effect models. Allele substitution effects were calculated using a regression model that included the genotype term as covariate. Single-SNP association analyses indicated that eight SNP within the genes IGF-1, IGF-1R, IGFBP5, PAPPA1, PMCH, PRLR, SOCS5, and SSTR2 were associated with SPC (P < 0.05), whereas four SNP in the genes GHR, PAPPA2, PRLR, and SOCS4 were associated with DO (P < 0.05). In conclusion, SNP within genes of the PRL and GH/IGF-1 pathways resulted as predictors of reproductive phenotypes in heat-stressed Holstein cows, and these SNP are proposed as candidates for a marker-assisted selection program intended to improve fertility of dairy cattle raised in warm climates.

Metabolic Reprogramming Regulates the Proliferative and Inflammatory Phenotype of Adventitial Fibroblasts in Pulmonary Hypertension Through the Transcriptional Corepressor C-Terminal Binding Protein-1.

BACKGROUND: Changes in metabolism have been suggested to contribute to the aberrant phenotype of vascular wall cells, including fibroblasts, in pulmonary hypertension (PH). Here, we test the hypothesis that metabolic reprogramming to aerobic glycolysis is a critical adaptation of fibroblasts in the hypertensive vessel wall that drives proliferative and proinflammatory activation through a mechanism involving increased activity of the NADH-sensitive transcriptional corepressor C-terminal binding protein 1 (CtBP1). METHODS: RNA sequencing, quantitative polymerase chain reaction,13C-nuclear magnetic resonance, fluorescence-lifetime imaging, mass spectrometry-based metabolomics, and tracing experiments with U-13C-glucose were used to assess glycolytic reprogramming and to measure the NADH/NAD+ ratio in bovine and human adventitial fibroblasts and mouse lung tissues. Immunohistochemistry was used to assess CtBP1 expression in the whole-lung tissues. CtBP1 siRNA and the pharmacological inhibitor 4-methylthio-2-oxobutyric acid (MTOB) were used to abrogate CtBP1 activity in cells and hypoxic mice. RESULTS: We found that adventitial fibroblasts from calves with severe hypoxia-induced PH and humans with idiopathic pulmonary arterial hypertension (PH-Fibs) displayed aerobic glycolysis when cultured under normoxia, accompanied by increased free NADH and NADH/NAD+ ratios. Expression of the NADH sensor CtBP1 was increased in vivo and in vitro in fibroblasts within the pulmonary adventitia of humans with idiopathic pulmonary arterial hypertension and animals with PH and cultured PH-Fibs, respectively. Decreasing NADH pharmacologically with MTOB or genetically blocking CtBP1 with siRNA upregulated the cyclin-dependent genes (p15 and p21) and proapoptotic regulators (NOXA and PERP), attenuated proliferation, corrected the glycolytic reprogramming phenotype of PH-Fibs, and augmented transcription of the anti-inflammatory gene HMOX1. Chromatin immunoprecipitation analysis demonstrated that CtBP1 directly binds the HMOX1 promoter. Treatment of hypoxic mice with MTOB decreased glycolysis and expression of inflammatory genes, attenuated proliferation, and suppressed macrophage numbers and remodeling in the distal pulmonary vasculature. CONCLUSIONS: CtBP1 is a critical factor linking changes in cell metabolism to cell phenotype in hypoxic and other forms of PH and a therapeutic target.

The evolution of the natural killer complex; a comparison between mammals using new high-quality genome assemblies and targeted annotation.

Natural killer (NK) cells are a diverse population of lymphocytes with a range of biological roles including essential immune functions. NK cell diversity is in part created by the differential expression of cell surface receptors which modulate activation and function, including multiple subfamilies of C-type lectin receptors encoded within the NK complex (NKC). Little is known about the gene content of the NKC beyond rodent and primate lineages, other than it appears to be extremely variable between mammalian groups. We compared the NKC structure between mammalian species using new high-quality draft genome assemblies for cattle and goat; re-annotated sheep, pig, and horse genome assemblies; and the published human, rat, and mouse lemur NKC. The major NKC genes are largely in the equivalent positions in all eight species, with significant independent expansions and deletions between species, allowing us to propose a model for NKC evolution during mammalian radiation. The ruminant species, cattle and goats, have independently evolved a second KLRC locus flanked by KLRA and KLRJ, and a novel KLRH-like gene has acquired an activating tail. This novel gene has duplicated several times within cattle, while other activating receptor genes have been selectively disrupted. Targeted genome enrichment in cattle identified varying levels of allelic polymorphism between the NKC genes concentrated in the predicted extracellular ligand-binding domains. This novel recombination and allelic polymorphism is consistent with NKC evolution under balancing selection, suggesting that this diversity influences individual immune responses and may impact on differential outcomes of pathogen infection and vaccination.

Parallel responses of human epidermal keratinocytes to inorganic SbIII and AsIII.

SbIII and AsIII are known to exhibit similar chemical properties, but the degree of similarity in their effects on biological systems merits further exploration. Present work compares the responses of human epidermal keratinocytes, a known target cell type for arsenite-induced carcinogenicity, to these metalloids after treatment for a week at environmentally relevant concentrations. Previous work with these cells has shown that arsenite and antimonite have parallel effects in suppressing differentiation, altering levels of several critical enzymes and maintaining colony forming ability. More globally, protein profiling now reveals parallels in SbIII and AsIII effects. The more sensitive technique of transcriptional profiling also shows considerable parallels. Thus, gene expression changes were almost entirely in the same directions for the two treatments, although the degree of change was sometimes significantly different. Inspection of the changes revealed that RYR1 and LRIG1 were among the genes strongly suppressed, consistent with reduced calcium-dependent differentiation and maintenance of EGF-dependent proliferative potential. Moreover, levels of miRNAs in the cells were altered in parallel, with nearly 90% of the 198 most highly expressed ones being suppressed. Among these was miR-203, which is known to decrease proliferative potential. Finally, both SbIII and AsIII were seen to attenuate bone morphogenetic protein 6 induction of dual specificity phosphatases 2 and 14, consistent with maintaining epidermal growth factor receptor signaling. These findings raise the question whether SbIII, like AsIII, could act as a human skin carcinogen.

Overlapping mouse subcongenic strains successfully separate two linked body fat QTL on distal MMU 2.

BACKGROUND: Mouse chromosome 2 is linked to growth and body fat phenotypes in many mouse crosses. With the goal to identify the underlying genes regulating growth and body fat on mouse chromosome 2, we developed five overlapping subcongenic strains that contained CAST/EiJ donor regions in a C57BL/6J (hg/hg) background (hg is a spontaneous deletion of 500 Kb on mouse chromosome 10). To fine map QTL on distal mouse chromosome 2 a total of 1,712 F2 mice from the five subcongenic strains, plus 278 F2 mice from the HG2D founder congenic strain were phenotyped and analyzed. Interval mapping (IM) and composite IM (CIM) were performed on body weight and body fat traits on a combination of SNP and microsatellite markers, which generated a high-density genotyping panel. RESULTS: Phenotypic analysis and interval mapping of total fat mass identified two QTL on distal mouse chromosome 2. One QTL between 150 and 161 Mb, Fatq2a, and the second between 173.3 and 175.6 Mb, Fatq2b. The two QTL reside in different congenic strains with significant total fat differences between homozygous cast/cast and b6/b6 littermates. Both of these QTL were previously identified only as a single QTL affecting body fat, Fatq2. Furthermore, through a novel approach referred here as replicated CIM, Fatq2b was mapped to the Gnas imprinted locus. CONCLUSIONS: The integration of subcongenic strains, high-density genotyping, and CIM succesfully partitioned two previously linked QTL 20 Mb apart, and the strongest QTL, Fatq2b, was fine mapped to a ~2.3 Mb region interval encompassing the Gnas imprinted locus.

What's in your next-generation sequence data? An exploration of unmapped DNA and RNA sequence reads from the bovine reference individual.

BACKGROUND: Next-generation sequencing projects commonly commence by aligning reads to a reference genome assembly. While improvements in alignment algorithms and computational hardware have greatly enhanced the efficiency and accuracy of alignments, a significant percentage of reads often remain unmapped. RESULTS: We generated de novo assemblies of unmapped reads from the DNA and RNA sequencing of the Bos taurus reference individual and identified the closest matching sequence to each contig by alignment to the NCBI non-redundant nucleotide database using BLAST. As expected, many of these contigs represent vertebrate sequence that is absent, incomplete, or misassembled in the UMD3.1 reference assembly. However, numerous additional contigs represent invertebrate species. Most prominent were several species of Spirurid nematodes and a blood-borne parasite, Babesia bigemina. These species are either not present in the US or are not known to infect taurine cattle and the reference animal appears to have been host to unsequenced sister species. CONCLUSIONS: We demonstrate the importance of exploring unmapped reads to ascertain sequences that are either absent or misassembled in the reference assembly and for detecting sequences indicative of parasitic or commensal organisms.

A clone-free, single molecule map of the domestic cow (Bos taurus) genome.

BACKGROUND: The cattle (Bos taurus) genome was originally selected for sequencing due to its economic importance and unique biology as a model organism for understanding other ruminants, or mammals. Currently, there are two cattle genome sequence assemblies (UMD3.1 and Btau4.6) from groups using dissimilar assembly algorithms, which were complemented by genetic and physical map resources. However, past comparisons between these assemblies revealed substantial differences. Consequently, such discordances have engendered ambiguities when using reference sequence data, impacting genomic studies in cattle and motivating construction of a new optical map resource--BtOM1.0--to guide comparisons and improvements to the current sequence builds. Accordingly, our comprehensive comparisons of BtOM1.0 against the UMD3.1 and Btau4.6 sequence builds tabulate large-to-immediate scale discordances requiring mediation. RESULTS: The optical map, BtOM1.0, spanning the B. taurus genome (Hereford breed, L1 Dominette 01449) was assembled from an optical map dataset consisting of 2,973,315 (439 X; raw dataset size before assembly) single molecule optical maps (Rmaps; 1 Rmap = 1 restriction mapped DNA molecule) generated by the Optical Mapping System. The BamHI map spans 2,575.30 Mb and comprises 78 optical contigs assembled by a combination of iterative (using the reference sequence: UMD3.1) and de novo assembly techniques. BtOM1.0 is a high-resolution physical map featuring an average restriction fragment size of 8.91 Kb. Comparisons of BtOM1.0 vs. UMD3.1, or Btau4.6, revealed that Btau4.6 presented far more discordances (7,463) vs. UMD3.1 (4,754). Overall, we found that Btau4.6 presented almost double the number of discordances than UMD3.1 across most of the 6 categories of sequence vs. map discrepancies, which are: COMPLEX (misassembly), DELs (extraneous sequences), INSs (missing sequences), ITs (Inverted/Translocated sequences), ECs (extra restriction cuts) and MCs (missing restriction cuts). CONCLUSION: Alignments of UMD3.1 and Btau4.6 to BtOM1.0 reveal discordances commensurate with previous reports, and affirm the NCBI's current designation of UMD3.1 sequence assembly as the "reference assembly" and the Btau4.6 as the "alternate assembly." The cattle genome optical map, BtOM1.0, when used as a comprehensive and largely independent guide, will greatly assist improvements to existing sequence builds, and later serve as an accurate physical scaffold for studies concerning the comparative genomics of cattle breeds.

Single nucleotide polymorphism (SNP) variation of wolves (Canis lupus) in Southeast Alaska and comparison with wolves, dogs, and coyotes in North America.

There is considerable interest in the genetics of wolves (Canis lupus) because of their close relationship to domestic dogs (C. familiaris) and the need for informed conservation and management. This includes wolf populations in Southeast Alaska for which we determined genotypes of 305 wolves at 173662 single nucleotide polymorphism (SNP) loci. After removal of invariant and linked SNP, 123801 SNP were used to quantify genetic differentiation of wolves in Southeast Alaska and wolves, coyotes (C. latrans), and dogs from other areas in North America. There is differentiation of SNP allele frequencies between the species (wolves, coyotes, and dogs), although differentiation is relatively low between some wolf and coyote populations. There are varying levels of differentiation among populations of wolves, including low differentiation of wolves in interior Alaska, British Columbia, and the northern US Rocky Mountains. There is considerable differentiation of SNP allele frequencies of wolves in Southeast Alaska from wolves in other areas. However, wolves in Southeast Alaska are not a genetically homogeneous group and there are comparable levels of genetic differentiation among areas within Southeast Alaska and between Southeast Alaska and other geographic areas. SNP variation and other genetic data are discussed regarding taxonomy and management.

Peptidomic analysis of healthy and subclinically mastitic bovine milk.

A variety of proteases release hundreds of endogenous peptide fragments from intact bovine milk proteins. Mass spectrometry-based peptidomics allows for high throughput sequence assignment of a large number of these peptides. Mastitis is known to result in increased protease activity in the mammary gland. Therefore, we hypothesized that subclinically mastitic milks would contain higher concentrations of released peptides. In this work, milks were sampled from three cows and, for each, one healthy and one subclinically mastitic teat were sampled for milk. Peptides were analyzed by nano-liquid chromatography quadrupole time of flight tandem mass spectrometry and identified with database searching. In total, 682 peptides were identified. The total number of released peptides increased 146% from healthy to subclinically mastitic milks (p < 0.05), and the total abundance of released peptides also increased significantly (p < 0.05). Bioinformatic analysis of enzyme cleavage revealed increases in activity of cathepsin D and elastase (p < 0.05) with subclinical mastitis.

Bayesian analysis of additive epistasis arising from new mutations in mice.

The continuous uploading of polygenic additive mutational variability has been reported in several studies in laboratory species with an inbred genetic background. These studies have focused on the direct contribution of new mutations without considering the possibility of epistatic effects derived from the interaction of new mutations with pre-existing polymorphisms. In this work we focused on this main topic and analysed the statistical and biological relevance of the epistatic variance for 9 week body weight in two populations of inbred mice. We developed a new linear mixed model parameterization where founder-related additive genetic variability, additive mutational variability and the interaction terms between both sources of variation were accounted for under a Bayesian design and without requiring the inversion of a matrix of epistatic genetic covariances. The analyses focused on a six-generations data set from C57BL/6J mice (n = 3736) and a five-generations data set from C57BL/6J(hg/hg) mice (n = 2843). The deviance information criterion (DIC) clearly favoured the model accounting for epistatic variability with reductions larger than 50 DIC units in both populations. Modal estimates for founder related, mutational and epistatic heritabilities were 0·068, 0·011 and 0·095 in C57BL/6J and 0·060, 0·010 and 0·113 in C57BL/6J(hg/hg), ruling out any doubt about the biological relevance of epistasis originating from new mutations in mice. These results contribute new insights on the relevance of epistasis in the genetic architecture of mammals and serve as an important component of an additional source of genetic heterogeneity for inbred strains of laboratory mice.

Molecular phylogeny and SNP variation of polar bears (Ursus maritimus), brown bears (U. arctos), and black bears (U. americanus) derived from genome sequences.

We assessed the relationships of polar bears (Ursus maritimus), brown bears (U. arctos), and black bears (U. americanus) with high throughput genomic sequencing data with an average coverage of 25× for each species. A total of 1.4 billion 100-bp paired-end reads were assembled using the polar bear and annotated giant panda (Ailuropoda melanoleuca) genome sequences as references. We identified 13.8 million single nucleotide polymorphisms (SNP) in the 3 species aligned to the polar bear genome. These data indicate that polar bears and brown bears share more SNP with each other than either does with black bears. Concatenation and coalescence-based analysis of consensus sequences of approximately 1 million base pairs of ultraconserved elements in the nuclear genome resulted in a phylogeny with black bears as the sister group to brown and polar bears, and all brown bears are in a separate clade from polar bears. Genotypes for 162 SNP loci of 336 bears from Alaska and Montana showed that the species are genetically differentiated and there is geographic population structure of brown and black bears but not polar bears.

Glycosylation of human milk lactoferrin exhibits dynamic changes during early lactation enhancing its role in pathogenic bacteria-host interactions.

Human milk lactoferrin (hmLF) is the most abundant glycoprotein present in human milk and displays a broad range of protective functions in the gut of newborn infants. hmLF is N-glycosylated, but little is known about the lactation stage-related development of the glycosylation phenotype. hmLF glycosylation from milk samples from five donors during the first 10 weeks of lactation was assessed and observed to be more diverse than previously reported. During this period dynamic changes in glycosylation were observed corresponding to a decrease in glycosylation in the second week followed by an increase in total glycosylation as well as higher order fucosylation thereafter. Gene expression analysis was performed in milk somatic cells from a sixth subject. It was found that fucosyltransferase expression increased during entire period, whereas expression of genes for the oligosaccharyl transferase complex decreased in the second week. The effect of hmLF glycosylation was examined for the protein's ability to affect bacterial binding to epithelial cells. hmLF significantly inhibited pathogen adhesion and purified hmLF glycans significantly reduced Salmonella invasion of colonic epithelial cells to levels associated with non-invasive deletion mutants. This study indicates that hmLF glycosylation is tightly regulated by gene expression and that glyco-variation is involved in modulating pathogen association.

Genomic variation and population structure detected by single nucleotide polymorphism arrays in Corriedale, Merino and Creole sheep.

THE AIM OF THIS STUDY WAS TO INVESTIGATE THE GENETIC DIVERSITY WITHIN AND AMONG THREE BREEDS OF SHEEP: Corriedale, Merino and Creole. Sheep from the three breeds (Merino n = 110, Corriedale n = 108 and Creole n = 10) were genotyped using the Illumina Ovine SNP50 beadchip(®). Genetic diversity was evaluated by comparing the minor allele frequency (MAF) among breeds. Population structure and genetic differentiation were assessed using STRUCTURE software, principal component analysis (PCA) and fixation index (FST). Fixed markers (MAF = 0) that were different among breeds were identified as specific breed markers. Using a subset of 18,181 single nucleotide polymorphisms (SNPs), PCA and STUCTURE analysis were able to explain population stratification within breeds. Merino and Corriedale divergent lines showed high levels of polymorphism (89.4% and 86% of polymorphic SNPs, respectively) and moderate genetic differentiation (FST = 0.08) between them. In contrast, Creole had only 69% polymorphic SNPs and showed greater genetic differentiation from the other two breeds (FST = 0.17 for both breeds). Hence, a subset of molecular markers present in the OvineSNP50 is informative enough for breed assignment and population structure analysis of commercial and Creole breeds.