Genetic Differentiation among Livestock Breeds-Values for Fst

Genome-wide re-sequencing reveals selection signatures for important economic traits in Taihang chickens

Taihang chickens is precious genetic resource with excellent adaptability and disease resistance, as well as high-quality eggs and meat. However, the genetic mechanism underlying important economic traits remain largely unknown. To address this gap, we conducted whole-genome resequencing of 66 Taihang and 15 White Plymouth rock chicken (Baiyu). The population structure analysis revealed that Taihang chickens and Baiyu are 2 independent populations. The genomic regions with strong selection signals and some candidate genes related to economic and appearance traits were identified. Additionally, we found a continuously selected 1.2 Mb region on chromosome 2 that is closely related to disease resistance. Therefore, our findings were helpful in further understanding the genetic architecture of the Taihang chickens and provided a worthy theoretical basis and technological support to improve high-quality Taihang chickens.

Identification and validation of novel breed-specific biomarker for the purpose of village chicken authentication using genomics approaches

Local village chicken, or "Ayam kampung" as it's known in Malaysia, is considered a premium chicken breed with a higher price than other chicken breeds. As a result of their comparable appearances and sizes, colored broiler chickens are often sold as village chickens, which is a form of food fraud that can result in a 3- to 4-fold rise in profit. Therefore, developing a breed-specific authentication method is crucial for preventing food fraud in the poultry industry. This study aims to investigate the genetic diversity of village chickens from other commercial chicken breed populations available in the market (broiler [Cobb], colored broiler [Hubbard], and layer [DeKalb]) to identify breed-specific DNA fragments as biomarkers for village chicken authentication. The Whole-genome sequencing and mutation calling of 12 chickens (3 chickens/breed) led to the identification of a total of 73,454,654 single nucleotide polymorphisms (SNP) and 8,762,338 insertion and deletions (InDel) variants, with more variants detected in the village chicken population (6,346,704 SNPs; 752,408 InDels) compared to commercial breeds. Therefore, this study revealed that village chickens were more genetically variable compared to other breeds in Malaysia. Furthermore, the breed-specific genomic region located on chromosome 1 (1:84,405,652) harboring SNP (C-T) with high discrimination power was discovered and validated which can be considered as a novel breed-specific biomarker to develop a method for accurate authentication of village chickens in Malaysia. This authentication method offers potentialw applications in the chicken industry and food safety.

Genetic characterization of cashmere goat (Capra hircus) populations in Mongolia

Objective

Characterization studies of the phenotypic and genetic diversity of Mongolian goats are limited, despite several goat breeds being registered in the country. This study aimed to evaluate the phenotypic and genetic diversity of 14 cashmere goat populations in Mongolia, consisting largely of identified goat breeds.

Methods

Body weight, cashmere quality, and coat color were the phenotypic traits considered in this study. A linear model was used to fit body weight and cashmere traits, and least squares means (LSMs) were estimated for the region and location classes. Genetic diversity and structure were assessed using a goat 50K SNP array.

Results

The studied populations exhibited greater phenotypic diversity at the regional level. A very small overall differentiation index (Fst: 0.017) was revealed by Wright's Fst and a very small overall inbreeding index (F ROH1 :0.019) was revealed based on runs of homozygosity. Genetic clustering of populations by principal components showed large variances for the two goat populations of the Russian admixture (Gobi Gurvan Saikhan and Uuliin Bor), and smaller but differentiated clusters for the remaining populations. Similar results were observed in the admixture analysis, which identified populations with the highest (Govi Gurvan Saikhan and Uuliin Bor) and lowest (Tsagaan Ovoo Khar) exotic admixtures. A genomewide association study (GWAS) of body weight and cashmere traits identified a few significant variants on chromosomes 2, 4, 5, 9, and 15, with the strongest variant for cashmere yield on chromosome 4. The GWAS on coat color yielded nine significant variants, with the strongest variants located on chromosomes 6, 13, and 18 and potential associations with KIT, ASIP, and MC1R genes. These signals were also found in other studies on coat color and patterns in goats.

Conclusion

Mongolian cashmere goats showed relatively low genetic differentiation and low inbreeding levels, possibly caused by the traditional pastoral livestock management system and the practice of trading breeding bucks across provinces, along with a recent increase in the goat population. Further investigation of cashmere traits using larger samples and alternative methods may help identify the genes or genomic regions underlying cashmere quality in goats.

Analysis of Genetic Diversity in Romanian Carpatina Goats Using SNP Genotyping Data

Animal husbandry is one of man's oldest occupations. It began with the domestication of animals and developed continuously, in parallel with the evolution of human society. The selection and improvement of goats in Romania was not a clearly defined objective until around 1980. In recent years, with the increasing economic value given to goats, breeding programs are becoming established. In Romania, a few goat genetic studies using microsatellites and mtDNA have been carried out; however, a systematic characterization of the country's goat genomic resources remains missing. In this study, we analyzed the genetic variability of Carpatina goats from four distinct geographical areas (northern, north-eastern, eastern and southern Romania), using the Illumina OvineSNP60 (RefSeq ARS1) high-density chip for 67 goats. Heterozygosity values, inbreeding coefficients and effective population size across all autosomes were calculated for those populations that inhabit high- and low-altitude and high- and low-temperature environments. Diversity, as measured by expected heterozygosity (HE), ranged from 0.413 in the group from a low-temperature environment to 0.420 in the group from a high-temperature environment. Within studied groups, the HT (high temperature) goats were the only group with a positive but low average inbreeding coefficient value, which was 0.009. After quality control (QC) analysis, 46,965 SNPs remained for analysis (MAF < 0.01). LD was calculated for each chromosome separately. The Ne has been declining since the time of domestication, having recently reached 123, 125, 185 and 92 for the HA (high altitude), LA (low altitude), HT (high temperature) and LT (low temperature) group, respectively. Our study revealed a low impact of inbreeding in the Carpatina population, and the Ne trend also indicated a steep decline in the last hundred years. These results will contribute to the genetic improvement of the Carpatina breed.

Persistence of autozygosity in crossbreds between autochthonous and cosmopolitan breeds of swine: a simulation study

Crossbreeding might be a valid strategy to valorize local pig breeds. Crossbreeding should reduce homozygosity and, as a consequence, yield hybrid vigor for fitness and production traits. This study aimed to quantify the persistence of autozygosity in terminal crossbred pigs compared with purebreds and, in turn, identify genomic regions where autozygosity's persistence would not be found. The study was based on genotyping data from 20 European local pig breeds and three cosmopolitan pig breeds used to simulate crossbred offspring. This study consisted of two steps. First, one hundred matings were simulated for each pairwise combination of the 23 considered breeds (for a total of 276 combinations), ignoring the sex of the parent individuals in order to generate purebred and crossbred matings leveraging all the germplasm available. Second, a few preselected terminal-maternal breed pairs were used to mimic a realistic terminal crossbreeding system: (i) Mora Romagnola (boars) or Cinta Senese (boars) crossed with Large White (sows) or Landrace (sows); (ii) Duroc (boars) crossed with Mora Romagnola (sows) or Cinta Senese (sows). Runs of homozygosity was used to estimate genome-wide autozygosity (FROH). Observed FROH was higher in purebreds than in crossbreds, although some crossbred combinations showed higher FROH than other purebred combinations. Among the purebreds, the highest FROH values were observed in Mora Romagnola and Turopolje (0.50 and 0.46, respectively). FROH ranged from 0.04 to 0.16 in the crossbreds Alentejana × Large White and Alentejana × Iberian, respectively. Persistence of autozygosity was found in several genomic segments harboring regions where quantitative trait loci (QTLs) were found in the literature. The regions were enriched in QTLs involved in fatty acid metabolism and associated with performance traits. This simulation shows that autozygosity persists in most breed combinations of terminal crosses. Results suggest that a strategy for crossbreeding is implemented when leveraging autochthonous and cosmopolitan breeds to obtain most of the hybrid vigor.

Geographic population structure and distinct intra-population dynamics of globally abundant freshwater bacteria

Implications of geographic separation and temporal dynamics on the evolution of free-living bacterial species are widely unclear. However, the vast amount of metagenome sequencing data generated during the last decades from various habitats around the world provides an unprecedented opportunity for such investigations. Here, we exploited publicly available and new freshwater metagenomes in combination with the genomes of abundant freshwater bacteria to reveal geographic and temporal population structure. We focused on species that were detected across broad geographic ranges at high enough sequence coverage for meaningful population genomic analyses, associated with the predominant freshwater taxa acI, LD12, Polynucleobacter, and Candidatus Methylopumilus. Despite the broad geographic ranges, each species appeared as a sequence-discrete cluster, in contrast to abundant marine taxa, for which continuous diversity structures were reported on a global scale. Population differentiation increased significantly with spatial distance in all species, but notable dispersal barriers (e.g. oceanic) were not apparent. Yet, the different species showed contrasting rates of geographic divergence and strikingly different intra-population dynamics in time series within individual habitats. The change in an LD12 population over 7 years was minor (FST = 0.04) compared to differentiation between lakes, whereas a Polynucleobacter population displayed strong changes within merely 2 months (FST up to 0.54), similar in scale to differentiation between populations separated by thousands of kilometers. The slowly and steadily evolving LD12 population showed high strain diversity, whereas the dynamic Polynucleobacter population exhibited alternating clonal expansions of mostly two strains only. Based on the contrasting population structures, we propose distinct models of speciation.

Genetic differentiation and structuration of the Gobra zebu cattle breeds reared in Senegal

The Gobra zebu genetic breeding program has resulted in the genetic improvement of a new population. This population gained genetic characteristics that set them apart from the other cattle populations reared in Senegal. The cause of these differences might be the reproductive isolation and selection to which this population of the "Centre de Recherches Zootechniques" of Dahra has been subjected since the 1950s. This study aimed to assess the genetic differentiation and structuration of this population compared to the main cattle breeds used in Senegal. A total of 180 individuals, selected from the Gobra selection nucleus and bovine populations from four main breeds in Senegal, were included in this study. We used a panel of 21 microsatellite markers among those recommended by the Food Agriculture Organization, to conduct the molecular genotyping of our sampled populations. The basic genetic parameters of differentiation and structuration were calculated using various bioinformatics software. The results of this study, particularly the degrees of genetic differentiation (Fst), the coefficient of genetic homogeneity (Gst), and the gene flow (Nm), show a significant genetic differentiation of the Gobra from the station compared to the other populations studied. Structuring and phylogeny analyses reveal a micro-structuring within the Gobra population as a novelty. This micro-structuring clearly identifies the Gobra individuals from Dahra's station among the other Gobra populations studied. The main causes of these observations would be reproductive isolation and the selection pressure exerted on this population for several decades.

Candidate Genes and Gene Networks Change with Age in Japanese Black Cattle by Blood Transcriptome Analysis

Age is an important physiological factor that affects the metabolism and immune function of beef cattle. While there have been many studies using the blood transcriptome to study the effects of age on gene expression, few have been reported on beef cattle. To this end, we used the blood transcriptomes of Japanese black cattle at different ages as the study subjects and screened 1055, 345, and 1058 differential expressed genes (DEGs) in the calf vs. adult, adult vs. old, and calf vs. old comparison groups, respectively. The weighted co-expression network consisted of 1731 genes. Finally, blue, brown, and yellow age-specific modules were obtained, in which genes were enriched in signaling pathways related to growth and development and immune metabolic dysfunction, respectively. Protein-protein interaction (PPI) analysis showed gene interactions in each specific module, and 20 of the highest connectivity genes were chosen as potential hub genes. Finally, we identified 495, 244, and 1007 genes by exon-wide selection signature (EWSS) analysis of different comparison groups. Combining the results of hub genes, we found that VWF, PARVB, PRKCA, and TGFB1I1 could be used as candidate genes for growth and development stages of beef cattle. CORO2B and SDK1 could be used as candidate marker genes associated with aging. In conclusion, by comparing the blood transcriptome of calves, adult cattle, and old cattle, the candidate genes related to immunity and metabolism affected by age were identified, and the gene co-expression network of different age stages was constructed. It provides a data basis for exploring the growth, development, and aging of beef cattle.