Whole genome sequencing analysis of alpaca suggests TRPV3 as a candidate gene for the suri phenotype.

BACKGROUND: Alpaca is a domestic South American camelid probably arising from the domestication of two wild camelids, the vicugna and the guanaco. Two phenotypes are described for alpaca, known as huacaya and suri. Huacaya fleece is characterized by compact, soft, and highly crimped fibers, while suri fleece is longer, straight, less crimped, and lustrous. The gene variants determining these phenotypes are still unknown, although previous studies suggested a dominant inheritance of the suri. Based on that, the aim of this study was the identification of the gene variants determining alpaca coat phenotypes through whole genome sequencing (WGS) analysis. RESULTS: The sample used includes two test-cross alpaca families, suri × huacaya, which produced two offspring, one with the suri phenotype and one with the huacaya phenotype. The analyzed sample was expanded through the addition of WGS data from six vicugnas and six guanacos; this because we assumed the absence of the gene variants linked to the suri phenotype in these wild species. The analysis of gene variant segregation with the suri phenotype, coupled with the filtering of gene variants present in the wild species, disclosed the presence in all the suri samples of a premature termination codon (PTC) in TRPV3 (transient receptor potential cation channel subfamily V member 3), a gene known to be involved in hair growth and cycling, thermal sensation, cold tolerance and adaptation in several species. Mutations in TRPV3 were previously associated with the alteration of hair structure leading to an impaired formation of the hair canal and the hair shaft in mouse. This PTC in TRPV3, due to a G > T substitution (p.Glu475*), results in a loss of 290 amino acids from the canonical translated protein, plausibly leading to a physiological dysfunction. CONCLUSION: The present results suggest that the suri phenotype may arise from a TRPV3 gene variant which may explain some of the suri features such as its longer hair fibre with lower number of cuticular scales compared to huacaya.

Sex- and APOE-specific genetic risk factors for late-onset Alzheimer's disease: Evidence from gene-gene interaction of longevity-related loci.

Advanced age is the largest risk factor for late-onset Alzheimer's disease (LOAD), a disease in which susceptibility correlates to almost all hallmarks of aging. Shared genetic signatures between LOAD and longevity were frequently hypothesized, likely characterized by distinctive epistatic and pleiotropic interactions. Here, we applied a multidimensional reduction approach to detect gene-gene interactions affecting LOAD in a large dataset of genomic variants harbored by genes in the insulin/IGF1 signaling, DNA repair, and oxidative stress pathways, previously investigated in human longevity. The dataset was generated from a collection of publicly available Genome Wide Association Studies, comprising a total of 2,469 gene variants genotyped in 20,766 subjects of Northwestern European ancestry (11,038 LOAD cases and 9,728 controls). The stratified analysis according to APOE*4 status and sex corroborated evidence that pathways leading to longevity also contribute to LOAD. Among the significantly interacting genes, PTPN1, TXNRD1, and IGF1R were already found enriched in gene-gene interactions affecting survival to old age. Furthermore, interacting variants associated with LOAD in a sex- and APOE-specific way. Indeed, while in APOE*4 female carriers we found several inter-pathway interactions, no significant epistasis was found in APOE*4 negative females; conversely, in males, significant intra- and inter-pathways epistasis emerged according to APOE*4 status. These findings suggest that interactions of risk factors may drive different trajectories of cognitive aging. Beyond helping to disentangle the genetic architecture of LOAD, such knowledge may improve precision in predicting the risk of dementia and enable effective sex- and APOE-stratified preventive and therapeutic interventions for LOAD.

Genome-wide scan for runs of homozygosity in South American Camelids.

BACKGROUND: Alpaca (Vicugna pacos), llama (Lama glama), vicugna (Vicugna vicugna) and guanaco (Lama guanicoe), are the camelid species distributed over the Andean high-altitude grasslands, the Altiplano, and the Patagonian arid steppes. Despite the wide interest on these animals, most of the loci under selection are still unknown. Using whole-genome sequencing (WGS) data we investigated the occurrence and the distribution of Runs Of Homozygosity (ROHs) across the South American Camelids (SACs) genome to identify the genetic relationship between the four species and the potential signatures of selection. RESULTS: A total of 37 WGS samples covering the four species was included in the final analysis. The multi-dimensional scaling approach showed a clear separation between the four species; however, admixture analysis suggested a strong genetic introgression from vicugna and llama to alpaca. Conversely, very low genetic admixture of the guanaco with the other SACs was found. The four species did not show significant differences in the number, length of ROHs (100-500 kb) and genomic inbreeding values. Longer ROHs (> 500 kb) were found almost exclusively in alpaca. Seven overlapping ROHs were shared by alpacas, encompassing nine loci (FGF5, LOC107034918, PRDM8, ANTXR2, LOC102534792, BSN, LOC116284892, DAG1 and RIC8B) while nine overlapping ROHs were found in llama with twenty-five loci annotated (ERC2, FZD9, BAZ1B, BCL7B, LOC116284208, TBL2, MLXIPL, PHF20, TRNAD-AUC, LOC116284365, RBM39, ARFGEF2, DCAF5, EXD2, HSPB11, LRRC42, LDLRAD1, TMEM59, LOC107033213, TCEANC2, LOC102545169, LOC116278408, SMIM15, NDUFAF2 and RCOR1). Four overlapping ROHs, with three annotated loci (DLG1, KAT6B and PDE4D) and three overlapping ROHs, with seven annotated genes (ATP6V1E1, BCL2L13, LOC116276952, BID, KAT6B, LOC116282667 and LOC107034552), were detected for vicugna and guanaco, respectively. CONCLUSIONS: The signatures of selection revealed genomic areas potentially selected for production traits as well as for natural adaptation to harsh environment. Alpaca and llama hint a selection driven by environment as well as by farming purpose while vicugna and guanaco showed selection signals for adaptation to harsh environment. Interesting, signatures of selection on KAT6B gene were identified for both vicugna and guanaco, suggesting a positive effect on wild populations fitness. Such information may be of interest to further ecological and animal production studies.

Phenotypic Correlation between Guard Hair and Down Hair in Chinese Alashan Left Banner White Cashmere Goat: A Preliminary Study.

In cashmere production studies, few trials have considered the guard hair features and their correlation with down fiber attributes. In this preliminary work, early observations on 158 one year old Chinese Alashan Left Banner White Cashmere goats were carried out. The aim was to describe the phenotypic correlation between the guard hair length and other fiber traits. The guard hair length was positively correlated with guard hair diameter and the down fiber length. Negative correlations were found between guard hair length and the coefficient of variation of guard hair diameter, between the guard hair diameter and its coefficient of variation, and between the down fiber diameter and the coefficient of variation of down fiber diameter. The body weight at first combing was not correlated with any of the other traits.

Dog-human translational genomics: state of the art and genomic resources.

Innovative models for medical research are strongly required nowadays. Convincing evidence supports dog as the most suitable spontaneous model for several human genetic diseases. Decades of studies on dog genome allowed the identification of hundreds of mutations causing genetic disorders, many of which are proposed as counterparts responsible for human diseases. Traditionally, the murine model is the most extensively used in human translational research. However, this species shows large physiological differences from humans, and it is kept under a controlled artificial environment. Conversely, canine genetic disorders often show pathophysiological and clinical features highly resembling the human counterpart. In addition, dogs share the same environment with humans; therefore, they are naturally exposed to many risk factors. Thus, different branches of translational medicine aim to study spontaneously occurring diseases in dogs to provide a more reliable model for human disorders. This review offers a comprehensive overview of the knowledge and resources available today for all the researchers involved in the field of dog-human translational medicine. Some of the main successful examples from dog-human translational genomics are reported, such as the canine association studies which helped to identify the causal mutation in the human counterpart. We also illustrated the ongoing projects aiming to create public canine big datasets. Finally, specific online databases are discussed along with several information resources that can speed up clinical translational research.

PAGEANT: personal access to genome and analysis of natural traits.

GWASs have identified numerous genetic variants associated with a wide variety of diseases, yet despite the wide availability of genetic testing the insights that would enhance the interpretability of these results are not widely available to members of the public. As a proof of concept and demonstration of technological feasibility, we developed PAGEANT (Personal Access to Genome & Analysis of Natural Traits), usable through Graphical User Interface or command line-based version, aiming to serve as a protocol and prototype that guides the overarching design of genetic reporting tools. PAGEANT is structured across five core modules, summarized by five Qs: (i) quality assurance of the genetic data; (ii) qualitative assessment of genetic characteristics; (iii) quantitative assessment of health risk susceptibility based on polygenic risk scores and population reference; (iv) query of third-party variant databases (e.g. ClinVAR and PharmGKB) and (v) quick Response code of genetic variants of interest. Literature review was conducted to compare PAGEANT with academic and industry tools. For 2504 genomes made publicly available through the 1000 Genomes Project, we derived their genomic characteristics for a suite of qualitative and quantitative traits. One exemplary trait is susceptibility to COVID-19, based on the most up-to-date scientific findings reported.

Genome-wide expression of the residual lung reacting to experimental Pneumonectomy.

BACKGROUND: Acute or chronic irreversible respiratory failure may occur in patients undergoing pneumonectomy. Aim of this study was to determine transcriptome expression changes after experimental pneumonectomy in swine model. Experimental left pneumonectomy was performed in five pigs under general anaesthesia. Both the resected and the remaining lung, after 60 post-operative completely uneventful days, underwent genome-wide bulk RNA-Sequencing (RNA-Seq). RESULTS: Histological analysis showed dilation of air spaces and rupture of interalveolar septa. In addition, mild inflammation, no fibrosis, radial stretch of the bronchus, strong enlargement of airspaces and thinning of the blood supply were observed. Bioinformatic analyses of bulk RNA-Seq data identified 553 Differentially Expressed Genes (DEGs) at adjusted P-value below 0.001, between pre- and post-pneumonectomy. The top 10 up-regulated DEGs were Edn1, Areg, Havcr2, Gadd45g, Depp1, Cldn4, Atf3, Myc, Gadd45b, Socs3; the top 10 down-regulated DEGs were Obscn, Cdkn2b, ENSSSCG00000015738, Prrt2, Amer1, Flrt3, Efnb2, Tox3, Znf793, Znf365. Leveraging digital cytometry tools, no difference in cellular abundance was found between the two experimental groups, while the analysis of cell type-specific gene expression patterns highlighted a striking predominance of macrophage-specific genes among the DEGs. DAVID-based gene ontology analysis showed a significant enrichment of "Extrinsic apoptotic signaling pathway" (FDR q = 7.60 × 10- 3) and "Response to insulin" (FDR q = 7.60 × 10- 3) genes, along with an enrichment of genes involved as "Negative regulators of DDX58/IFIH1 signaling" (FDR q = 7.50 × 10- 4) found by querying the REACTOME pathway database. Gene network analyses indicated a general dysregulation of gene inter-connections. CONCLUSION: This translational genomics study highlighted the existence both of individual genes, mostly dysregulated in certain cellular populations (e.g., macrophages), and gene-networks involved in pulmonary reaction after left pneumonectomy. Their involvement in lung homeostasis is largely supported by previous studies, carried out both in humans and in other animal models (under homeostatic or disease-related conditions), that adopted candidate-gene approaches. Overall, the present findings represent a preliminary assessment for future, more focused, studies on compensatory lung adaptation, pulmonary regeneration and functional reload.

Exon-1 skipping and intron-1 retaining by alternative splicing of the c-KIT gene encodes a novel splice variant in the skin of Merino sheep (Ovis aries).

c-KIT, a type III receptor protein tyrosine kinase, plays an essential role in melanocyte development, migration, and survival. Mutations within the c-KIT gene are previously shown to cause the white coat color phenotypes in pigs, mice, goats, and humans. However, up so far, the splicing isoform(s), genomic architecture of c-KIT have not been characterized well in merino sheep. Reverse transcriptase (RT)-PCR analysis with molecular prediction identified two basic splice variants: Transcript Variant-1, 2 for 12 bp insertion coding sequences (CDS) corresponding to the four amino acids 'GNSK', respectively. Using 5' RACE, here we report for the first time a novel c-KIT 'Transcript Variant-3' from the skin of merino sheep by comparative genome analyses at exon(1)-intron(1)-exon(2) boundaries. In contrast, a single product of 795 bp was characterized by 3' RACE. We also demonstrated that the c-KIT gene expression at the transcript level is not mediated via an intron-9 splicing event. Overall, beyond what was observed in other mammals, our data provide novel insights into the molecular structure of the c-KIT gene in sheep.

PERHAPS: Paired-End short Reads-based HAPlotyping from next-generation Sequencing data.

The identification of rare haplotypes may greatly expand our knowledge in the genetic architecture of both complex and monogenic traits. To this aim, we developed PERHAPS (Paired-End short Reads-based HAPlotyping from next-generation Sequencing data), a new and simple approach to directly call haplotypes from short-read, paired-end Next Generation Sequencing (NGS) data. To benchmark this method, we considered the APOE classic polymorphism (*1/*2/*3/*4), since it represents one of the best examples of functional polymorphism arising from the haplotype combination of two Single Nucleotide Polymorphisms (SNPs). We leveraged the big Whole Exome Sequencing (WES) and SNP-array data obtained from the multi-ethnic UK BioBank (UKBB, N=48,855). By applying PERHAPS, based on piecing together the paired-end reads according to their FASTQ-labels, we extracted the haplotype data, along with their frequencies and the individual diplotype. Concordance rates between WES directly called diplotypes and the ones generated through statistical pre-phasing and imputation of SNP-array data are extremely high (>99%), either when stratifying the sample by SNP-array genotyping batch or self-reported ethnic group. Hardy-Weinberg Equilibrium tests and the comparison of obtained haplotype frequencies with the ones available from the 1000 Genome Project further supported the reliability of PERHAPS. Notably, we were able to determine the existence of the rare APOE*1 haplotype in two unrelated African subjects from UKBB, supporting its presence at appreciable frequency (approximatively 0.5%) in the African Yoruba population. Despite acknowledging some technical shortcomings, PERHAPS represents a novel and simple approach that will partly overcome the limitations in direct haplotype calling from short read-based sequencing.

PDGFA in Cashmere Goat: A Motivation for the Hair Follicle Stem Cells to Activate.

The cashmere hair follicle (HF) perpetually goes through cycles of growth, involution and rest. The photoperiod is the main factor in the control of seasonal coat change in cashmere goats while stem cells play a crucial role in the HF growth. Several factors, including Platelet-Derived Growth Factor A (PDGFA), Bone Morphogenetic Protein 2 (BMP2) and Lim-Homeobox gene 2 (LHX2) are implicated in HF morphogenesis and cycle. In this work, the mentioned molecules were investigated to evaluate their role in follicular cycle activation. The study was performed on skin samples collected at different periods of HF cycle and the molecular expression of PDGFA, BMP2 and LHX2 was evaluated by Real-Time PCR (qPCR) at each time point. Since PDGFA showed the most variation, the goat PDGFA gene was sequenced and the protein localization was investigated by immunohistochemistry together with PDGF receptor α (PDGFRα). PDGFA immunostaining was observed in the basal layer of the HF outer root sheath and the immunoreaction appeared stronger in the regressive HFs compared to those in the anagen phase according to qPCR analysis. PDGFRα was observed in the HF epithelium, proving the effect of PDGFA on the follicular structure. The data obtained suggest that PDGFA and BMP2 are both implicated in HF cycle in goat. In particular, PDGFA secreted by the HF is involved in the anagen activation.

Evidence of post-transcriptional readthrough regulation in FGF5 gene of alpaca.

Two different phenotypes are described in alpaca, identified as suri and huacaya, which differ in the type of fleece. The huacaya fleece is characterized by compact, soft and highly crimped fibers, while the suri fleece is longer, straight, less-crimped and lustrous. In our study, the Fibroblast growth factor 5 (FGF5) was investigated as a possible candidate gene for hair length in alpaca (Vicugna pacos). As previously identified in other mammals, our results show that the alpaca FGF5 gene gives rise to a short (FGF5S) and a long (FGF5) isoform. Interestingly, in the long isoform, we observed a point mutation (i.e., a transition C>T at position 499 downstream of the ATG codon) that is able to generate a premature termination codon (PTC). The highly conserved nucleotide and amino acid sequence after PTC suggested a readthrough event (RT) that was confirmed by western blot analysis. The analysis of cDNA sequence revealed motifs and structures of mRNA undergoing RT. In fact, the event is positively influenced by particular signals harbored by the transcript. To the best of our knowledge, this is the first case of a readthrough event on PTC reported for the FGF5 gene and the first case of this translational mechanism in alpaca.