Next Article in Journal
Existing and Emerging Metabolomic Tools for ALS Research
Next Article in Special Issue
Regulation of AMH, AMHR-II, and BMPs (2,6) Genes of Bovine Granulosa Cells Treated with Exogenous FSH and Their Association with Protein Hormones
Previous Article in Journal
Transcriptome Meta-Analysis Deciphers a Dysregulation in Immune Response-Associated Gene Signatures during Sepsis
Previous Article in Special Issue
DDB1 Regulates Sertoli Cell Proliferation and Testis Cord Remodeling by TGFβ Pathway
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 10
Issue 12
10.3390/genes10121010
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Sex Determination Using RNA-Sequencing Analyses in Early Prenatal Pig Development

by
Susana A. Teixeira
1,
Adriana M. G. Ibelli
2,
Maurício E. Cantão
2,
Haniel C. de Oliveira
1,
Mônica C. Ledur
2,
Jane de O. Peixoto
2,
Daniele B. D. Marques
1,
Karine A. Costa
1,
Luiz. L. Coutinho
3 and
Simone E. F. Guimarães
1,*
1
Animal Science Department, Universidade Federal de Viçosa, Viçosa 36570-900, Brazil
2
Animal Genetics Laboratory, Embrapa Swine and Poultry National Research Center, Concórdia 89715-899, Brazil
3
Functional Genomics Center, ESALQ, Universidade de São Paulo, Piracicaba 89715-899, Brazil
*
Author to whom correspondence should be addressed.
Genes 2019, 10(12), 1010; https://0-doi-org.brum.beds.ac.uk/10.3390/genes10121010
Submission received: 30 October 2019 / Revised: 20 November 2019 / Accepted: 27 November 2019 / Published: 5 December 2019
(This article belongs to the Special Issue Genomics of Sexual Development and Reproduction in Mammals)
Browse Figure
Review Reports Versions Notes

Abstract

:
Sexual dimorphism is a relevant factor in animal science, since it can affect the gene expression of economically important traits. Eventually, the interest in the prenatal phase in a transcriptome study may not comprise the period of development in which male and female conceptuses are phenotypically divergent. Therefore, it would be interesting if sex differentiation could be performed using transcriptome data, with no need for extra techniques. In this study, the sex of pig conceptuses (embryos at 25 days-old and fetuses at 35 days-old) was determined by reads counts per million (CPM) of Y chromosome-linked genes that were discrepant among samples. Thus, ten genes were used: DDX3Y, KDM5D, ZFY, EIF2S3Y, EIF1AY, LOC110255320, LOC110257894, LOC396706, LOC100625207, and LOC110255257. Conceptuses that presented reads CPM sum for these genes (ΣCPMchrY) greater than 400 were classified as males and those with ΣCPMchrY below 2 were classified as females. It was demonstrated that the sex identification can be performed at early stages of pig development from RNA-sequencing analysis of genes mapped on Y chromosome. Additionally, these results reinforce that sex determination is a mechanism conserved across mammals, highlighting the importance of using pigs as an animal model to study sex determination during human prenatal development.

1. Introduction

Maternal nutrition may affect progeny anatomy, physiology and metabolism during critical periods of prenatal life, which is described as fetal programming [1]. The effects of maternal nutrition on fetal programming have been recently addressed in several studies using pigs as a model organism [2,3,4,5,6], since there are striking similarities between pigs and humans regarding the anatomy, physiology, metabolism, and nutrition, which provide the basis for the use of this animal in several studies. In addition, Sus scrofa is a prolific specie with a relatively short developmental period, besides being a cheaper and easier mammalian model organism to maintain [7,8]. In this context, transcriptome studies using RNA-sequencing (RNA-seq) performed in pigs highlight the relevance of the data obtained in this specie for human-related researches.
The RNA sequencing (RNA-seq) approach generates a comprehensive picture of gene expression levels at different developmental stages and physiological conditions [9,10,11]. However, these analyses may be influenced by sexual dimorphism, since males and females are marked by significant biological differences that can modify the gene expression of economically important traits [12]. Therefore, researchers should take into account the conceptuses sex and sex ratios in their research before drawing any conclusion [13].
In mammals, sex determination is basically organized in four stages: (1) Determination of the chromosomal sex, which is established at fertilization when a Y- or a X-chromosome sperm fuses with the oocyte to determine the zygote genetic sex; (2) the differentiation of gonads into testicles or ovaries from the presence or absence of sex-determining region Y gene (SRY) at the critical time window during embryonic development; (3) the differentiation of male and female internal and external genitalia from the non-sexually dimorphic structures present in the embryo, and (4) anatomical and physiological differences, i.e. sexual differentiation. Thus, all sex determination stages are mainly related to general gonadal function [14,15].
At early pig development, the absence of gonadal phenotype may impair the accurate visual sex determination, since the tunica albuginea, a marker of testis formation, is histologically identified only at 27 days [16] and the beginning of testicular descent occurs around 60 days [17]. Therefore, if research projects are designed at early development stages in which male and female conceptuses are not phenotypically divergent and there is response to sexual dimorphism, the sex determination could be performed using alternative tools rather than phenotypic differences.
Several studies have demonstrated how to identify the conceptuses sex at early stages of development using molecular techniques, as PCR (polymerase chain reaction) from Y chromosome-linked genes [13,18,19,20]. The PCR provides sensitive, precise, rapid, and reliable results [20]. However, if a RNA-seq study has already been performed on conceptuses to obtain a broader knowledge of the transcriptome under experimental conditions (e.g. inclusion of additives on maternal diets), these data can also be applied to determine the conceptuses sex with no need for extra techniques. A similar approach has been used by Petropoulos et al. [21] to determine the sex of human embryos through single-cell sequencing.
Therefore, considering that in early prenatal pig development male and female conceptuses are not phenotypically divergent and the sexual dimorphism may affect the transcriptome analyses, we aimed to demonstrate that the sex of 25 and 35 days-old conceptuses can be determined by read counting of genes mapped on Y chromosome using already available RNA-seq data.

2. Materials and Methods

2.1. Experimental Animals and Design

The experimental protocols used in this study have followed ethical principles in animal research (CONCEA, 2016) and were approved by the Ethical Committee on Animal Use of the Universidade Federal de Viçosa (UFV), MG, Brazil [protocol # 06/2017].
All experimental protocols were performed according to the experimental design described in Costa et al. [22]. Briefly, 24 hours after the second insemination, 11 gilts received a basal diet for pregnant animals without supplementation (CONT) and 12 gilts received a CONT diet supplemented with 1.0% L-arginine (ARG), and two gestational ages were considered (25 and 35 days). From these gilts, 20 became pregnant. Additional details regarding reproductive management and nutritional information of the diets have been previously described in Costa et al. [22].
At 25 days of gestation, five females of CONT (n = 5) and five females of ARG (n = 5) were rendered unconscious using head-only electrical stunning (240V, 1.3A) and immediately exsanguinated. The same procedure was followed for four females of CONT (n = 4) and six females of ARG (n = 6) at 35 days of gestation. After slaughter, four conceptuses were collected per female at each gestational age, totaling in average 20 conceptuses per treatment.
The conceptuses were quickly washed with PBS (Phosphate Buffered Saline) solution, individually identified, stored in liquid nitrogen and transported to the Animal Biotechnology Laboratory (LABTEC) at the Department of Animal Science, UFV. At LABTEC, each conceptus (embryos from gilts slaughtered at 25 days and fetuses from gilts slaughtered at 35 days) was entirely and separately macerated in liquid nitrogen. Although an average of five gilts was initially used per treatment in the original trial described in Costa et al. [22], subsequent analyses were performed on a subset of gilts (n = 3) per treatment, since at least three biological replicates are recommended for describing results in RNA-seq experiments [9,23]. Moreover, three embryos from each CONT gilt (n = 3 gilts) and three embryos from each ARG gilt (n = 3 gilts) at 25 days of gestation (25DC and 25DA, respectively) and three fetuses from each CONT gilt (n = 3 gilts) and three fetuses from each ARG gilt (n = 3 gilts) at 35 days of gestation (35DC and 35DA, respectively), totaling 36 samples, were randomly chosen for RNA-seq analysis. Afterwards, these samples were transported in liquid nitrogen to the Animal Genetics Laboratory at the Embrapa Swine and Poultry National Research Center, Concordia, SC, Brazil, for further RNA extraction and library preparation.

2.2. RNA Extraction and Library Preparation

Total RNA extraction of the 25- and 35-day-old conceptuses from the ARG and CONT females was performed with TRIzol (Invitrogen, San Diego, CA, USA). The macerated conceptuses (100 mg) and TRIzol (1 mL) were mixed with vortex and then incubated for 5 minutes at room temperature (RT, 25 °C). Then, 200 µL of chloroform were added, shaking vigorously for 15 seconds and incubated at RT for 5 minutes. Centrifugation was performed at 11,000× g at 4 °C for 15 minutes. Approximately 600 μL of the clear upper aqueous phase containing only RNA were carefully removed and transferred to a new tube, and 600 μL of 70% ethanol were added and homogenized by inversion. This volume was added to the silica column RNeasy mini kit (Qiagen, Hilden, Germany) and centrifuged for 15 seconds at 8000× g. The eluate was discarded and 700 μL RW1 buffer were added, followed by centrifugation for 15 seconds at 8000× g. Two washes with 500 μL RPE buffer were done and, finally, RNAs were eluted in 50 μL of RNase free water.
After RNA extraction, the quantification was performed in a QUBIT fluorimeter (Thermo Scientific, Waltham, MA, USA) and the integrity was determined in 1.0% agarose gel. In addition, the Agilent 2100 BioAnalyzer (Agilent Technologies, Santa Clara, CA, USA) was used for integrity measurement, in which samples with RNA integrity number (RIN) higher than eight were used for library preparation. In this sense, one fetus sample was excluded from subsequent analysis due to low RIN score (RIN < 8). Therefore, the 35 remaining conceptuses samples (nine ARG and nine CONT at 25 days and eight ARG and nine CONT at 35 days) were submitted to RNA-seq library preparation using the TruSeq Stranded mRNA Library Prep Kit (Illumina, Inc., San Diego, CA, USA), followed by purification of the poly-A tail using 2 μg of total RNA, according to the manufacturer´s recommendations.

2.3. Sequencing, Quality Control and Mapping

The libraries were sequenced in Illumina HiSeq2500 (Illumina, Inc.; San Diego, CA, USA), following the 2 × 100 bp paired-end protocol, at the Functional Genomics Center, ESALQ, Universidade de São Paulo, Piracicaba, SP, Brazil. The FASTQ files were deposited in the SRA database, with Bioproject number PRJNA576701 and Biosample numbers SAMN13003023, SAMN13003024, SAMN13003025, SAMN13003026, SAMN13003027, SAMN13003028, SAMN13003029, SAMN13003030, SAMN13003031, SAMN13003032, SAMN13003033, SAMN13003034, SAMN13003035, SAMN13003036, SAMN13003037, SAMN13003038, SAMN13003039, SAMN13003040, SAMN13003041, SAMN13003042, SAMN13003043, SAMN13003044, SAMN13003045, SAMN13003046, SAMN13003047, SAMN13003048, SAMN13003049, SAMN13003050, SAMN13003051, SAMN13003052, SAMN13003053, SAMN13003054, SAMN13003055, SAMN13003056 and SAMN13003057.
The BAQCOM pipeline () [24] was used to perform quality control (QC), mapping and reads counting. BAQCOM uses a set of software to analyze RNA-seq data, as Trimmomatic [25], version 0.38, to identify and remove adapter and low quality sequences, HISAT2 [26], version 0.11.2, to map reads against the genome and HTseq-count [27], version 0.11.2, to count reads in features. Using BAQCOM pipeline, only reads with Phred quality ≥ 20 and length ≥ 70 pb were mapped against the pig reference genome (Sus scrofa, v. 11.1) and the reads counting for each feature was based on Ensembl annotation release 95.

2.4. Sex Identification

The conceptuses sex determination (n = 18 embryos and n = 17 fetuses) was performed using discrepant reads counts per million (CPM) of Y chromosome-linked genes from the conceptuses RNA-seq dataset. After the samples classification as males or females, the HTseq files containing all samples transcripts were renamed and used to build a multi-dimensional scale (MDS) plot using plotMDS function of EdgeR package [28] of R software [29] to evaluate conceptuses dispersion.

3. Results

Averages of 15.3 and 14.4 million reads/samples were generated for embryos and fetuses, respectively. After the data QC, averages of 13.6 million reads/embryos and 12.9 million reads/fetuses remained for further analyses. More than 98% of the reads were mapped against the pig reference genome (Sus scrofa, v.11.1), with 83.8% and 79.8% of the reads counted into genes for embryos and fetuses, respectively. An average of 0.05% of the reads was uniquely mapped on Y chromosome. At 25 and 35 days of age, male and female conceptuses could not be phenotypically differentiated, since gonadal phenotype could not be observed (Figure S1).
From 79 Y chromosome-linked genes, ten showed discrepant reads CPM among samples, i.e, high and low reads CPM values (Table 1 and Table 2), and were selected for sex determination: DDX3Y (ATP-dependent RNA helicase DDX3X), KDM5D (Lysine demethylase 5D), ZFY (Zinc Finger Protein, Y-Linked), EIF2S3Y (eukaryotic translation initiation factor 2 subunit 3), EIF1AY (eukaryotic translation initiation factor 1A, Y-linked), LOC110255320 (lysine-specific demethylase 6A-like), LOC110257894 (gamma-taxilin-like), LOC396706 (U2 small nuclear ribonucleoprotein auxiliary factor 35 kDa subunit-related protein 2), LOC100625207 (probable ubiquitin carboxyl-terminal hydrolase FAF-X) and LOC110255257 (oral-facial-digital syndrome 1 protein-like). The CPM sum for these genes (ΣCPMchrY) was calculated for each conceptus in order to determine its sex. In this way, conceptuses that presented ΣCPMchrY greater than 400 were classified as males and those with ΣCPMchrY below 2 were classified as females. Male conceptuses showed mean (standard deviation), minimum and maximum ΣCPMchrY of 499.08 (44.74), 427.29, and 573.09, respectively (Table 1), while females showed mean (standard deviation), minimum and maximum ΣCPMchrY of 0.38 (0.42), 0.00 and 1.75, respectively (Table 2). In the MDS plot, a segregation of conceptuses into two distinct sex groups was observed (Figure 1).

4. Discussion

The control of animals sex ratio is desirable in livestock production [20], since the expression of genes that affect a wide range of economically important traits may be influenced by sexual dimorphism [12]. Considering that the sex determination is equivalent to testis determination [30], the testis could be easily used to determine male and female individual rates. However, in pig conceptuses, the tunica albuginea, a marker of testis formation, is histologically identified only at 27 days [16] and the beginning of testicular descent occurs around 60 days [17]. In this context, we demonstrated that it is possible to determine the sex of pig embryos (25 days) and fetuses (35 days) using reads CPM of ten Y chromosome-linked genes through RNA-seq analysis, with no need to use other molecular tools for this purpose, such as PCR. Although most often considered a cheaper and faster technique than RNA-seq, PCR would represent another cost source in this study, since RNA-seq analyses had been previously performed aiming to test a hypothesis on a maternal nutrition experiment [22].
Several studies have also used Y chromosome-linked genes to determine the sex of pigs at early development stages [13,18,19], however, all of them have applied DNA-based approaches. The Y chromosome harbors genes that are essential for testis development and function, such as the master gene for testis determination (SRY), which is expressed during a critical period of embryonic development, and other genes that are important for spermatogenic function, as the ZFY, EIF1AY and KDM5D. The importance of these Y-linked genes during male development has been elucidated for mammals [31,32,33,34]. These genes are located in the non-recombining region of the Y chromosome (NRY), which is also known as the male-specific Y (MSY) region, transmitted from fathers to sons without recombination with the X chromosome [35]. The MSY region contains only 27 genes coding for distinct proteins, as the SRY, EIF1AY, KDM5D, ZFY, EIF2S3Y, and DDX3Y [31,35,36].
In this study, among all abovementioned genes, SRY was not used to differentiate male and female conceptuses (data not shown). This gene seems to follow a strict expression pattern during mammalian embryonic development. In mouse, its expression starts at 10.5 days post coitum (dpc), reaches a peak at 11.5 dpc and then wanes by 12.5 dpc [32,37]. In pigs, the SRY expression is detectable in the embryos genital ridge around 21 and 23 days, with faint expression at 26 days and absence of expression at 31 days [38]. Consequently, both 25 and 35 days of pig prenatal periods do not comprise the main SRY expression phase, reinforcing the critical time window during embryonic development in which SRY is expressed.
Subsequently to the sex determination by the SRY gene, other male-specific genes are important to maintain testicular development and spermatogenesis [39]. According to Bellott et al. [40], during spermatogenesis, each stage of the molecular central dogma is regulated by the following Y-linked genes: KDM5D, ZFY, and EIF1AY. The KDM5D gene encodes a histone demethylase enzyme that acts removing tri- and di-methylations of lysine 4 of histone H3 at the start site of transcription in actively transcribed genes [41]. The transcription factor ZFY, the first coding gene identified in the human Y chromosome, regulates the transcription of a number of Y-linked genes [31] and mediates multiple aspects of spermatogenesis and reproduction, such as morphology, motility, capacitation, acrosome reaction, and oocyte activation, as well as chromosomal aberrations [42]. The translation initiation factor EIF1AY is requested for a high rate of protein biosynthesis and the absence of this gene may contribute to failures in the spermatogenic process [43].
Another important translation initiation factor during spermatogenesis is the EIF2S3Y gene, located in the short arm of the Y chromosome. This gene has been identified as a key mouse-specific regulator of spermatogonia proliferation and differentiation [44]. Finally, the DDX3Y gene, along with KDM5D and EIF1AY, are azoospermia factors (AZFa) [45]. Azoospermia is defined as the absence of sperm in at least two different ejaculates [31]. The DDX3Y is the major AZFa gene and belongs to a highly conserved subfamily of the DEAD-box RNA helicase family [33]. This gene family is involved in oogenesis and spermatogenesis, silencing mobile elements and other repetitive genomic regions in germinal tissues, and is important for gonads formation during embryo development [46]. Therefore, the KDM5D, ZFY, EIF1AY, EIF2S3Y, and DDX3Y genes showed discrepant reads CPM among samples and have important spermatogenic and gonadal functions, being crucial to differentiate male and female pig conceptuses in the present study. Similar pattern has already been observed in humans blastocyst [21], which highlights the importance of using pigs as an animal model to study sex determination throughout human prenatal development.
The reads CPM observed for the Y-linked genes in female conceptuses may be considered technical artifacts from RNA-seq, since mammalian X and Y chromosomes share regions of high sequence similarity, causing mismapping of short reads to a reference genome [47]. In this context, Ballouz et al. [48] found reads of female samples mapped in Y chromosome, which was also characterized as reads mismapping. Although considered technical artifacts, these Y chromosome reads CPM found in females were not a major issue in the present study, since the discrepant reads CPM among samples for the ten Y chromosome-linked genes were sufficient to differentiate males and females.
The LOC110255257, LOC100625207, LOC110255320, LOC110257894, and LOC396706 genes are located in the Y chromosome, however, they have not yet been fully characterized in Sus scrofa genome. According to Skinner et al. [49], there are few data available on the porcine Y chromosome. Therefore, although these genes are located on this chromosome and have shown contribution to sex determination of pig conceptuses due to the discrepant reads CPM among samples, their functions will not be addressed in this study. Further studies might be performed to characterize these Y-linked genes in pigs.

5. Conclusions

We have demonstrated that sex identification can be performed from read counts per million of genes mapped in the Y chromosome at early stages of pig development using already available conceptuses RNA-seq data, with no need for extra techniques. Additionally, the results reinforce that molecular sex determination is a mechanism conserved across mammalian species, since the role of these genes is also important during male development in humans, highlighting the applicability of pigs as an animal model to study sex determination during human prenatal development.

Supplementary Materials

The following is available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4425/10/12/1010/s1, Figure S1: 25 days-old pig embryo (left) and 35 days-old pig fetus (right) do not present evident gonadal phenotype.

Author Contributions

S.A.T., A.M.G.I., M.C.L., J.d.O.P., and S.E.F.G.; methodology, S.A.T., A.M.G.I., M.E.C., H.C.d.O., L.L.C.; software, S.A.T., A.M.G.I., M.E.C., H.C.d.O.; formal Analysis, S.A.T., A.M.G.I., M.E.C., H.C.d.O. and D.B.D.M.; writing—Original draft preparation S.A.T.; writing—Review and editing, A.M.G.I., M.C.L., J.d.O.P., M.E.C., H.C.d.O., D.B.D.M., K.A.C. and S.E.F.G.; supervision, S.E.F.G.; project administration, S.E.F.G.; funding acquisition, M.C.L., J.d.O.P., S.E.F.G.

Funding

This research was funded by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq, Grant # 423558/2018-8), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES, Grant numbers PROEX 0577/2018 and 320.020.170.11pq), Instituto Nacional de Ciência e Tecnologia- Ciência Animal (INCT-CA, Grant # 465377/2014-9) and Fundação de Amparo à Pesquisa do Estado de Minas Gerais (FAPEMIG, Grant # CVZ APQ 03434/16).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Godfrey, K.M.; Barker, D.J.P. Fetal programming and adult health. Public Health Nutr. 2001, 4, 611–624. [Google Scholar] [CrossRef] [PubMed]
  2. Li, J.; Xia, H.; Yao, W.; Wang, T.; Li, J.; Piao, X.; Thacker, P.; Wu, G.; Wang, F. Effects of arginine supplementation during early gestation (day 1 to 30) on litter size and plasma metabolites in gilts and sows 1. J. Anim. Sci. 2015, 93, 5291–5303. [Google Scholar] [CrossRef] [PubMed]
  3. Wu, G.; Bazer, F.W.; Cudd, T.A.; Meininger, C.J.; Spencer, T.E. Recent Advances in Nutritional Sciences Maternal Nutrition and Fetal Development. J. Nutr. 2004, 134, 2169–2172. [Google Scholar] [CrossRef] [PubMed]
  4. Zhang, S.; Heng, J.; Song, H.; Zhang, Y.; Lin, X.; Tian, M.; Chen, F.; Guan, W. Role of Maternal Dietary Protein and Amino Acids on Fetal Programming, Early Neonatal Development, and Lactation in Swine. Animals 2019, 9, 19. [Google Scholar] [CrossRef] [Green Version]
  5. Assis, K.; Botelho, D.; Marques, D.; Filardi, C.; Campos, D.; Saraiva, A.; Domingos, J.; Eliza, S.; Guimarães, F. Nutrition in fluence on sow reproductive performance and conceptuses development and survival: A review about L-arginine supplementation. Livest. Sci. 2019, 228, 97–103. [Google Scholar]
  6. Peñagaricano, F.; Wang, X.; Rosa, G.J.M.; Radunz, A.E.; Khatib, H. Maternal nutrition induces gene expression changes in fetal muscle and adipose tissues in sheep. BMC Genom. 2014, 15, 1034. [Google Scholar] [CrossRef] [Green Version]
  7. Verma, N.; Rettenmeier, A.W.; Schmitz-Spanke, S. Recent advances in the use of Sus scrofa (pig) as a model system for proteomic studies. Proteomics 2011, 11, 776–793. [Google Scholar] [CrossRef]
  8. Papatheodorou, I.; Fonseca, N.A.; Keays, M.; Tang, Y.A.; Barrera, E.; Bazant, W.; Burke, M.; Füllgrabe, A.; Fuentes, A.M.P.; George, N.; et al. Expression Atlas: Gene and protein expression across multiple studies and organisms. Nucleic Acids Res. 2018, 46, D246–D251. [Google Scholar] [CrossRef]
  9. Conesa, A.; Madrigal, P.; Tarazona, S.; Gomez-cabrero, D.; Cervera, A.; Mcpherson, A.; Szcze, W.; Gaffney, D.J.; Elo, L.L.; Zhang, X. A survey of best practices for RNA-seq data analysis. Genome Biol. 2016, 17, 13. [Google Scholar] [CrossRef] [Green Version]
  10. Wang, Z.; Gerstein, M.; Snyder, M. RNA-Seq: A revolutionary tool for transcriptomics. Nat. Rev. Genet. 2009, 10, 57–63. [Google Scholar] [CrossRef]
  11. Sanglard, L.P.; Nascimento, M.; Moriel, P.; Sommer, J.; Ashwell, M.; Poore, M.H.; Duarte, M.D.S.; Serão, N.V.L. Impact of energy restriction during late gestation on the muscle and blood transcriptome of beef calves after preconditioning. BMC Genom. 2018, 19, 702. [Google Scholar] [CrossRef] [PubMed]
  12. Cui, L.; Zhang, J.; Ma, J.; Guo, Y.; Li, L.; Xiao, S.; Ren, J.; Yang, B.; Huang, L. Sexually dimorphic genetic architecture of complex traits in a large-scale F2 cross in pigs. Genet. Sel. Evol. 2014, 46, 76. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Blanes, M.S.; Tsoi, S.C.M.; Dyck, M.K. Accurate and Phenol Free DNA Sexing of Day 30 Porcine Embryos by PCR. J. Vis. Exp. 2016, 108, e53301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Miura, K.; Tomita, A.; Kanai, Y. Sex Determination and Differentiation in Mammals. In Reproductive and Developmental Strategies; Kazuya, K., Takeshi, K., Yasuhiro, I., Mariko, K., Eds.; Springer: Tokyo, Japan, 2018; pp. 407–433. ISBN 9784431566090. [Google Scholar]
  15. Goodfellow, P.N.; Darling, S.M. Genetics of sex determination in man and mouse. Development 1988, 102, 251–258. [Google Scholar] [PubMed]
  16. Pelliniemi, L.J. Ultrastructure of the early ovary and testis in pig embryos. Am. J. Anat. 1975, 144, 89–111. [Google Scholar] [CrossRef]
  17. Hyttel, P.; Sinowatz, F.; Vejlsted, M.; Betteridge, K. Essentials of Domestic Animal Embryology; Edwards, R., Rodenhuis, J., Betteridge, K., Eds.; Saunders: Philadelphia, PA, USA, 2010; ISBN 978-0-7020-2899-1. [Google Scholar]
  18. Tsoi, S.; Blanes, M.; Chen, T.Y.; Langendijk, P.; Athorn, R.; Foxcroft, G.; Dyck, M. Identification of differentially expressed genes in sexed pig embryos during post-hatching development in primiparous sows exposed to differing intermittent suckling and breeding strategies. Genom. Data 2016, 9, 30–34. [Google Scholar] [CrossRef] [Green Version]
  19. Bick, J.T.; Flöter, V.L.; Robinson, M.D.; Bauersachs, S.; Ulbrich, S.E. Small RNA-seq analysis of single porcine blastocysts revealed that maternal estradiol-17beta exposure does not affect miRNA isoform (isomiR) expression. BMC Genom. 2018, 19, 590. [Google Scholar] [CrossRef] [Green Version]
  20. Han, S.H.; Yang, B.C.; Ko, M.S.; Oh, H.S.; Lee, S.S. Length difference between equine ZFX and ZFY genes and its application for molecular sex determination. J. Assist. Reprod. Genet. 2010, 27, 725–728. [Google Scholar] [CrossRef] [Green Version]
  21. Petropoulos, S.; Edsgärd, D.; Reinius, B.; Deng, Q.; Panula, S.P.; Codeluppi, S.; Reyes, A.P.; Linnarsson, S.; Sandberg, R.; Lanner, F. Single-Cell RNA-Seq Reveals Lineage and X Chromosome Dynamics in Human Preimplantation Embryos. Cell 2016, 165, 1012–1026. [Google Scholar] [CrossRef] [Green Version]
  22. Costa, K.A.; Saraiva, A.; Guimaraes, J.D.; Marques, D.B.D.; Machado-neves, M.; Barbosa, L.M.R.; Villadiego, F.A.C.; Veroneze, R.; de Oliveira, L.F.; Garcia, I.S.; et al. Dietary L-arginine supplementation during early gestation of gilts affects conceptuses development. Theriogenology 2019, 140, 62–71. [Google Scholar] [CrossRef]
  23. Nie, H.; Wang, Y.; Su, Y.; Hua, J. Exploration of miRNAs and target genes of cytoplasmic male sterility line in cotton during flower bud development. Funct. Integr. Genom. 2018, 18, 457–476. [Google Scholar] [CrossRef] [PubMed]
  24. BAQCOM: Bioinformatics Analysis for Quality Control and Mapping. Available online: https://github.com/hanielcedraz/BAQCOM (accessed on 15 April 2019).
  25. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Anders, S.; Pyl, P.T.; Huber, W. HTSeq—A Python framework to work with high-throughput sequencing data. Bioinformatics 2015, 31, 166–169. [Google Scholar] [CrossRef] [PubMed]
  28. Robinson, M.D.; Mccarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. R Core Team. A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  30. Koopman, P.; Gubbay, J.; Vivian, N.; Goodfellow, P.; Lovell-Badge, R. Male development of chromosomally female mice transgenic for Sry. Nature 1991, 351, 117–121. [Google Scholar] [CrossRef]
  31. Krausz, C.; Casamonti, E. Spermatogenic failure and the Y chromosome. Hum. Genet. 2017, 136, 637–655. [Google Scholar] [CrossRef]
  32. Kashimada, K.; Koopman, P. Sry: The master switch in mammalian sex determination. Development 2010, 137, 3921–3930. [Google Scholar] [CrossRef] [Green Version]
  33. Rauschendorf, M.A.; Zimmer, J.; Ohnmacht, C.; Vogt, P.H. DDX3X, the X homologue of AZFa gene DDX3Y, expresses a complex pattern of transcript variants only in the male germ line. Mol. Hum. Reprod. 2014, 20, 1208–1222. [Google Scholar] [CrossRef]
  34. Peneder, P.; Wallner, B.; Vogl, C. Exchange of genetic information between therian X and Y chromosome gametologs in old evolutionary strata. Ecol. Evol. 2017, 7, 8478–8487. [Google Scholar] [CrossRef]
  35. Fan, H.; Hu, Y.; Shan, L.; Yu, L.; Wang, B.; Li, M.; Wu, Q.; Wei, F. Synteny search identifies Carnivore Y chromosome for evolution of male specific genes. Integr. Zool. 2018, 14, 224–234. [Google Scholar] [CrossRef] [PubMed]
  36. Maan, A.A.; Eales, J.; Akbarov, A.; Rowland, J.; Xu, X.; Jobling, M.A.; Charchar, F.J.; Tomaszewski, M. The y chromosome: A blueprint for men’s health? Eur. J. Hum. Genet. 2017, 25, 1181–1188. [Google Scholar] [CrossRef] [PubMed]
  37. Hacker, A.; Capel, B.; Goodfellow, P.; Lovell-Badge, R. Expression of Sry, the mouse sex determining gene. Development 1995, 121, 1603–1614. [Google Scholar] [PubMed]
  38. Daneau, I.; Ethier, J.-F.; Lussier, J.G.; Silversides, D.W. Porcine SRY Gene Locus and Genital Ridge Expression. Biol. Reprod. 2005, 55, 47–53. [Google Scholar] [CrossRef] [Green Version]
  39. Akimoto, C.; Kitagawa, H.; Matsumoto, T.; Kato, S. Spermatogenesis-specific association of SMCY and MSH5. Genes Cells 2008, 13, 623–633. [Google Scholar] [CrossRef]
  40. Bellott, D.W.; Hughes, J.F.; Skaletsky, H.; Brown, L.G.; Pyntikova, T.; Cho, T.; Koutseva, N.; Zaghlul, S.; Graves, T.; Rock, S.; et al. Mammalian Y chromosomes retain widely expressed dosage-sensitive regulators. Nature 2014, 508, 494–499. [Google Scholar] [CrossRef]
  41. Rasmussen, P.B.; Staller, P. The KDM5 family of histone demethylases as targets in oncology drug discovery. Epigenomics 2014, 6, 277–286. [Google Scholar] [CrossRef]
  42. Nakasuji, T.; Ogonuki, N.; Chiba, T.; Kato, T.; Shiozawa, K.; Yamatoya, K.; Tanaka, H.; Kondo, T.; Miyado, K.; Miyasaka, N.; et al. Complementary Critical Functions of Zfy1 and Zfy2 in Mouse Spermatogenesis and Reproduction. PLoS Genet. 2017, 13, e1006578. [Google Scholar] [CrossRef]
  43. Kleiman, S.E.; Yogev, L.; Hauser, R.; Botchan, A.; Maymon, B.B.S.; Paz, G.; Yavetz, H. Expression profile of AZF genes in testicular biopsies of azoospermic men. Hum. Reprod. 2007, 22, 151–158. [Google Scholar] [CrossRef] [Green Version]
  44. Mazeyrat, S.; Saut, N.; Grigoriev, V.; Mahadevaiah, S.K.; Ojarikre, O.A.; Rattigan, Á.; Bishop, C.; Eicher, E.M.; Mitchell, M.J.; Burgoyne, P.S. A Y-encoded subunit of the translation initiation factor Eif2 is essential for mouse spermatogenesis. Nat. Genet. 2002, 29, 49–53. [Google Scholar] [CrossRef]
  45. Mitchell, M.J.; Metzler-Guillemain, C.; Toure, A.; Coutton, C.; Arnoult, C.; Ray, P.F. Single gene defects leading to sperm quantitative anomalies. Clin. Genet. 2017, 91, 208–216. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Kotov, A.A.; Akulenko, N.V.; Kibanov, M.V.; Olenina, L.V. DEAD-Box RNA Helicases in Animal Gametogenesis. Mol. Biol. 2014, 48, 16–28. [Google Scholar] [CrossRef]
  47. Webster, T.H.; Couse, M.; Grande, B.M.; Karlins, E.; Phung, T.N.; Richmond, P.A.; Whitford, W.; Wilson, M.A. Identifying, understanding, and correcting technical artifacts on the sex chromosomes in next-generation sequencing data. GigaScience 2019, 8, giz074. [Google Scholar] [CrossRef] [Green Version]
  48. Ballouz, S.; Dobin, A.; Gingeras, T.R.; Gillis, J. The fractured landscape of RNA-seq alignment: The default in our STARs. Nucleic Acids Res. 2018, 46, 5125–5138. [Google Scholar] [CrossRef]
  49. Skinner, B.M.; Sargent, C.A.; Churcher, C.; Hunt, T.; Herrero, J.; Loveland, J.E.; Dunn, M.; Louzada, S.; Fu, B.; Chow, W.; et al. The pig X and Y chromosomes: Structure, sequence and evolution. Genome Res. 2016, 26, 130–139. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Genes 10 01010 g001 550
Figure 1. Multi-dimensional scale (MDS) plot showing segregation of pig conceptuses (25D: embryos at 25 days-old and 35D: fetuses at 35 days-old) into two distinct sex groups. Red: Female conceptuses; Black: Male conceptuses; Samples identification: male embryos from supplemented (E-25DA-M) or non-supplemented (E-25DC-M) gilts; male fetuses from supplemented (F-35DA-M) or non-supplemented (F-35DC-M) gilts; female embryos from supplemented (E-25DA-F) or non-supplemented (E-25DC-F) gilts; female fetuses from supplemented (F-35DA-F) or non-supplemented (F-35DC-F) gilts.
Figure 1. Multi-dimensional scale (MDS) plot showing segregation of pig conceptuses (25D: embryos at 25 days-old and 35D: fetuses at 35 days-old) into two distinct sex groups. Red: Female conceptuses; Black: Male conceptuses; Samples identification: male embryos from supplemented (E-25DA-M) or non-supplemented (E-25DC-M) gilts; male fetuses from supplemented (F-35DA-M) or non-supplemented (F-35DC-M) gilts; female embryos from supplemented (E-25DA-F) or non-supplemented (E-25DC-F) gilts; female fetuses from supplemented (F-35DA-F) or non-supplemented (F-35DC-F) gilts.
Genes 10 01010 g001
Table
Table 1. Reads counts per million of 10 genes located in the Y chromosome of male pig conceptuses.
Table 1. Reads counts per million of 10 genes located in the Y chromosome of male pig conceptuses.
Conceptuses Sample_ID 1Counts Per Million (CPM)ΣCPMchrY 3
DDX3Y 2KDM5DZFYEIF2S3YEIF1AYLOC110255320LOC110257894LOC396706LOC100625207LOC110255257
E1-25DA-M151.467.5224.81162.2340.6412.049.711.5361.7315.17486.84
E2-25DA-M166.547.6431.68197.7345.1612.3011.811.9561.0017.72553.53
E3-25DC-M167.219.5628.67211.9745.8613.2312.792.2166.5215.07573.09
E4-25DC-M170.7610.1529.78189.7946.3911.4312.412.1171.8117.90562.53
E5-25DC-M149.286.5624.94180.7836.5513.4710.471.4355.2916.46495.23
E6-25DC-M156.649.1930.92192.2340.7612.1410.992.3068.5616.24539.97
E7-25DC-M166.109.4333.57192.2344.9212.3111.431.8471.7817.03560.64
F1-35DC-M193.415.3526.07130.6744.659.796.251.2352.4613.16483.04
F2-35DC-M167.287.1128.43135.3839.8410.257.191.4156.2015.04468.13
F3-35DC-M168.097.0531.32141.5337.8010.967.261.7858.0913.03476.91
F4-35DA-M157.518.1726.77119.4533.459.306.462.0651.2712.85427.29
F5-35DA-M200.338.0424.02126.1345.0912.206.562.1345.8312.20482.53
F6-35DA-M187.016.4724.02128.6141.5810.358.012.0245.5413.43467.04
F7-35DA-M180.787.0629.27133.4241.1411.746.572.1050.7813.20476.06
F8-35DC-M178.1510.6833.04155.5743.4911.299.232.0662.9417.32523.77
F9-35DA-M143.3310.2827.83135.3331.4212.338.651.5553.7912.73437.24
F10-35DA-M173.5110.4026.69128.3236.7112.626.881.9158.1915.37470.60
1 Conceptuses Sample Identification: Seven male embryos samples from supplemented (E-25DA-M) and non-supplemented (E-25DC-M) gilts; ten male fetuses samples from supplemented (F-35DA-M) and non-supplemented (F-35DC-M) gilts; 2 DDX3Y: ATP-dependent RNA helicase DDX3X; KDM5D: Lysine demethylase 5D; ZFY: Zinc Finger Protein Y-Linked; EIF2S3Y: eukaryotic translation initiation factor 2 subunit 3; EIF1AY: eukaryotic translation initiation factor 1A, Y-linked; LOC110255320: Lysine-specific demethylase 6A-like; LOC110257894: Gamma-taxilin-like; LOC396706: U2 small nuclear ribonucleoprotein auxiliary factor 35 kDa subunit-related protein 2; LOC100625207: Probable ubiquitin carboxyl-terminal hydrolase FAF-X; LOC110255257: Oral-facial-digital syndrome 1 protein-like; 3 ΣCPMchrY: sum of reads CPM of 10 Y chromosome-linked genes in male conceptuses.
Table
Table 2. Reads counts per million of 10 Y chromosome-linked genes that showed reads counting in female pig conceptuses.
Table 2. Reads counts per million of 10 Y chromosome-linked genes that showed reads counting in female pig conceptuses.
Conceptuses Sample_ID 1Counts Per Million (CPM)ΣCPMchrY 3
DDX3Y2KDM5DZFYEIF2S3YEIF1AYLOC110255320LOC110257894LOC396706LOC100625207LOC110255257
E1-25DA-F0.000.070.000.000.070.000.000.000.000.000.14
E2-25DA-F0.000.000.000.000.070.000.000.000.000.070.14
E3-25DC-F0.000.000.070.070.000.000.000.000.000.000.14
E4-25DC-F0.280.000.070.070.070.000.000.000.000.000.49
E5-25DC-F0.000.000.000.150.000.000.000.000.070.150.37
E6-25DC-F0.000.000.140.000.070.000.000.000.000.000.21
E7-25DC-F0.000.000.150.080.000.000.000.000.000.000.23
E8-25DC-F0.000.000.000.000.000.000.000.000.000.000.00
E9-25DC-F0.000.000.000.000.070.000.000.000.000.070.14
E10-25DC-F0.000.000.080.080.080.000.000.000.000.000.24
E11-25DA-F0.130.000.070.130.070.000.000.000.200.000.60
F1-35DA-F0.080.000.080.000.000.000.000.000.000.080.24
F2-35DA-F0.000.080.230.000.080.000.000.000.000.000.39
F3-35DA-F0.160.000.000.250.080.000.000.000.080.000.57
F4-35DC-F0.000.000.070.000.000.000.000.000.000.000.07
F5-35DA-F0.710.000.320.480.080.000.000.000.160.001.75
F6-35DA-F0.220.070.140.290.070.000.000.000.220.001.01
F7-35DA-F0.000.000.070.000.000.000.000.000.000.000.07
1 Conceptuses Sample Identification: Eleven female embryos samples from supplemented (E-25DA-F) and non-supplemented (E-25DC-F) gilts; seven female fetuses samples from supplemented (F-35DA-F) and non-supplemented (F-35DC-F) gilts; 2 DDX3Y: ATP-dependent RNA helicase DDX3X; KDM5D: Lysine-specific Demethylase 5D; ZFY: Zinc Finger Protein Y-Linked; EIF2S3Y: eukaryotic translation initiation factor 2 subunit 3; EIF1AY: eukaryotic translation initiation factor 1A, Y-linked; LOC110255320: Lysine-specific demethylase 6A-like; LOC110257894: Gamma-taxilin-like; LOC396706: U2 small nuclear ribonucleoprotein auxiliary factor 35 kDa subunit-related protein 2; LOC100625207: Probable ubiquitin carboxyl-terminal hydrolase FAF-X; LOC110255257: Oral-facial-digital syndrome 1 protein-like; 3 ΣCPMchrY: sum of reads CPM of 10 Y chromosome-linked genes in females conceptuses.

Share and Cite

MDPI and ACS Style

Teixeira, S.A.; Ibelli, A.M.G.; Cantão, M.E.; Oliveira, H.C.d.; Ledur, M.C.; Peixoto, J.d.O.; Marques, D.B.D.; Costa, K.A.; Coutinho, L.L.; Guimarães, S.E.F. Sex Determination Using RNA-Sequencing Analyses in Early Prenatal Pig Development. Genes 2019, 10, 1010. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10121010

AMA Style

Teixeira SA, Ibelli AMG, Cantão ME, Oliveira HCd, Ledur MC, Peixoto JdO, Marques DBD, Costa KA, Coutinho LL, Guimarães SEF. Sex Determination Using RNA-Sequencing Analyses in Early Prenatal Pig Development. Genes. 2019; 10(12):1010. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10121010

Chicago/Turabian Style

Teixeira, Susana A., Adriana M. G. Ibelli, Maurício E. Cantão, Haniel C. de Oliveira, Mônica C. Ledur, Jane de O. Peixoto, Daniele B. D. Marques, Karine A. Costa, Luiz. L. Coutinho, and Simone E. F. Guimarães. 2019. "Sex Determination Using RNA-Sequencing Analyses in Early Prenatal Pig Development" Genes 10, no. 12: 1010. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10121010

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop