Assessing the power of principal components and wright's fixation index analyzes applied to reveal the genome-wide genetic differences between herds of Holstein cows

Genetic diversity of Plasmodium falciparum reticulocyte binding protein homologue-5, which is a potential malaria vaccine candidate: baseline data from areas of varying malaria endemicity in Mainland Tanzania

BACKGROUND

The limited efficacy of the two recently approved malaria vaccines, RTS,S/AS01 and R21/Matrix- M™, highlights the need for alternative vaccine candidate genes. Plasmodium falciparum Reticulocyte Binding Protein Homologue 5 (Pfrh5) is a promising malaria vaccine candidate, given its limited polymorphism, its essential role in parasite survival, a lack of immune selection pressure and higher efficacy against multiple parasites strains. This study evaluated the genetic diversity of Pfrh5 gene among parasites from regions with varying malaria transmission intensities in Mainland Tanzania, to generate baseline data for this potential malaria vaccine candidate.

METHODS

This study utilized secondary data of 697 whole-genome sequences which were generated by the MalariaGEN Community Network. The samples which were sequenced to generated the data were collected between 2010 and 2015 from five districts within five regions of Mainland Tanzania, with varying endemicities (Morogoro-urban district in Morogoro region, Muheza in Tanga, Kigoma-Ujiji in Kigoma, Muleba in Kagera, and Nachingwea district in Lindi region). Wright's fixation index (FST), Wright's inbreeding coefficient (Fws), Principal component analysis (PCA), nucleotide diversity (π), haplotype network, haplotype diversity (Hd), Tajima's D, and Linkage disequilibrium (LD) were used to assess the diversity of the gene.

RESULTS

Of the sequences used in this study, 84.5% (n = 589/697) passed quality control and 313 (53.1%) were monoclonal (contained infections from a single strain of P. falciparum) and were used for haplotype diversity and haplotype network analysis. High within-host diversity (Fws < 0.95) was reported in Kigoma-Ujiji (60.7%), Morogoro-urban (53.1%), and Nachingwea (50.8%), while Muleba (53.9%) and Muheza (61.6%) had low within-host diversity (Fws ≥ 0.95). PCA did not show any population structure and the mean FST value was 0.015. Low nucleotide diversity values were observed across the study sites (mean π = 0.00056). A total of 27 haplotypes were observed among the 313 monoclonal samples and under-fives exhibited higher haplotype counts. The Pf3D7 was detected as Hap_1, which occurred in 16/313 (5.1%) monoclonal sequences. Negative Tajima's D values were observed among the parasite populations in all the study sites.

CONCLUSION

Low levels of polymorphism in the pfrh5 gene were observed based on low nucleotide and haplotype diversity, a lack of population structure and negative Tajima's D values. This study provides essential data on the diversity of the Pfrh5 gene indicating that it can be considered in the development of the next generation malaria vaccines. Robust and intensive studies of this and other candidate genes are crucial to support the prioritization of the Pfrh5 gene for potential inclusion in a broadly cross-protective malaria vaccine.

Genetic diversity of United States Rambouillet, Katahdin and Dorper sheep

BACKGROUND

Managing genetic diversity is critically important for maintaining species fitness. Excessive homozygosity caused by the loss of genetic diversity can have detrimental effects on the reproduction and production performance of a breed. Analysis of genetic diversity can facilitate the identification of signatures of selection which may contribute to the specific characteristics regarding the health, production and physical appearance of a breed or population. In this study, breeds with well-characterized traits such as fine wool production (Rambouillet, N = 745), parasite resistance (Katahdin, N = 581) and environmental hardiness (Dorper, N = 265) were evaluated for inbreeding, effective population size (Ne), runs of homozygosity (ROH) and Wright's fixation index (FST) outlier approach to identify differential signatures of selection at 36,113 autosomal single nucleotide polymorphisms (SNPs).

RESULTS

Katahdin sheep had the largest current Ne at the most recent generation estimated with both the GONe and NeEstimator software. The most highly conserved ROH Island was identified in Rambouillet with a signature of selection on chromosome 6 containing 202 SNPs called in an ROH in 50 to 94% of the individuals. This region contained the DCAF16, LCORL and NCAPG genes that have been previously reported to be under selection and have biological roles related to milk production and growth traits. The outlier regions identified through the FST comparisons of Katahdin with Rambouillet and Dorper contained genes with known roles in milk production and mastitis resistance or susceptibility, and the FST comparisons of Rambouillet with Katahdin and Dorper identified genes related to wool growth, suggesting these traits have been under natural or artificial selection pressure in these populations. Genes involved in the cytokine-cytokine receptor interaction pathways were identified in all FST breed comparisons, which indicates the presence of allelic diversity between these breeds in genomic regions controlling cytokine signaling mechanisms.

CONCLUSIONS

In this paper, we describe signatures of selection within diverse and economically important U.S. sheep breeds. The genes contained within these signatures are proposed for further study to understand their relevance to biological traits and improve understanding of breed diversity.

Identification of homozygosity-rich regions in the Holstein genome

In this study, 371 Holstein cows from six herds and 26 Holstein bulls, which were used in these herds, were genotyped by the Illumina BovineSNP50 array. For runs of homozygosity (ROH) identification, consecutive and sliding runs were performed by the detectRUNS and Plink software. The missing calls did not significantly affect the ROH data. The mean number of ROH identified by consecutive runs was 95.4 ± 2.7, and that by sliding runs was 86.0 ± 2.6 in cows, while this number for Holstein bulls was lower 58.9 ± 1.9. The length of the ROH segments varied from 1 Mb to over 16 Mb, with the largest number of ROH having a length of 1-2 Mb. Of the 29 chromosomes, BTA 14, BTA 16, and BTA 7 were the most covered by ROH. The mean coefficient of inbreeding across the herds was 0.111 ± 0.003 and 0.104 ± 0.004 based on consecutive and sliding runs, respectively, and 0.078 ± 0.005 for bulls based on consecutive runs. These values do not exceed those for Holstein cattle in North America. The results of this study confirmed the more accurate identification of ROH by consecutive runs, and also that the number of allowed heterozygous SNPs may have a significant effect on ROH data.

Estimation of genetic variability and identification of regions under selection based on runs of homozygosity in Beijing-You Chickens

The genetic composition of populations is the result of a long-term process of selection and adaptation to specific environments and ecosystems. Runs of homozygosity (ROHs) are homozygous segments of the genome where the 2 haplotypes inherited from the parents are identical. The detection of ROH can be used to describe the genetic variability and quantify the level of inbreeding in an individual. Here, we investigated the occurrence and distribution of ROHs in 40 Beijing-You Chickens from the random breeding population (BJY_C) and 40 Beijing-You Chickens from the intramuscular fat (IMF) selection population (BJY_S). Principal component analysis (PCA) and maximum likelihood (ML) analyses showed that BJY_C was completely separated from the BJY_S. The nucleotide diversity of BJY_C was higher than that of BJY_S, and the decay rate of LD of BJY_C was faster. The ROHs were identified for a total of 7,101 in BJY_C and 9,273 in BJY_S, respectively. The ROH-based inbreeding estimate (F_ROH) of BJY_C was 0.079, which was significantly lower than that of BJY_S (F_ROH = 0.114). The results were the same as the estimates of the inbreeding coefficients calculated based on homozygosity (F_HOM), the correlation between uniting gametes (F_UNI), and the genomic relationship matrix (F_GRM). Additionally, the distribution and number of ROH islands in chromosomes of BJY_C and BJY_S were significantly different. The ROH islands of BJY_S that included genes associated with lipid metabolism and fat deposition, such as CIDEA and S1PR1, were absent in BJY_C. However, GPR161 was detected in both populations, which is a candidate gene for the formation of the unique five-finger trait in Beijing-You chickens. Our findings contributed to the understanding of the genetic diversity of random or artificially selected populations, and allowed the accurate monitoring of population inbreeding using genomic information, as well as the detection of genomic regions that affect traits under selection.

Genetic Architecture and Signatures of Selection in the Caqueteño Creole (Colombian Native Cattle)

Evolutionary mechanisms have shaped the genomic architecture of Colombian Creole cattle breeds. The mating and selection processes have impacted several traits, promoting differences within and between populations. Studies of population structure and selection signatures in Colombian Creole breeds are scarce, and need more attention to better understand genetic differentiation, gene flow, and genetic distance. This study aimed to analyze the population structure and identify selection imprints in the Criollo Caqueteño (CAQ) population. It used 127 CAQ animals genotyped with Chip HD 777,000 SNPs. The population structure analyses used discriminant principal component analysis (DAPC), integrated haplotype scoring (iHS), and index-fixing (Fst) methodologies to detect selection signals. We can highlight SNP regions on the genes TMPRSS15, PGAM2, and EGFR, identified by the Fst method. Additionally, the iHS regions for cluster 1 identified candidate genes on BTA 3 (CMPK1 and FOXD2), BTA 11 (RCAN1), and BTA 22 (ARPP21). In group 2, we can highlight the genes on BTA 4 (SLC13A4, BRAF), BTA 9 (ULBP), BTA 14 (CSMD3) and BTA 19 (KRTAP9-2). These candidate genes have been associated with fertility traits, precocity, growth, and environmental and disease resistance, indicating a genetic potential in CAQ animals. All this promotes a better understanding of the diversity and genetic structure in the CAQ population. Based on that, our study can significantly assist the sustainable development and conservation of the breed in the Colombian Amazon.