A missense variant in the NSDHL gene in a Chihuahua with a congenital cornification disorder resembling inflammatory linear verrucous epidermal nevi

Congenital hemidysplasia with ichthyosiform nevus and limb defects syndrome in humans is a genodermatosis characterized by inflammatory linear verrucous epidermal nevi (ILVEN), often showing a striking lateralization pattern. It is caused by variants in the NSDHL gene encoding a 3β‐hydroxysteroid dehydrogenase involved in the cholesterol biosynthesis pathway. In the present study, we investigated a female Chihuahua, which showed clinical and histological signs of ILVEN. We performed a candidate gene analysis in the affected animal. This analysis revealed a single missense variant in the NSDHL gene in the affected dog (XM_014111859.2:c.700G>A). The variant is predicted to cause a non‐conservative amino acid change from glycine to arginine, XP_013967334.1:p.(Gly234Arg). The mutant allele was absent from WGS data of 594 genetically diverse dogs and eight wolves. Sanger sequencing confirmed that the variant was heterozygous in the affected dog and absent from 22 control Chihuahuas. Based on the knowledge about the functional impact of NSDHL variants in dogs and other species, c.700G>A is probably pathogenic and a convincing candidate causative variant for the observed skin lesions in the affected Chihuahua.     

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.     

A second KRT71 allele in curly coated dogs

Major characteristics of coat variation in dogs can be explained by variants in only a few genes. Until now, only one missense variant in the KRT71 gene, p.Arg151Trp, has been reported to cause curly hair in dogs. However, this variant does not explain the curly coat in all breeds as the mutant ¹⁵¹Trp allele, for example, is absent in Curly Coated Retrievers. We sequenced the genome of a Curly Coated Retriever at 22× coverage and searched for variants in the KRT71 gene. Only one protein‐changing variant was present in a homozygous state in the Curly Coated Retriever and absent or present in a heterozygous state in 221 control dogs from different dog breeds. This variant, NM_001197029.1:c.1266_1273delinsACA, was an indel variant in exon 7 that caused a frameshift and an altered and probably extended C‐terminus of the KRT71 protein NP_001183958.1:p.(Ser422ArgfsTer?). Using Sanger sequencing, we found that the variant was fixed in a cohort of 125 Curly Coated Retrievers and segregating in five of 14 additionally tested breeds with a curly or wavy coat. KRT71 variants cause curly hair in humans, mice, rats, cats and dogs. Specific KRT71 variants were further shown to cause alopecia. Based on this knowledge from other species and the predicted molecular consequence of the newly identified canine KRT71 variant, it is a compelling candidate causing a second curly hair allele in dogs. It might cause a slightly different coat phenotype than the previously published p.Arg151Trp variant and could potentially be associated with follicular dysplasia in dogs.     

Identification of compound heterozygous DNAH11 variants in a Han-Chinese family with primary ciliary dyskinesia

Primary ciliary dyskinesia (PCD) is a group of genetically and clinically heterogeneous disorders with motile cilia dysfunction. It is clinically characterized by oto-sino-pulmonary diseases and subfertility, and half of the patients have situs inversus (Kartagener syndrome). To identify the genetic cause in a Han-Chinese pedigree, whole-exome sequencing was conducted in the 37-year-old proband, and then, Sanger sequencing was performed on available family members. Minigene splicing assay was applied to verify the impact of the splice-site variant. Compound heterozygous variants including a splice-site variant (c.1974-1G>C, rs1359107415) and a missense variant (c.7787G>A, p.(Arg2596Gln), rs780492669), in the dynein axonemal heavy chain 11 gene (DNAH11) were identified and confirmed as the disease-associated variants of this lineage. The minigene expression in vitro revealed that the c.1974-1G>C variant could cause skipping over exon 12, predicted to result in a truncated protein. This discovery may enlarge the DNAH11 variant spectrum of PCD, promote accurate genetic counselling and contribute to PCD diagnosis.     

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.     

Novel variants in the KIT and PAX3 genes in horses with white‐spotted coat colour phenotypes

Variants in the EDNRB, KIT, MITF, PAX3 and TRPM1 genes are known to cause white spotting phenotypes in horses, which can range from the common white markings up to completely white horses. In this study, we investigated these candidate genes in 169 horses with white spotting phenotypes not explained by the previously described variants. We identified a novel missense variant, PAX3:p.Pro32Arg, in Appaloosa horses with a splashed white phenotype in addition to their leopard complex spotting patterns. We also found three novel variants in the KIT gene. The splice site variant c.1346+1G>A occurred in a Swiss Warmblood horse with a pronounced depigmentation phenotype. The missense variant p.Tyr441Cys was present in several part‐bred Arabians with sabino‐like depigmentation phenotypes. Finally, we provide evidence suggesting that the common and widely distributed KIT:p.Arg682His variant has a very subtle white‐increasing effect, which is much less pronounced than the effect of the other described KIT variants. We termed the new KIT variants W18–W20 to provide a simple and unambiguous nomenclature for future genetic testing applications.     

Identification of two TYRP1 loss‐of‐function alleles in Valais Red sheep

The Valais Red sheep breed is a local breed of the Swiss canton Valais. Although the breed is characterised by its brown colour, black animals occasionally occur and the objective of this study was to identify the causative genetic variants responsible for the obvious difference. A GWAS using high‐density SNP data to compare 51 brown and 38 black sheep showed a strong signal on chromosome 2 at the TYRP1 locus. Haplotype analyses revealed three different brown‐associated alleles. The WGS of three sheep revealed four protein‐changing variants within the TYRP1 gene. Three of these variants were associated with the recessively inherited brown coat colour. This includes the known missense variant TYRP1:c.869G>T designated as b^Soay and two novel loss‐of‐function variants. We propose to designate the frame‐shift variant TYRP1:c.86_87delGA as b^VS1 and the nonsense variant TYRP1:c.1066C>T as b^VS2. Interestingly, the b^VS1 allele occurs only in local breeds of Switzerland whereas the b^VS2 allele seems to be more widespread across Europe.     

A novel KRT71 variant in curly‐coated dogs

Curly fur is a common phenotype in many dog breeds, known to result from a missense variant (c.451C>T) in exon 2 of the keratin 71 (KRT71) gene. During screening for this variant across various breeds, we found that Curly Coated Retrievers (CCRs) fixed with the trait did not carry the known variant. By analysis of whole‐genome sequencing data of one CCR we identified a novel genetic cause for curly fur. We found a novel structural variant in exon 7 of the KRT71 gene (c.1266_1273delinsACA) that was predicted to result in a frameshift and stop loss, therefore significantly affecting the structure of the protein, if translated. The variant was also found at lower frequencies in five other breeds, including Lagotto Romagnolo, Bichon Frise, Spanish Water Dog, Chesapeake Bay Retriever and Irish Terrier. One curly‐coated Lagotto carried neither of the two KRT71 variants. These results identify a second variant for curly coat in KRT71 and suggest the existence of additional alleles. This study enables the development of an additional KRT71 gene test for breeders to understand and manage coat types.     

A Triglyceride‐Raising APOA5 Genetic Variant Is Negatively Associated With Obesity and BMI in the Chinese Population

Apolipoprotein A‐V (apo A‐V) exerts a potent triglyceride (TG)‐lowering effect through enhanced intravascular TG‐hydrolysis with increased uptake of TG‐derived free fatty acids into muscle and adipose tissue. Genetic variants in the APOA5 gene were strongly associated with plasma TG concentrations. The aim of this study was to examine whether APOA5 genetic variation was associated with obesity. We genotyped the missense c.553 G>T polymorphism (p.G185C) in the APOA5 gene in 1,085 Chinese (333 obese subjects and 752 nonobese controls). We analyzed the association between the c.553 G>T polymorphism and obesity and related metabolic phenotypes. The T allele at the c.553 G>T polymorphism was associated with higher plasma TG concentrations. Each additional T allele was associated with an increased TG concentration of 53.5 mg/dl (95% confidence interval (CI) 29.6–76.0, P < 0.0001). However, the T allele was associated lower risk of obesity (odds ratio (OR), 0.48; 95% CI 0.32–0.73, P = 0.0004). Each additional copy of the T allele was associated with a BMI decrease of 0.73 kg/m² (95% CI 0.26–1.16, P = 0.002), equivalent to 2.11 kg in a person 1.7 m tall. We may then conclude that the TG‐raising APOA5 genetic variant was associated with a decrease in BMI and reduced risk of obesity in the Chinese population.     

The prion‐related protein (testis‐specific) gene (PRNT) is highly polymorphic in Portuguese sheep

The objective of this study was to search for polymorphisms in the ovine prion‐related protein (testis‐specific) gene (PRNT). Sampling included 567 sheep from eight Portuguese breeds. The PRNT gene‐coding region was analyzed by single‐strand conformation polymorphism and sequencing, allowing the identification of the first ovine PRNT polymorphisms, in codons 6, 38, 43 and 48: c.17C>T (p.Ser6Phe, which disrupts a consensus arginine‐X‐X‐serine/threonine motif); c.112G>C (p.Gly38>Arg); c.129T>C and c.144A>G (synonymous) respectively. Polymorphisms in codons 6, 38 and 48 occur simultaneously in 50.6% of the animals, 38.8% presenting as heterozygous. To study the distribution of the polymorphism in codon 43, a restriction fragment length polymorphism analysis was performed. Polymorphic variant c.129C, identified in 89.8% of the animals with 32.8% presented as heterozygous, was considered the wild genotype in Portuguese sheep. Eight different haplotypes which have comparable distribution in all breeds were identified for the PRNT gene. In conclusion, the PRNT coding region is highly polymorphic in sheep, unlike the prion protein 2 dublet gene (PRND), in which we previously found only one synonymous substitution (c.78G>A), in codon 26. The absence or reduced number of PRND heterozygotes (c.78G>A) was significantly associated with three PRNT haplotypes (17C‐112G‐129T‐144A,17CT‐112GC‐129CT‐144AG and 17T‐112C‐129C‐144G), and the only three animals found homozygous at c.78A had the 17C‐112G‐129C‐144A PRNT haplotype. These results constitute evidence of an association between polymorphic variation in PRND and PRNT genes, as has already been observed for PRND and prion protein gene (PRNP).     

Non‐synonymous SNPs in MC1R gene are associated with the extended black variant in domestic ducks (Anas platyrhynchos)

In this study, we performed a sequence characterization of the duck melanocortin 1 receptor (alpha‐melanocyte stimulating hormone receptor) (MC1R) gene to analyze the relationship between MC1R polymorphism and the extended black variant in domestic ducks based on the extended black (E) and non‐extended black (e⁺) allele hypothesis of the duck MC1R gene. Both c.52G>A and c.376G>A substitutions are highly associated with the duck extended black variant (P < 0.01), but the novel c.52G>A substitution is more of a fit for the allele hypothesis of the duck MC1R gene.     

Substantial differences in bias between single‐digest and double‐digest RAD‐seq libraries: A case study

The trade‐offs of using single‐digest vs. double‐digest restriction site‐associated DNA sequencing (RAD‐seq) protocols have been widely discussed. However, no direct empirical comparisons of the two methods have been conducted. Here, we sampled a single population of Gulf pipefish (Syngnathus scovelli) and genotyped 444 individuals using RAD‐seq. Sixty individuals were subjected to single‐digest RAD‐seq (sdRAD‐seq), and the remaining 384 individuals were genotyped using a double‐digest RAD‐seq (ddRAD‐seq) protocol. We analysed the resulting Illumina sequencing data and compared the two genotyping methods when reads were analysed either together or separately. Coverage statistics, observed heterozygosity, and allele frequencies differed significantly between the two protocols, as did the results of selection components analysis. We also performed an in silico digestion of the Gulf pipefish genome and modelled five major sources of bias: PCR duplicates, polymorphic restriction sites, shearing bias, asymmetric sampling (i.e., genotyping fewer individuals with sdRAD‐seq than with ddRAD‐seq) and higher major allele frequencies. This combination of approaches allowed us to determine that polymorphic restriction sites, an asymmetric sampling scheme, mean allele frequencies and to some extent PCR duplicates all contribute to different estimates of allele frequencies between samples genotyped using sdRAD‐seq versus ddRAD‐seq. Our finding that sdRAD‐seq and ddRAD‐seq can result in different allele frequencies has implications for comparisons across studies and techniques that endeavour to identify genomewide signatures of evolutionary processes in natural populations.     

The CDKN2A/p16INK4a 5′UTR sequence and translational regulation: impact of novel variants predisposing to melanoma

Many variants of uncertain functional significance in cancer susceptibility genes lie in regulatory regions, and clarifying their association with disease risk poses significant challenges. We studied 17 germline variants (nine of which were novel) in the CDKN2A 5′UTR with independent approaches, which included mono and bicistronic reporter assays, Western blot of endogenous protein, and allelic representation after polysomal profiling to investigate their impact on CDKN2A mRNA translation regulation. Two of the novel variants (c.‐27del23, c.‐93‐91delAGG) were classified as causal mutations (score ≥3), along with the c.‐21C>T, c.‐34G>T, and c.‐56G>T, which had already been studied by a subset of assays. The novel c.‐42T>A as well as the previously described c.‐67G>C were classified as potential mutations (score 1 or 2). The remaining variants (c.‐14C>T, c.‐20A>G, c.‐25C>T+c.‐180G>A, c.‐30G>A, c.‐40C>T, c.‐45G>A, c.‐59C>G, c.‐87T>A, c.‐252A>T) were classified as neutral (score 0). In conclusion, we found evidence that nearly half of the variants found in this region had a negative impact on CDKN2A mRNA translation, supporting the hypothesis that 5′UTR can act as a cellular Internal Ribosome Entry Site (IRES) to modulate p16^INK4a translation.     

Interaction of MC1R and PMEL alleles on solid coat colors in Highland cattle

Six solid colors occur in Highland cattle: black, dun, silver dun and red, yellow, and white. These six coat colors are explained by a non‐epistatic interaction of the genotypes at the MC1R and PMEL genes. A three base pair deletion in the PMEL gene leading to the deletion of a leucine from the signal peptide is observed in dilute‐colored Highland cattle (c.50_52delTTC, p.Leu18del). The mutant PMEL allele acts in a semi‐dominant manner. Dun Galloway cattle also have one copy of the deletion allele, and silver dun Galloway cattle have two copies. The presence of two adjacent leucine residues at the site of this deletion is highly conserved in human, horse, mouse and chicken as well as in cattle with undiluted coat colors. Highland and Galloway cattle thus exhibit a similar dose‐dependent dilution effect based on the number of PMEL :c.50_51delTTC alleles, as Charolais cattle with PMEL :c.64G>A alleles. The PMEL :c.64G>A allele was not found in Highland or Galloway cattle.     

Population genomic analyses from low‐coverage RAD‐Seq data: a case study on the non‐model cucurbit bottle gourd

Restriction site‐associated DNA sequencing (RAD‐Seq), a next‐generation sequencing‐based genome ‘complexity reduction’ protocol, has been useful in population genomics in species with a reference genome. However, the application of this protocol to natural populations of genomically underinvestigated species, particularly under low‐to‐medium sequencing depth, has not been well justified. In this study, a Bayesian method was developed for calling genotypes from an F₂ population of bottle gourd [Lagenaria siceraria (Mol.) Standl.] to construct a high‐density genetic map. Low‐depth genome shotgun sequencing allowed the assembly of scaffolds/contigs comprising approximately 50% of the estimated genome, of which 922 were anchored for identifying syntenic regions between species. RAD‐Seq genotyping of a natural population comprising 80 accessions identified 3226 single nuclear polymorphisms (SNPs), based on which two sub‐gene pools were suggested for association with fruit shape. The two sub‐gene pools were moderately differentiated, as reflected by the Hudson's F_ST value of 0.14, and they represent regions on LG7 with strikingly elevated F_ST values. Seven‐fold reduction in heterozygosity and two times increase in LD (r²) were observed in the same region for the round‐fruited sub‐gene pool. Outlier test suggested the locus LX3405 on LG7 to be a candidate site under selection. Comparative genomic analysis revealed that the cucumber genome region syntenic to the high F_ST island on LG7 harbors an ortholog of the tomato fruit shape gene OVATE. Our results point to a bright future of applying RAD‐Seq to population genomic studies for non‐model species even under low‐to‐medium sequencing efforts. The genomic resources provide valuable information for cucurbit genome research.     

Development of an SNP‐based parentage verification panel for lovebirds

The genus Agapornis, or lovebirds, are popular pet parrots worldwide. Currently, breeders are dependent on pedigree records as a selection tool as no molecular parentage verification test is available for any of the nine species. The A. roseicollis reference genome was recently assembled. This was followed by the sequencing of the whole genomes of the parents of the reference genome individual at 30× coverage. The parents’ reads were mapped against the reference genome to identify SNPs. Over 1.6 million SNPs, shared between the parents, were discovered using the Genome Analysis Toolkit pipeline. SNPs were filtered to a panel of 480 SNPs based on Genome Analysis Toolkit parameters. The panel of 480 SNPs was genotyped in a population of 960 lovebirds across seven species. A panel of 262 SNPs was compiled that included SNPs successfully amplified across all species. The 262‐SNP panel was reduced based on the observed heterozygosity (H_O) and minor allele frequency (MAF) values per SNP to include the lowest number of SNPs with the highest exclusion power for parentage verification. Two smaller panels consisting of 195 SNPs with MAF and H_O values >0.1 and 40 SNPs with MAF and H_O values >0.3, were constructed. The panels were verified using 43 families from different species with known relationships to evaluate the exclusion power of each panel. The 195 SNP panel with an average exclusion probability of 99.9% and MAF and H_O values >0.1 was proposed as the routine Agapornis parentage verification panel.