Porcine SOX9 Gene Expression Is Influenced by an 18 bp Indel in the 5'-Untranslated Region

Whole-genome de novo sequencing reveals genomic variants associated with differences of sex development in SRY negative pigs

BACKGROUND

Differences of sex development (DSD) are congenital conditions in which chromosomal, gonadal, or phenotypic sex is atypical. In more than 50% of human DSD cases, a molecular diagnosis is not available. In intensively farmed pig populations, the incidence of XX DSD pigs is relatively high, leading to economic losses for pig breeders. Interestingly, in the majority of 38, XX DSD pigs, gonads still develop into testis-like structures or ovotestes despite the absence of the testis-determining gene (SRY). However, the current understanding of the molecular background of XX DSD pigs remains limited.

METHODS

Anatomical and histological characteristics of XX DSD pigs were analysed using necropsy and HE staining. We employed whole-genome sequencing (WGS) with 10× Genomics technology and used de novo assembly methodology to study normal female and XX DSD pigs. Finally, the identified variants were validated in 32 XX DSD pigs, and the expression levels of the candidate variants in the gonads of XX DSD pigs were further examined.

RESULTS

XX DSD pigs are characterised by the intersex reproductive organs and the absence of germ cells in the seminiferous tubules of the gonads. We identified 4,950 single-nucleotide polymorphisms (SNPs) from non-synonymous mutations in XX DSD pigs. Cohort validation results highlighted two specific SNPs, "c.218T > C" in the "Interferon-induced transmembrane protein 1 gene (IFITM1)" and "c.1043C > G" in the "Newborn ovary homeobox gene (NOBOX)", which were found exclusively in XX DSD pigs. Moreover, we verified 14 candidate structural variants (SVs) from 1,474 SVs, identifying a 70 bp deletion fragment in intron 5 of the WW domain-containing oxidoreductase gene (WWOX) in 62.5% of XX DSD pigs. The expression levels of these three candidate genes in the gonads of XX DSD pigs were significantly different from those of normal female pigs.

CONCLUSION

The nucleotide changes of IFITM1 (c.218T > C), NOBOX (c.1043 C > G), and a 70 bp deletion fragment of the WWOX were the most dominant variants among XX DSD pigs. This study provides a theoretical basis for better understanding the molecular background of XX DSD pigs. DSD are conditions affecting development of the gonads or genitalia. These disorders can happen in many different types of animals, including pigs, goats, dogs, and people. In people, DSD happens in about 0.02-0.13% of births, and in pigs, the rate is between 0.08% and 0.75%. Pigs have a common type of DSD where the animal has female chromosomes (38, XX) but no SRY gene, which is usually found on the Y chromosome in males. XX DSD pigs may look like both males and females on the outside and have testis-like or ovotestis (a mix of ovary and testis) gonads inside. XX DSD pigs often lead to not being able to have piglets, slower growth, lower chance of survival, and poorer meat quality. Here, we used a method called whole-genome de novo sequencing to look for variants in the DNA of XX DSD pigs. We then checked these differences in a larger group of pigs. Our results reveal the nucleotide changes in IFITM1 (c.218T > C), NOBOX (c.1043 C > G), and a 70 bp deletion fragment in intron 5 of the WWOX, all linked to XX DSD pigs. The expression levels of these three genes were also different in the gonads of XX DSD pigs compared to normal female pigs. These variants are expected to serve as valuable molecular markers for XX DSD pigs. Because pigs are a lot like humans in their genes, physiology, and body structure, this research could help us learn more about what causes DSD in people.

Transcriptome Analyses Show Changes in Gene Expression Triggered by a 31-bp InDel within OsSUT3 5'UTR in Rice Panicle

Pollen development and its fertility are obligatory conditions for the reproductive success of flowing plants. Sucrose transporter 3 (OsSUT3) is known to be preferentially expressed and may play critical role in developing pollen. A 31-bp InDel was identified as a unique variation and was shown to be responsible for the expression of downstream gene in our previous study. In this study, to analyze the changes of gene expression triggered by 31-bp InDel during pollen development, two vectors (p385-In/Del::OsSUT3-GUS) were constructed and then stably introduced into rice. Histochemical and quantitative real-time PCR (qRT-PCR) analysis of transgenic plants showed that 31-bp deletion drastically reduced the expressions of downstream genes, including both OsSUT3 and GUS in rice panicle at booting stage, especially that of OsSUT3. The transcriptome profile of two types of panicles at booting stage revealed a total of 1028 differentially expressed genes (DEGs) between 31-bp In and 31-bp Del transgenic plants. Further analyses showed that 397 of these genes were significantly enriched for the 'metabolic process' and 'binding'. Among them, nineteen genes had a strong relationship with starch and sucrose metabolism and were identified as candidate genes potentially associated with the starch accumulation in rice pollen, which that was also verified via qRT-PCR. In summary, 31-bp InDel plays a crucial role not only in the regulation of downstream genes but in the expression of sucrose-starch metabolizing genes in multiple biological pathways, and provides a different regulation mechanism for sucrose metabolism in pollen.

Screening Indels from the whole genome to identify the candidates and their association with economic traits in several goat breeds

In the present study, the re-sequencing of our previous whole-genome sequencing (WGS) for selected individuals of Dazu-black goat (DBG) and Inner-Mongolia Cashmere goat (IMCG) breeds was used to detect and compare the differentiation in Indels depending on the reference genome of goat. Then, three selected candidate Indels rs668795676, rs657996810, and rs669452874 of the three genes SUFU, SYCP2L and GLIPR1L1, respectively, have been chosen, based on the results of prior GWAS across the genome, and examined for their association with body weight and dimensions (body height, body length, heart girth, chest width, cannon circumference, and chest depth) by kompetitive allele specific PCR assay for 342 goats from the three studied goat breeds (DBG, n = 203; ♂99, ♀104), IMCG (n = 65; 15♂, 50♀), and Hechuan white goat (HWG, n = 74; 34♂, 40♀) breeds. The analysis of 192.747 Gb WGS revealed an average 334,151 Indels in the whole genome of DBG and IMCG breeds. Chromosome 1 had a maximum number of mutations (Indels) of 58,497 and 55,527 for IMCG and DBG, respectively, while chromosome 25 had the least number of mutations of 15,680 and 16,103 for IMCG and DBG, respectively. The majority of Indels were either Ins or Del of short fragments of 1-5 bp, which covered 79.06 and 71.78% of the total number of Indels mutations in IMCG and DBG, respectively. Comparing the differences of Indels between the studied goat breeds revealed 100 and 110 unique Indels for IMCG and DBG, respectively. The Indels loci in the intron region were unique for both studied goat breeds which were related to 30 and 38 candidate genes in IMCG and DBG, respectively, including SUFU, SYCP2L, and GLIPR1L1 genes. Concerning rs669452874 locus, body height and body length of Del/Del genotype in DBG were significantly higher (P < 0.05) than that of Ins/Del genotype, while body height and body length of Del/Del genotype in IMCG were significantly higher (P < 0.01) than those of Ins/Ins genotype, whereas body height and body length and heart girth of Del/Del genotype in HWG were significantly higher (P < 0.01) than those of the Ins/Del and Ins/Ins genotypes. Thus, Del/Del genotype of rs669452874 locus can be used as a candidate molecular marker related to the body dimensions in the studied goat breeds.

Alternative Splicing Isoforms of Porcine CREB Are Differentially Involved in Transcriptional Transactivation

The cAMP response element-binding protein (CREB), a basic leucine zipper transcription factor, is involved in the activation of numerous genes in a variety of cell types. The CREB gene is rich in alternative splicing (AS) events. However, studies on the AS of CREB genes in pigs are limited, and few reports have compared the roles of isoforms in activating gene expression. Here, five AS transcripts, V1–5, were characterized by RT-PCR and two, V3 and V5, were new identifications. Both V1 and V2 have all the functional domains of the CREB protein, with similar tissue expression profiles and mRNA stability, suggesting that they have similar roles. The transcriptional transactivation activities of four isoforms encoding complete polypeptides were analyzed on the expression of the B-cell CLL/lymphoma 2-like protein 2 and the poly (A)-binding protein, nuclear 1 genes with a dual-luciferase reporter system, and differential activities were observed. Both V1 and V2 have promoting effects, but their roles are gene-specific. V3 has no effect on the promoter of the two genes, while V4 functions as a repressor. The mechanisms underlying the differential roles of V1 and V2 were analyzed with RNA-seq, and the genes specifically regulated by V1 and V2 were identified. These results will contribute to further revealing the role of CREB and to analyzing the significance of AS in genes.

A 104-bp Structural Variation of the ADPRHL1 Gene Is Associated With Growth Traits in Chickens

Analyzing marker-assisted breeding is an important method utilized in modern molecular breeding. Recent studies have determined that a large number of molecular markers appear to explain the impact of "lost heritability" on human height. Therefore, it is necessary to locate molecular marker sites in poultry and investigate the possible molecular mechanisms governing their effects. In this study, we found a 104-bp insertion/deletion polymorphism in the 5'UTR of the ADPRHL1 gene through resequencing. In cross-designed F₂ resource groups, the indel was significantly associated with weight at 0, 2, 4, 6, and 10 weeks and a number of other traits [carcass weight (CW), semi-evisceration weight (SEW), evisceration weight (EW), claw weight (CLW), wings weight (DWW), gizzard weight (GW), pancreas weight (PW), chest muscle weight (CMW), leg weight (LW), leg muscle weight (LMW), shedding Weight (SW), liver rate (LR), and leg muscle rate (LMR)] (P < 0.05). In brief, the insertion-insertion (II) genotype was significantly associated with the greatest growth traits and meat quality traits, whereas the values associated with the insertion-deletion (ID) genotype were the lowest in the F₂ reciprocal cross chickens. The mutation sites were genotyped in 4,526 individuals from 12 different chicken breeds and cross-designed F₂ resource groups. The II genotype is the most important genotype in commercial broilers, and the I allele frequency observed in these breeds is relatively high. Deletion mutations tend to be fixed in commercial broilers. However, there is still considerable great potential for breeding in dual-purpose chickens and commercial laying hens. A luciferase reporter assay showed that the II genotype of the ADPRHL1 gene possessed 2.49-fold higher promoter activity than the DD genotype (P < 0.05). We hypothesized that this indel might affect the transcriptional activity of ADPRHL1, thereby affecting the growth traits of chickens. These findings may help to elucidate the function of the ADPRHL1 gene and facilitate enhanced reproduction in the chicken industry.

Comprehensive transcriptome analysis of hypothalamus reveals genes associated with disorders of sex development in pigs

XX sex reversal, also called XX disorders of sex development (XX-DSD), is a condition affecting the development of the gonads or genitalia, and is relatively common in pigs. However, its genetic etiology and transcriptional regulation mechanism in the hypothalamic-pituitary-gonadal axis (HPGA) remain mostly unknown. XX-DSD (SRY-negative) pigs and normal sows were selected by external genitalia observation. The hypothalamus, which is the integrated center of the HPGA was sampled for whole-transcriptome RNA-seq. The role of DEmiRNA was validated by its overexpression and knockdown in vitro. A total of 1,258 lncRNAs, 1,086 mRNAs, and 61 microRNAs differentially expressed in XX-DSD pigs compared with normal female pigs. Genes in the hormone biosynthesis and secretion pathway significantly up-regulated, and the up-regulation of GNRH1, KISS1 and AVP may associate with the abnormal secretion of GnRH. We also predicted the lncRNA-miRNA-mRNA co-expression triplets and constructed three competing endogenous RNA (ceRNA) potentially associated with XX-DSD. Functional enrichment studies suggested that TCONS_00340886, TCONS_00000204 and miR-181a related to GnRH secretion. Further, miR-181a inhibitor up-regulated GNRH1, PAK6, and CAMK4 in the GT1-7 cells. Conversely, transfection of miR-181a mimics obtained the opposite trends. The expression levels of FSHR, LHR, ESR1 and ESR2 were significantly higher in XX-DSD gondas than those in normal sows. Taken together, we proposed that the balance of endocrine had broken in XX-DSD pigs. The current study is the first to examine the transcriptomic profile in the hypothalamus of XX-DSD pigs. It provides new insight into coding and non-coding RNAs that may be associated with DSD in pigs.

Goat DNMT3B: An indel mutation detection, association analysis with litter size and mRNA expression in gonads

DNA methyltransferase 3β (DNMT3B) is a gene encoding a de novo methylation enzyme that is required for DNA methylation during mammalian embryo development. Previous genome-wide association analysis suggested DNMT3B is a candidate gene for goat fertility, but there is no study on the effect of DNMT3B on litter size in goat. The aim of this study was to identify possible insertion/deletion (indel) mutations associated with litter size. Seven putative indels were designed to study their association with litter size, but just one 11-bp insertion variant of intron 22 (the last intron) was found in healthy female Shaanbei white cashmere goats (SBWC goats) (n = 1534). Statistical analysis showed that the 11-bp insertion was related to the first-born litter size (P < 0.01) and the goats with the deletion/deletion genotype had a higher average first-born litter size (P < 0.01). In addition, the expression profile of the DNMT3B mRNA in goat was detected, which revealed significant differences in DNMT3B mRNA expression in the gonads. Additionally, the results of western blotting revealed that the ovaries of mothers of multi-lamb (MML) had a higher level of DNMT3B protein than the ovaries of mothers of single-lamb (MSL). Furthermore, the mRNA of DNMT3B was widely expressed in male goats. Differences in mRNA expression levels were observed in the ovaries of MSL and MML. These findings indicated that the 11-bp indel in DNMT3B was significantly associated with first-born litter size, which can be used for marker-assisted selection (MAS) of goats for breeding.

Pig SOX9: Expression profiles of Sertoli cell (SCs) and a functional 18 bp indel affecting testis weight

Sex determining region Y-box 9 (SOX9), an important member of the SRY- type HMGbox (SOX) gene family, plays an important role in the regulation of mammalian reproduction, including sex differentiation during the embryonic development stage and spermatogenesis after birth. To explore the roles of polymorphism and expression of the SOX9 gene in the development of testes, we analyzed the indel of SOX9 in pigs and the corresponding expression level of the SOX9 gene in 7-day and 5-month-old porcine Sertoli cells. Results revealed that the DD haplotype of SOX9 gene as well as the ID genotype were significantly associated with larger testicular weight, while the II haplotype was closely related to the smaller testicular weight. More importantly, the SOX9 gene expression of ID genotyped group was significantly higher than that in II genotyped group. Our results first revealed that the indel polymorphism and expression of SOX9 were significantly associated with pig reproduction traits indicating the critical roles of SOX9 gene in testes development. The study provides a new clue for understanding the regulation of animal reproductive activities.

Two insertion/deletion variants in the promoter region of the QPCTL gene are significantly associated with body weight and carcass traits in chickens

Glutaminyl-peptide cyclotransferase-like (QPCTL) is an isoenzyme of glutaminyl-peptide cyclotransferase (QPCT). QPCTL and QPCT catalyze the formation of N-terminal modified pyroglutamate-fractalkine and the chemokine CCL2. The objective of this study was to investigate the association between insertions/deletions in the chicken QPCTL promoter region with growth traits in chickens. We first detected two insertion/deletion variants of QPCTL via whole-genome resequencing analysis of DNA samples from Xichuan chickens. A total of 1896 individuals from 12 breeds were genotyped for 52- and 224-bp insertions/deletions. We found two novel insertions/deletions in the promoter region of the chicken QPCTL gene and studied their association with chicken body weight and carcass traits. Our findings show that QPCTL can be a molecular marker for chicken genetics and breeding programs.

Identification of a novel 12-bp insertion/deletion (indel) of iPS-related Oct4 gene and its association with reproductive traits in male piglets

As a key factor of cellular reprogramming, Oct4 is one of vital transcription factors for induced pluripotent stem cells (iPSCs). Loss of its function or deletion causes apoptosis in primordial germ cells (PGCs), which affect reproductive traits in mammals. In this study, a novel 12-bp insertion/deletion (indel) polymorphism (NC_010449:g.2759-2760insGGTTTTTGTCTA) within the Oct4 gene was identified in 442 pigs of Large White (LW) and Landrace (LD) breeds, showing three genotypes designated as II, ID, and DD. The frequencies of allele "I" in LW and LD pigs were 0.587 and 0.648, respectively. The male piglets with homozygous II or DD genotypes of Oct4 gene exhibited better reproductive traits than those with heterozygous ID genotype. Moreover, there were two significant associations between this 12-bp indel polymorphism and testis long circumference (TLC) (P=0.005) and testis short girth (TSG) (P=0.003) as well as 15-day testis weight (TW) (P=0.013) in the LW male piglets. These findings suggest that the 12-bp indel polymorphism of the Oct4 gene might be a potential DNA marker for selecting preferred individuals in relation to reproductive traits in pig marker-assisted selection (MAS) breeding, which could contribute to the breeding and genetics in male piglets.

Sex Reversal in Non-Human Placental Mammals

Gonads are very peculiar organs given their bipotential competence. Indeed, early differentiating genital ridges evolve into either of 2 very distinct organs: the testis or the ovary. Accumulating evidence now demonstrates that both genetic pathways must repress the other in order for the organs to differentiate properly, meaning that if this repression is disrupted or attenuated, the other pathway may completely or partially be expressed, leading to disorders of sex development. Among these disorders are the cases of XY male-to-female and XX female-to-male sex reversals as well as true hermaphrodites, in which there is a discrepancy between the chromosomal and gonadal sex. Here, we review known cases of XY and XX sex reversals described in mammals, focusing mostly on domestic animals where sex reversal pathologies occur and on wild species in which deviations from the usual XX/XY system have been documented.