FABP4 gene has a very large effect on feed efficiency in lactating Israeli Holstein cows.

Improving feed efficiency (FE) is a major goal for the livestock industry. Previously, we have identified 48 SNP markers distributed over 32 genes significantly associated with residual feed intake (RFI) in Israeli Holstein male calves, the most significant of which are located in the bovine FABP4 gene. In the present study, we tested associations of eight of the FABP4 markers with RFI and feed conversion ratio (FCR), along with milk composition and feeding behavioral traits, in 114 lactating Israeli Holstein cows. Large allele effects were found, along with large contributions of FABP4 markers to the phenotypic variation [mean contribution of all significant markers (P < 0.05), 15.4 and 12.0% for RFI and FCR, respectively] and genotypic variation [means of all significant markers (P < 0.05), 75.7 and 32.4% in RFI and FCR, respectively]. However, the association of all significant FABP4 markers with FE and milk content traits was found in opposite directions, such that improved FE was accompanied by decreased milk content. Hence, before inclusion in breeding programs, the gain in FE must be economically balanced with the loss in milk contents. On the other hand, these findings imply that in any current improvement program concentrated on milk traits alone, without taking into account the effect on FE, the progress in milk composition is probably accompanied by deterioration of FE. These results, if confirmed in other populations and breeds, set FABP4 as a prime candidate in any marker-assisted selection program targeting FE as a whole and RFI in particular.

FABP4 is a leading candidate gene associated with residual feed intake in growing Holstein calves.

Ecological and economic concerns drive the need to improve feed utilization by domestic animals. Residual feed intake (RFI) is one of the most acceptable measures for feed efficiency (FE). However, phenotyping RFI-related traits is complex and expensive and requires special equipment. Advances in marker technology allow the development of various DNA-based selection tools. To assimilate these technologies for the benefit of RFI-based selection, reliable phenotypic measures are prerequisite. In the current study, we identified single nucleotide polymorphisms (SNPs) associated with RFI phenotypic consistency across different ages and diets (named RFI 1-3), using DNA samples of high or low RFI ranked Holstein calves. Using targeted sequencing of chromosomal regions associated with FE- and RFI-related traits, we identified 48 top SNPs significantly associated with at least one of three defined RFIs. Eleven of these SNPs were harbored by the fatty acid binding protein 4 (FABP4). While 10 significant SNPs found in FABP4 were common for RFI 1 and RFI 3, one SNP (FABP4_5; A

Whole-genome resequencing of two elite sires for the detection of haplotypes under selection in dairy cattle.

Using a combination of whole-genome resequencing and high-density genotyping arrays, genome-wide haplotypes were reconstructed for two of the most important bulls in the history of the dairy cattle industry, Pawnee Farm Arlinda Chief ("Chief") and his son Walkway Chief Mark ("Mark"), each accounting for ∼7% of all current genomes. We aligned 20.5 Gbp (∼7.3× coverage) and 37.9 Gbp (∼13.5× coverage) of the Chief and Mark genomic sequences, respectively. More than 1.3 million high-quality SNPs were detected in Chief and Mark sequences. The genome-wide haplotypes inherited by Mark from Chief were reconstructed using ∼1 million informative SNPs. Comparison of a set of 15,826 SNPs that overlapped in the sequence-based and BovineSNP50 SNPs showed the accuracy of the sequence-based haplotype reconstruction to be as high as 97%. By using the BovineSNP50 genotypes, the frequencies of Chief alleles on his two haplotypes then were determined in 1,149 of his descendants, and the distribution was compared with the frequencies that would be expected assuming no selection. We identified 49 chromosomal segments in which Chief alleles showed strong evidence of selection. Candidate polymorphisms for traits that have been under selection in the dairy cattle population then were identified by referencing Chief's DNA sequence within these selected chromosome blocks. Eleven candidate genes were identified with functions related to milk-production, fertility, and disease-resistance traits. These data demonstrate that haplotype reconstruction of an ancestral proband by whole-genome resequencing in combination with high-density SNP genotyping of descendants can be used for rapid, genome-wide identification of the ancestor's alleles that have been subjected to artificial selection.

Use of proximity loggers to Estimate Affiliative and Agonistic Relationships among Group-housed Holstein calves.

We used proximity loggers to map the social network of 15 group-housed suckling Holstein calves over a 70-day period divided into five trial periods. Tag ID, encountered tag ID, and encounter start time and length were collected from the proximity loggers. Use of 3D matrices and trendlines distinguished three encounter types: "rapid decline" (up to ~ 120 s), "moderate decline" (~120 to ~ 600 s), and "exponential decline" (longer than ~ 600 s). We hypothesized that a random encounter between two calves with an agonistic relationship would terminate quickly, whereas affiliative encounters would linger. Quadratic assignment with 5000 random permutations revealed a significant negative correlation between short encounters ("agonistic") and long encounters ("affiliative") matrices in periods 1 and 5 (r = -0.490 and -0.767, respectively), tendency toward negative correlation in period 3 (r = -0.141, p = 0.104), and no significant correlations in periods 2 and 4. These results were in accordance with the social network analyses revealing agonistic (short encounters) relationships in those periods. We suggest using this method in further studies on the effects of various factors on social network dynamics and its effect on health, intake, growth, and efficiency.

Re-Evaluation of Genotyping Methodologies in Cattle: The Proficiency of Imputation.

In dairy cattle, identifying polymorphisms that contribute to complex economical traits such as residual feed intake (RFI) is challenging and demands accurate genotyping. In this study, we compared imputed genotypes (n = 192 cows) to those obtained using the TaqMan and high-resolution melting (HRM) methods (n = 114 cows), for mutations in the FABP4 gene that had been suggested to have a large effect on RFI. Combining the whole genome sequence (n = 19 bulls) and the cows' BovineHD BeadChip allowed imputing genotypes for these mutations that were verified by Sanger sequencing, whereas, an error rate of 11.6% and 10.7% were encountered for HRM and TaqMan, respectively. We show that this error rate seriously affected the linkage-disequilibrium analysis that supported this gene candidacy over other BTA14 gene candidates. Thus, imputation produced superior genotypes and should also be regarded as a method of choice to validate the reliability of the genotypes obtained by other methodologies that are prone to genotyping errors due to technical conditions. These results support the view that RFI is a complex trait and that searching for the causative sequence variation underlying cattle RFI should await the development of statistical methods suitable to handle additive and epistatic interactions.

Multisite haplotype on cattle chromosome 3 is associated with quantitative trait locus effects on lactation traits.

The goal of this study was to identify candidate genes and DNA polymorphisms for quantitative trait loci (QTL) affecting milk yield (MY), fat yield (FY), and protein yield (PY) previously mapped to bovine chromosome 3 (BTA3). To accomplish this, 373 half-siblings sired by three bulls previously shown to be segregating for lactation trait QTL, and 263 additional sires in the U.S. Dairy Bull DNA Repository (DBDR) were genotyped for 2,500 SNPs within a 16.3 Mbp QTL critical region on BTA3. Targeted resequencing of ∼1.8 Mbp within the QTL critical region of one of the QTL heterozygous sires identified additional polymorphisms useful for association studies. Twenty-three single nucleotide polymorphisms (SNPs) within a fine-mapped region were associated with effects on breeding values for MY, FY, or PY in DBDR sires, of which five SNPs were in strong linkage disequilibrium in the population. This multisite haplotype included SNPs located within exons or promoters of four tightly linked genes: RAP1A, ADORA3, OVGP1, and C3H1orf88. An SNP within RAP1A showed strong evidence of a recent selective sweep based on integrated haplotype score and was also associated with breeding value for PY. Because of its known function in alveolar lumen formation in the mammary gland, RAP1A is thus a strong candidate gene for QTL effects on lactation traits. Our results provide a detailed assessment of a QTL region that will be a useful guide for complex traits analysis in humans and other noninbred species.

The Meat Quality Characteristics of Holstein Calves: The Story of Israeli 'Dairy Beef'.

Global animal production systems are often criticized for their lack of sustainability and insufficient resilience to ensure food security. The 'farm-to-fork' approach aims at orienting food systems towards the creation of a positive environmental impact, nutritious, healthy, safe and sufficient foods, and fairer economic returns for primary producers. Many countries rely on an imported supply of live animals to fulfill their needs for fresh meat. In Israel, ~60% of the sources of fresh beef come from the import of live animals. In order to encourage sustainable beef production in Israel, the proportion of local beef should be raised at the expense of imported animals. However, for this to be achieved, the superior performance of local beef should be justified. The current study was conducted to compare between the meat quality characteristics of local (Israeli Holstein; N = 205) vs. imported (Australian; N = 169) animals. Generally, while the imported calves presented a higher dressing percentage (p < 0.0001), the local animals were characterized by tenderer meat (p < 0.0001), longer sarcomeres (p < 0.0001), higher a* color attributes and pH (p < 0.001), superior cooking (p = 0.002) and thawing loss (p < 0.0001), higher intra-muscular fat (IMF) content, and a higher PUFA proportion (p < 0.01 and p < 0.0001, respectively) and PUFA:SFA ratio. The findings shown herein may provide sound arguments for stakeholders and policy makers to facilitate sustainable local beef production in Israel.

Dietary pomegranate peel improves milk quality of lactating ewes: Emphasis on milk fat globule membrane properties and antioxidative traits.

Concentrated pomegranate peel extract (CPE) was supplemented to ewes, and milk yield and fat content-fatty acid (FA) and phospholipid (PL) composition-were monitored. CPE-fed ewes had higher milk yield, and fat, protein and lactose contents than controls. Milk PL content-20% higher in the CPE-supplemented group-was regulated by treatment and not by total fat content; milk phosphatidylethanolamine and phosphatidylcholine increased by 22 and 26%, respectively, in CPE-supplemented vs. control ewes. Milk saturated FA concentration was higher, and total polyunsaturated and monounsaturated FA content lower in the CPE vs. control group, regardless of milk total fat content. CPE supplementation increased milk antioxidant capacity, suggesting antioxidant transfer from dietary source to milk, increasing stability and nutritive value. Our study provides first evidence for milk quality improvement in terms of antioxidants and PL enrichment without compromising total milk fat, suggesting strategies to improve dairy animals' milk composition without compromising total production.

Efficient inference of haplotypes from genotypes on a large animal pedigree.

We present a simple algorithm for reconstruction of haplotypes from a sample of multilocus genotypes. The algorithm is aimed specifically for analysis of very large pedigrees for small chromosomal segments, where recombination frequency within the chromosomal segment can be assumed to be zero. The algorithm was tested both on simulated pedigrees of 155 individuals in a family structure of three generations and on real data of 1149 animals from the Israeli Holstein dairy cattle population, including 406 bulls with genotypes, but no females with genotypes. The rate of haplotype resolution for the simulated data was >91% with a standard deviation of 2%. With 20% missing data, the rate of haplotype resolution was 67.5% with a standard deviation of 1.3%. In both cases all recovered haplotypes were correct. In the real data, allele origin was resolved for 22% of the heterozygous genotypes, even though 70% of the genotypes were missing. Haplotypes were resolved for 36% of the males. Computing time was insignificant for both data sets. Despite the intricacy of large-scale real pedigree genotypes, the proposed algorithm provides a practical rule-based solution for resolving haplotypes for small chromosomal segments in commercial animal populations.

Structural maintenance of chromosome complexes and bone development: the beginning of a wonderful relationship?

Bone development depends on environmental, nutritional and hormonal factors. Yet, an ordered and timed activation of genes and their associated molecular pathways are central for the growth and development of healthy bones. The correct expression of genes depends on both cis- and trans-regulatory elements. Of these, the elusive role of chromatin ultrastructure is just beginning to become appreciated. Changes in the higher-order structure of chromatin are affecting the expression of genes in response to intrinsic and environmental signals. Cohesin and condensin are members of the structural maintenance of chromosome (SMC) family of protein complexes, which mediate higher-order chromatin structure by tethering distinct regions of chromatin either inter- or intra-molecularly. In recent years, SMCs had been identified for their function in the regulation of gene expression and developmental processes, whereas malfunction of cohesin or condensin has an impact on human health. However, little is known about the specific roles of SMC complexes in bone development and their possible effect on bone health. Here, we review studies that suggest an intimate link between SMCs and bone development, as well as a plausible effect, direct or indirect, on the bone health. We describe genetic syndromes associated with SMCs with distinctive bone phenotypes and identify links between SMCs and bone-related molecular pathways. Future studies of the relationship between SMCs and bone development will reveal new understandings of both the cellular and molecular roles of SMC complexes and provide new insights into the growth and developmental processes in the bone.

Osteoporosis genetics: year 2011 in review.

Increased rates of osteoporotic fractures represent a worldwide phenomenon, which result from a progressing aging in the population around the world and creating socioeconomic problems. This review will focus mostly on human genetic studies identifying genomic regions, genes and mutations associated with osteoporosis (bone mineral density (BMD) and bone loss) and related fractures, which were published during 2011. Although multiple genome-wide association studies (GWAS) were performed to date, the genetic cause of osteoporosis and fractures has not yet been found, and only a small fraction of high heritability of bone mass was successfully explained. GWAS is a successful tool to initially define and prioritize specific chromosomal regions showing associations with the desired traits or diseases. Following the initial discovery and replication, targeted sequencing is needed in order to detect those rare variants which GWAS does not reveal by design. Recent GWAS findings for BMD included WNT16 and MEF2C. The role of bone morphogenetic proteins in fracture healing has been explored by several groups, and new single-nucleotide polymorphisms present in genes such as NOGGIN and SMAD6 were found to be associated with a greater risk of fracture non-union. Finding new candidate genes, and mutations associated with BMD and fractures, also provided new biological connections. Thus, candidates for molecular link between bone metabolism and lactation (for example, RAP1A gene), as well as possible pleiotropic effects for bone and muscle (ACTN3 gene) were suggested. The focus of contemporary studies seems to move toward whole-genome sequencing, epigenetic and functional genomics strategies to find causal variants for osteoporosis.

The genetic pleiotropy of musculoskeletal aging.

Musculoskeletal aging is detrimental to multiple bodily functions and starts early, probably in the fourth decade of an individual's life. Sarcopenia is a health problem that is expected to only increase as a greater portion of the population lives longer; prevalence of the related musculoskeletal diseases is similarly expected to increase. Unraveling the biological and biomechanical associations and molecular mechanisms underlying these diseases represents a formidable challenge. There are two major problems making disentangling the biological complexity of musculoskeletal aging difficult: (a) it is a systemic, rather than "compartmental," problem, which should be approached accordingly, and (b) the aging per se is neither well defined nor reliably measurable. A unique challenge of studying any age-related condition is a need of distinguishing between the "norm" and "pathology," which are interwoven throughout the aging organism. We argue that detecting genes with pleiotropic functions in musculoskeletal aging is needed to provide insights into the potential biological mechanisms underlying inter-individual differences insusceptibility to the musculoskeletal diseases. However, exploring pleiotropic relationships among the system's components is challenging both methodologically and conceptually. We aimed to focus on genetic aspects of the cross-talk between muscle and its "neighboring" tissues and organs (tendon, bone, and cartilage), and to explore the role of genetics to find the new molecular links between skeletal muscle and other parts of the "musculoskeleton." Identification of significant genetic variants underlying the musculoskeletal system's aging is now possible more than ever due to the currently available advanced genomic technologies. In summary, a "holistic" genetic approach is needed to study the systems's normal functioning and the disease predisposition in order to improve musculoskeletal health.

Accelerated expansion of group IID-like phospholipase A2 genes in Bos taurus.

Low-molecular-weight, calcium-dependent phospholipase A2 genes (PLA2s) that belong to the secreted type of PLA2s are clustered within a syntenic group on human 1p35-p36 and mouse 4qD3. We reassembled trace files available from the Whole Genome Sequencing (WGS) Project, obtaining an 86-kb contig with three tandem PLA2G2D duplications in the Hereford strain. We used mate-pair data to monitor the assembly and to exclude chimeric clones, demonstrating that the current WGS data may be assembled even in a highly repetitive region with a coverage exceeding fivefold. The genomic structure indicated that most of the PLA2G2D transcripts are formed by four exons. Two alternative first exons were present in all duplications. In two duplications insertions of satellite DNA in the third intron created a novel exon that gave rise to a two-exon product. Linkage and comparative mapping placed the bovine PLA2G2 locus on BTA2, indicating that it evolved from an ancestral PLA2G2D locus common to human, cattle, and rodents. Bovine PLA2G2D variants were capable of encoding 147-amino-acid polypeptides that consisted of putative signal peptide and metal-binding domains. Cysteine residues were conserved in positions analogous to those forming the seven disulfide bonds characteristic of PLA2G2 genes. Quantitative PCR analysis of bovine PLA2G2D transcripts indicated that their expression levels varied between the dry period and lactation in the mammary gland samples and that their expression was polymorphic in liver tissue. The recent burst of duplication and divergence of the bovine PLA2G2D genes and their polymorphic nature are typical of innate immune response genes.

Identification of a missense mutation in the bovine ABCG2 gene with a major effect on the QTL on chromosome 6 affecting milk yield and composition in Holstein cattle.

We previously localized a quantitative trait locus (QTL) on chromosome 6 affecting milk fat and protein concentration to a 4-cM confidence interval, centered on the microsatellite BM143. We characterized the genes and sequence variation in this region and identified common haplotypes spanning five polymorphic sites in the genes IBSP, SPP1, PKD2, and ABCG2 for two sires heterozygous for this QTL. Expression of SPP1 and ABCG2 in the bovine mammary gland increased from parturition through lactation. SPP1 and all the coding exons of ABCG2 and PKD2 were sequenced for these two sires. The single nucleotide change capable of encoding a substitution of tyrosine-581 to serine (Y581S) in the ABCG2 transporter was the only polymorphism corresponding to the segregation status of all 3 heterozygous and 15 homozygous sires for the QTL in the Israeli and U.S. Holstein populations. The allele substitution fixed effects on the genetic evaluations of 335 Israeli sires were -341 kg milk, +0.16% fat, and +0.13% protein (F-value = 200). No other polymorphism gave significant effect for fat and protein concentration in models that also included Y581S. The allele substitution effects on the genetic evaluations of 670 cows, daughters of two heterozygous sires, were -226 kg milk, 0.09% fat, and 0.08% protein (F-value = 394), with partial dominance towards the 581S homozygotes. We therefore propose that Y581S in ABCG2 is the causative site for this QTL.