Genome-wide association for testis weight in the diversity outbred mouse population

Genetic Diversity, Selection Signatures, and Genome-Wide Association Study Identify Candidate Genes Related to Litter Size in Hu Sheep

Hu sheep is a renowned prolific local sheep breed in China, widely distributed across the country due to its excellent reproductive performance. Deciphering the molecular mechanisms underlying the high fecundity of Hu sheep is crucial for improving the litter size of ewes. In this study, we genotyped 830 female Hu sheep using the Illumina OvineSNP50 BeadChip and performed genetic diversity analysis, selection signature detection, and a genome-wide association study (GWAS) for litter size. Our results revealed that the Hu sheep population exhibits relatively high genetic diversity. A total of 4927 runs of homozygosity (ROH) segments were detected, with the majority (74.73%) being short in length. Different genomic inbreeding coefficients (FROH, FHOM, FGRM, and FUNI) ranged from -0.0060 to 0.0126, showing low levels of inbreeding in this population. Additionally, we identified 91 candidate genomic regions through three complementary selection signature methods, including ROH, composite likelihood ratio (CLR), and integrated haplotype score (iHS), and annotated 189 protein-coding genes. Moreover, we observed two significant SNPs related to the litter size of Hu sheep using GWAS analysis based on a repeatability model. Integrating the selection signatures and the GWAS results, we identified 15 candidate genes associated with litter size, among which BMPR1B and UNC5C were particularly noteworthy. These findings provide valuable insights for improving the reproductive performance and breeding of high-fecundity lines of Hu sheep.

Molecular characterization of SPATA6 and association of its SNPs with testicular size in sheep

The testis is an important organ for maintaining fertility in males, and testis size is positively correlated with ejaculate volume, sperm motility, thus fertility. Spermatogenesis-associated 6 (SPATA6) is an evolutionarily conserved testis-specific gene reported in many species. However, the effect of SPATA6 expression levels on testicular development and the effect of single nucleotide polymorphisms (SNPs) on testis and epididymis phenotype in sheep have not been studied. The purpose of the research was to investigate the expression profile of SPATA6 and its effect on testicular development and to confirm the effect of SNPs on the testis and epididymis phenotype. In this study, we detected a 1245bp coding sequence (CDS) of SPATA6 and encoded 414 amino acids. The expression levels of SPATA6 were significantly higher in the testis than in other tissues and gradually increased with testis development. Moreover, the expression level in the large testis was significantly higher than that in the small testis at six months. A total of 11 SNPs were detected in the coding region of SPATA6 by cDNA-pooling sequencing and improved multiplex ligation detection reaction (iMLDR) methods. Correlation analysis showed that SNP2 (c. 3631C > G) significantly affected left epididymis weight (LEW) and right epididymis weight (REW), and SNP10 (c. 937 A > G) significantly affected REW. And the combined genotype of SNP1 (c. 4245 G > A) and SNP2 significantly affected REW. The current study concluded that SPATA6 plays an important role in testicular development and the SNPs significantly associated with the epididymis phenotype can provide molecular markers for the early selection of high-fertility Hu sheep.

17α-estradiol does not adversely affect sperm parameters or fertility in male mice: implications for reproduction-longevity trade-offs

17α-estradiol (17α-E2) is referred to as a nonfeminizing estrogen that was recently found to extend healthspan and lifespan in male, but not female, mice. Despite an abundance of data indicating that 17α-E2 attenuates several hallmarks of aging in male rodents, very little is known with regard to its effects on feminization and fertility. In these studies, we evaluated the effects of 17α-E2 on several markers of male reproductive health in two independent cohorts of mice. In alignment with our previous reports, chronic 17α-E2 treatment prevented gains in body mass, but did not adversely affect testes mass or seminiferous tubule morphology. We subsequently determined that chronic 17α-E2 treatment also did not alter plasma 17β-estradiol or estrone concentrations, while mildly increasing plasma testosterone levels. We also determined that chronic 17α-E2 treatment did not alter plasma follicle-stimulating hormone or luteinizing hormone concentrations, which suggests 17α-E2 treatment does not alter gonadotropin-releasing hormone neuronal function. Sperm quantity, morphology, membrane integrity, and various motility measures were also unaffected by chronic 17α-E2 treatment in our studies. Lastly, two different approaches were used to evaluate male fertility in these studies. We found that chronic 17α-E2 treatment did not diminish the ability of male mice to impregnate female mice, or to generate successfully implanted embryos in the uterus. We conclude that chronic treatment with 17α-E2 at the dose most commonly employed in aging research does not adversely affect reproductive fitness in male mice, which suggests 17α-E2 does not extend lifespan or curtail disease parameters through tradeoff effects with reproduction.

Connecting the Dots: N6-Methyladenosine (m6 A) Modification in Spermatogenesis

N6-methyladenosine (m6 A) is the most common RNA modification found in eukaryotes and is involved in multiple biological processes, including neuronal development, tumorigenesis, and gametogenesis. It is well known that methylation-modifying enzymes (classified into writers, erasers, and readers) mediate catalysis, clearance, and recognition of m6 A. Recent studies suggest that these genes may be associated with spermatogenesis. Numerous studies have revealed the m6 A role during spermatogenesis. However, the expression patterns and relationships of these m6 A enzymes during various stages of spermatogenesis remain unknown. In this review, it is aimed to provide an overview of m6 A enzyme functions and elucidate their potential mechanisms and regulatory relationships at a specific phase during spermatogenesis, providing new insights into the m6 A modification underlying the spermatogenesis process.

Genome-wide association study reveals candidate markers related to field fertility and semen quality traits in Holstein-Friesian bulls

In vitro assessment of bull semen quality is routinely used in bull semen processing centres in order to ensure that semen destined to be used in the field has passed minimum standards. Despite these stringent quality control checks, individual bulls that pass the quality control checks can still vary in field fertility by up to 25%. A genome-wide association study was undertaken to determine genetic markers associated with prefreeze and post-thaw bull sperm quality traits as well as field fertility. Genome-wide association analysis was performed using a single nucleotide polymorphism (SNP) regression mixed linear model in WOMBAT. Genes within a 250 Kb span of a suggestive (P ≤ 1 × 10^-5) SNP were considered as candidate genes. One SNP was associated with adjusted pregnancy rate, and 21 SNPs were associated across the seven semen quality traits (P ≤ 1 × 10^-5). Functional candidate genes include SIPA1L2 which was associated with adjusted pregnancy rate. This encodes a Rap GTPase-activating protein involved in Rap1 signalling pathway and was previously found to play a role in the process of sperm differentiation. Gene ontology (GO) analysis also identified significantly enriched biological processes involved protein tyrosine kinase activity including genes such as DYRK1A, TEC and TXK that were associated with sperm motility prior to freezing. Another candidate gene associated with post-thaw sperm motility was FHDC1 which coordinates actin filament and microtubule dynamics. The induced 11 GO terms in the ejaculates rejected after freezing trait were related to ATPase, phosphatase and hydrolase activity. These results reveal novel specific genomic regions and candidate genes associated with economically important phenotypes such as field fertility and semen quality traits.

Reproductive genomics of the mouse: implications for human fertility and infertility

Genetic analyses of mammalian gametogenesis and fertility have the potential to inform about two important and interrelated clinical areas: infertility and contraception. Here, we address the genetics and genomics underlying gamete formation, productivity and function in the context of reproductive success in mammalian systems, primarily mouse and human. Although much is known about the specific genes and proteins required for meiotic processes and sperm function, we know relatively little about other gametic determinants of overall fertility, such as regulation of gamete numbers, duration of gamete production, and gamete selection and function in fertilization. As fertility is not a binary trait, attention is now appropriately focused on the oligogenic, quantitative aspects of reproduction. Multiparent mouse populations, created by complex crossing strategies, exhibit genetic diversity similar to human populations and will be valuable resources for genetic discovery, helping to overcome current limitations to our knowledge of mammalian reproductive genetics. Finally, we discuss how what we know about the genomics of reproduction can ultimately be brought to the clinic, informing our concepts of human fertility and infertility, and improving assisted reproductive technologies.

A mutation modulating DDX3Y gene expression cosegregates with the major Y-chromosomal haplogroups and with testis size in Hu sheep

Variations in the Y-chromosome are usually correlated with male-specific traits. However, this condition has been described only sporadically, even in human genetics. The present study was conducted to clone the full-length gene sequence of ovine DEAD-box helicase 3, Y-linked (DDX3Y), and investigate the effect of the expression and variation of DDX3Y on the reproductive traits of Hu sheep. Consequently, we identified the full coding sequence and genomic sequence of ovine DDX3Y. Quantitative PCR (qPCR) analysis showed that ovine DDX3Y was highly expressed in testis, and the expression level increased during testicular development. Furthermore, individuals with larger testis at 6 months expressed significantly more DDX3Y mRNA in the testis than individuals with smaller testis. Notably, a novel SNP (g. 12657 C>A) in the 3' untranslated region was identified in Hu sheep and Tan sheep according to the investigation of the full DDX3Y genomic sequence of 1069 individuals from nine sheep breeds. Association analysis revealed that the SNP was significantly related to testis size in Hu sheep. Meanwhile, Hu rams with the derived C allele showed significantly higher expression levels of DDX3Y in testis than those with the ancestral A allele. In addition, data mining in a previous study showed that the C allele cosegregated with the globally major Y-chromosomal haplogroups y-HA and y-HC, and the A allele is found in all rams with haplogroups y-HB1, y-HB2 and y-HD. This study suggests that the association of the Y-chromosomal haplogroups with testis size in Hu sheep can be extrapolated to the sheep population worldwide.

Genes underlying the evolution of tetrapod testes size

BACKGROUND

Testes vary widely in mass relative to body mass across species, but we know very little about which genes underlie and contribute to such variation. This is partly because evidence for which genes are implicated in testis size variation tends to come from investigations involving just one or a few species. Contemporary comparative phylogenetic methods provide an opportunity to test candidate genes for their role in phenotypic change at a macro-evolutionary scale-across species and over millions of years. Previous attempts to detect genotype-phenotype associations across species have been limited in that they can only detect where genes have driven directional selection (e.g. brain size increase).

RESULTS

Here, we introduce an approach that uses rates of evolutionary change to overcome this limitation to test whether any of twelve candidate genes have driven testis size evolution across tetrapod vertebrates-regardless of directionality. We do this by seeking a relationship between the rates of genetic and phenotypic evolution. Our results reveal five genes (Alkbh5, Dmrtb1, Pld6, Nlrp3, Sp4) that each have played unique and complex roles in tetrapod testis size diversity. In all five genes, we find strong significant associations between the rate of protein-coding substitutions and the rate of testis size evolution. Such an association has never, to our knowledge, been tested before for any gene or phenotype.

CONCLUSIONS

We describe a new approach to tackle one of the most fundamental questions in biology: how do individual genes give rise to biological diversity? The ability to detect genotype-phenotype associations that have acted across species has the potential to build a picture of how natural selection has sculpted phenotypic change over millions of years.

Heritability of fat distributions in male mice from the founder strains of the Diversity Outbred mouse population

Specific fat distributions are risk factors for complex diseases, including coronary heart disease and obstructive sleep apnea. To demonstrate the utility of high-diversity mouse models for elucidating genetic associations, we describe the phenotyping and heritability of fat distributions within the five classical inbred and three wild-derived founder mouse strains of the Collaborative Cross and Diversity Outbred mice. Measurements of subcutaneous and internal fat volumes in the abdomen, thorax and neck, and fat volumes in the tongue and pericardium were obtained using magnetic resonance imaging in male mice from the A/J (n = 12), C57BL/6J (n = 17), 129S1/SvlmJ (n = 12), NOD/LtJ (n = 14), NZO/HILtJ (n = 12), CAST/EiJ (n = 14), PWK/PhJ (n = 12), and WSB/EiJ (n = 15) strains. Phenotypes were compared across strains using analysis of variance and heritability estimated as the proportion of phenotypic variability attributable to strain. Heritability ranged from 44 to 91% across traits, including >70% heritability of tongue fat. A majority of heritability estimates remained significant controlling for body weight, suggesting genetic influences independent of general obesity. Principal components analysis supports genetic influences on overall obesity and specific to increased pericardial and intra-neck fat. Thus, among the founder strains of the Collaborative Cross and Diversity Outbred mice, we observed significant heritability of subcutaneous and internal fat volumes in the neck, thorax and abdomen, pericardial fat volume and tongue fat volume, consistent with genetic architecture playing an important role in explaining trait variability. Findings pave the way for studies utilizing high-diversity mouse models to identify genes affecting fat distributions and, in turn, influencing risk for associated complex disorders.

High-throughput sleep phenotyping produces robust and heritable traits in Diversity Outbred mice and their founder strains

STUDY OBJECTIVES

This study describes high-throughput phenotyping strategies for sleep and circadian behavior in mice, including examinations of robustness, reliability, and heritability among Diversity Outbred (DO) mice and their eight founder strains.

METHODS

We performed high-throughput sleep and circadian phenotyping in male mice from the DO population (n = 338) and their eight founder strains: A/J (n = 6), C57BL/6J (n = 14), 129S1/SvlmJ (n = 6), NOD/LtJ (n = 6), NZO/H1LtJ (n = 6), CAST/EiJ (n = 8), PWK/PhJ (n = 8), and WSB/EiJ (n = 6). Using infrared beam break systems, we defined sleep as at least 40 s of continuous inactivity and quantified sleep-wake amounts and bout characteristics. We developed assays to measure sleep latency in a new environment and during a modified Murine Multiple Sleep Latency Test, and estimated circadian period from wheel-running experiments. For each trait, broad-sense heritability (proportion of variability explained by all genetic factors) was derived in founder strains, while narrow-sense heritability (proportion of variability explained by additive genetic effects) was calculated in DO mice.

RESULTS

Phenotypes were robust to different inactivity durations to define sleep. Differences across founder strains and moderate/high broad-sense heritability were observed for most traits. There was large phenotypic variability among DO mice, and phenotypes were reliable, although estimates of heritability were lower than in founder mice. This likely reflects important nonadditive genetic effects.

CONCLUSIONS

A high-throughput phenotyping strategy in mice, based primarily on monitoring of activity patterns, provides reliable and heritable estimates of sleep and circadian traits. This approach is suitable for discovery analyses in DO mice, where genetic factors explain some proportion of phenotypic variation.