Whole-Genome Analyses of Korean Native and Holstein Cattle Breeds by Massively Parallel Sequencing

A main goal of cattle genomics is to identify DNA differences that account for variations in economically important traits. In this study, we performed whole-genome analyses of three important cattle breeds in Korea—Hanwoo, Jeju Heugu, and Korean Holstein—using the Illumina HiSeq 2000 sequencing platform. We achieved 25.5-, 29.6-, and 29.5-fold coverage of the Hanwoo, Jeju Heugu, and Korean Holstein genomes, respectively, and identified a total of 10.4 million single nucleotide polymorphisms (SNPs), of which 54.12% were found to be novel. We also detected 1,063,267 insertions–deletions (InDels) across the genomes (78.92% novel). Annotations of the datasets identified a total of 31,503 nonsynonymous SNPs and 859 frameshift InDels that could affect phenotypic variations in traits of interest. Furthermore, genome-wide copy number variation regions (CNVRs) were detected by comparing the Hanwoo, Jeju Heugu, and previously published Chikso genomes against that of Korean Holstein. A total of 992, 284, and 1881 CNVRs, respectively, were detected throughout the genome. Moreover, 53, 65, 45, and 82 putative regions of homozygosity (ROH) were identified in Hanwoo, Jeju Heugu, Chikso, and Korean Holstein respectively. The results of this study provide a valuable foundation for further investigations to dissect the molecular mechanisms underlying variation in economically important traits in cattle and to develop genetic markers for use in cattle breeding.     

Genome-Wide Study of Structural Variants in Bovine Holstein,Montbéliarde and Normande Dairy Breeds

High-throughput sequencing technologies have offered in recent years new opportunities to study genome variations. These studies have mostly focused on single nucleotide polymorphisms, small insertions or deletions and on copy number variants. Other structural variants, such as large insertions or deletions, tandem duplications, translocations, and inversions are less well-studied, despite that some have an important impact on phenotypes. In the present study, we performed a large-scale survey of structural variants in cattle. We report the identification of 6,426 putative structural variants in cattle extracted from whole-genome sequence data of 62 bulls representing the three major French dairy breeds. These genomic variants affect DNA segments greater than 50 base pairs and correspond to deletions, inversions and tandem duplications. Out of these, we identified a total of 547 deletions and 410 tandem duplications which could potentially code for CNVs. Experimental validation was carried out on 331 structural variants using a novel high-throughput genotyping method. Out of these, 255 structural variants (77%) generated good quality genotypes and 191 (75%) of them were validated. Gene content analyses in structural variant regions revealed 941 large deletions removing completely one or several genes, including 10 single-copy genes. In addition, some of the structural variants are located within quantitative trait loci for dairy traits. This study is a pan-genome assessment of genomic variations in cattle and may provide a new glimpse into the bovine genome architecture. Our results may also help to study the effects of structural variants on gene expression and consequently their effect on certain phenotypes of interest.     

Novel polymorphisms of SIX4 gene and their association with body measurement traits in Qinchuan cattle

Sine oculis homeobox homolog 4 (SIX4) gene belongs to the sine oculis/SIX gene family, which includes six members in vertebrates. SIX4 gene plays a crucial role in skeletal myogenesis, and its genetic variations or deficiency may cause hypopituitarism, suggesting that SIX4 gene is a potential candidate gene affecting body measurement traits (BMTs) in animals. Herein, the objectives of this study were to identify genetic polymorphisms of bovine SIX4 gene and to analyze potential association between single nucleotide polymorphisms (SNPs) and body measurement traits in Qinchuan cattle. In the present study, we investigated polymorphisms of SIX4 gene in 426 Qinchuan cattle using DNA sequencing and polymerase chain reaction–restriction fragment length polymorphisms. Three novel SNPs were identified within bovine SIX4 gene. Associations between body measurement traits and SIX4 gene polymorphisms were investigated, and significant statistical associations were found between polymorphisms of these three SNPs and body measurement traits (P < 0.05). Hence, based on results obtained from this study, we conjectured that SIX4 gene may have potential effects on body measurement traits in Qinchuan cattle population and could be used for marker-assisted selection.     

Rice-Map: a new-generation rice genome browser

Background

Because the Japanese native cattle Kuchinoshima-Ushi have been isolated in a small island and their lineage has been intensely protected, it has been assumed to date that numerous and valuable genomic variations are conserved in this cattle breed.

Results

In this study, we evaluated genetic features of this breed, including single nucleotide polymorphism (SNP) information, by whole-genome sequencing using a Genome Analyzer II. A total of 64.2 Gb of sequence was generated, of which 86% of the obtained reads were successfully mapped to the reference sequence (Btau 4.0) with BWA. On an average, 93% of the genome was covered by the reads and the number of mapped reads corresponded to 15.8-fold coverage across the covered region. From these data, we identified 6.3 million SNPs, of which more than 5.5 million (87%) were found to be new. Out of the SNPs annotated in the bovine sequence assembly, 20,432 were found in protein-coding regions containing 11,713 nonsynonymous SNPs in 4,643 genes. Furthermore, phylogenetic analysis using sequence data from 10 genes (more than 10 kbp) showed that Kuchinoshima-Ushi is clearly distinct from European domestic breeds of cattle.

Conclusions

These results provide a framework for further genetic studies in the Kuchinoshima-Ushi population and research on functions of SNP-containing genes, which would aid in understanding the molecular basis underlying phenotypic variation of economically important traits in cattle and in improving intrinsic defects in domestic cattle breeds.     

Whole-Genome Resequencing Analysis of Hanwoo and Yanbian Cattle to Identify Genome-Wide SNPs and Signatures of Selection

Over the last 30 years, Hanwoo has been selectively bred to improve economically important traits. Hanwoo is currently the representative Korean native beef cattle breed, and it is believed that it shared an ancestor with a Chinese breed, Yanbian cattle, until the last century. However, these two breeds have experienced different selection pressures during recent decades. Here, we whole-genome sequenced 10 animals each of Hanwoo and Yanbian cattle (20 total) using the Illumina HiSeq 2000 sequencer. A total of approximately 3.12 and 3.07 billion sequence reads were mapped to the bovine reference sequence assembly (UMD 3.1) at an average of approximately 10.71- and 10.53-fold coverage for Hanwoo and Yanbian cattle, respectively. A total of 17,936,399 single nucleotide polymorphisms (SNPs) were yielded, of which 22.3% were found to be novel. By annotating the SNPs, we further retrieved numerous nonsynonymous SNPs that may be associated with traits of interest in cattle. Furthermore, we performed whole-genome screening to detect signatures of selection throughout the genome. We located several promising selective sweeps that are potentially responsible for economically important traits in cattle; the PPP1R12A gene is an example of a gene that potentially affects intramuscular fat content. These discoveries provide valuable genomic information regarding potential genomic markers that could predict traits of interest for breeding programs of these cattle breeds.     

Next-generation sequencing-based genome-wide mutation analysis of l-lysine-producing Corynebacterium glutamicum ATCC 21300 strain

In order to identify single nucleotide polymorphism and insertion/deletion mutations, we performed whole-genome re-sequencing of the enhanced l-lysine-producing Corynebacterium glutamicum ATCC 21300 strain. In total, 142 single nucleotide polymorphisms and 477 insertion/deletion mutations were identified in the ATCC 21300 strain when compared to 3,434 predicted genes of the wild-type C. glutamicum ATCC 13032 strain. Among them, 110 transitions and 29 transversions of single nucleotide polymorphisms were found from genes of the ATCC 21300 strain. In addition, 11 genes, involved in the L-lysine biosynthetic pathway and central carbohydrate metabolism, contained mutations including single nucleotide polymorphisms and insertions/deletions. Interestingly, RT-PCR analysis of these 11 genes indicated that they were normally expressed in the ATCC 21300 strain. This information of genome-wide gene-associated variations will be useful for genome breeding of C. glutamicum in order to develop an industrial amino acid-producing strain with minimal mutation.     

Genome-wide analysis of Hanwoo and Chikso populations using the BovineSNP50 genotyping array

Hanwoo and Chikso are classified as Korean native cattle breeds that are currently registered with the Food and Agriculture Organization. However, there is still a lack of genomic studies to compare Hanwoo to Chikso populations. The objective of this study was to perform genome-wide analysis of Hanwoo and Chikso populations, investigating the genetic relationships between these two populations. We genotyped a total of 319 cattle including 214 Hanwoo and 105 Chikso sampled from Gangwon Province Livestock Technology Research Institute, using the Illumina Bovine SNP50K Beadchip. After performing quality control on the initially generated datasets, we assessed linkage disequilibrium patterns for all the possible SNP pairs within 1 Mb apart. Overall, average r² values in Hanwoo (0.048) were lower than Chikso (0.074) population. The genetic relationship between the populations was further assured by the principal component analysis, exhibiting clear clusters in each of the Hanwoo and Chikso populations, respectively. Overall heterozygosity for Hanwoo (0.359) was slightly higher than Chikso (0.345) and inbreeding coefficient was also a bit higher in Hanwoo (??0.015) than Chikso (??0.035). The average F_ST value was 0.036 between Hanwoo and Chikso, indicating little genetic differentiation between those two breeds. Furthermore, we found potential selection signatures including LRP1B and NTRK2 genes that might be implicated with meat and reproductive traits in cattle. In this study, the results showed that both Hanwoo and Chikso populations were not under severe level of inbreeding. Although the principal component analysis exhibited clear clusters in each of the populations, we did not see any clear evidence that those two populations are highly differentiated each other.     

SNP identification in FBXO32 gene and their associations with growth traits in cattle

The F-box protein 32 (FBXO32), also known as Atrogin-1, is one of the four subunits of the ubiquitin protein ligase complex. FBXO32 has been previously shown to be involved in regulation of initiation and development of muscle mass. In the present study, we investigated the polymorphism of FBXO32 gene in 1313 cattle from seven bovine breeds using DNA sequencing, polymerase chain reaction–restriction fragment length polymorphism (PCR–RFLP) and PCR-based amplification-created restriction site (PCR-ACRS) methods. Four novel single nucleotide polymorphisms (SNPs) were identified within bovine FBXO32, and were deposited in the GenBank database. The association studies between these four SNPs and growth traits were performed in NanYang cattle. Notably, the SNPs ss411628932 and ss411628936 were shown to be significantly associated with body length of 24-month-old NanYang cattle. Based on the above four SNPs, 16 haplotypes were identified. The main haplotype was AATA, which occurred at a frequency of more than 40%. Additionally, phylogenetic analysis showed that geographical distance was essential to gene flow among seven cattle breeds. Indigenous bovine breeds displayed genetic difference in comparison to hybrid bovine breeds that have foreign origins. We herein describe for the first time a comprehensive study on the variability of bovine FBXO32 gene that is predictive of genetic potential for body length phenotype.     

Genome-wide discovery of DNA polymorphism in Brassica rapa

Single nucleotide polymorphisms (SNPs) and/or insertion/deletions (InDels) are frequent sequence variations in the plant genome, which can be developed as molecular markers for genetic studies on crop improvement. The ongoing Brassica rapa genome sequencing project has generated vast amounts of sequence data useful in genetic research. Here, we report a genome-wide survey of DNA polymorphisms in the B. rapa genome based on the 557 bacterial artificial clone sequences of B. rapa ssp. pekinensis cv. Chiifu. We identified and characterized 21,311 SNPs and 6,753 InDels in the gene space of the B. rapa genome by re-sequencing 1,398 sequence-tagged sites (STSs) in eight genotypes. Comparison of our findings with a B. rapa genetic linkage map confirmed that STS loci were distributed randomly over the B. rapa whole genome. In the 1.4 Mb of aligned sequences, mean nucleotide polymorphism and diversity were θ = 0.00890 and π = 0.00917, respectively. Additionally, the nucleotide diversity in introns was almost three times greater than that in exons, and the frequency of observed InDel was almost 17 times higher in introns than in exons. Information regarding SNPs/InDels obtained here will provide an important resource for genetic studies and breeding programs of B. rapa.     

Identification of genetic polymorphisms in FABP3 and FABP4 and putative association with back fat thickness in Korean native cattle

The aim of this study was to determine whether single nucleotide polymorphisms (SNP) in the beef cattle adipocyte fatty-acid binding protein 3 and 4 (FABP3 and FABP4) genes are associated with carcass weight (CW) and back fat thickness (BF) of beef cattle. By direct DNA sequencing in 24 unrelated Korean native cattle, we identified 20 SNPs in FABP3 and FABP4. Among them, 10 polymorphic sites were selected for genotyping in our beef cattle. We performed SNP, haplotype and linkage disequilibrium studies on 419 Korean native cattle with the 10 SNPs in the FABP genes. Statistical analysis revealed that 220AG (I74V) and 348+303TC polymorphisms in FABP4 showed putative associations with BF traits (P=0.02 and 0.01, respectively). Our findings suggest that the polymorphisms in FABP4 may play a role in determining one of the important genetic factors that influence BF in beef cattle.     

A novel polymorphism of the myogenin gene is associated with body measurement traits in native Chinese breeds

Using PCR-SSCP and DNA sequencing technology, we examined the association of single nucleotide polymorphisms (SNPs) in the bovine MyoG gene with body measurement traits in 779 individuals of six native Chinese cattle breeds, namely Luxi, Luxi × Simmental crossbred, Nanyang, Xia'nan, Jiaxian red, and Qinchuan. A novel SNP, T314C, was detected. Allelic frequencies of MyoG-T/C in the six breeds were 0.8308/0.1692, 0.8774/0.1226, 0.8021/0.1979, 0.8209/0.1791, 0.8630/0.1370, 0.8044/0.1956, respectively. Least squares analysis revealed a significant (P < 0.05) association of the MyoG SNP with rump length in four breeds (Luxi, Xia'nan, Jiaxian red, and Qinchuan), with hucklebone width in three breeds (Luxi × Simmental crossbred, Nanyang and Xia'nan), with waist height in two breeds (Luxi × Simmental crossbred and Nanyang) and with body length in the Luxi breed. We conclude that the MyoG SNP has potential as a genetic marker for economically relevant body measurement traits in native Chinese cattle breeds.     

Genome‐wide survey of artificial mutations induced by ethyl methanesulfonate and gamma rays in tomato

Genome‐wide mutations induced by ethyl methanesulfonate (EMS) and gamma irradiation in the tomato Micro‐Tom genome were identified by a whole‐genome shotgun sequencing analysis to estimate the spectrum and distribution of whole‐genome DNA mutations and the frequency of deleterious mutations. A total of ~370 Gb of paired‐end reads for four EMS‐induced mutants and three gamma‐ray‐irradiated lines as well as a wild‐type line were obtained by next‐generation sequencing technology. Using bioinformatics analyses, we identified 5920 induced single nucleotide variations and insertion/deletion (indel) mutations. The predominant mutations in the EMS mutants were C/G to T/A transitions, while in the gamma‐ray mutants, C/G to T/A transitions, A/T to T/A transversions, A/T to G/C transitions and deletion mutations were equally common. Biases in the base composition flanking mutations differed between the mutagenesis types. Regarding the effects of the mutations on gene function, >90% of the mutations were located in intergenic regions, and only 0.2% were deleterious. In addition, we detected 1 140 687 spontaneous single nucleotide polymorphisms and indel polymorphisms in wild‐type Micro‐Tom lines. We also found copy number variation, deletions and insertions of chromosomal segments in both the mutant and wild‐type lines. The results provide helpful information not only for mutation research, but also for mutant screening methodology with reverse‐genetic approaches.     

ZBTB38 gene polymorphism associated with body measurement traits in native Chinese cattle breeds

Body measurement traits, influenced by genes and environmental factors, play numerous important roles in the value assessment of productivity and economy. In this study, we investigated the association between genetic polymorphisms of the zinc finger and BTB domain containing 38 gene (ZBTB38) and body measurement traits in native Chinese cattle. Using direct DNA sequencing in 423 individuals of 8 different cattle subpopulations, 9 novel polymorphisms were identified for genotyping within 647 bp region of exon 1 of ZBTB38. Linkage disequilibrium and association analysis revealed that two coding exon polymorphisms (2323 G>A and 2325 C>T polymorphisms), missense mutations valine GTC(T) to isoleucine ATC(T), were associated with body length (BL), withers height (WH) and rump length (RL). Furthermore, the analysis of 2323 G>A and 2325 C>T SNP markers shows that there are significant effects on the BL (P = 0.0488), WH (P = 0.0044) and RL (P = 0.0314) in the total population. These results clearly suggest that the ZBTB38 gene is among the target genes for body measurement traits in bovine breeding, and provide data for establishment of an animal model using cattle to study big animal body type.     

Whole-genome sequencing of the efficient industrial fuel-ethanol fermentative Saccharomyces cerevisiae strain CAT-1

The Saccharomyces cerevisiae strains widely used for industrial fuel-ethanol production have been developed by selection, but their underlying beneficial genetic polymorphisms remain unknown. Here, we report the draft whole-genome sequence of the S. cerevisiae strain CAT-1, which is a dominant fuel-ethanol fermentative strain from the sugarcane industry in Brazil. Our results indicate that strain CAT-1 is a highly heterozygous diploid yeast strain, and the ~12-Mb genome of CAT-1, when compared with the reference S228c genome, contains ~36,000 homozygous and ~30,000 heterozygous single nucleotide polymorphisms, exhibiting an uneven distribution among chromosomes due to large genomic regions of loss of heterozygosity (LOH). In total, 58 % of the 6,652 predicted protein-coding genes of the CAT-1 genome constitute different alleles when compared with the genes present in the reference S288c genome. The CAT-1 genome contains a reduced number of transposable elements, as well as several gene deletions and duplications, especially at telomeric regions, some correlated with several of the physiological characteristics of this industrial fuel-ethanol strain. Phylogenetic analyses revealed that some genes were likely associated with traits important for bioethanol production. Identifying and characterizing the allelic variations controlling traits relevant to industrial fermentation should provide the basis for a forward genetics approach for developing better fermenting yeast strains.     

Phylogenetic classification of Escherichia coli O157:H7 strains of human and bovine origin using a novel set of nucleotide polymorphisms

Background

Cattle are a reservoir of Shiga toxin-producing Escherichia coli O157:H7 (STEC O157), and are known to harbor subtypes not typically found in clinically ill humans. Consequently, nucleotide polymorphisms previously discovered via strains originating from human outbreaks may be restricted in their ability to distinguish STEC O157 genetic subtypes present in cattle. The objectives of this study were firstly to identify nucleotide polymorphisms in a diverse sampling of human and bovine STEC O157 strains, secondly to classify strains of either bovine or human origin by polymorphism-derived genotypes, and finally to compare the genotype diversity with pulsed-field gel electrophoresis (PFGE), a method currently used for assessing STEC O157 diversity.

Results

High-throughput 454 sequencing of pooled STEC O157 strain DNAs from human clinical cases (n = 91) and cattle (n = 102) identified 16,218 putative polymorphisms. From those, 178 were selected primarily within genomic regions conserved across E. coli serotypes and genotyped in 261 STEC O157 strains. Forty-two unique genotypes were observed that are tagged by a minimal set of 32 polymorphisms. Phylogenetic trees of the genotypes are divided into clades that represent strains of cattle origin, or cattle and human origin. Although PFGE diversity surpassed genotype diversity overall, ten PFGE patterns each occurred with multiple strains having different genotypes.

Conclusions

Deep sequencing of pooled STEC O157 DNAs proved highly effective in polymorphism discovery. A polymorphism set has been identified that characterizes genetic diversity within STEC O157 strains of bovine origin, and a subset observed in human strains. The set may complement current techniques used to classify strains implicated in disease outbreaks.     

Copy number variation and variant discovery in Bullmastiff dogs

Identification of genomic variants within dogs is important for understanding genetic factors contributing to breed diversity and phenotypic traits. This study aimed to identify sources of variation in the Bullmastiff using high‐density signal intensity and whole‐genome sequence data. Close to 3000 copy number variants (CNVs) were identified in Bullmastiff dogs using Canine HD BeadChip data. When CNVs were collated, 82 CNV regions (CNVRs) were detected, 50% in transcribed regions encompassing 432 genes. Fifty of the CNVRs detected have not been reported in other breeds and represent potential breed‐specific variants. A proportion of the CNVR variants with predicted modifying effects on gene pathways may contribute to breed traits. Approximately 5 million putative variants per dog, inclusive of single nucleotide polymorphisms (SNPs), multi‐nucleotide polymorphisms (MNPs) and insertion and deletions (INDELs), were identified from DNA sequence data on a small number of animals. Identification of genetic variants in the Bullmastiff highlights sources of variation in the breed and molecular markers that will assist in future trait and disease investigations in dogs.     

Investigation of transferability of BovineSNP50 BeadChip from cattle to water buffalo for genome wide association study

Cattle and water buffalo belong to the same subfamily Bovinae and share chromosome banding and gene order homology. In this study, we used genome-wide Illumina BovineSNP50 BeadChip to analyze 91 DNA samples from three breeds of water buffalo (Nili-Ravi, Murrah and their crossbred with local GuangXi buffalos in China), to demonstrate the genetic divergence between cattle and water buffalo through a large single nucleotide polymorphism (SNP) transferability study at the whole genome level, and performed association analysis of functional traits in water buffalo as well. A total of 40,766 (75.5 %) bovine SNPs were found in the water buffalo genome, but 49,936 (92.5 %) were with only one allele, and finally 935 were identified to be polymorphic and useful for association analysis in water buffalo. Therefore, the genome sequences of water buffalo and cattle shared a high level of homology but the polymorphic status of the bovine SNPs varied between these two species. The different patterns of mutations between species may associate with their phenotypic divergence due to genome evolution. Among 935 bovine SNPs, we identified a total of 9 and 7 SNPs significantly associated to fertility and milk production traits in water buffalo, respectively. However, more works in larger sample size are needed in future to verify these candidate SNPs for water buffalo.     

SNP and haplotype analysis reveal IGF2 variants associated with growth traits in Chinese Qinchuan cattle

Insulin-like growth factor 2 (IGF2) is a potent cell growth and differentiation factor and is implicated in mammals’ growth and development. The objective of this study was to evaluate the effects of the mutations in the bovine IGF2 with growth traits in Chinese Qinchuan cattle. Four single nucleotide polymorphisms (SNPs) were detected of the bovine IGF2 by DNA pool sequencing and forced polymerase chain reaction–restriction fragment length polymorphism (forced PCR–RFLP) methods. We also investigated haplotype structure and linkage disequilibrium (LD) coefficients for four SNPs in 817 individuals representing two main cattle breeds from China. The result of haplotype analysis showed eight different haplotypes and 27 combined genotypes within the study population. The statistical analyses indicated that the four SNPs, combined genotypes and haplotypes are associated with the withers height, body length, chest breadth, chest depth and body weight in Qinchuan cattle population (P < 0.05 or <0.01). The mutant-type variants and mutant haplotype (Hap 8: ATGG; likely to be the beneficial QTN allele) was superior for growth traits; the heterozygote diplotype was associated with higher growth traits compared to wild-type homozygote. Our results provide evidence that polymorphisms in the IGF2 gene are associated with growth traits, and may be used for marker-assisted selection in beef cattle breeding program.