Deleterious protein-coding variants in diverse cattle breeds of the world

Growth resilience to weather variation in commercial free-ranging chickens in Ethiopia

BACKGROUND

The poultry industry in sub-Saharan Africa is a rapidly developing sector mostly based on smallholder farming. Increased demand for poultry-derived products, driven by the growing economy and population, has intensified importations of highly productive exotic breeds and crossbreeding with local ecotypes. However, commercial chickens with exotic genes often struggle to adapt to the local climate under smallholder farmers management. Understanding the chicken response to weather changes is crucial for developing selection schemes that ensure proper adaptation. In the present study, we derived individual phenotypes for growth resilience of commercial free-ranging chickens to changing weather conditions in Ethiopia. In addition, we performed genomic association analyses to assess the genetic background of these phenotypes and identify potential candidate genes of interest.

RESULTS

Novel resilience phenotypes describing changes in chicken growth profiles in response to weather fluctuation were developed. Variations in daily air temperature, relative humidity and amount of precipitation had the strongest impact on growth. Significant genomic variance was detected for growth resilience to changes in air temperature measurements and a temperature-humidity index. Genomic markers correlated with these resilience traits were mostly located within or near candidate genes associated with lipid metabolism and adipocyte homeostasis. Some of these genes have been previously linked to animal responses to environmental stressors in other species.

CONCLUSIONS

The phenotypes of growth resilience of chickens to changing weather conditions exhibited significant genomic variation. The outcomes of this study may facilitate the genomic selection of commercial chickens that are not only highly productive, but also capable of maintaining their production levels under varying weather conditions.

ZFAT (isoform-specific) and its antisense RNA 1 (ZFAT-AS1) are two allele-specific monoallelically expressed genes in cattle

In mammals, imprinted genes are characterised by a monoallelic expression, which is based on parental origin and is essential for both foetal and placental development. The ZFAT gene encodes a transcriptional factor, and its non-coding antisense RNA, ZFAT-AS1, overlaps with the ZFAT locus. Both ZFAT and ZFAT-AS1 are maternally imprinted in human placentas. In bovines, the imprinting status of the ZFAT and ZFAT-AS1 genes has yet to be reported. In this study, we analysed the allelic expression of three transcript variants (X1-X3) of the bovine ZFAT and ZFAT-AS1 genes in somatic tissues and placentas using a single nucleotide polymorphism-based method. The results showed that bovine ZFAT exhibited isoform-specific paternal expression. The ZFAT X2 variant exhibited monoallelic expression in the bovine placentas and biallelic expression in the six bovine somatic tissues (heart, liver, spleen, lung, kidney and brain). However, the ZFAT X1 and X3 variants were biallelically expressed in both bovine tissues and placentas. A 311 bp bovine ZFAT-AS1 complementary DNA (cDNA) sequence was obtained by aligning the human ZFAT-AS1 cDNA sequence with the bovine genome and conducting reverse transcription polymerase chain reaction amplification. Bovine ZFAT-AS1 have monoallelic expression in bovine placentas and somatic tissues. In addition, the DNA methylation of two regions was characterised, including the partial promoter, and exon 1 and intron 1 regions of ZFAT, and there were no differentially methylated regions.

Purifying Selection Influences the Comparison of Heterozygosities between Populations

Heterozygosity is a fundamental measure routinely used to compare between populations to infer the level of genetic variation and their relative effective population sizes. However, such comparison is highly influenced by the magnitude of selection pressure on the genomic regions used. Using over 2 million Single Nucleotide Variants (SNVs) from chimpanzee and mouse populations, this study shows that the heterozygosities estimated using neutrally evolving sites of large populations were two times higher than those of small populations. However, this difference was only ~1.6 times for the heterozygosities estimated using nonsynonymous sites. This suggests an excess in the nonsynonymous heterozygosities due to the segregation of deleterious variants in small populations. This excess in the nonsynonymous heterozygosities of the small populations was estimated to be 23-31%. Further analysis revealed that the magnitude of the excess is modulated by effective population size (Ne) and selection intensity (s). Using chimpanzee populations, this investigation found that the excess in nonsynonymous diversity in the small population was little (6%) when the difference between the Ne values of large and small populations was small (2.4 times). Conversely, this was high (23%) when the difference in Ne was large (5.9 times). Analysis using mouse populations showed that the excess in the nonsynonymous diversity of highly constrained genes of the small population was much higher (38%) than that observed for the genes under relaxed selective constraints (21%). Similar results were observed when the expression levels of genes were used as a proxy for selection intensity. These results emphasize the use of neutral regions, less constrained genes, or lowly expressed genes when comparing the heterozygosities between populations.

Deleterious mutation load in the admixed mice population

Deleterious mutation loads are known to correlate negatively with effective population size (Ne). Due to this reason, previous studies observed a higher proportion of harmful mutations in small populations than that in large populations. However, the mutational load in an admixed population that derived from introgression between individuals from two populations with vastly different Ne is not known. We investigated this using the whole genome data from two subspecies of the mouse (Mus musculus castaneus and Mus musculus musculus) with significantly different Ne. We used the ratio of diversities at nonsynonymous and synonymous sites (dN/dS) to measure the harmful mutation load. Our results showed that this ratio observed for the admixed population was intermediate between those of the parental populations. The dN/dS ratio of the hybrid population was significantly higher than that of M. m. castaneus but lower than that of M. m. musculus. Our analysis revealed a significant positive correlation between the proportion of M. m. musculus ancestry in admixed individuals and their dN/dS ratio. This suggests that the admixed individuals with high proportions of M. m. musculus ancestry have large dN/dS ratios. We also used the proportion of deleterious nonsynonymous SNVs as a proxy for deleterious mutation load, which also produced similar results. The observed results were in concordance with those expected by theory. We also show a shift in the distribution of fitness effects of nonsynonymous SNVs in the admixed genomes compared to the parental populations. These findings suggest that the deleterious mutation load of the admixed population is determined by the proportion of the ancestries of the subspecies. Therefore, it is important to consider the status and the level of genetic admixture of the populations whilst estimating the mutation loads.

Comparative transcriptome in large-scale human and cattle populations

BACKGROUND

Cross-species comparison of transcriptomes is important for elucidating evolutionary molecular mechanisms underpinning phenotypic variation between and within species, yet to date it has been essentially limited to model organisms with relatively small sample sizes.

RESULTS

Here, we systematically analyze and compare 10,830 and 4866 publicly available RNA-seq samples in humans and cattle, respectively, representing 20 common tissues. Focusing on 17,315 orthologous genes, we demonstrate that mean/median gene expression, inter-individual variation of expression, expression quantitative trait loci, and gene co-expression networks are generally conserved between humans and cattle. By examining large-scale genome-wide association studies for 46 human traits (average n = 327,973) and 45 cattle traits (average n = 24,635), we reveal that the heritability of complex traits in both species is significantly more enriched in transcriptionally conserved than diverged genes across tissues.

CONCLUSIONS

In summary, our study provides a comprehensive comparison of transcriptomes between humans and cattle, which might help decipher the genetic and evolutionary basis of complex traits in both species.