Animal genetics最新文献

Variants in RPGRIP1 and MAP9, termed RPGRIP1ins44 and MAP9del respectively, are both associated with a form of canine progressive retinal atrophy referred to as RPGRIP1-CRD and have both been demonstrated to modify the development and progression of this disease. In the current study both variants were genotyped in at least 50 dogs of 132 diverse breeds and the data reveal that both segregate in multiple breeds. Individually, each variant is common within largely non-overlapping subsets of breed, and there is a negative correlation between their frequencies within breeds that segregate both variants. The frequency of both variants exceeds 0.05 in a single breed only, the Miniature Longhaired Dachshund. These data indicate that both variants are likely to be ancient and predate the development and genetic isolation of modern dog breeds. That both variants are present individually at high frequency in multiple breeds is consistent with the hypothesis that homozygosity of either variant alone is not associated with a clinically relevant phenotype, whereas the negative correlation between the two variants is consistent with the application of selective pressure, from dog breeders, against homozygosity at both loci, probably due to the more severe phenotype associated with homozygosity at both loci.

Genetic improvement of complex traits in animal and plant breeding depends on the efficient and accurate estimation of breeding values. Deep learning methods have been shown to be not superior over traditional genomic selection (GS) methods, partially due to the degradation problem (i.e. with the increase of the model depth, the performance of the deeper model deteriorates). Since the deep learning method residual network (ResNet) is designed to solve gradient degradation, we examined its performance and factors related to its prediction accuracy in GS. Here we compared the prediction accuracy of conventional genomic best linear unbiased prediction, Bayesian methods (BayesA, BayesB, BayesC, and Bayesian Lasso), and two deep learning methods, convolutional neural network and ResNet, on three datasets (wheat, simulated and real pig data). ResNet outperformed other methods in both Pearson's correlation coefficient (PCC) and mean squared error (MSE) on the wheat and simulated data. For the pig backfat depth trait, ResNet still had the lowest MSE, whereas Bayesian Lasso had the highest PCC. We further clustered the pig data into four groups and, on one separated group, ResNet had the highest prediction accuracy (both PCC and MSE). Transfer learning was adopted and capable of enhancing the performance of both convolutional neural network and ResNet. Taken together, our findings indicate that ResNet could improve GS prediction accuracy, affected potentially by factors such as the genetic architecture of complex traits, data volume, and heterogeneity.

Split paw pad disease is a scarcely defined phenotype characterized by skin lesions on the paw pads of dogs. We studied a family of German Shepherd dogs, in which four dogs developed intermittent paw pad lesions and lameness. The paw pads of two of the affected dogs were biopsied and demonstrated cleft formation in the stratum spinosum and stratum corneum, the outermost layers of the epidermis. Whole genome sequencing data from an affected dog revealed a private heterozygous 18 bp in frame deletion in the KRT5 gene. The deletion NM_001346035.1:c.988_1005del or NP_001332964.1:p.(Asn330_Asp335del) is predicted to lead to a loss of six amino acids in the L12 linker domain of the encoded keratin 5. KRT5 variants in human patients lead to various subtypes of epidermolysis bullosa simplex (EBS). Localized EBS is the mildest of the KRT5-related human diseases and may be caused by variants affecting the L12 linker domain of keratin 5. We therefore think that the detected KRT5 deletion in dogs represents a candidate causal variant for the observed skin lesions in dogs. However, while the clinical phenotype of KRT5-mutant dogs of this study closely resembles human patients with localized EBS, there are differences in the histopathology. EBS is defined by cleft formation within the basal layer of the epidermis while the cleft formation in the dogs described herein occurred in the outermost layers, a hallmark of split paw pad disease. Our study provides a basis for further studies into the exact relation of split paw pad disease and EBS.

Kashmir cattle, which were kept by local pastoralists for centuries, are exceptionally resilient and adaptive to harsh environments. Despite its significance, the genomic characteristics of this cattle breed remain elusive. This study utilized whole genome sequences of Kashmir cattle (n = 20; newly sequenced) alongside published whole genomes of 32 distinct breeds and seven core cattle populations (n = 135). The analysis identified ~25.87 million biallelic single nucleotide polymorphisms in Kashmir cattle, predominantly in intergenic and intron regions. Population structure analyses revealed distinct clustering patterns of Kashmir cattle with proximity to the South Asian, African and Chinese indicine cattle populations. Genetic diversity analysis of Kashmir cattle demonstrated lower inbreeding and greater nucleotide diversity than analyzed global breeds. Homozygosity runs indicated less consanguineous mating in Kashmir cattle compared with European taurine breeds. Furthermore, six selection sweep detection methods were used within Kashmir cattle and other cattle populations to identify genes associated with vital traits, including immunity (BOLA-DQA5, BOLA-DQB, TNFAIP8L, FCRL4, AOAH, HIF1AN, FBXL3, MPEG1, CDC40, etc.), reproduction (GOLGA4, BRWD1, OSBP2, LEO1 ADCY5, etc.), growth (ADPRHL1, NRG2, TCF12, TMOD4, GBP4, IGF2, RSPO3, SCD, etc.), milk composition (MRPS30 and CSF1) and high-altitude adaptation (EDNRA, ITPR2, AGBL4 and SCG3). These findings provide essential genetic insights into the characteristics and establish the foundation for the scientific conservation and utilization of Kashmir cattle breed.

The Yarkand hare (Lepus yarkandensis) inhabits arid desert areas and is endemic to China. It has evolved various adaptations to survive in hot arid environments, including stress responses, the ability to maintain water homeostasis and heat tolerance. Here, we performed a selective sweep analysis to identify the candidate genes for adaptation to hot arid environments in the Yarkand hare. A total of 397 237 single-nucleotide polymorphisms were obtained from 80 Yarkand hares, which inhabit hot arid environments, and 36 Tolai hares (Lepus tolai), which inhabit environments with a mild climate, via specific-locus amplified fragment sequencing. We identified several candidate genes that were associated with the heat stress response (HSPE1), oxidative stress response (SLC23A and GLRX2), immune response (IL1R1 and IRG1), central nervous system development (FGF13, THOC2, FMR1 and MECP2) and regulation of water homeostasis (CDK1) according to fixation index values and θ_π ratios in the selective sweep analysis, and six of these genes (GLRX2, IRG1, FGF13, FMR1, MECP2 and CDK1) are newly discovered genes. To the best of our knowledge, this is the first study to identify candidate genes for adaptation to hot arid environments in the Yarkand hare. The results of this study enhance our understanding of the adaptation of the Yarkand hare to hot arid environments and will aid future studies aiming to functionally verify these candidate genes.

Cats with a distinctive white hair pattern of unknown molecular cause have been discovered in the Finnish domestic cat population. Based on the unique appearance of these cats, we have named this phenotype salmiak (“salty licorice”). The use of a commercially available panel test to genotype four salmiak-colored cats revealed the absence of all known variants associated with white-haired phenotypic loci: full White (W), Spotting (W^s) and the Birman white Gloves associated (w^g) allele of the KIT proto-oncogene (KIT) gene. Whole-genome sequencing on two salmiak-colored cats was conducted to search for candidate causal variants in the KIT gene. Despite a lack of coding variants, visual inspection of the short read alignments revealed a large ~95 kb deletion located ~65 kb downstream of the KIT gene in the salmiak cats. Additional PCR genotyping of 180 domestic cats and three salmiak-colored cats confirmed the homozygous derived variant genotype fully concordant with the salmiak phenotype. We suggest the newly identified variant be designated as w^sal for “w salmiak”.

The conservation of animal genetic resources refers to measures taken to prevent the loss of genetic diversity in livestock populations, including the protection of breeds from extinction. Creole cattle populations have suffered a drastic reduction in recent decades owing to absorbent crosses or replacement with commercial breeds of European or Indian origin. Genetic characterization can serve as a source of information for conservation strategies to maintain genetic variation. The objective of this work was to evaluate the levels of inbreeding and kinship through the use of genomic information. A total of 903 DNAs from 13 cattle populations from Argentina, Bolivia and Uruguay were genotyped using an SNP panel of 48 K. Also, a dataset of 76 K SNPs from Peruvian Creole was included. Two inbreeding indices (F_ROH and Fhat2) and kinship relationships were calculated. In addition, effective population size (N_e), linkage disequilibrium, population composition and phylogenetic relationships were estimated. In Creole cattle, F_ROH ranged from 0.14 to 0.03, and Fhat2 was close to zero. The inferred N_e trends exhibited a decline toward the present for all populations, whereas Creole cattle presented a lower magnitude of N_e than foreign breeds. Cluster analysis clearly differentiated the taurine and Zebu components (K2) and showed that Bolivian Creole cattle presented Zebu gene introgression. Despite the population reduction, Creole populations did not present extreme values of consanguinity and kinship and maintain high levels of genetic diversity. The information obtained in this work may be useful for planning conservation programmes for these valuable local animal genetic resources.

Neuronal ceroid lipofuscinoses (NCL) are among the most prevalent neurodegenerative disorders of early life in humans. Disease-causing variants have been described for 13 different NCL genes. In this study, a refined pathological characterization of a female rabbit with progressive neurological signs reminiscent of NCL was performed. Cytoplasmic pigment present in neurons was weakly positive with Sudan black B and autofluorescent. Immunohistology revealed astrogliosis, microgliosis and axonal degeneration. During the subsequent genetic investigation, the genome of the affected rabbit was sequenced and examined for private variants in NCL candidate genes. The analysis revealed a homozygous ~10.7 kb genomic duplication on chromosome 15 comprising parts of the MFSD8 gene, NC_013683.1:g.103,727,963_103,738,667dup. The duplication harbors two internal protein coding exons and is predicted to introduce a premature stop codon into the transcript, truncating ~50% of the wild-type MFSD8 open reading frame encoding the major facilitator superfamily domain containing protein 8, XP_002717309.2:p.(Glu235Leufs*23). Biallelic loss-of-function variants in MFSD8 have been described to cause NCL7 in human patients, dogs and a single cat. The available clinical and pathological data, together with current knowledge about MFSD8 variants and their functional impact in other species, point to the MFSD8 duplication as a likely causative defect for the observed phenotype in the affected rabbit.

Using seven indicator traits, we investigated the genetic basis of bull fertility and predicted gene interactions from SNP associations. We used percent normal sperm as the key phenotype for the association weight matrix–partial correlation information theory (AWM-PCIT) approach. Beyond a simple list of candidate genes, AWM-PCIT predicts significant gene interactions and associations for the selected traits. These interactions formed a network of 537 genes: 38 genes were transcription cofactors, and 41 genes were transcription factors. The network displayed two distinct clusters, one with 294 genes and another with 243 genes. The network is enriched in fertility-associated pathways: steroid biosynthesis, p53 signalling, and the pentose phosphate pathway. Enrichment analysis also highlighted gene ontology terms associated with ‘regulation of neurotransmitter secretion’ and ‘chromatin formation’. Our network recapitulates some genes previously implicated in another network built with lower-density genotypes. Sequence-level data also highlights additional candidate genes relevant to bull fertility, such as FOXO4, FOXP3, GATA1, CYP27B1, and EBP. A trio of regulatory genes—KDM5C, LRRK2, and PME—was deemed core to the network because of their overarching connections. This trio probably influences bull fertility through their interaction with genes, both known and unknown as to their role in male fertility. Future studies may target the trio and their target genes to enrich our understanding of male fertility further.