Reduced representation approach for identification of genome-wide SNPs and their annotation for economically important traits in Indian Tharparkar cattle

Genome-wide selection signatures address trait specific candidate genes in cattle indigenous to arid regions of India

The peculiarity of Indian cattle lies in milk quality, resistance to diseases and stressors as well as adaptability. The investigation addressed selection signatures in Gir and Tharparkar cattle, belonging to arid ecotypes of India. Double digest restriction-site associated DNA sequencing (ddRAD-seq) yielded nearly 26 million high-quality reads from unrelated seven Gir and seven Tharparkar cows. In all, 19,127 high-quality SNPs were processed for selection signature analysis. An approach involving within-population composite likelihood ratio (CLR) statistics and between-population FST statistics was used to capture selection signatures within and between the breeds, respectively. A total of 191 selection signatures were addressed using CLR and FST approaches. Selection signatures overlapping 86 and 73 genes were detected as Gir- and Tharparkar-specific, respectively. Notably, genes related to production (CACNA1D, GHRHR), reproduction (ESR1, RBMS3), immunity (NOSTRIN, IL12B) and adaptation (ADAM22, ASL) were annotated to selection signatures. Gene pathway analysis revealed genes in insulin/IGF pathway for milk production, gonadotropin releasing hormone pathway for reproduction, Wnt signalling pathway and chemokine and cytokine signalling pathway for adaptation. This is the first study where selection signatures are identified using ddRAD-seq in indicine cattle breeds. The study shall help in conservation and leveraging genetic improvements in Gir and Tharparkar cattle.

Phenotypic and genetic characterization of unexplored, potential cattle population of Madhya Pradesh

Bawri or Garri, a non-descript cattle population managed under an extensive system in Madhya Pradesh state of India, was identified and characterized both genetically and phenotypically to check whether or not it can be recognised as a breed. The cattle have white and gray colour and are medium sized with 122.5 ± 7.5 cm and 109.45 ± 0.39 cm height at withers in male and female, respectively. Double-digest restriction site associated DNA (ddRAD) sequencing was employed to identify ascertainment bias free SNPs representing the entire genome cost effectively; resulting in calling 1,156,650 high quality SNPs. Observed homozygosity was 0.76, indicating Bawri as a quite unique population. However, the inbreeding coefficient was 0.025, indicating lack of selection. SNPs found here can be used in GWAS and genetic evaluation programs. Considering the uniqueness of Bawri cattle, it can be registered as a breed for its better genetic management.

Genomic clues of association between clinical mastitis and SNPs identified by ddRAD sequencing in Murrah buffaloes

The total milk production of India is 209.96 MT out of which 45% is contributed by the indigenous buffalo and due to their high producing virtue, the prevalence of mastitis is 5-20%. Despite the increasing level of technological advancement, mastitis is still an issue of concern for dairy industry in India as well as across the world. Therefore, the present study aimed to identify the SNPs and associate them with the incidence of clinical mastitis in Murrah buffalo using the ddRAD sequencing approach taking mastitis incidence data of 96 Murrah buffaloes. A total of 246 million quality controlled reads were obtained with an average alignment rate of 99.01% and at a read depth of 10, quality controlled SNPs obtained were 18,056. The logistic regression model was used and a total of seven SNPs were found significantly associated (p < 0.001) with mastitis incidence and seven genes were identified viz., NCBP1, FOXN3, TPK1, XYLT2, CPXM2, HERC1, and OPCML. The majority of them were having tumor suppressing action, related to immunogenetics or glycolytic and energy production. Conclusively, the SNPs identified in this study may be useful for future studies on mastitis incidence in Murrah buffalo and the SNP associations can be further validated.

Genome-wide genetic variation and population structure of native and cosmopolitan cattle breeds reared in Türkiye

This is the first comprehensive study to reveal genetic variation and population structure at genome level in six Anatolian (Anatolian Black, East Anatolian Red, South Anatolian Red, South Anatolian Yellow, Turkish Grey Steppe, and Zavot) and two cosmopolitan (Brown Swiss and Holstein Friesian) cattle breeds reared in Türkiye. Being 20 samples from each population, a total of 160 blood samples retrieved from representative herds were utilized to generate genomic libraries by ddRADseq method. Genomic libraries sequenced by Illumina HiSeq X Ten instrument revealed a total of 211,119 bi-allelic SNPs with high call rate. Compared to cosmopolitan cattle breeds, a higher genetic variation was observed in native Turkish cattle with an average of 0.380 observed heterozygosity. Genetic distances were comparatively low between native cattle breeds, whereas the highest genetic distance (0.064) was detected between South Anatolian Yellow and Brown Swiss. Population structure analyses showed that the native Turkish and cosmopolitan cattle breeds were clearly different from each other according to their phylogenetic origin. Besides, a high level of genetic admixture was detected among Anatolian Black, Turkish Grey Steppe, South Anatolian Red, and South Anatolian Yellow, whereas East Anatolian Red and Zavot were distinct from the other native and cosmopolitan cattle breeds. TreeMix algorithm under the assumption of one and two migration events revealed a migration route from Anatolian clade to Anatolian Black, while a second migration edge was drawn from Brown Swiss to East Anatolian Red. This study demonstrates the importance of national conservation studies in the native breeds whose population size has dramatically decreased. In addition, SNP arrays and next-generation sequencing platforms are recommended for future studies to reveal the genetic variation of other local Turkish livestock species to arrange effective conservation programs.

Genome-wide identification and annotation of SNPs for economically important traits in Frieswal™, newly evolved crossbred cattle of India

The Frieswal™ is a crossbred cattle evolved by ICAR-Central Institute for Research on Cattle utilizing more than 15,000 cattle maintained at more than 37 military farms spread all over the agro-climatic regions of the country. The ddRAD sequencing method was used to identify and annotate the SNPs and INDELs. The results of variant calling revealed 1,487,851 SNPs and 128,175 INDELs at a read depth of 10. A total of 3,775,079 effects were identified, and majority (66.41%) of the effects were in the intron region of the genome followed by intergenic (21.87%). Majority (99.18%) of the variants had the modifier effect. The results revealed a higher magnitude of transitions as compared to the transversion. The classification of SNPs by functional class revealed a majority of missense (43%) and silent (56%) effects. Out of 26,278 genes identified, 1841 SNPs were annotated in 207 candidate genes responsible for various milk production and reproduction traits. The observed heterozygosity was 0.2804 against the expected heterozygosity value of 0.2978. The overall average inbreeding coefficient (FIS) was 0.0604. The pathway analysis revealed that the prolactin signaling pathway (GO:0038161) was significant biological process complete for both milk production and reproduction traits. The SNP variations can be effectively used as markers for early and accurate identification of the QTLs and for formulating an efficient and effective breed improvement program in Frieswal™ cattle.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03701-0.

Exploring molecular dynamic indicators associated with reproductive performance of Bos indicus cattle in blood plasma samples through data-independent acquisition mass spectrometry

Improving reproductive performance of cattle is of paramount importance for sustainable dairy farming. Poor reproduction performance (RP) hinders the genetic improvement of important Bos indicus cattle breeds. It is well known that incorporation of molecular information along with conventional breeding method is far better than use of conventional method alone for the genetic improvement of reproductive performance traits in cattle. Therefore, the present study sought to investigate the plasma proteome of the Deoni cows in cyclical (n = 6) and pregnant (n = 6) reproductive phases with varying reproductive performance (high and low). High-throughput data independent acquisition (DIA) based proteomics was performed to understand corresponding proteome. We identified a total of 430 plasma proteins. Among cyclic cows, twenty proteins were differentially regulated in low RP as compared to high RP. BARD1 and AFP proteins were observed upregulated in cyclical cows whose upregulation reported to affect reproductive performance in cattle. Among the pregnant cows, thirty-five proteins were differentially regulated, including the downregulation of FGL2 and ZNFX1 that modulates the maternal immune response mechanism which is required for successful implantation of the embryo. Also, proteins such as AHSG, CLU and SERPINA6 were upregulated in the pregnant cows whose upregulation reported to reduced reproductive performance. The results of this study will be helpful in establishing a framework for future research on the aspect of improving reproductive performance in Bos indicus cattle breeds. SIGNIFICANCE: The Indian subcontinent is the center of domestication for Bos indicus cattle breeds and they are known for their disease resistance, heat tolerance, ability to survive in low input regime and harsh climatic conditions. In recent times, population of many important Bos indicus breeds including Deoni cattle is declining due to various factors, especially due to reproductive performance. Traditional breeding methods are not sufficient enough to understand and improve the reproductive performance traits in important Bos indicus cattle breeds. Proteomics approach is a promising technology to understand the complex biological factors which leads to poor reproductive performance in cattle. The present study utilized DIA based LC- MS/MS analysis to identify the plasma proteins associated with reproductive performance in cyclical and pregnant cows. This study if improved further, can be used to develop potential protein markers associated with reproductive performance which is useful for the selection and genetic improvement of important Bos indicus breeds.

ddRAD sequencing based genotyping of six indigenous dairy cattle breeds of India to infer existing genetic diversity and population structure

The present investigation aimed to identify genome wide SNPs and to carry out diversity and population structure study using ddRAD-seq based genotyping of 58 individuals of six indigenous milch cattle breeds (Bos indicus) such as Sahiwal, Gir, Rathi, Tharparkar, Red Sindhi and Kankrej of India. A high percentage of reads (94.53%) were mapped to the Bos taurus (ARS-UCD1.2) reference genome assembly. Following filtration criteria, a total of 84,027 high quality SNPs were identified across the genome of 6 cattle breeds with the highest number of SNPs observed in Gir (34,743), followed by Red Sindhi (13,092), Kankrej (12,812), Sahiwal (8956), Tharparkar (7356) and Rathi (7068). Most of these SNPs were distributed in the intronic regions (53.87%) followed by intergenic regions (34.94%) while only 1.23% were located in the exonic regions. Together with analysis of nucleotide diversity (π = 0.373), Tajima's D (D value ranging from - 0.295 to 0.214), observed heterozygosity (HO ranging from 0.464 to 0.551), inbreeding coefficient (FIS ranging from - 0.253 to 0.0513) suggested for the presence of sufficient within breed diversity in the 6 major milch breeds of India. The phylogenetic based structuring, principal component and admixture analysis revealed genetic distinctness as well as purity of almost all of the 6 cattle breeds. Overall, our strategy has successfully identified thousands of high-quality genome wide SNPs that will further enrich the Bos indicus representation basic information about genetic diversity and structure of 6 major Indian milch cattle breeds which should have implications for better management and conservation of valuable indicine cattle diversity.

Genome-wide identification and annotation of SNPs and their mapping in candidate genes related to milk production and fertility traits in Badri cattle

This study was conducted in Badri cattle using a double digest restriction-site associated DNA sequencing approach. The study aimed to identify and annotate high confidence single nucleotide polymorphisms (SNPs) and their mapping in candidate genes related to production and fertility in dairy cattle. A total of 7,168,552 genome-wide SNPs were initially identified in Badri cattle by alignment with the Bos indicus reference genome. After filtration of SNPs, 65,483 high confidence SNPs were retained and further used for downstream analysis. Annotation of high confidence SNPs revealed 99.197% SNPs had modifier impact, 0.326% SNPs were low impact, 0.036% were high impact, and 0.441% were moderate impact SNPs. Most SNPs in Badri cattle were found in intergenic, transcript and intronic regions. The candidate genes for milk production PRKCE, ABCG2, GHR, EPS8, CAST and NRXN1 were found to harbour maximum high confidence variants. Among candidate genes for fertility in cattle, ATP2B1, SOX5, WDR27, ARHGAP12, CACNA1D, ANKRD6, GRIA3, ZNF521 and CAST822 have maximum high confidence variants mapped in them. The SNPs found mapped in the candidate genes will be important genetic tools in the search for phenotype-modifying nucleotide changes and will aid in formulating relevant genetic improvement programmes for dairy cattle.