Luciana Correia de Almeida Regitano

Analysis of copy number variations among diverse cattle breeds

Genomic structural variation is an important and abundant source of genetic and phenotypic variation. Here, we describe the first systematic and genome-wide analysis of copy number variations (CNVs) in modern domesticated cattle using array comparative genomic hybridization (array CGH), quantitative PCR (qPCR), and fluorescent in situ hybridization (FISH). The array CGH panel included 90 animals from 11 Bos taurus, three Bos indicus, and three composite breeds for beef, dairy, or dual purpose. We identified over 200 candidate CNV regions (CNVRs) in total and 177 within known chromosomes, which harbor or are adjacent to gains or losses. These 177 high-confidence CNVRs cover 28.1 megabases or approximately 1.07% of the genome. Over 50% of the CNVRs (89/177) were found in multiple animals or breeds and analysis revealed breed-specific frequency differences and reflected aspects of the known ancestry of these cattle breeds. Selected CNVs were further validated by independent methods using qPCR and FISH. Approximately 67% of the CNVRs (119/177) completely or partially span cattle genes and 61% of the CNVRs (108/177) directly overlap with segmental duplications. The CNVRs span about 400 annotated cattle genes that are significantly enriched for specific biological functions, such as immunity, lactation, reproduction, and rumination. Multiple gene families, including ULBP, have gone through ruminant lineage-specific gene amplification. We detected and confirmed marked differences in their CNV frequencies across diverse breeds, indicating that some cattle CNVs are likely to arise independently in breeds and contribute to breed differences. Our results provide a valuable resource beyond microsatellites and single nucleotide polymorphisms to explore the full dimension of genetic variability for future cattle genomic research.

Worldwide Patterns of Ancestry, Divergence, and Admixture in Domesticated Cattle

The domestication and development of cattle has considerably impacted human societies, but the histories of cattle breeds and populations have been poorly understood especially for African, Asian, and American breeds. Using genotypes from 43,043 autosomal single nucleotide polymorphism markers scored in 1,543 animals, we evaluate the population structure of 134 domesticated bovid breeds. Regardless of the analytical method or sample subset, the three major groups of Asian indicine, Eurasian taurine, and African taurine were consistently observed. Patterns of geographic dispersal resulting from co-migration with humans and exportation are recognizable in phylogenetic networks. All analytical methods reveal patterns of hybridization which occurred after divergence. Using 19 breeds, we map the cline of indicine introgression into Africa. We infer that African taurine possess a large portion of wild African auroch ancestry, causing their divergence from Eurasian taurine. We detect exportation patterns in Asia and identify a cline of Eurasian taurine/indicine hybridization in Asia. We also identify the influence of species other than Bos taurus taurus and B. t. indicus in the formation of Asian breeds. We detect the pronounced influence of Shorthorn cattle in the formation of European breeds. Iberian and Italian cattle possess introgression from African taurine. American Criollo cattle originate from Iberia, and not directly from Africa with African ancestry inherited via Iberian ancestors. Indicine introgression into American cattle occurred in the Americas, and not Europe. We argue that cattle migration, movement and trading followed by admixture have been important forces in shaping modern bovine genomic variation.

Genome-Wide Association for Growth Traits in Canchim Beef Cattle

Studies are being conducted on the applicability of genomic data to improve the accuracy of the selection process in livestock, and genome-wide association studies (GWAS) provide valuable information to enhance the understanding on the genetics of complex traits. The aim of this study was to identify genomic regions and genes that play roles in birth weight (BW), weaning weight adjusted for 210 days of age (WW), and long-yearling weight adjusted for 420 days of age (LYW) in Canchim cattle. GWAS were performed by means of the Generalized Quasi-Likelihood Score (GQLS) method using genotypes from the BovineHD BeadChip and estimated breeding values for BW, WW, and LYW. Data consisted of 285 animals from the Canchim breed and 114 from the MA genetic group (derived from crossings between Charolais sires and ½ Canchim + ½ Zebu dams). After applying a false discovery rate correction at a 10% significance level, a total of 4, 12, and 10 SNPs were significantly associated with BW, WW, and LYW, respectively. These SNPs were surveyed to their corresponding genes or to surrounding genes within a distance of 250 kb. The genes DPP6 (dipeptidyl-peptidase 6) and CLEC3B (C-type lectin domain family 3 member B) were highlighted, considering its functions on the development of the brain and skeletal system, respectively. The GQLS method identified regions on chromosome associated with birth weight, weaning weight, and long-yearling weight in Canchim and MA animals. New candidate regions for body weight traits were detected and some of them have interesting biological functions, of which most have not been previously reported. The observation of QTL reports for body weight traits, covering areas surrounding the genes (SNPs) herein identified provides more evidence for these associations. Future studies targeting these areas could provide further knowledge to uncover the genetic architecture underlying growth traits in Canchim cattle.

Association of GH and IGF-1 polymorphisms with growth traits in a synthetic beef cattle breed

The Canchim beef cattle (5/8 Charolais + 3/8 Zebu) has been selected for meat production in Brazil since late 1950. In the present work the effects of growth hormone (GH) and insulin-like growth factor 1 (IGF-1) polymorphisms were investigated in 688 animals born between 1998 and 2000. These animals belonged to two genetic groups, i.e., traditional and new lineages. Genotype effects on expected breeding values for birth weight (BW), weaning weight (WW) and yearling weight (YW) were investigated by the least square method. Significant effects were found for GH genotype on YW (p ≤ 0.05), with positive effects associated with the LV (leucine/valine) genotype. For IGF-1 genotypes, significant effects were found on BW (p ≤ 0.01) and YW (p ≤ 0.01). Average substitution effects for IGF-1 alleles estimated by regression analysis suggested a positive effect of the IGF-1 225 bp allele on BW and of the 229 bp allele on YW.

Genome-wide association study for intramuscular fat deposition and composition in Nellore cattle

BackgroundMeat from Bos taurus and Bos indicus breeds are an important source of nutrients for humans and intramuscular fat (IMF) influences its flavor, nutritional value and impacts human health. Human consumption of fat that contains high levels of monounsaturated fatty acids (MUFA) can reduce the concentration of undesirable cholesterol (LDL) in circulating blood. Different feeding practices and genetic variation within and between breeds influences the amount of IMF and fatty acid (FA) composition in meat. However, it is difficult and costly to determine fatty acid composition, which has precluded beef cattle breeding programs from selecting for a healthier fatty acid profile. In this study, we employed a high-density single nucleotide polymorphism (SNP) chip to genotype 386 Nellore steers, a Bos indicus breed and, a Bayesian approach to identify genomic regions and putative candidate genes that could be involved with deposition and composition of IMF.ResultsTwenty-three genomic regions (1-Mb SNP windows) associated with IMF deposition and FA composition that each explain ≥ 1% of the genetic variance were identified on chromosomes 2, 3, 6, 7, 8, 9, 10, 11, 12, 17, 26 and 27. Many of these regions were not previously detected in other breeds. The genes present in these regions were identified and some can help explain the genetic basis of deposition and composition of fat in cattle.ConclusionsThe genomic regions and genes identified contribute to a better understanding of the genetic control of fatty acid deposition and can lead to DNA-based selection strategies to improve meat quality for human consumption.

Global liver gene expression differences in Nelore steers with divergent residual feed intake phenotypes

BackgroundEfficiency of feed utilization is important for animal production because it can reduce greenhouse gas emissions and improve industry profitability. However, the genetic basis of feed utilization in livestock remains poorly understood. Recent developments in molecular genetics, such as platforms for genome-wide genotyping and sequencing, provide an opportunity to identify genes and pathways that influence production traits. It is known that transcriptional networks influence feed efficiency-related traits such as growth and energy balance. This study sought to identify differentially expressed genes in animals genetically divergent for Residual Feed Intake (RFI), using RNA sequencing methodology (RNA-seq) to obtain information from genome-wide expression profiles in the liver tissues of Nelore cattle.ResultsDifferential gene expression analysis between high Residual Feed Intake (HRFI, inefficient) and low Residual Feed Intake (LRFI, efficient) groups was performed to provide insights into the molecular mechanisms that underlie feed efficiency-related traits in beef cattle. A total of 112 annotated genes were identified as being differentially expressed between animals with divergent RFI phenotypes. These genes are involved in ion transport and metal ion binding; act as membrane or transmembrane proteins; and belong to gene clusters that are likely related to the transport and catalysis of molecules through the cell membrane and essential mechanisms of nutrient absorption. Genes with functions in cellular signaling, growth and proliferation, cell death and survival were also differentially expressed. Among the over-represented pathways were drug or xenobiotic metabolism, complement and coagulation cascades, NRF2-mediated oxidative stress, melatonin degradation and glutathione metabolism.ConclusionsOur data provide new insights and perspectives on the genetic basis of feed efficiency in cattle. Some previously identified mechanisms were supported and new pathways controlling feed efficiency in Nelore cattle were discovered. We potentially identified genes and pathways that play key roles in hepatic metabolic adaptations to oxidative stress such as those involved in antioxidant mechanisms. These results improve our understanding of the metabolic mechanisms underlying feed efficiency in beef cattle and will help develop strategies for selection towards the desired phenotype.

Linkage disequilibrium and haplotype block structure in a composite beef cattle breed.

BackgroundThe development of linkage disequilibrium (LD) maps and the characterization of haplotype block structure at the population level are useful parameters for guiding genome wide association (GWA) studies, and for understanding the nature of non-linear association between phenotypes and genes. The elucidation of haplotype block structure can reduce the information of several single nucleotide polymorphisms (SNP) into the information of a haplotype block, reducing the number of SNPs in a coherent way for consideration in GWA and genomic selection studies.ResultsThe maximum average LD, measured by r2 varied between 0.33 to 0.40 at a distance of < 2.5 kb, and the minimum average values of r2 varied between 0.05 to 0.07 at distances ranging from 400 to 500 kb, clearly showing that the average r2 reduced with the increase in SNP pair distances. The persistence of LD phase showed higher values at shorter genomic distances, decreasing with the increase in physical distance, varying from 0.96 at a distance of < 2.5 kb to 0.66 at a distance from 400 to 500 kb. A total of 78% of all SNPs were clustered into haplotype blocks, covering 1,57 Mb of the total autosomal genome size.ConclusionsThis study presented the first high density linkage disequilibrium map and haplotype block structure for a composite beef cattle population, and indicates that the high density SNP panel over 700 k can be used for genomic selection implementation and GWA studies for Canchim beef cattle.

QTL affecting body weight in a candidate region of cattle chromosome 5

The objective of this work was to identify QTLs for liveweight in a candidate region of bovine chromosome 5. Half-sib families from two lines, one traditional and the other new, of Canchim beef cattle (5/8 Charolais + 3/8 Zebu) were genotyped for four microsatellite markers, including the microsatellite in the IGF-1 (insulin-like growth factor-1) promoter region. Significant differences in allele distribution between the two lines were found for three markers. Interval mapping analyses in this region indicated the presence of a QTL controlling birth weight (p < 0.05) and of a QTL influencing breeding value for yearling weight (p < 0.01) in the newer line of the breed. The previously identified interaction between the IGF-1 genotype and genetic group strengthens the hypothesis of a linked QTL rather than an IGF-1 effect on growth traits in the Canchim cattle.

k -Casein, b -lactoglobulin and growth hormone allele frequencies and genetic distances in Nelore, Gyr, Guzerá, Caracu, Charolais, Canchim and Santa Gertrudis cattle

The genotypes for k-casein (k-CN), b-lactoglobulin (b-LG) and growth hormone (GH) were determined by polymerase chain reaction (PCR) and restriction enzyme digestion in seven breeds of cattle (Nelore, Gyr, Guzera, Caracu, Charolais, Canchim and Santa Gertrudis). k-Casein had two alleles with the A allele occurring at a higher frequency in Bos indicus breeds (0.93, 0.92 and 0.91% for Gyr, Guzera and Nelore, respectively). The b-lactoglobulin locus had two alleles in all of the breeds. European breeds had a higher frequency of the b-LG A allele than Zebu breeds. The GH locus had two alleles (L and V) in Bos taurus and was monomorphic (L allele only) in all of the Bos indicus breeds evaluated. The highest frequency for the V allele was observed in Charolais cattle. The markers used revealed a considerable similarity among breeds, with two main groups being discernible. One group consisted of Zebu and Santa Gertrudis breeds and the other consisted of European and Canchim breeds.

Genetic characterization of Aberdeen Angus cattle using molecular markers

Aberdeen Angus beef cattle from the Brazilian herd were studied genetically using restriction fragment length polymorphism (RFLP) of the k-casein - HinfI (CSN3 - HinfI), b-lactoglobulin - HaeIII (LGB - HaeIII) and growth hormone AluI (GH- AluI) genes, as well as four microsatellites (TEXAN15, CSFM50, BM1224 and BM7160). The RFLP genotypes were determined using the polymerase chain reaction (PCR) followed by digestion with restriction endonucleases and electrophoresis in agarose gels. With the exception of the microsatellite BM7160, which was analyzed in an automatic sequencer, the PCR products were genotyped by silver staining. The allele and genotype frequencies, heterozygosities and gene diversity were estimated. The values for these parameters of variability were comparable to other cattle breeds. The genetic relationship of the Aberdeen Angus to other breeds (Caracu, Canchim, Charolais, Guzerath, Gyr, Nelore, Santa Gertrudis and Simmental) was investigated using Neis genetic distance. Cluster analysis placed the Aberdeen Angus in an isolated group in the Bos taurus breeds branch. This fact is in agreement with the geographic origin of this breed.