Genome-wide assessment and mapping of inbreeding depression identifies candidate genes associated with semen traits in Holstein bulls

BACKGROUND

The reduction in phenotypic performance of a population due to mating between close relatives is called inbreeding depression. The genetic background of inbreeding depression for semen traits is poorly understood. Thus, the objectives were to estimate the effect of inbreeding and to identify genomic regions underlying inbreeding depression of semen traits including ejaculate volume (EV), sperm concentration (SC), and sperm motility (SM). The dataset comprised ~ 330 K semen records from ~ 1.5 K Holstein bulls genotyped with 50 K single nucleotide polymorphism (SNP) BeadChip. Genomic inbreeding coefficients were estimated using runs of homozygosity (i.e., F_ROH > 1 Mb) and excess of SNP homozygosity (F_SNP). The effect of inbreeding was estimated by regressing phenotypes of semen traits on inbreeding coefficients. Associated variants with inbreeding depression were also detected by regressing phenotypes on ROH state of the variants.

RESULTS

Significant inbreeding depression was observed for SC and SM (p < 0.01). A 1% increase in F_ROH reduced SM and SC by 0.28% and 0.42% of the population mean, respectively. By splitting F_ROH into different lengths, we found significant reduction in SC and SM due to longer ROH, which is indicative of more recent inbreeding. A genome-wide association study revealed two signals positioned on BTA 8 associated with inbreeding depression of SC (p < 0.00001; FDR < 0.02). Three candidate genes of GALNTL6, HMGB2, and ADAM29, located in these regions, have established and conserved connections with reproduction and/or male fertility. Moreover, six genomic regions on BTA 3, 9, 21 and 28 were associated with SM (p < 0.0001; FDR < 0.08). These genomic regions contained genes including PRMT6, SCAPER, EDC3, and LIN28B with established connections to spermatogenesis or fertility.

CONCLUSIONS

Inbreeding depression adversely affects SC and SM, with evidence that longer ROH, or more recent inbreeding, being especially detrimental. There are genomic regions associated with semen traits that seems to be especially sensitive to homozygosity, and evidence to support some from other studies. Breeding companies may wish to consider avoiding homozygosity in these regions for potential artificial insemination sires.

A genome-wide perspective on the diversity and selection signatures in indigenous goats using 53 K single nucleotide polymorphism array

Tibetan goats, Taihang goats, Jining grey goats, and Meigu goats are the representative indigenous goats in China, found in Qinghai-Tibet Plateau, Western pastoral area, Northern and Southern agricultural regions. Very few studies have conducted a comprehensive analysis of the genomic diversity and selection of these breeds. We genotyped 96 unrelated individuals, using goat 53 K Illumina BeadChip array, of the following goat breeds: Tibetan (TG), Taihang (THG), Jining grey (JGG), and Meigu (MGG). A total of 45 951 single nucleotide polymorphisms were filtered to estimate the genetic diversity and selection signatures. All breeds had a high proportion (over 95%) of polymorphic loci. The observed and excepted heterozygosity ranged from 0.338 (MGG) to 0.402 (JGG) and 0.339 (MGG) to 0.395 (JGG), respectively. Clustering analysis displayed a genetically distinct lineage for each breed, and their Fst were greater than 0.25, indicating that they had a higher genetic differentiation between groups. Furthermore, effective population size reduced in all four populations, indicating a loss of genetic diversity. In addition, runs of homozygosity were mainly distributed in 5-10 Mb. Lastly, we identified signature genes, which were closely related to high-altitude adaptation (ADIRF) and prolificity (CNTROB, SMC3, and PTEN). This study provides a valuable resource for future studies on genome-wide perspectives on the diversity and selection signatures of Chinese indigenous goats.

Interchromosomal template-switching as a novel molecular mechanism for imprinting perturbations associated with Temple syndrome

Background

Intrachromosomal triplications (TRP) can contribute to disease etiology via gene dosage effects, gene disruption, position effects, or fusion gene formation. Recently, post-zygotic de novo triplications adjacent to copy-number neutral genomic intervals with runs of homozygosity (ROH) have been shown to result in uniparental isodisomy (UPD). The genomic structure of these complex genomic rearrangements (CGRs) shows a consistent pattern of an inverted triplication flanked by duplications (DUP-TRP/INV-DUP) formed by an iterative DNA replisome template-switching mechanism during replicative repair of a single-ended, double-stranded DNA (seDNA), the ROH results from an interhomolog or nonsister chromatid template switch. It has been postulated that these CGRs may lead to genetic abnormalities in carriers due to dosage-sensitive genes mapping within the copy-number variant regions, homozygosity for alleles at a locus causing an autosomal recessive (AR) disease trait within the ROH region, or imprinting-associated diseases.

Methods

Here, we report a family wherein the affected subject carries a de novo 2.2-Mb TRP followed by 42.2 Mb of ROH and manifests clinical features overlapping with those observed in association with chromosome 14 maternal UPD (UPD(14)mat). UPD(14)mat can cause clinical phenotypic features enabling a diagnosis of Temple syndrome. This CGR was then molecularly characterized by high-density custom aCGH, genome-wide single-nucleotide polymorphism (SNP) and methylation arrays, exome sequencing (ES), and the Oxford Nanopore long-read sequencing technology.

Results

We confirmed the postulated DUP-TRP/INV-DUP structure by multiple orthogonal genomic technologies in the proband. The methylation status of known differentially methylated regions (DMRs) on chromosome 14 revealed that the subject shows the typical methylation pattern of UPD(14)mat. Consistent with these molecular findings, the clinical features overlap with those observed in Temple syndrome, including speech delay.

Conclusions

These data provide experimental evidence that, in humans, triplication can lead to segmental UPD and imprinting disease. Importantly, genotype/phenotype analyses further reveal how a post-zygotically generated complex structural variant, resulting from a replication-based mutational mechanism, contributes to expanding the clinical phenotype of known genetic syndromes. Mechanistically, such events can distort transmission genetics resulting in homozygosity at a locus for which only one parent is a carrier as well as cause imprinting diseases.

Electronic supplementary material

The online version of this article (10.1186/s13073-019-0633-y) contains supplementary material, which is available to authorized users.

High-resolution population structure and runs of homozygosity reveal the genetic architecture of complex traits in the Lipizzan horse

BACKGROUND

The sample ascertainment bias due to complex population structures remains a major challenge in genome-wide investigations of complex traits. In this study we derived the high-resolution population structure and levels of autozygosity of 377 Lipizzan horses originating from five different European stud farms utilizing the SNP genotype information of the high density 700 k Affymetrix Axiom™ Equine genotyping array. Scanning the genome for overlapping runs of homozygosity (ROH) shared by more than 50% of horses, we identified homozygous regions (ROH islands) in order to investigate the gene content of those candidate regions by gene ontology and enrichment analyses.

RESULTS

The high-resolution population network approach revealed well-defined substructures according to the origin of the horses (Austria, Slovakia, Croatia and Hungary). The highest mean genome coverage of ROH (SROH) was identified in the Austrian (SROH = 342.9), followed by Croatian (SROH = 214.7), Slovakian (SROH = 205.1) and Hungarian (SROH = 171.5) subpopulations. ROH island analysis revealed five common islands on ECA11 and ECA14, hereby confirming a closer genetic relationship between the Hungarian and Croatian as well as between the Austrian and Slovakian samples. Private islands were detected for the Hungarian and the Austrian Lipizzan subpopulations. All subpopulations shared a homozygous region on ECA11, nearly identical in position and length containing among other genes the homeobox-B cluster, which was also significantly (p < 0.001) highlighted by enrichment analysis. Gene ontology terms were mostly related to biological processes involved in embryonic morphogenesis and anterior/posterior specification. Around the STX17 gene (causative for greying), we identified a ROH island harbouring the genes NR4A3, STX17, ERP44 and INVS. Within further islands on ECA14, ECA16 and ECA20 we detected the genes SPRY4, NDFIP1, IMPDH2, HSP90AB1, whereas SPRY4 and HSP90AB1 are involved in melanoma metastasis and survival rate of melanoma patients in humans.

CONCLUSIONS

We demonstrated that the assessment of high-resolution population structures within one single breed supports the downstream genetic analyses (e.g. the identification of ROH islands). By means of ROH island analyses, we identified the genes SPRY4, NDFIP1, IMPDH2, HSP90AB1, which might play an important role for further studies on equine melanoma. Furthermore, our results highlighted the impact of the homeobox-A and B cluster involved in morphogenesis of Lipizzan horses.

A genetic assessment of the English bulldog

Background

This study examines genetic diversity among 102 registered English Bulldogs used for breeding based on maternal and paternal haplotypes, allele frequencies in 33 highly polymorphic short tandem repeat (STR) loci on 25 chromosomes, STR-linked dog leukocyte antigen (DLA) class I and II haplotypes, and the number and size of genome-wide runs of homozygosity (ROH) determined from high density SNP arrays. The objective was to assess whether the breed retains enough genetic diversity to correct the genotypic and phenotypic abnormalities associated with poor health, to allow for the elimination of deleterious recessive mutations, or to make further phenotypic changes in body structure or coat. An additional 37 English bulldogs presented to the UC Davis Veterinary Clinical Services for health problems were also genetically compared with the 102 registered dogs based on the perception that sickly English bulldogs are products of commercial breeders or puppy-mills and genetically different and inferior.

Results

Four paternal haplotypes, with one occurring in 93 % of dogs, were identified using six Y-short tandem repeat (STR) markers. Three major and two minor matrilines were identified by mitochondrial D-loop sequencing. Heterozygosity was determined from allele frequencies at genomic loci; the average number of alleles per locus was 6.45, with only 2.7 accounting for a majority of the diversity. However, observed and expected heterozygosity values were nearly identical, indicating that the population as a whole was in Hardy-Weinberg equilibrium (HWE). However, internal relatedness (IR) and adjusted IR (IRVD) values demonstrated that a number of individuals were the offspring of parents that were either more inbred or outbred than the population as a whole. The diversity of DLA class I and II haplotypes was low, with only 11 identified DLA class I and nine class II haplotypes. Forty one percent of the breed shared a single DLA class I and 62 % a single class II haplotype. Nineteen percent of the dogs were homozygous for the dominant DLA class I haplotype and 42 % for the dominant DLA class II haplotype. The extensive loss of genetic diversity is most likely the result of a small founder population and artificial genetic bottlenecks occurring in the past. The prominent phenotypic changes characteristic of the breed have also resulted in numerous large runs of homozygosity (ROH) throughout the genome compared to Standard Poodles, which were phenotypically more similar to indigenous-type dogs.

Conclusions

English bulldogs have very low genetic diversity resulting from a small founder population and artificial genetic bottlenecks. Although some phenotypic and genotypic diversity still exists within the breed, whether it is sufficient to use reverse selection to improve health, select against simple recessive deleterious traits, and/or to accommodate further genotypic/phenotypic manipulations without further decreasing existing genetic diversity is questionable.