Studies on the Relationships between Growth and Gonad Development during First Sexual Maturation of Macrobrachium nipponense and Associated SNPs Screening.

In this study, we used full-sib families to investigate the association between growth and gonad development during first sexual maturation of M. nipponense. We found that male GSI was significantly negatively correlated with growth traits (p < 0.01) and there were no significant correlations between female GSI (Gonadosomatic index) and growth traits (p > 0.05). HSI (Hepatopancreas index) in both males and females showed no significant correlations with growth traits (p > 0.05). We furthermore investigated the association between the specific allele of Mn-CTS L1 polymorphism and gonad development and growth traits. In total, 35 mutation loci were screened and 16 high-quality single-nucleotide polymorphisms (SNPs) loci were obtained after validation. Four and two SNPs proved to be strongly associated with all growth traits in female and male M. nipponense separately, among which A+118T might be a candidate SNP positively associated with large growth traits. Two and one SNPs were screened, respectively, in males and females to associate with GSI, while three SNPs were detected to associate with female HSI, among which A+1379C may be applied as a potential molecular marker for gene-assisted selection to improve both reproduction speed and growth traits in M. nipponense.

Sex-specific overdominance at the maturation vgll3 gene for reproductive fitness in wild Atlantic salmon.

Linking reproductive fitness with adaptive traits at the genomic level can shed light on the mechanisms that produce and maintain sex-specific selection. Here, we construct a multigenerational pedigree to investigate sex-specific selection on a maturation gene, vgll3, in a wild Atlantic salmon population. The vgll3 locus is responsible for ~40% of the variation in maturation (sea age at first reproduction). Genetic parentage analysis was conducted on 18,265 juveniles (parr) and 685 adults collected at the same spawning ground over eight consecutive years. A high proportion of females (26%) were iteroparous and reproduced two to four times in their lifetime. A smaller proportion of males (9%) spawned at least twice in their lifetime. Sex-specific patterns of reproductive fitness were related to vgll3 genotype. Females showed a pattern of overdominance where vgll3*EL genotypes had three-fold more total offspring than homozygous females. In contrast, males demonstrated that late-maturing vgll3*LL individuals had two-fold more offspring than either vgll3*EE or vgll3*EL males. Taken together, these data suggest that balancing selection in females contributes to the maintenance of variation at this locus via increased fitness of iteroparous vgll3*EL females. This study demonstrates the utility of multigenerational pedigrees for uncovering complex patterns of reproduction, sex-specific selection and the maintenance of genetic variation.

Deregulation of oxidative phosphorylation pathways in embryos derived in vitro from prepubertal and pubertal heifers based on whole-transcriptome sequencing.

BACKGROUND: Although, oocytes from prepubertal donors are known to be less developmentally competent than those from adult donors it does not restrain their ability to produce full-term pregnancies. The transcriptomic profile of embryos could be used as a predictor for embryo's individual developmental competence. The aim of the study was to compare transcriptomic profile of blastocysts derived from prepubertal and pubertal heifers oocytes. Bovine cumulus-oocyte complexes (COCs) were obtained by ovum pick- up method from prepubertal and pubertal heifers. After in vitro maturation COCs were fertilized and cultured to the blastocyst stage. Total RNA was isolated from both groups of blastocysts and RNA-seq was performed. Gene ontology analysis was performed by DAVID (Database for Annotation, Visualization and Integrated Discovery). RESULTS: A higher average blastocyst rate was obtained in the pubertal than in the pre-pubertal group. There were no differences in the quality of blastocysts between the examined groups. We identified 436 differentially expressed genes (DEGs) between blastocysts derived from researched groups, of which 247 DEGs were downregulated in blastocysts derived from pubertal compared to prepubertal heifers oocytes, and 189 DEGs were upregulated. The genes involved in mitochondrial function, including oxidative phosphorylation (OXPHOS) were found to be different in studied groups using Kyoto Encyclopedia of Genes (KEGG) pathway analysis and 8 of those DEGs were upregulated and 1 was downregulated in blastocysts derived from pubertal compared to prepubertal heifers oocytes. DEGs associated with mitochondrial function were found: ATP synthases (ATP5MF-ATP synthase membrane subunit f, ATP5PD- ATP synthase peripheral stalk subunit d, ATP12A- ATPase H+/K + transporting non-gastric alpha2 subunit), NADH dehydrogenases (NDUFS3- NADH: ubiquinone oxidoreductase subunit core subunit S3, NDUFA13- NADH: ubiquinone oxidoreductase subunit A13, NDUFA3- NADH: ubiquinone oxidoreductase subunit A3), cytochrome c oxidase (COX17), cytochrome c somatic (CYCS) and ubiquinol cytochrome c reductase core protein 1 (UQCRC1). We found lower number of apoptotic cells in blastocysts derived from oocytes collected from prepubertal than those obtained from pubertal donors. CONCLUSIONS: Despite decreased expression of genes associated with OXPHOS pathway in blastocysts from prepubertal heifers oocytes, the increased level of ATP12A together with the lower number of apoptotic cells in these blastocysts might support their survival after transfer.

Single-Cell Transcriptional Profile Construction of Rat Pituitary Glands before and after Sexual Maturation and Identification of Novel Marker Spp1 in Gonadotropes.

The mammalian pituitary gland drives highly conserved physiological processes such as somatic cell growth, pubertal transformation, fertility, and metabolism by secreting a variety of hormones. Recently, single-cell transcriptomics techniques have been used in pituitary gland research. However, more studies have focused on adult pituitary gland tissues from different species or different sexes, and no research has yet resolved cellular differences in pituitary gland tissue before and after sexual maturation. Here, we identified a total of 15 cell clusters and constructed single-cell transcriptional profiles of rats before and after sexual maturation. Furthermore, focusing on the gonadotrope cluster, 106 genes were found to be differentially expressed before and after sexual maturation. It was verified that Spp1, which is specifically expressed in gonadotrope cells, could serve as a novel marker for this cell cluster and has a promotional effect on the synthesis and secretion of follicle-stimulating hormone. The results provide a new resource for further resolving the regulatory mechanism of pituitary gland development and pituitary hormone synthesis and secretion.

Characterization of sexual maturity-associated N6-methyladenosine in boar testes.

BACKGROUND: The health and size of the testes are crucial for boar fertility. Testicular development is tightly regulated by epigenetics. N6-methyladenosine (m6A) modification is a prevalent internal modification on mRNA and plays an important role in development. The mRNA m6A methylation in boar testicular development still needs to be investigated. RESULTS: Using the MeRIP-seq technique, we identify and profile m6A modification in boar testes between piglets and adults. The results showed 7783 distinct m6A peaks in piglets and 6590 distinct m6A peaks in adults, with 2,471 peaks shared between the two groups. Enrichment of GO and KEGG analysis reveal dynamic m6A methylation in various biological processes and signalling pathways. Meanwhile, we conjointly analyzed differentially methylated and expressed genes in boar testes before and after sexual maturity, and reproductive related genes (TLE4, TSSK3, TSSK6, C11ORF94, PATZ1, PHLPP1 and PAQR7) were identified. Functional enrichment analysis showed that differential genes are associated with important biological functions, including regulation of growth and development, regulation of metabolic processes and protein catabolic processes. CONCLUSION: The results demonstrate that m6A methylation, differential expression and the related signalling pathways are crucial for boar testicular development. These results suggest a role for m6A modification in boar testicular development and provided a resource for future studies on m6A function in boar testicular development.

Genetic parameters for cloacal gland, sexual maturity, reproductive organs weight, and body weight in meat-type quail.

Selection to increase body weight in poultry can hamper reproduction traits and compromise production efficiency. Thus, attention to reproduction traits is essential to improving the sustainability of breeding programs. Data from a domestic quail breeding program for meat production were used to estimate genetic parameters. We analyzed five traits: 4-week body weight, age at sexual maturity for males and females, cloacal gland area, female, and male reproductive organs weights. A multi-trait mixed model analysis with fixed effects of generation/hatch was performed, assuming environmental covariance equals zero for sex-limited traits. Heritability estimates range from low to moderate for male sexual maturity and cloacal gland area, and high for other traits. Intersexual genetic correlation for age at sexual maturity is positive, which can lead to correlated responses in the other sex. Reproductive organs weights are genetically correlated with body weight, but not significantly between sexes and nor with sexual maturity. Genetic correlations for the cloacal gland area were positive with body weight and negative with age at sexual maturity of males and females, demonstrating a potential use of this trait for selection with favorable outcomes in reproduction. The use of the cloacal gland area can be used in the same way as the scrotal circumference in mammals, improving female reproduction traits by selecting a trait recorded in males.

Loss of Fshr Prevents Testicular Maturation in Atlantic Salmon (Salmo salar L.).

Early puberty poses a significant challenge for male Atlantic salmon in aquaculture due to its negative impact on growth and welfare. The regulation of puberty in vertebrates involves 2 key reproductive hormones: follicle-stimulating hormone (FSH) and luteinizing hormone (LH) and their gonadal receptors. In male mice lacking FSH receptor, testes size is reduced, but fertility is maintained, while medaka and zebrafish with a disrupted fshr gene exhibit near normal testis size and fertility. In these fishes both Fsh and Lh are present during puberty and Lh may rescue fertility, while in salmonid fish only Fsh is present in the circulation during puberty. Using CRISPR-Cas9, we produced crispants with a high prevalence of fshr mutations at the target site, which remained fertile, although more than half showed a testis development deviating from wild-type (wt) males. Crossing out these F0 crispants to each other produced a viable F1 generation showing frameshift (fshr-/-) or in-frame mutations (fshrif/if). Nearly all wt males matured while all fshr-/- males remained immature with small testes containing A spermatogonia as the furthest developed germ cell type and prepubertal plasma androgen levels. Also, the pituitary transcript levels of gnrhr2bba and lhb, but not for fshb, were reduced in the fshr-/- males compared with maturing males. More than half of the fshrif/if mutant males showed no or a delayed maturation. In conclusion, Atlantic salmon show the unique characteristic that loss of Fshr function alone results in male infertility, offering new opportunities to control precocious puberty or fertility in salmon.

Clustering of multi-tissue transcriptomes in gilts with normal cyclicity or delayed puberty reveals genes related to pubertal development†.

In gilts, puberty is marked by standing estrus in the presence of a boar. Delayed puberty (DP; failure to display pubertal estrus) is a major reason for gilt removal. To investigate the physiological determinants underlying DP in gilts, transcriptomic data from tissues relevant to estrus and puberty, such as mediobasal hypothalamus, anterior pituitary gland, ovarian cortex, olfactory bulb, amygdala, and hippocampus, were obtained from age-matched DP (n = 8) and cyclic control gilts at follicular phase (n = 8) and luteal phase (n = 8) of the estrous cycle. A gene expression module analysis via three-way gene × individual × tissue clustering using tensor decomposition identified pituitary and ovary gene modules contributing to regulation of pubertal development. Analysis of gene expression in the hypothalamic-pituitary-ovary axis identified reduced expression of hypothalamic genes critical for stimulating gonadotropin secretion (KISS1 and TAC3) and reduced expression of LHB in the anterior pituitary of DP gilts compared with their cyclic counterparts. Consequently, luteinizing hormone-induced genes in the ovary important for folliculogenesis (OXTR, RUNX2, and PTX3) were less expressed in DP gilts. Other intrafollicular genes (AHR, PTGS2, PTGFR, and IGFBP7) and genes in the steroidogenesis pathways (STAR and CYP11A1) necessary to complete the ovulatory cascade were also less expressed in DP gilts. This is the first clustering of multi-tissue expression data from DP and cyclic gilts to identify genes differentially expressed in gilts of similar ages but at different levels of sexual development. A critical lack of gonadotropin support and reduced ovarian responsiveness underlie DP in gilts.

Probabilistic reaction norm reveals family-related variation in the association between size, condition, and sexual maturation onset in Atlantic salmon (Salmo salar).

This study investigated the relationship between the size, condition, year class, family, and sexual maturity of Atlantic salmon (Salmo salar) using data collected in an aquaculture selective breeding programme. Males that were sexually mature at 2 years of age (maiden spawn) have, on average, greater fork length and condition factor (K) at 1 year of age than their immature counterparts. For every 10-mm increase in fork length or 0.1 increase in K at 1 year of age, the odds of sexual maturity at 2 years of age increased by 1.48 or 1.22 times, respectively. Females that were sexually mature at 3 years of age (maiden spawn) have, on average, greater fork length and K at 2 years of age than their immature counterparts. For every 10-mm increase in fork length or 0.1 increase in K at 2 years of age, the odds of sexual maturity at 3 years of age increased by 1.06 or 1.44 times, respectively. The family explained 34.93% of the variation in sexual maturity among 2-year-old males that was not attributable to the average effects of fork length and K at 1 year of age and year class. The proportion of variation in sexual maturity among 3-year-old females explained by the family could not be investigated. These findings suggest that the onset of sexual maturation in Atlantic salmon is conditional on performance (with respect to energy availability) surpassing a threshold, the magnitude of which can vary between families and is determined by a genetic component. This could support the application of genetic selection to promote or inhibit the onset of sexual maturation in farmed stocks.

Genome-wide association study of age at puberty and its (co)variances with fertility and stature in growing and lactating Holstein-Friesian dairy cattle.

Reproductive performance is a key determinant of cow longevity in a pasture-based, seasonal dairy system. Unfortunately, direct fertility phenotypes such as intercalving interval or pregnancy rate tend to have low heritabilities and occur relatively late in an animal's life. In contrast, age at puberty (AGEP) is a moderately heritable, early-in-life trait that may be estimated using an animal's age at first measured elevation in blood plasma progesterone (AGEP4) concentrations. Understanding the genetic architecture of AGEP4 in addition to genetic relationships between AGEP4 and fertility traits in lactating cows is important, as is its relationship with body size in the growing animal. Thus, the objectives of this research were 3-fold. First, to estimate the genetic and phenotypic (co)variances between AGEP4 and subsequent fertility during first and second lactations. Second, to quantify the associations between AGEP4 and height, length, and BW measured when animals were approximately 11 mo old (standard deviation = 0.5). Third, to identify genomic regions that are likely to be associated with variation in AGEP4. We measured AGEP4, height, length, and BW in approximately 5,000 Holstein-Friesian or Holstein-Friesian × Jersey crossbred yearling heifers across 54 pasture-based herds managed in seasonal calving farm systems. We also obtained calving rate (CR42, success or failure to calve within the first 42 d of the seasonal calving period), breeding rate (PB21, success or failure to be presented for breeding within the first 21 d of the seasonal breeding period) and pregnancy rate (PR42, success or failure to become pregnant within the first 42 d of the seasonal breeding period) phenotypes from their first and second lactations. The animals were genotyped using the Weatherby's Versa 50K SNP array (Illumina, San Diego, CA). The estimated heritabilities of AGEP4, height, length, and BW were 0.34 (90% credibility interval [CRI]: 0.30, 0.37), 0.28 (90% CRI: 0.25, 0.31), 0.21 (90% CRI: 0.18, 0.23), and 0.33 (90% CRI: 0.30, 0.36), respectively. In contrast, the heritabilities of CR42, PB21 and PR42 were all <0.05 in both first and second lactations. The genetic correlations between AGEP4 and these fertility traits were generally moderate, ranging from 0.11 to 0.60, whereas genetic correlations between AGEP4 and yearling body-conformation traits ranged from 0.02 to 0.28. Our GWAS highlighted a genomic window on chromosome 5 that was strongly associated with variation in AGEP4. We also identified 4 regions, located on chromosomes 14, 6, 1, and 11 (in order of decreasing importance), that exhibited suggestive associations with AGEP4. Our results show that AGEP4 is a reasonable predictor of estimated breeding values for fertility traits in lactating cows. Although the GWAS provided insights into genetic mechanisms underpinning AGEP4, further work is required to test genomic predictions of fertility that use this information.

Regulatory Effects of the Kiss1 Gene in the Testis on Puberty and Reproduction in Hezuo and Landrance Boars.

Kisspeptin, a neuropeptide encoded by the Kiss1 gene, combines with its receptor Kiss1R to regulate the onset of puberty and male fertility by the hypothalamic-pituitary-gonadal axis. However, little is known regarding the expression signatures and molecular functions of Kiss1 in the testis. H&E staining revealed that well-arranged spermatogonia, spermatocytes, round and elongated spermatids, and spermatozoa, were observed in 4-, 6-, and 8-month-old testes compared to 1- and 3-month-old testes of Hezuo pigs; however, these were not observed in Landrance until 6 months. The diameter, perimeter, and cross-sectional area of seminiferous tubules and the perimeter and area of the tubular lumen increased gradually with age in both pigs. Still, Hezuo pigs grew faster than Landrance. The cloning results suggested that the Hezuo pigs' Kiss1 CDS region is 417 bp in length, encodes 138 amino acids, and is highly conserved in the kisspeptin-10 region. qRT-PCR and Western blot indicated that the expression trends of Kiss1 mRNA and protein were essentially identical, with higher expression levels at post-pubertal stages. Immunohistochemistry demonstrated that the Kiss1 protein was mainly located in Leydig cells and post-pubertal spermatogenic cells, ranging from round spermatids to spermatozoa. These studies suggest that Kiss1 is an essential regulator in the onset of puberty and spermatogenesis of boars.

Molecular characterization of Smtdc-1 and Smddc-1 discloses roles as male-competence factors for the sexual maturation of Schistosoma mansoni females.

Introduction: Schistosomes are the only mammalian flatworms that have evolved separate sexes. A key question of schistosome research is the male-dependent sexual maturation of the female since a constant pairing contact with a male is required for the onset of gonad development in the female. Although this phenomenon is long known, only recently a first peptide-based pheromone of males was identified that contributes to the control of female sexual development. Beyond this, our understanding of the molecular principles inducing the substantial developmental changes in a paired female is still rudimentary. Objectives: Previous transcriptomic studies have consistently pointed to neuronal genes being differentially expressed and upregulated in paired males. These genes included Smp_135230 and Smp_171580, both annotated as aromatic-L-amino-acid decarboxylases (DOPA decarboxylases). Here, we characterized both genes and investigated their roles in male-female interaction of S. mansoni. Methodologies/findings: Sequence analyses indicated that Smp_135230 represents an L-tyrosine decarboxylase (Smtdc-1), whereas Smp_171580 represents a DOPA decarboxylase (Smddc-1). By qRT-PCR, we confirmed the male-specific and pairing-dependent expression of both genes with a significant bias toward paired males. RNA-interference experiments showed a strong influence of each gene on gonad differentiation in paired females, which was enhanced by double knockdown. Accordingly, egg production was significantly reduced. By confocal laser scanning microscopy, a failure of oocyte maturation was found in paired knockdown females. Whole-mount in situ hybridization patterns exhibited the tissue-specific occurrence of both genes in particular cells at the ventral surface of the male, the gynecophoral canal, which represents the physical interface of both genders. These cells probably belong to the predicted neuronal cluster 2 of S. mansoni. Conclusion: Our results suggest that Smtdc-1 and Smddc-2 are male-competence factors that are expressed in neuronal cells at the contact zone between the genders as a response of pairing to subsequently control processes of female sexual maturation.

Genotype-specific variation in seasonal body condition at a large-effect maturation locus.

Organisms use resource allocation strategies to survive seasonal environmental changes and life-history stage transitions. Earlier studies found a transcription cofactor, vgll3, associating with maturation timing that inhibits adipogenesis in mice and affects body condition in juvenile salmon. Owing to a lack of temporal studies examining seasonality effects on phenotypes such as vgll3 genotype, body condition, maturation and different life stages, we investigated the influence of different larval and juvenile temperatures, vgll3 genotype and interactions with body condition and maturation rate. We reared Atlantic salmon for 2 years in four larval-juvenile phase temperature groups until the occurrence of mature males. We found no effect of larval temperature on the measured phenotypes or maturation rate. However, we observed an increased maturation rate in individuals of the warm juvenile temperature treatment and differences in body condition associated with vgll3 genotype. Early maturation genotype individuals had a less variable body condition across seasons compared with late maturation genotype individuals. This result suggests a vgll3 influence on resource allocation strategies; possibly linked with the early maturation process, with early maturation genotype individuals having a higher maturation rate and a higher body condition in the spring.

Enduring sex-dependent effects of lipopolysaccharide treatment on the hypothalamic-pituitary-gonadal axis in mice.

Pubertal stress causes enduring sexual behavior dysfunction in males and females, but the underlying mechanism remains unknown. These changes may arise from pubertal programming of the hypothalamic-pituitary-gonadal axis. Previous findings show that stress exposure downregulates the hypothalamic-pituitary-gonadal axis, particularly through the reduction of the neuropeptide kisspeptin (Kiss1) and its receptor (Kiss1R). Although acute changes in kiss1 and Kiss1r genes have been observed following pubertal immune stress, it is unclear whether immune stress-induced downregulation of kiss1 and kiss1r persists beyond puberty. The current study investigated the enduring sex-specific consequences of lipopolysaccharide on the expression of Kiss1 and Kiss1r in 160 pubertal or adult mice at multiple time points. Six-week and 10-week-old male and female mice were treated with either saline or with lipopolysaccharide. Mice were euthanized either 8 h or 4 weeks following treatment. Although we did not identify any sex differences, our results revealed that lipopolysaccharide treatment decreases hypothalamic Kiss1 and Kiss1r in both pubertal and adult mice within 8 h of treatment. The decreased hypothalamic Kiss1 expression persists 4 weeks later only in mice treated with lipopolysaccharide during puberty. Our findings highlight the age-dependent vulnerability of the hypothalamic-pituitary-gonadal axis to immune stress, providing a better understanding of the mechanisms implicated in allostatic shift during immune stress. Finally, our findings also show the effects of immune stress on various components of the hypothalamic-pituitary-gonadal axis, which could have implications for sexual and fertility-related dysfunctions.

gBLUP-GWAS identifies candidate genes, signaling pathways, and putative functional polymorphisms for age at puberty in gilts.

Successful development of replacement gilts determines their reproductive longevity and lifetime productivity. Selection for reproductive longevity is challenging due to low heritability and expression late in life. In pigs, age at puberty is the earliest known indicator for reproductive longevity and gilts that reach puberty earlier have a greater probability of producing more lifetime litters. Failure of gilts to reach puberty and display a pubertal estrus is a major reason for early removal of replacement gilts. To identify genomic sources of variation in age at puberty for improving genetic selection for early age at puberty and related traits, gilts (n = 4,986) from a multigeneration population representing commercially available maternal genetic lines were used for a genomic best linear unbiased prediction-based genome-wide association. Twenty-one genome-wide significant single nucleotide polymorphisms (SNP) located on Sus scrofa chromosomes (SSC) 1, 2, 9, and 14 were identified with additive effects ranging from -1.61 to 1.92 d (P < 0.0001 to 0.0671). Novel candidate genes and signaling pathways were identified for age at puberty. The locus on SSC9 (83.7 to 86.7 Mb) was characterized by long range linkage disequilibrium and harbors the AHR transcription factor gene. A second candidate gene on SSC2 (82.7 Mb), ANKRA2, is a corepressor for AHR, suggesting a possible involvement of AHR signaling in regulating pubertal onset in pigs. Putative functional SNP associated with age at puberty in the AHR and ANKRA2 genes were identified. Combined analysis of these SNP showed that an increase in the number of favorable alleles reduced pubertal age by 5.84 ± 1.65 d (P < 0.001). Candidate genes for age at puberty showed pleiotropic effects with other fertility functions such as gonadotropin secretion (FOXD1), follicular development (BMP4), pregnancy (LIF), and litter size (MEF2C). Several candidate genes and signaling pathways identified in this study play a physiological role in the hypothalamic-pituitary-gonadal axis and mechanisms permitting puberty onset. Variants located in or near these genes require further characterization to identify their impact on pubertal onset in gilts. Because age at puberty is an indicator of future reproductive success, these SNP are expected to improve genomic predictions for component traits of sow fertility and lifetime productivity expressed later in life.

Selecting for replacement gilts is challenging because sow reproductive traits are lowly heritable and expressed late in life. Age at puberty is the earliest indicator of future reproductive success of gilts. Genetic selection for early onset of puberty could be feasible with the availability of molecular genetic predictors for age at puberty. To identify genomic sources associated with variation in age at puberty in gilts, a large-scale genome-wide association study was conducted at the U.S Meat Animal Research Center, Clay Center, Nebraska. Novel genomic associations for age at puberty were identified. Several candidate genes identified for age at puberty are involved in signaling pathways that regulate ovarian functions and pubertal onset. Potential causative genetic variants for age at puberty were identified within the candidate genes. These novel SNP are important new markers for use in genomic selection of replacement gilts with early puberty and provide critical new insight into biological mechanisms important for pubertal development in gilts.

A pituitary gene network linking vgll3 to regulators of sexual maturation in male Atlantic salmon.

Age at maturity is a key life history trait and a significant contributor to life history strategy variation. The maturation process is influenced by genetic and environmental factors, but specific causes of variation in maturation timing remain elusive. In many species, the increase in the regulatory gonadotropin-releasing hormone 1 (GnRH1) marks the onset of puberty. Atlantic salmon, however, lacks the gnrh1 gene, suggesting gnrh3 and/or other regulatory factors are involved in the maturation process. Earlier research in Atlantic salmon has found a strong association between alternative alleles of vgll3 and maturation timing. Recently we reported strong induction of gonadotropin genes (fshb and lhb) in the pituitary of Atlantic salmon homozygous for the early maturation allele (E) of vgll3. The induction of gonadotropins was accompanied by increased expression of their direct upstream regulators, c-jun and sf1 (nr5a1b) but the regulatory connection between vgll3 and these regulators has never been investigated in any organism. In this study, we investigated the potential regulatory connection between vgll3 genotypes and these regulators through a stepwise approach of identifying a gene regulatory network (GRN) containing c-jun and sf1, and transcription factor motif enrichment analysis. We found a GRN containing c-jun with predicted upstream regulators, e2f1, egr1, foxj1 and klf4, to be differentially expressed in the pituitary. Finally, we suggest a vgll3 and Hippo pathway -dependent model for transcriptional regulation of c-jun and sf1 in the pituitary, which may have broader implications across vertebrates.

Correlation scan: identifying genomic regions that affect genetic correlations applied to fertility traits.

Although the genetic correlations between complex traits have been estimated for more than a century, only recently we have started to map and understand the precise localization of the genomic region(s) that underpin these correlations. Reproductive traits are often genetically correlated. Yet, we don't fully understand the complexities, synergism, or trade-offs between male and female fertility. In this study, we used reproductive traits in two cattle populations (Brahman; BB, Tropical Composite; TC) to develop a novel framework termed correlation scan (CS). This framework was used to identify local regions associated with the genetic correlations between male and female fertility traits. Animals were genotyped with bovine high-density single nucleotide polymorphisms (SNPs) chip assay. The data used consisted of ~1000 individual records measured through frequent ovarian scanning for age at first corpus luteum (AGECL) and a laboratory assay for serum levels of insulin growth hormone (IGF1 measured in bulls, IGF1b, or cows, IGF1c). The methodology developed herein used correlations of 500-SNP effects in a 100-SNPs sliding window in each chromosome to identify local genomic regions that either drive or antagonize the genetic correlations between traits. We used Fisher's Z-statistics through a permutation method to confirm which regions of the genome harboured significant correlations. About 30% of the total genomic regions were identified as driving and antagonizing genetic correlations between male and female fertility traits in the two populations. These regions confirmed the polygenic nature of the traits being studied and pointed to genes of interest. For BB, the most important chromosome in terms of local regions is often located on bovine chromosome (BTA) 14. However, the important regions are spread across few different BTA's in TC. Quantitative trait loci (QTLs) and functional enrichment analysis revealed many significant windows co-localized with known QTLs related to milk production and fertility traits, especially puberty. In general, the enriched reproductive QTLs driving the genetic correlations between male and female fertility are the same for both cattle populations, while the antagonizing regions were population specific. Moreover, most of the antagonizing regions were mapped to chromosome X. These results suggest regions of chromosome X for further investigation into the trade-offs between male and female fertility. We compared the CS with two other recently proposed methods that map local genomic correlations. Some genomic regions were significant across methods. Yet, many significant regions identified with the CS were overlooked by other methods.

Genetic Basis of Sexual Maturation Heterosis: Insights From Ovary lncRNA and mRNA Repertoire in Chicken.

Sexual maturation is fundamental to the reproduction and production performance, heterosis of which has been widely used in animal crossbreeding. However, the underlying mechanism have long remained elusive, despite its profound biological and agricultural significance. In the current study, the reciprocal crossing between White Leghorns and Beijing You chickens were performed to measure the sexual maturation heterosis, and the ovary lncRNAs and mRNAs of purebreds and crossbreeds were profiled to illustrate molecular mechanism of heterosis. Heterosis larger than 20% was found for pubic space and oviduct length, whereas age at first egg showed negative heterosis in both crossbreeds. We identified 1170 known lncRNAs and 1994 putative lncRNAs in chicken ovary using a stringent pipeline. Gene expression pattern showed that nonadditivity was predominant, and the proportion of nonadditive lncRNAs and genes was similar between two crossbreeds, ranging from 44.24% to 49.15%. A total of 200 lncRNAs and 682 genes were shared by two crossbreeds, respectively. GO and KEGG analysis showed that the common genes were significantly enriched in the cell cycle, animal organ development, gonad development, ECM-receptor interaction, calcium signaling pathway and GnRH signaling pathway. Weighted gene co-expression network analysis (WGCNA) identified that 7 out of 20 co-expressed lncRNA-mRNA modules significantly correlated with oviduct length and pubic space. Interestingly, genes harbored in seven modules were also enriched in the similar biological process and pathways, in which nonadditive lncRNAs, such as MSTRG.17017.1 and MSTRG.6475.20, were strongly associated with nonadditive genes, such as CACNA1C and TGFB1 to affect gonad development and GnRH signaling pathway, respectively. Moreover, the results of real-time quantitative PCR (RT-qPCR) correlated well with the transcriptome data. Integrated with positive heterosis of serum GnRH and melatonin content detected in crossbreeds, we speculated that nonadditive genes involved in the GnRH signaling pathway elevated the gonad development, leading to the sexual maturation heterosis. We characterized a systematic landscape of ovary lncRNAs and mRNAs related to sexual maturation heterosis in chicken. The quantitative exploration of hybrid transcriptome changes lays foundation for genetic improvement of sexual maturation traits and provides insights into endocrine control of sexual maturation.

Hypothalamic Overexpression of Makorin Ring Finger Protein 3 Results in Delayed Puberty in Female Mice.

Makorin ring finger protein 3 (MKRN3) is an important neuroendocrine player in the control of pubertal timing and upstream inhibitor of gonadotropin-releasing hormone secretion. In mice, expression of Mkrn3 in the hypothalamic arcuate and anteroventral periventricular nucleus is high early in life and declines before the onset of puberty. Therefore, we aimed to explore if the persistence of hypothalamic Mkrn3 expression peripubertally would result in delayed puberty. Female mice that received neonatal bilateral intracerebroventricular injections of a recombinant adeno-associated virus expressing Mkrn3 had delayed vaginal opening and first estrus compared with animals injected with control virus. Subsequent estrous cycles and fertility were normal. Interestingly, male mice treated similarly did not exhibit delayed puberty onset. Kiss1, Tac2, and Pdyn mRNA levels were increased in the mediobasal hypothalamus in females at postnatal day 28, whereas kisspeptin and neurokinin B protein levels in the arcuate nucleus were decreased, following Mkrn3 overexpression, compared to controls. Cumulatively, these data suggest that Mkrn3 may directly or indirectly target neuropeptides of Kiss1 neurons to degradation pathways. This mouse model suggests that MKRN3 may be a potential contributor to delayed onset of puberty, in addition to its well-established roles in central precocious puberty and the timing of menarche.

Dicer ablation in Kiss1 neurons impairs puberty and fertility preferentially in female mice.

Kiss1 neurons, producing kisspeptins, are essential for puberty and fertility, but their molecular regulatory mechanisms remain unfolded. Here, we report that congenital ablation of the microRNA-synthesizing enzyme, Dicer, in Kiss1 cells, causes late-onset hypogonadotropic hypogonadism in both sexes, but is compatible with pubertal initiation and preserved Kiss1 neuronal populations at the infantile/juvenile period. Yet, failure to complete puberty and attain fertility is observed only in females. Kiss1-specific ablation of Dicer evokes disparate changes of Kiss1-cell numbers and Kiss1/kisspeptin expression between hypothalamic subpopulations during the pubertal-transition, with a predominant decline in arcuate-nucleus Kiss1 levels, linked to enhanced expression of its repressors, Mkrn3, Cbx7 and Eap1. Our data unveil that miRNA-biosynthesis in Kiss1 neurons is essential for pubertal completion and fertility, especially in females, but dispensable for initial reproductive maturation and neuronal survival in both sexes. Our results disclose a predominant miRNA-mediated inhibitory program of repressive signals that is key for precise regulation of Kiss1 expression and, thereby, reproductive function.