Development and validation of the Axiom-MaruPri SNP chip for genetic analyses of domesticated old world camelids.

Camels play a crucial socio-economic role in sustaining the livelihoods of millions in arid and semi-arid regions. They possess remarkable physiological attributes which enable them to thrive in extreme environments, and provide a source of meat, milk and transportation. With their unique traits, camels embody an irreplaceable source of untapped genomic knowledge. This study introduces Axiom-MaruPri, a medium-density SNP chip meticulously designed and validated for both Camelus bactrianus and Camelus dromedarius. Comprising of 182,122 SNP markers, derived from the re-sequenced data of nine Indian dromedary breeds and the double-humped Bactrian camel, this SNP chip offers 34,894 markers that display polymorphism in both species. It achieves an estimated inter-marker distance of 14 Kb, significantly enhancing the coverage of the camel genome. The medium-density chip has been successfully genotyped using 480 camel samples, achieving an impressive 99 % call rate, with 96 % of the 182,122 SNPs being highly reliable for genotyping. Phylogenetic analysis and Discriminant Analysis of Principal Components yield clear distinctions between Bactrian camels and dromedaries. Moreover, the discriminant functions substantially enhance the classification of dromedary camels into different breeds. The clustering of various camel breeds reveals an apparent correlation between geographical and genetic distances. The results affirm the efficacy of this SNP array, demonstrating high genotyping precision and clear differentiation between Bactrian and dromedary camels. With an enhanced genome coverage, accuracy and economic efficiency the Axiom_MaruPri SNP chip is poised to advance genomic breeding research in camels. It holds the potential to serve as an invaluable genetic resource for investigating population structure, genome-wide association studies and implementing genomic selection in domesticated camelid species.

Whole genome resequencing revealed genomic variants and functional pathways related to adaptation in Indian yak populations.

The present study aims to identify genomic variants through a whole genome sequencing (WGS) approach and uncover biological pathways associated with adaptation and fitness in Indian yak populations. A total of 30 samples (10 from each population) were included from Arunachali, Himachali and Ladakhi yak populations. WGS analysis revealed a total of 32171644, 27260825, and 32632460 SNPs and 4865254, 4429941, and 4847513 Indels in the Arunachali, Himachali, and Ladakhi yaks, respectively. Genes such as RYR2, SYNE2, BOLA, HF1, and the novel transcript ENSBGRG00000011079 were found to have the maximum number of high impact variants in all three yak populations, and might play a major role in local adaptation. Functional enrichment analysis of genes harboring high impact SNPs revealed overrepresented pathways related to response to stress, immune system regulation, and high-altitude adaptation. This study provides comprehensive information about genomic variants and their annotation in Indian yak populations, thus would serve as a data resource for researchers working on the yaks. Furthermore, it could be well exploited for better yak conservation strategies by estimating population genetics parameters viz., effective population size, inbreeding, and observed and expected heterozygosity.

Design and validation of high-density SNP array of goats and population stratification of Indian goat breeds.

Goats are the supporting pillars of rural economy contributing significantly to meat and milk production in India. It is a species targeted for fulfilling the interdependent goals of poverty reduction and creation of employment for supporting the rural income. The increased demand for goat products necessitates their genetic characterization and improvement to augment the production of native breeds. Bi-allelic, genome wide, densely placed single nucleotide polymorphism (SNP) markers are most suitable for this purpose. This paper describes the design and validation of an Affymetrix Axiom-based high-density (HD) SNP chip for goats. The array was designed using a panel of 225 samples from 15 diverse goat breeds of India. In total, more than 38 million high quality SNPs were subjected to stringent filtering and 626,975 SNPs were finally tiled on the array. The average coverage of SNPs in our chip is one SNP per four kilobase (kb), providing a denser coverage of the goat genome than previously available arrays. The HD chip (Axiom_Cahi) was validated by genotyping 443 samples from 26 indigenous goat breeds/populations. The results revealed 95.83% markers to be highly informative and polymorphic in Indian goats. Multivariate analysis indicated population structuring, as 15 breeds could be segregated using the designed array. Phylogenetic analysis suggested stratification of breeds by geographic proximity. This HD SNP chip for goats is a valuable resource for genomic selection, genome wide association as well as population genetic studies in goats.

Transcriptomics reveals key genes responsible for functional diversity in pectoralis major muscles of native black Kadaknath and broiler chicken.

RNA sequencing-based expression profiles from pectoralis major muscles of black meat (Kadaknath) and white meat (broiler) chicken were compared to identify differentially expressed genes. A total of 156 genes with log2 fold change ≥ ± 2.0 showed higher expression in Kadaknath and 68 genes were expressed at a lower level in comparison to broiler. Significantly enriched biological functions of up-regulated genes in Kadaknath were skeletal muscle cell differentiation, regulation of response to reactive oxygen, positive regulation of fat cell differentiation and melanosome. Significant ontology terms up-regulated in broiler included DNA replication origin binding, G-protein coupled receptor signaling pathway and chemokine activity. Highly inter-connected differentially expressed genes in Kadaknath (ATFs, C/EPDs) were observed to be important regulators of cellular adaptive functions, while in broiler, the hub genes were involved in cell cycle progression and DNA replication. The study is an attempt to get an insight into the transcript diversity of pectoralis major muscles of Kadaknath and broiler chicken. Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-023-03682-0.

Validation of stable reference genes in peripheral blood mononuclear cells for expression studies involving vector-borne haemoparasitic diseases in bovines.

Normalization of gene expression data using appropriate reference genes is critical to diminish any technical bias in an experiment involving quantitative real-time PCR (qPCR). To the best of our knowledge, this is the first report offering a systematic assessment of 14 potential reference genes (RPLP0, ACTB, RPS28, YWHAZ, SDHA, PPIA, RPS9, RPS15, UXT, GAPDH, B2M, BACH1, HMBS, and PPIB) for the identification of the most stable normalizers for qPCR of target genes in peripheral blood mononuclear cells (PBMCs) of bovines for vector-borne haemoparasitic diseases such as anaplasmosis, babesiosis, theileriosis, and trypanosomiasis. A total of 38 blood samples were collected from healthy as well as diseased cattle and buffaloes representing different haemoparasitic diseases. RNA isolated from the PBMCs was subjected to qPCR for the 14 prospective internal control genes. The comprehensive ranking of the genes was accomplished by the RefFinder tool that integrates the results of three algorithms (geNorm, NormFinder, and BestKeeper) and the comparative CT method. RPS15, B2M, and GAPDH were ranked to be the most stable genes, whereas, PPIA and HMBS emerged to be the least suitable genes. Validation of the selected reference genes by the qPCR analysis of two immunity genes, ISG15 and GPX7 was congruent with the observations of this study. We recommend that a panel of three reference genes including RPS15, B2M, and GAPDH could prove useful in delineating the transcriptional landscape of PBMCs for vector-borne haemoparasitic diseases in bovines.

Comparative expression profiling of cytokine genes in Theileria annulata-infected and healthy cattle.

Theileriosis is one of the top ten economically important diseases in cattle in India. Cytokines are considered important mediators and regulators of the immune response to an infection. In the present study, the gene expression profiles of fourteen cytokines (IL1A, IL1B, IL2, IL4, IL6, IL8, IL10, IL12A, IL12B, IL16, TGFB1, TNFA, IFNA and IFNB) were compared in Theileria annulata-infected and healthy crossbred cattle. Blood samples were obtained from the District Disease Diagnostic Laboratory, Karnal. The presence/absence of T. annulata infection in the animals was determined on the basis of blood smear examination and molecular detection through PCR using the genus-specific primers. Total RNA was extracted from peripheral blood mononuclear cells, which was further reverse transcribed to cDNA. Primer3 software was employed to design the primers for Real-Time qPCR. The results were examined using 2-∆∆Ct method with RPS15 and GAPDH as the reference genes. The expression of IL1B, IL6, IL8, IL10, IL12A, IL12B, TNFA, IFNA and IFNB was significantly higher, whereas the expression of IL2 was lower in the infected animals. The transcript levels of IL1A and TGFB1 were also higher in the diseased animals, but the results were non-significant. This study profiles the expression kinetics of various pro-inflammatory and anti-inflammatory cytokine genes in response to bovine theileriosis.

Expression profiling of cytokine genes in peripheral blood mononuclear cells from Anaplasma marginale infected and healthy cattle.

In this study, changes in expression profiles of genes encoding 14 cytokines (IL1A, IL1B, IL2, IL4, IL6, IL8, IL10, IL12A, IL12B, IL16, IFNA, IFNB, TGFB1, and TNFA) were investigated amongst six Anaplasma marginale infected and six healthy crossbred cattle. Health status of the animals was determined based on clinical signs, blood smear examination and molecular detection using A. marginale-specific primers. Total RNA was isolated from the peripheral blood mononuclear cells of the infected animals as well as the healthy controls, which was further reverse transcribed to cDNA. Primers for real time PCR were designed using Primer3 software and the results were analyzed by the 2-ΔΔCt method with RPS15 and GAPDH as the reference genes. The expression levels of IL1A, IL1B, IL6, IL10, IL12A, IL12B, and TNFA varied significantly between the two groups, with higher expression in the infected cattle. The transcript abundance of IL4, IL16, and TGFB1 did not vary between the diseased and healthy animals. The expression of IL2 and IL8 was higher in the healthy animals, but the results were non-significant. Taken together, this study provides evidence for difference in expression of cytokine genes in response to anaplasmosis in crossbred cattle.

Muscle transcriptome provides the first insight into the dynamics of gene expression with progression of age in sheep.

The dynamic synergy of genes and pathways in muscles in relation to age affects the muscle characteristics. Investigating the temporal changes in gene expression will help illustrate the molecular mechanisms underlying muscle development. Here we report the gene expression changes in skeletal muscles through successive age groups in Bandur, a meat type sheep of India. RNA sequencing data was generated from the longissimus thoracis muscles from four age groups, ranging from lamb to adult. Analysis of 20 highest expressed genes common across the groups revealed muscle protein, phosphorylation, acetylation, metal binding and transport as significant functions. Maximum differentiation was observed after 2.5-3 years on transition from lambs to adult. Transcriptional regulation by the TFAP2 transcription factors, IL-6 signaling and PI3K/AKT signaling pathways were enriched in younger animals. The gene-protein network demarcated key interactive genes involved in muscle development and proliferation that can be used as candidates for future research on improvement of muscle characteristics.

Transcriptomic diversity in longissimus thoracis muscles of Barbari and Changthangi goat breeds of India.

The present study is an attempt to examine the differential expression of genes in longissimus thoracis muscles between meat and wool type Indian goat breeds. Barbari goat is considered the best meat breed while Changthangi is famous for its fine fibre quality. RNA sequencing data was generated from four biological replicates of longissimus thoracis muscles of Barbari and Changthangi goats. A clear demarcation could be observed between the breeds in terms of expression of genes associated with lipid metabolism (FASN, SCD, THRSP, DGAT2 and FABP3). Most significant genes with high connectivity identified by gene co-expression network analysis were associated with triacylglycerol biosynthesis pathway in Barbari goat. Highly interactive genes identified in Changthangi goat were mainly associated with muscle fibre type. This study provides an insight into the differential expression of genes in longissimus thoracis muscles between Barbari and Changthangi goats that are adapted to and reared in different agro-climatic regions.

Comparative gene expression profiling of milk somatic cells of Sahiwal cattle and Murrah buffaloes.

India is the world's largest milk producing country because of massive contribution made by cattle and buffaloes. In the present investigation, comprehensive comparative profiling of transcriptomic landscape of milk somatic cells of Sahiwal cattle and Murrah buffaloes was carried out. Genes with highest transcript abundance in both species were enriched for biological processes such as lactation, immune response, cellular oxidant detoxification and response to hormones. Analysis of differential expression identified 377 significantly up-regulated and 847 significantly down-regulated genes with fold change >1.5 in Murrah buffaloes as compared to Sahiwal cattle (padj <0.05). Marked enrichment of innate and adaptive immune response related GO terms and higher expression of genes for various host defense peptides such as lysozyme, defensin ß and granzymes were evident in buffaloes. Genes related to ECM-receptor interaction, complement and coagulation cascades, cytokine-cytokine receptor interaction and keratinization pathway showed more abundant expression in cattle. Network analysis of the up-regulated genes delineated highly connected genes representing immunity and haematopoietic cell lineage (CBL, CD28, CD247, PECAM1 and ITGA4). For the down-regulated dataset, genes with highest interactions were KRT18, FGFR1, GPR183, ITGB3 and DKK3. Our results lend support to more robust immune mechanisms in buffaloes, possibly explaining lower susceptibility to mammary infections as compared to cattle.

Differential expression of miRNAs in skeletal muscles of Indian sheep with diverse carcass and muscle traits.

The study presents the miRNA profiles of two Indian sheep populations with divergent carcass and muscle traits. The RNA sequencing of longissimus thoracis muscles from the two populations revealed a total of 400 known miRNAs. Myomirs or miRNAs specific to skeletal muscles identified in our data included oar-miR-1, oar-miR-133b, oar-miR-206 and oar-miR-486. Comparison of the two populations led to identification of 100 differentially expressed miRNAs (p < 0.05). A total of 45 miRNAs exhibited a log2 fold change of ≥ ( ±) 3.0. Gene Ontology analysis revealed cell proliferation, epithelial to mesenchymal transition, apoptosis, immune response and cell differentiation as the most significant functions of the differentially expressed miRNAs. The differential expression of some miRNAs was validated by qRT-PCR analysis. Enriched pathways included metabolism of proteins and lipids, PI3K-Akt, EGFR and cellular response to stress. The microRNA-gene interaction network revealed miR-21, miR-155, miR-143, miR-221 and miR-23a as the nodal miRNAs, with multiple targets. MicroRNA-21 formed the focal point of the network with 42 interactions. The hub miRNAs identified in our study form putative regulatory candidates for future research on meat quality traits in Indian sheep. Our results provide insight into the biological pathways and regulatory molecules implicated in muscling traits of sheep.

Milk somatic cell derived transcriptome analysis identifies regulatory genes and pathways during lactation in Indian Sahiwal cattle (Bos indicus).

BACKGROUND: The present study is an effort to understand the genomic drivers of lactation in Sahiwal (Bos indicus), the best milch cattle breed of the tropics. METHODS: RNA sequencing of four animals from early, mid and late lactation stages was performed using milk somatic cells as source of RNA. RESULTS: The genes encoding the milk casein and whey proteins showed highest expression in early and mid lactation, with a declining trend towards the late stage. The enhanced expression of PLIN2, FABP5 and FABP3 genes in mid lactation suggests enrichment of the PPARα pathway which is linked to fatty acid metabolism. A gradual decline in the percentage of genes involved in metabolism of proteins, mRNA and insulin synthesis from early to late lactation reflected transition from lactogenesis to involution. Major biological pathways maintained throughout lactation were adaptive immune system, FGF signaling, EGFR signaling, activated TLR4 signaling, NFkB and MAP kinases activation mediated by TLR4 signaling repertoire. Differential expression analysis revealed 547, 1010 and 1313 differentially expressed genes (p < 0.05) between early-late, early-mid and mid-late stages, respectively. The topmost regulatory genes identified by network analysis from the differentially expressed genes, were involved in Chemokine receptor, GPCR and EGFR1 pathways. CONCLUSION: The genes and pathways delineated in this study have regulatory implications in cell morphogenesis, lipid droplet formation and protein synthesis in the course of lactation. The study provides an insight into the expression profile of genes influencing milk properties and lactation in Sahiwal cattle.

Skin transcriptome profiling of Changthangi goats highlights the relevance of genes involved in Pashmina production.

Pashmina, the world's finest natural fiber is derived from secondary hair follicles of Changthangi goats which are domesticated in Ladakh region of Jammu and Kashmir by nomadic pastoralists. Complex epithelial-mesenchymal interactions involving numerous signal molecules and signaling pathways govern hair follicle morphogenesis and mitosis across different species. The present study involved transcriptome profiling of skin from fiber type Changthangi goats and meat type Barbari goats to unravel gene networks and metabolic pathways that might contribute to Pashmina development. In Changthangi goats, 525 genes were expressed at significantly higher levels and 54 at significantly lower levels with fold change >2 (padj < 0.05). Functional annotation and enrichment analysis identified significantly enriched pathways to be formation of the cornified envelope, keratinization and developmental biology. Expression of genes for keratins (KRTs) and keratin-associated proteins (KRTAPs) was observed to be much higher in Changthangi goats. A host of transcriptional regulator genes for hair follicle keratin synthesis such as GPRC5D, PADI3, HOXC13, FOXN1, LEF1 and ELF5 showed higher transcript abundance in Pashmina producing goats. Positive regulation of Wnt signaling pathway and negative regulation of Oncostatin M signaling pathway may be speculated to be important contributors to hair follicle development and hair shaft differentiation in Changthangi goats.

Transcriptome profiling of longissimus thoracis muscles identifies highly connected differentially expressed genes in meat type sheep of India.

This study describes the muscle transcriptome profile of Bandur breed, a consumer favoured, meat type sheep of India. The transcriptome was compared to the less desirable, unregistered local sheep population, in order to understand the molecular factors related to muscle traits in Indian sheep breeds. Bandur sheep have tender muscles and higher backfat thickness than local sheep. The longissimus thoracis transcriptome profiles of Bandur and local sheep were obtained using RNA sequencing (RNA Seq). The animals were male, non-castrated, with uniform age and reared under similar environment, as well as management conditions. We could identify 568 significantly up-regulated and 538 significantly down-regulated genes in Bandur sheep (p≤0.05). Among these, 181 up-regulated and 142 down-regulated genes in Bandur sheep, with a fold change ≥1.5, were considered for further analysis. Significant Gene Ontology terms for the up-regulated dataset in Bandur sheep included transporter activity, substrate specific transmembrane, lipid and fatty acid binding. The down-regulated activities in Bandur sheep were mainly related to RNA degradation, regulation of ERK1 and ERK2 cascades and innate immune response. The MAPK signaling pathway, Adipocytokine signaling pathway and PPAR signaling pathway were enriched for Bandur sheep. The highly connected genes identified by network analysis were CNOT2, CNOT6, HSPB1, HSPA6, MAP3K14 and PPARD, which may be important regulators of energy metabolism, cellular stress and fatty acid metabolism in the skeletal muscles. These key genes affect the CCR4-NOT complex, PPAR and MAPK signaling pathways. The highly connected genes identified in this study, form interesting candidates for further research on muscle traits in Bandur sheep.

Buffalo milk transcriptome: A comparative analysis of early, mid and late lactation.

The expression of genes and their regulation during lactation in buffaloes remains less understood. To understand the interplay of various genes and pathways, the milk transcriptome from three lactation stages of Murrah buffalo was analyzed by RNA sequencing. The filtered reads were mapped to the Bubalus bubalis as well as Bos taurus reference assemblies. The average mapping rate to water buffalo and Btau 4.6 reference sequence, was 75.5% and 75.7% respectively. Highly expressed genes (RPKM > 3000), throughout lactation included CSN2, CSN1S1, CSN3, LALBA, SPP1 and TPT1. A total of 12833 transcripts were common across all the stages, while 271, 205 and 418 were unique to early, mid and late lactation respectively. Majority of the genes throughout lactation were linked to biological functions like protein metabolism, transport and immune response. A discernible shift from metabolism in early stage to metabolism and immune response in mid stage, and an increase in immune response functions in late lactation was observed. The results provide information of candidate genes and pathways involved in the different stages of lactation in buffalo. The study also identified 14 differentially expressed and highly connected genes across the three lactation stages, which can be used as candidates for future research.

Understanding the cryptic introgression and mixed ancestry of Red Junglefowl in India.

Red Junglefowls (RJFs), the wild progenitor of modern day chickens (DCs), are believed to be in genetic endangerment due to introgression of domestic genes through opportunistic matings with domestic or feral chickens. Previous studies from India reported rare hybridization of RJFs in the wild. However, RJF population genetic structure, pattern of gene flow and their admixture with DC populations are poorly understood at the landscape level. We conducted this study with a large sample size, covering the predicted natural distribution range of RJFs in India. We documented strong evidence of directional gene flow from DCs to free-ranging wild RJFs, with the Northeastern RJF population exhibiting the most genetic variants in their nuclear and mitochondrial genomes, indicating it to be the ancestral population from which early radiation may have occurred. The results provide evidence that landscape features do not act as a barrier to gene flow and the distribution pattern could not be explored due to physical sharing or exchange of wild birds in the past when forests were continuous across RJF range in India.

Mitochondrial D-loop analysis for uncovering the population structure and genetic diversity among the indigenous duck (Anas platyrhynchos) populations of India.

The indigenous domestic duck (Anas platyrhynchos domestica) which is domesticated from Mallard (Anas platyrhynchos) contributes significantly to poor farming community in coastal and North Eastern regions of India. For conservation and maintenance of indigenous duck populations it is very important to know the existing genetic diversity and population structure. To unravel the population structure and genetic diversity among the five indigenous duck populations of India, the mitochondrial D-loop sequences of 120 ducks were analyzed. The sequence analysis by comparison of mtDNA D-loop region (470 bp) of five Indian duck populations revealed 25 mitochondrial haplotypes. Pairwise FST value among populations was 0.4243 (p < .01) and the range of nucleotide substitution per site (Dxy) between the five Indian duck populations was 0.00034-0.00555, and the net divergence (Da) was 0-0.00355. The phylogenetic analysis in the present study unveiled three clades. The analysis revealed genetic continuity among ducks of coastal region of the country which formed a separate group from the ducks of the inland area. Both coastal as well as the land birds revealed introgression of the out group breed Khaki Campbell, which is used for breed improvement programs in India. The observations revealed very less selection and a single matrilineal lineage of indigenous domestic ducks.

Molecular characterization of RNA binding motif protein 3 (RBM3) gene from Pashmina goat.

Pashmina goat inhabits the high altitude cold arid desert of Ladakh, India. This goat is known for its finest and costliest under fiber. Though the under fiber may be a part of its complex thermoregulation mechanism, the genetics of its adaptability under cold conditions is not known. As an attempt to understand its adaptive genetics, and the role of RNA-binding proteins at the cellular response, this study was conducted to characterize the RBM3 gene in Pashmina goat and its expression during hypothermia. The ORF of Pashmina RBM3 gene was 273 bp. Phylogenetic analysis revealed that Pashmina RBM3 is closely related to Bos taurus RBM3. Pashmina RBM3 was characterized by comparative modeling studies. The final 3-D model contained two α-helices and four ß-sheets. qRT-PCR data showed that Pashmina RBM3 gene expression was significantly higher (P < 0.05) at moderate (30 °C) hypothermic stress conditions as compared with deep (15 °C) hypothermia.

Development of parthenote following in vivo transfer of embryos in Capra hircus.

The aim of this study is to generate parthenogenetic embryos from chemically activated in vitro matured caprine oocytes and to study the in vivo developmental potency of such embryos. The parthenogenetic embryos (2-8 and 16 cells to morula stage) were surgically transferred in 26 recipients. Pregnancy in recipients following embryo transfer was monitored by ultrasonography. The recipient aborted a foetus on day 34 post transfer. Sexing of parthenogenetic foetus showed a single band of amelogenin gene indicating female cell DNA. Microsatellite analysis revealed that the recipient has not contributed genetically to the parthenogenetic foetus confirming the identity of aborted foetus of parthenogenetic origin. The authors believe that this is the first authentic report on in vivo development of parthenogenetic foetus in Capra hircus.