The Pashmina and Barbari are two famous goat breeds found in the wide areas of the Indo-Pak region. Pashmina is famous for its long hair-fiber (Cashmere) production while Barbari is not-selected for this trait. So, the mRNA expression profiling in the skin samples of both breeds would be an attractive and judicious approach for detecting putative genes involved in this valued trait. Here, we performed differential gene expression analysis on publicly available RNA-Seq data from both breeds. Out of 44,617,994 filtered reads of Pashmina and 55,995,999 of Barbari which are 76.48% and 73.69% mapped to the ARS1 reference transcriptome assembly respectively. A pairwise comparison of both breeds resulted in 47,159 normalized expressed transcripts while 8,414 transcripts are differentially expressed above the significant threshold. Among these, 4,788 are upregulated in Pashmina while 3,626 transcripts are upregulated in Barbari. Fifty-nine transcripts harbor 57 genes including 32 LOC genes and 24 are annotated genes which were selected on the basis of TMM counts > 500. Genes with ectopic expressions other than uncharacterized and LOC symbol genes are Keratins (KRT) and Keratin Associated Proteins (KRTAPs), CystatinA&6, TCHH, SPRR4, PPIA, SLC25A4, S100A11, DMKN, LOR, ANXA2, PRR9 and SFN. All of these genes are likely to be involved in keratinocyte differentiation, sulfur matrix proteins, dermal papilla cells, hair follicles proliferation, hair curvature, wool fiber diameter, hair transition, hair shaft differentiation and its keratinization. These differentially expressed reported genes are critically valuable for enhancing the quality and quantity of the pashmina fiber and overall breed improvement. This study will also provide important information on hair follicle differentiation for further enrichment analyses and introducing this valued trait to other goat breeds as well.

An in-silico WES approach using the Galaxy platform was adopted in the current study to predict the genetic basis of Premature Ovarian Failure (POF), where three affected patients in a Saudi Arabian family of seven, found associated with X-linked recessive mutations. The current analysis discovered 518,054 variants using FreeBayes variant caller that had 1,461,864 effects on variable sites in the genome revealed by SnpEff software. The causal genetic mutations were filtered and annotated with the ClinVar database using the GEMINI tool. This tool retained 369 pathogenic mutations harboring 130 genes. Among the total, 268 variants positioned on 69 genes are shared with three affected individuals, 61 variants on 23 genes are shared by any two of the affected individuals, and 40 of the variants on 38 genes are present in any one of the affected sample. Two mutations in one of the already POF-associated, POF1B gene were also observed e.g. (i) g.84563135T>A; p.M349L and (ii) g.84563194C>T; p.R329Q in the two affected individuals i.e. IV-I-C & IV-6 in the current data. This gene consists of 17 exons that span the region of >100 kb. The putative function of this gene in regulating the actin cytoskeleton due to homology with myosin tail and maintains a number of oocytes during fetal ovary development. In a nutshell, this Galaxy pipeline facilitates all-in-one to pinpoint not only the known pathogenic gene mutations for this disorder but few other novel genetic variants as well, whose gene-disease association may be validated by further experimental studies.

Mitochondrial Encephalohepatopathy (MEH) is an autosomal recessive neurodevelopmental disorder usually accompanied by microcephaly, white matter changes, cardiac and hepatic failure. Here, we applied the whole-exome sequencing (WES) framework on a trio family data with unaffected non-consanguineous parents and proband (neonate girl) with this inherited disorder. A total of 2,928,402 variants were observed with 2,613,746 SNPs, 112,336 multiple nucleotide polymorphisms (MNPs), 72,610 insertions, 113,207 deletions and 16,503 mixed variants. These variations are responsible for 82,813,631 effects on various genomic regions. Our pipeline uncovered candidate gene mutations from these variants and retained a handful of 5,277 variants harboring 3,598 genes, out of which, 8 genes codes for non-coding RNA while 178 genes are those with high impact severity. Among these 178 variants, 125 are de-novo variants that are not previously reported in the ClinVar database. Consistent to previous studies, the leftover high impact severity genes are involved in encephalopathy, Leigh syndrome, Charcot–Marie–Tooth disease, global developmental disorder, seizures, spastic paraplegia, premature ovarian failure, mitochondrial myopathy-cerebellar, ataxia-pigmentary, retinopathy syndrome, ocular and retinal degeneration, deafness, intellectual disability, cardiofacioneurodevelopmental syndrome etc. All these clinical features were also observed in the patient studied. The current analysis highlights and expands the genetic architecture of the MEH phenotype. Furthermore, this pipeline on trio family data significantly broadens the concept of its usefulness as a first-tier diagnostic method in the detection of complex multisystem phenotypic disorders.

Whole genome pooled sequence data of 12 Pakistani Teddy goats is analyzed for positive selection signatures as their breed defining characteristics. Selection imprints left in the Teddy genome are unveiled by genomic differentiation after the successful paired-end alignment of 635,357,043 reads with (ARS1) reference genome assembly. Pooled-heterozygosity ( ) and Tajima’s D (TD) are applied for validation and getting better hits of selection signals, while pairwise F_ST statistics is conducted on Teddy vs. Bezoar (wild goat ancestor) for genomic differentiation. Annotation of regions under positive selection reveals 59 genes underlying production and adaptive traits. score ≥ 5 detected six windows having highest scores on Chr. 29, 9, 25, 15 and 14 that harbor HRASLS5, LACE1 and AXIN1 genes which are candidate for embryonic development, lactation and body height. Secondly, TD value of ≤ -2.2 showed 4 windows with very strong hits on Chr.5 & 9 harbor STIM1 and ADM genes related to body mass and weight. Lastly, F_ST analysis generated three strong signals with threshold ≤ 0.42 on Chr.12 & 5 harbor ITGB1 gene associated with milk production & lactation traits. Other significant selection signatures encompass genes associated with wool production, prolificacy, immunity and coat colors. In brief, this study identified the genes under selection in this Pakistani goat breed that will be helpful to refining future breeding policies and converging required productive traits within and across other goat breeds and to explore full genetic potential of this valued livestock species.