The genetics of brown coat color and white spotting in domestic yaks (Bos grunniens)

Domestic yaks (Bos grunniens) exhibit two major coat color variations: a brown vs. wild‐type black pigmentation and a white spotting vs. wild‐type solid color pattern. The genetic basis for these variations in color and distribution remains largely unknown and may be complicated by a breeding history involving hybridization between yaks and cattle. Here, we investigated 92 domestic yaks from China using a candidate gene approach. Sequence variations in MC1R, PMEL and TYRP1 were surveyed in brown yaks; TYRP1 was unassociated with the coloration and excluded. Recessive mutations from MC1R, or p.Gln34*, p.Met73Leu and possibly p.Arg142Pro, are reported in bovids for the first time and accounted for approximately 40% of the brown yaks in this study. The remaining 60% of brown individuals correlated with a cattle‐derived deletion mutation from PMEL (p.Leu18del) in a dominant manner. Degrees of white spotting found in yaks vary from color sidedness and white face, to completely white. After examining the candidate gene KIT, we suggest that color‐sided and all‐white yaks are caused by the serial translations of KIT (Cs₆ or Cs₂₉) as reported for cattle. The white‐faced phenotype in yaks is associated with the KIT haplotype S^wf. All KIT mutations underlying the serial phenotypes of white spotting in yaks are identical to those in cattle, indicating that cattle are the likely source of white spotting in yaks. Our results reveal the complex genetic origins of domestic yak coat color as either native in yaks through evolution and domestication or as introduced from cattle through interspecific hybridization.     

Genomic regions influencing coat color saturation and facial markings in Fleckvieh cattle

Genomic regions associated with coat color and pigmented areas of the head were identified for Fleckvieh (dual‐purpose Simmental), a red‐spotted and white‐headed cattle breed. Coat color was measured with a chromameter, implementing the CIELAB color space and resulting in numerical representation of lightness, color intensity, red/green and blue/yellow color components, rather than subjective classification. Single marker regression analyses with fixed effects of the sex and barn were applied, and significant regions were determined with the local false discovery rate methodology. The PMEL and ERBB3 genes on chromosome 5 were in the most significant region for the color measurements. In addition to the blue/yellow color component and color intensity, the AP3B2 gene on chromosome 21 was identified. Its function was confirmed for similar traits in a range of model species. The KIT gene on chromosome 6 was found to be strongly associated with the inhibition of circum‐ocular pigmentation and pigmented spots on the cheek.     

Temporal variation in coat colour (genotypes) supports major changes in the Nordic cattle population after Iron Age

Variation in coat colour genotypes of archaeological cattle samples from Finland was studied by sequencing 69 base pairs of the extension locus (melanocortin 1‐receptor, MC1R) targeting both a transition and a deletion defining the three main alleles, such as dominant black (E^D), wild type (E⁺) and recessive red (e). The 69‐bp MC1R sequence was successfully analysed from 23 ancient (1000–1800 AD) samples. All three main alleles and genotype combinations were detected with allele frequencies of 0.26, 0.17 and 0.57 for E^D, E⁺ and e respectively. Recessive red and dominant black alleles were detected in both sexes. According to the best of our knowledge, this is the first ancient DNA study defining all three main MC1R alleles. Observed MC1R alleles are in agreement with calculated phenotype frequencies from historical sources. The division of ancient Finnish cattle population into modern Finnish breeds with settled colours was dated to the 20th century. From the existing genotyped populations in Europe (43 breeds, n = 2360), the closest match to ancient MC1R genotype frequencies was with the Norwegian native multicoloured breeds. In combined published genotype data of ancient (n = 147) and genotypes and phenotypes of modern Nordic cattle (n = 738), MC1R allele frequencies showed temporal changes similar to neutral mitochondrial DNA and Y‐chromosomal haplotypes analysed earlier. All three markers indicate major change in genotypes in Nordic cattle from the Late Iron Age to the Medieval period followed by slower change through the historical periods until the present.     

Whole genome sequencing reveals a novel deletion variant in the KIT gene in horses with white spotted coat colour phenotypes

White spotting phenotypes in horses can range in severity from the common white markings up to completely white horses. EDNRB, KIT, MITF, PAX3 and TRPM1 represent known candidate genes for such phenotypes in horses. For the present study, we re‐investigated a large horse family segregating a variable white spotting phenotype, for which conventional Sanger sequencing of the candidate genes’ individual exons had failed to reveal the causative variant. We obtained whole genome sequence data from an affected horse and specifically searched for structural variants in the known candidate genes. This analysis revealed a heterozygous ~1.9‐kb deletion spanning exons 10–13 of the KIT gene (chr3:77,740,239_77,742,136del1898insTATAT). In continuity with previously named equine KIT variants we propose to designate the newly identified deletion variant W22. We had access to 21 horses carrying the W22 allele. Four of them were compound heterozygous W20/W22 and had a completely white phenotype. Our data suggest that W22 represents a true null allele of the KIT gene, whereas the previously identified W20 leads to a partial loss of function. These findings will enable more precise genetic testing for depigmentation phenotypes in horses.     

Molecular variation in pigmentation genes contributing to coat colour in native Korean Hanwoo cattle

Pigmentation genes such as TYR (tyrosinase), TYRP1 (tyrosinase-related protein 1), DCT (previously TYRP2, or tyrosinase-related protein 2), ASIP (agouti) and MC1R (melanocortin receptor 1) play a major role in cattle coat colour. To understand the genotypic profile underlying coat colour in native Korean Hanwoo cattle and Angus black cattle, portions of the above-mentioned genes were amplified. Sequence analysis revealed variation in the TYRP1 (exon 5) and MC1R genes. Restriction enzyme analysis of these two genes could distinguish between different colours of Hanwoo cattle. Quantitative estimates of melanin and eumelanin in hair from three different-coloured Hanwoo phenotypes and Angus black showed significant differences at the breed and phenotypic levels. Finally, sequence variants in MC1R were associated with total melanin and eumelanin in breeds as well as in Hanwoo phenotypes.     

A novel MLPH variant in dogs with coat colour dilution

Coat colour dilution may be the result of altered melanosome transport in melanocytes. Loss‐of‐function variants in the melanophilin gene (MLPH) cause a recessively inherited form of coat colour dilution in many mammalian and avian species including the dog. MLPH corresponds to the D locus in many domestic animals, and recessive alleles at this locus are frequently denoted with d. In this study, we investigated dilute coloured Chow Chows whose coat colour could not be explained by their genotype at the previously known MLPH:c.–22G>A variant. Whole genome sequencing of such a dilute Chow Chow revealed another variant in the MLPH gene: MLPH:c.705G>C. We propose to designate the corresponding mutant alleles at these two variants d¹ and d². We performed an association study in a cohort of 15 dilute and 28 non‐dilute Chow Chows. The dilute dogs were all either compound heterozygous d¹/d² or homozygous d²/d², whereas the non‐dilute dogs carried at least one wildtype allele D. The d² allele did not occur in 417 dogs from diverse other breeds. However, when we genotyped a Sloughi family, in which a dilute coloured puppy had been born out of non‐dilute parents, we again observed perfect co‐segregation of the newly discovered d² allele with coat colour dilution. Finally, we identified a blue Thai Ridgeback with the d¹/d² genotype. Thus, our data identify the MLPH:c.705G>C as a variant explaining a second canine dilution allele. Although relatively rare overall, this d² allele is segregating in at least three dog breeds, Chow Chows, Sloughis and Thai Ridgebacks.     

Analysis of ROH patterns in the Noriker horse breed reveals signatures of selection for coat color and body size

Overlapping runs of homozygosity (ROH islands) shared by the majority of a population are hypothesized to be the result of selection around a target locus. In this study we investigated the impact of selection for coat color within the Noriker horse on autozygosity and ROH patterns. We analyzed overlapping homozygous regions (ROH islands) for gene content in fragments shared by more than 50% of horses. Long‐term assortative mating of chestnut horses and the small effective population size of leopard spotted and tobiano horses resulted in higher mean genome‐wide ROH coverage (S_ROH) within the range of 237.4–284.2 Mb, whereas for bay, black and roan horses, where rotation mating is commonly applied, lower autozygosity (S_ROH from 176.4–180.0 Mb) was determined. We identified seven common ROH islands considering all Noriker horses from our dataset. Specific islands were documented for chestnut, leopard spotted, roan and bay horses. The ROH islands contained, among others, genes associated with body size (ZFAT, LASP1 and LCORL/NCAPG), coat color (MC1R in chestnut and the factor PATN1 in leopard spotted horses) and morphogenesis (HOXB cluster in all color strains except leopard spotted horses). This study demonstrates that within a closed population sharing the same founders and ancestors, selection on a single phenotypic trait, in this case coat color, can result in genetic fragmentation affecting levels of autozygosity and distribution of ROH islands and enclosed gene content.     

Reverse genetic screen for loss‐of‐function mutations uncovers a frameshifting deletion in the melanophilin gene accountable for a distinctive coat color in Belgian Blue cattle

In the course of a reverse genetic screen in the Belgian Blue cattle breed, we uncovered a 10‐bp deletion (c.87_96del) in the first coding exon of the melanophilin gene (MLPH), which introduces a premature stop codon (p.Glu32Aspfs*1) in the same exon, truncating 94% of the protein. Recessive damaging mutations in the MLPH gene are well known to cause skin, hair, coat or plumage color dilution phenotypes in numerous species, including human, mice, dog, cat, mink, rabbit, chicken and quail. Large‐scale array genotyping undertaken to identify p.Glu32Aspfs*1 homozygous mutant animals revealed a mutation frequency of 5% in the breed and allowed for the identification of 10 homozygous mutants. As expression of a colored coat requires at least one wild‐type allele at the co‐dominant Roan locus encoded by the KIT ligand gene (KITLG), homozygous mutants for p.Ala227Asp corresponding with the missense mutation were excluded. The six remaining colored calves displayed a distinctive dilution phenotype as anticipated. This new coat color was named ‘cool gray’. It is the first damaging mutation in the MLPH gene described in cattle and extends the already long list of species with diluted color due to recessive mutations in MLPH and broadens the color palette of gray in this breed.     

Fine‐mapping the POLL locus in Brahman cattle yields the diagnostic marker CSAFG29

The POLL locus has been mapped to the centromeric region of bovine chromosome 1 (BTA1) in both taurine breeds and taurine–indicine crosses in an interval of approximately 1 Mb. It has not yet been mapped in pure‐bred zebu cattle. Despite several efforts, neither causative mutations in candidate genes nor a singular diagnostic DNA marker has been identified. In this study, we genotyped a total of 68 Brahman cattle and 20 Hereford cattle informative for the POLL locus for 33 DNA microsatellites, 16 of which we identified de novo from the bovine genome sequence, mapping the POLL locus to the region of the genes IFNAR2 and SYNJ1. The 303‐bp allele of the new microsatellite, CSAFG29, showed strong association with the POLL allele. We then genotyped 855 Brahman cattle for CSAFG29 and confirmed the association between the 303‐bp allele and POLL. To determine whether the same association was found in taurine breeds, we genotyped 334 animals of the Angus, Hereford and Limousin breeds and 376 animals of the Brangus, Droughtmaster and Santa Gertrudis composite taurine–zebu breeds. The association between the 303‐bp allele and POLL was confirmed in these breeds; however, an additional allele (305 bp) was also associated but not fully predictive of POLL. Across the data, CSAFG29 was in sufficient linkage disequilibrium to the POLL allele in Australian Brahman cattle that it could potentially be used as a diagnostic marker in that breed, but this may not be the case in other breeds. Further, we provide confirmatory evidence that the scur phenotype generally occurs in animals that are heterozygous for the POLL allele.     

Exploring Symbiodinium diversity and host specificity in Acropora corals from geographical extremes of Western Australia with 454 amplicon pyrosequencing

Scleractinian corals have demonstrated the ability to shuffle their endosymbiotic dinoflagellate communities (genus Symbiodinium) during periods of acute environmental stress. This has been proposed as a mechanism of acclimation, which would be increased by a diverse and flexible association with Symbiodinium. Conventional molecular techniques used to evaluate Symbiodinium diversity are unable to identify genetic lineages present at background levels below 10%. Next generation sequencing (NGS) offers a solution to this problem and can resolve microorganism diversity at much finer scales. Here we apply NGS to evaluate Symbiodinium diversity and host specificity in Acropora corals from contrasting regions of Western Australia. The application of 454 pyrosequencing allowed for detection of Symbiodinium operational taxonomic units (OTUs) occurring at frequencies as low as 0.001%, offering a 10 000‐fold increase in sensitivity compared to traditional methods. All coral species from both regions were overwhelmingly dominated by a single clade C OTU (accounting for 98% of all recovered sequences). Only 8.5% of colonies associated with multiple clades (clades C and D, or C and G), suggesting a high level of symbiont specificity in Acropora assemblages in Western Australia. While only 40% of the OTUs were shared between regions, the dominance of a single OTU resulted in no significant difference in Symbiodinium community structure, demonstrating that the coral‐algal symbiosis can remain stable across more than 15° of latitude and a range of sea surface temperature profiles. This study validates the use of NGS platforms as tools for providing fine‐scale estimates of Symbiodinium diversity and can offer critical insight into the flexibility of the coral‐algal symbiosis.     

Development and characterization of EST‐SSR markers in an East Asian temperate plant genus Diabelia (Caprifoliaceae)

A set of expressed sequence tag (EST) simple sequence repeat (SSR) markers were developed and characterized using next‐generation sequencing technology for the genus Diabelia (Caprifoliaceae). De novo assembly of RNA‐seq reads resulted in 58 669 contigs with the N50 length of 1211 bp. A total of 2746 contigs were identified to harbor SSR motifs, of which 48 primer pairs were designed and 11 were shown to be polymorphic across three morphospecies of Diabelia. When evaluated with 30 individuals, the number of alleles per locus ranged from 2 to 11 and the expected heterozygosity varied from 0.399 to 0.873, respectively. Distance‐based clustering indicated that the EST‐SSR markers can provide sufficient power to distinguish the three species (or populations). These markers will be useful for evaluating the range‐wide genetic diversity of each species and examining genetic divergence and gene flow between the three species.     

The de novo genome assembly and annotation of a female domestic dromedary of North African origin

The single‐humped dromedary (Camelus dromedarius) is the most numerous and widespread of domestic camel species and is a significant source of meat, milk, wool, transportation and sport for millions of people. Dromedaries are particularly well adapted to hot, desert conditions and harbour a variety of biological and physiological characteristics with evolutionary, economic and medical importance. To understand the genetic basis of these traits, an extensive resource of genomic variation is required. In this study, we assembled at 65× coverage, a 2.06 Gb draft genome of a female dromedary whose ancestry can be traced to an isolated population from the Canary Islands. We annotated 21 167 protein‐coding genes and estimated ~33.7% of the genome to be repetitive. A comparison with the recently published draft genome of an Arabian dromedary resulted in 1.91 Gb of aligned sequence with a divergence of 0.095%. An evaluation of our genome with the reference revealed that our assembly contains more error‐free bases (91.2%) and fewer scaffolding errors. We identified ~1.4 million single‐nucleotide polymorphisms with a mean density of 0.71 × 10^?3 per base. An analysis of demographic history indicated that changes in effective population size corresponded with recent glacial epochs. Our de novo assembly provides a useful resource of genomic variation for future studies of the camel's adaptations to arid environments and economically important traits. Furthermore, these results suggest that draft genome assemblies constructed with only two differently sized sequencing libraries can be comparable to those sequenced using additional library sizes, highlighting that additional resources might be better placed in technologies alternative to short‐read sequencing to physically anchor scaffolds to genome maps.     

Development of diagnostic microsatellite markers from whole‐genome sequences of Ammodramus sparrows for assessing admixture in a hybrid zone

Studies of hybridization and introgression and, in particular, the identification of admixed individuals in natural populations benefit from the use of diagnostic genetic markers that reliably differentiate pure species from each other and their hybrid forms. Such diagnostic markers are often infrequent in the genomes of closely related species, and genomewide data facilitate their discovery. We used whole‐genome data from Illumina HiSeqS2000 sequencing of two recently diverged (600,000 years) and hybridizing, avian, sister species, the Saltmarsh (Ammodramus caudacutus) and Nelson's (A. nelsoni) Sparrow, to develop a suite of diagnostic markers for high‐resolution identification of pure and admixed individuals. We compared the microsatellite repeat regions identified in the genomes of the two species and selected a subset of 37 loci that differed between the species in repeat number. We screened these loci on 12 pure individuals of each species and report on the 34 that successfully amplified. From these, we developed a panel of the 12 most diagnostic loci, which we evaluated on 96 individuals, including individuals from both allopatric populations and sympatric individuals from the hybrid zone. Using simulations, we evaluated the power of the marker panel for accurate assignments of individuals to their appropriate pure species and hybrid genotypic classes (F1, F2, and backcrosses). The markers proved highly informative for species discrimination and had high accuracy for classifying admixed individuals into their genotypic classes. These markers will aid future investigations of introgressive hybridization in this system and aid conservation efforts aimed at monitoring and preserving pure species. Our approach is transferable to other study systems consisting of closely related and incipient species.     

Diet assessment of the Atlantic Sea Nettle Chrysaora quinquecirrha in Barnegat Bay,New Jersey,using next‐generation sequencing

Next‐generation sequencing (NGS) methodologies have proven useful in deciphering the food items of generalist predators, but have yet to be applied to gelatinous animal gut and tentacle content. NGS can potentially supplement traditional methods of visual identification. Chrysaora quinquecirrha (Atlantic sea nettle) has progressively become more abundant in Mid‐Atlantic United States’ estuaries including Barnegat Bay (New Jersey), potentially having detrimental effects on both marine organisms and human enterprises. Full characterization of this predator's diet is essential for a comprehensive understanding of its impact on the food web and its management. Here, we tested the efficacy of NGS for prey item determination in the Atlantic sea nettle. We implemented a NGS ‘shotgun’ approach to randomly sequence DNA fragments isolated from gut lavages and gastric pouch/tentacle picks of eight and 84 sea nettles, respectively. These results were verified by visual identification and co‐occurring plankton tows. Over 550 000 contigs were assembled from ~110 million paired‐end reads. Of these, 100 contigs were confidently assigned to 23 different taxa, including soft‐bodied organisms previously undocumented as prey species, including copepods, fish, ctenophores, anemones, amphipods, barnacles, shrimp, polychaete worms, flukes, flatworms, echinoderms, gastropods, bivalves and hemichordates. Our results not only indicate that a ‘shotgun’ NGS approach can supplement visual identification methods, but targeted enrichment of a specific amplicon/gene is not a prerequisite for identifying Atlantic sea nettle prey items.     

Morphological and Molecular Characterization of Fusarium Isolated From Maize in Syria

Maize is the third most important cereal after wheat and barley in Syria. Maize plants are attacked by several Fusarium species causing mainly stalk and ear rot of maize which poses a major impact worldwide. Identification of Fusarium species is important for disease control and for assessment of exposure risk to mycotoxines. To identify Fusarium species attacking maize in Syria, a total of 32 Fusarium isolates were recovered from maize ears collected from four different geographical regions, mainly from Ghouta surrounding Damascus. Fusarium isolates were identified based on morphology and on partial DNA sequencing of the TEF1‐α and rDNA/ITS genes. The majority (26 of 32) of these isolates was identified as F. verticillioides (subdivided into four groups), whereas three isolates turned out to be F. thapsinum, F. equiseti and F. andiyazi. The remaining three isolates were close to F. andiyazi, although further investigation is needed to confirm whether they represent a yet undescribed species. Furthermore, our results showed that sequencing the TEF1‐α gene is much more informative than sequencing of the rDNA/ITS region for Fusarium identification at the species level. PCR analysis showed that only F. verticillioides isolates were potentially fumonisin producers and that only the F. equiseti isolate was potentially trichotecene producer. This is the first report on Fusarium thapsinum, F. equiseti and F. andiyazi attacking maize in Syria.     

Genomewide single nucleotide polymorphism discovery in Atlantic salmon (Salmo salar): validation in wild and farmed American and European populations

A considerable number of single nucleotide polymorphisms (SNPs) are required to elucidate genotype–phenotype associations and determine the molecular basis of important traits. In this work, we carried out de novo SNP discovery accounting for both genome duplication and genetic variation from American and European salmon populations. A total of 9 736 473 nonredundant SNPs were identified across a set of 20 fish by whole‐genome sequencing. After applying six bioinformatic filtering steps, 200 K SNPs were selected to develop an Affymetrix Axiom^® myDesign Custom Array. This array was used to genotype 480 fish representing wild and farmed salmon from Europe, North America and Chile. A total of 159 099 (79.6%) SNPs were validated as high quality based on clustering properties. A total of 151 509 validated SNPs showed a unique position in the genome. When comparing these SNPs against 238 572 markers currently available in two other Atlantic salmon arrays, only 4.6% of the SNP overlapped with the panel developed in this study. This novel high‐density SNP panel will be very useful for the dissection of economically and ecologically relevant traits, enhancing breeding programmes through genomic selection as well as supporting genetic studies in both wild and farmed populations of Atlantic salmon using high‐resolution genomewide information.