Charactering the ZFAND3 gene mapped in the sex-determining locus in hybrid tilapia (Oreochromis spp.)

Scrutinizing the human TEX genes in the context of human male infertility

BACKGROUND

Infertility affects around 15% of all couples worldwide and is increasingly linked to variants in genes specifically expressed in the testis. Well-established causes of male infertility include pathogenic variants in the genes TEX11, TEX14, and TEX15, while few studies have recently reported variants in TEX13B, TEX13C, FAM9A (TEX39A), and FAM9B (TEX39B).

OBJECTIVES

We aimed at screening for novel potential candidate genes among the human TEX ("testis expressed") genes as well as verifying previously described disease associations in this set of genes.

MATERIALS AND METHODS

To this end, we screened the exome sequencing data of 1305 men, including 1056 crypto- and azoospermic individuals, and determined cell-specific expression by analyzing testis-specific single-cell RNA sequencing data for genes with identified variants. To investigate the overarching role in male fertility, we generated testis-specific knockdown (KD) models of all 10 orthologous TEX genes in Drosophila melanogaster.

RESULTS

We detected rare potential disease-causing variants in TEX10, TEX13A, TEX13B, TEX13C, TEX13D, ZFAND3 (TEX27), TEX33, FAM9A (TEX39A), and FAM9B (TEX39B), in 28 infertile men, of which 15 men carried variants in TEX10, TEX27, and TEX33. The KD of TEX2, TEX9, TEX10, TEX13, ZFAND3 (TEX27), TEX28, TEX30, NFX1 (TEX42), TEX261, and UTP4 (TEX292) in Drosophila resulted in normal fertility.

DISCUSSION

Based on our findings, the autosomal dominant predicted genes TEX10 and ZFAND3 (TEX27) and the autosomal recessive predicted gene TEX33, which all three are conceivably required for germ cell maturation, were identified as novel potential candidate genes for human non-obstructive azoospermia. We additionally identified hemizygous loss-of-function (LoF) variants in TEX13B, TEX13C, and FAM9A (TEX39A) as unlikely monogenic culprits of male infertility as LoF variants were also found in control men.

CONCLUSION

Our findings concerning the X-linked genes TEX13B, TEX13C, and FAM9A (TEX39A) contradict previous reports and will decrease false-positive reports in genetic diagnostics of azoospermic men.

Absence of Figla-like Gene Is Concordant with Femaleness in Cichlids Harboring the LG1 Sex-Determination System

Oreochromis niloticus has been used as a reference genome for studies of tilapia sex determination (SD) revealing segregating genetic loci on linkage groups (LGs) 1, 3, and 23. The master key regulator genes (MKR) underlying the SD regions on LGs 3 and 23 have been already found. To identify the MKR in fish that segregate for the LG1 XX/XY SD-system, we applied short variant discovery within the sequence reads of the genomic libraries of the Amherst hybrid stock, Coptodon zillii and Sarotherodon galilaeus, which were aligned to a 3-Mbp-region of the O. aureus genome. We obtained 66,372 variants of which six were concordant with the XX/XY model of SD and were conserved across these species, disclosing the male specific figla-like gene. We further validated this observation in O. mossambicus and in the Chitralada hybrid stock. Genome alignment of the 1252-bp transcript showed that the figla-like gene’s size was 2664 bp, and that its three exons were capable of encoding 99 amino acids including a 45-amino-acid basic helix–loop–helix domain that is typical of the ovary development regulator—factor-in-the-germline-alpha (FIGLA). In Amherst gonads, the figla-like gene was exclusively expressed in testes. Thus, the figla-like genomic presence determines male fate by interrupting the female developmental program. This indicates that the figla-like gene is the long-sought SD MKR on LG1.

Detecting Local Adaptation between North and South European Atlantic Salmon Populations

Pollution and other anthropogenic effects have driven a decrease in Atlantic salmon (Salmo salar) in the Iberian Peninsula. The restocking effort carried out in the 1980s, with salmon from northern latitudes with the aim of mitigating the decline of native populations, failed, probably due to the deficiency in adaptation of foreign salmon from northern Europe to the warm waters of the Iberian Peninsula. This result would imply that the Iberian populations of Atlantic salmon have experienced local adaptation in their past evolutionary history, as has been described for other populations of this species and other salmonids. Local adaptation can occur by divergent selections between environments, favoring the fixation of alleles that increase the fitness of a population in the environment it inhabits relative to other alleles favored in another population. In this work, we compared the genomes of different populations from the Iberian Peninsula (Atlantic and Cantabric basins) and Scotland in order to provide tentative evidence of candidate SNPs responsible for the adaptive differences between populations, which may explain the failures of restocking carried out during the 1980s. For this purpose, the samples were genotyped with a 220,000 high-density SNP array (Affymetrix) specific to Atlantic salmon. Our results revealed potential evidence of local adaptation for North Spanish and Scottish populations. As expected, most differences concerned the comparison of the Iberian Peninsula with Scotland, although there were also differences between Atlantic and Cantabric populations. A high proportion of the genes identified are related to development and cellular metabolism, DNA transcription and anatomical structure. A particular SNP was identified within the NADP-dependent malic enzyme-2 (_mMEP-2*), previously reported by independent studies as a candidate for local adaptation in salmon from the Iberian Peninsula. Interestingly, the corresponding SNP within the _mMEP-2* region was consistent with a genomic pattern of divergent selection.

Genome-Wide Association Study and Pathway Analysis for Female Fertility Traits in Iranian Holstein Cattle

Female fertility is an important trait that contributes to cow’s profitability and it can be improved by genomic information. The objective of this study was to detect genomic regions and variants affecting fertility traits in Iranian Holstein cattle. A data set comprised of female fertility records and 3,452,730 pedigree information from Iranian Holstein cattle were used to predict the breeding values, which were then employed to estimate the de-regressed proofs (DRP) of genotyped animals. A total of 878 animals with DRP records and 54k SNP markers were utilized in the genome-wide association study (GWAS). The GWAS was performed using a linear regression model with SNP genotype as a linear covariate. The results showed that an SNP on BTA19, ARS-BFGL-NGS-33473, was the most significant SNP associated with days from calving to first service. In total, [69] significant SNPs were located within 27 candidate genes. Novel potential candidate genes include OSTN, DPP6, EphA5, CADPS2, Rfc1, ADGRB3, Myo3a, C10H14orf93, KIAA1217, RBPJL, SLC18A2, GARNL3, NCALD, ASPH, ASIC2, OR3A1, CHRNB4, CACNA2D2, DLGAP1, GRIN2A and ME3. These genes are involved in different pathways relevant to female fertility and other characteristics in mammals. Gene set enrichment analysis showed that thirteen GO terms had significant overrepresentation of genes statistically associated with female fertility traits. The results of network analysis identified CCNB1 gene as a hub gene in the progesterone-mediated oocyte maturation pathway, significantly associated with age at first calving. The candidate genes identified in this study can be utilized in genomic tests to improve reproductive performance in Holstein cattle.

Fine Mapping Using Whole-Genome Sequencing Confirms Anti-Müllerian Hormone as a Major Gene for Sex Determination in Farmed Nile Tilapia (Oreochromis niloticus L.)

Nile tilapia (Oreochromis niloticus) is one of the most cultivated and economically important species in world aquaculture. Intensive production promotes the use of monosex animals, due to an important dimorphism that favors male growth. Currently, the main mechanism to obtain all-male populations is the use of hormones in feeding during larval and fry phases. Identifying genomic regions associated with sex determination in Nile tilapia is a research topic of great interest. The objective of this study was to identify genomic variants associated with sex determination in three commercial populations of Nile tilapia. Whole-genome sequencing of 326 individuals was performed, and a total of 2.4 million high-quality bi-allelic single nucleotide polymorphisms (SNPs) were identified after quality control. A genome-wide association study (GWAS) was conducted to identify markers associated with the binary sex trait (males = 1; females = 0). A mixed logistic regression GWAS model was fitted and a genome-wide significant signal comprising 36 SNPs, spanning a genomic region of 536 kb in chromosome 23 was identified. Ten out of these 36 genetic variants intercept the anti-Müllerian (Amh) hormone gene. Other significant SNPs were located in the neighboring Amh gene region. This gene has been strongly associated with sex determination in several vertebrate species, playing an essential role in the differentiation of male and female reproductive tissue in early stages of development. This finding provides useful information to better understand the genetic mechanisms underlying sex determination in Nile tilapia.

Characterization of a male specific region containing a candidate sex determining gene in Atlantic cod

The genetic mechanisms determining sex in teleost fishes are highly variable and the master sex determining gene has only been identified in few species. Here we characterize a male-specific region of 9 kb on linkage group 11 in Atlantic cod (Gadus morhua) harboring a single gene named zkY for zinc knuckle on the Y chromosome. Diagnostic PCR test of phenotypically sexed males and females confirm the sex-specific nature of the Y-sequence. We identified twelve highly similar autosomal gene copies of zkY, of which eight code for proteins containing the zinc knuckle motif. 3D modeling suggests that the amino acid changes observed in six copies might influence the putative RNA-binding specificity. Cod zkY and the autosomal proteins zk1 and zk2 possess an identical zinc knuckle structure, but only the Y-specific gene zkY was expressed at high levels in the developing larvae before the onset of sex differentiation. Collectively these data suggest zkY as a candidate master masculinization gene in Atlantic cod. PCR amplification of Y-sequences in Arctic cod (Arctogadus glacialis) and Greenland cod (Gadus macrocephalus ogac) suggests that the male-specific region emerged in codfishes more than 7.5 million years ago.

SOX2 participates in spermatogenesis of Zhikong scallop Chlamys farreri

As an important transcription factor, SOX2 involves in embryogenesis, maintenance of stem cells and proliferation of primordial germ cell (PGC). However, little was known about its function in mature gonads. Herein, we investigated the SOX2 gene profiles in testis of scallop, Chlamys farreri. The level of C. farreri SOX2 (Cf-SOX2) mRNA increased gradually along with gonadal development and reached the peak at mature stage, and was located in all germ cells, including spermatogonia, spermatocytes, spermatids and spermatozoa. Knockdown of Cf-SOX2 using RNAi leaded to a mass of germ cells lost, and only a few spermatogonia retained in the nearly empty testicular acini after 21 days. TUNEL assay showed that apoptosis occurred in spermatocytes. Furthermore, transcriptome profiles of the testes were compared between Cf-SOX2 knockdown and normal scallops, 131,340 unigenes were obtained and 2,067 differential expression genes (DEGs) were identified. GO and KEGG analysis showed that most DEGs were related to cell apoptosis (casp2, casp3, casp8), cell proliferation (samd9, crebzf, iqsec1) and spermatogenesis (htt, tusc3, zmynd10, nipbl, mfge8), and enriched in p53, TNF and apoptosis pathways. Our study revealed Cf-SOX2 is essential in spermatogenesis and testis development of C. farreri and provided important clues for better understanding of SOX2 regulatory mechanisms in bivalve testis.

ddRADseq reveals determinants for temperature-dependent sex reversal in Nile tilapia on LG23

BACKGROUND

In Nile tilapia sex determination is governed by a male heterogametic system XX/XY either on LG1 or LG23. The latter carries a Y-specific duplicate of the amh gene, which is a testis-determining factor. Allelic variants in the amh gene demonstrated to be major triggers for autosomal and temperature-dependent sex reversal. Further, QTL on LG23 and LG20 show a temperature-responsiveness with influence on the phenotypic sex relative to the sex chromosomes. Here we present a ddRADseq based approach to identify genomic regions that show unusual large differentiation in terms of fixation index (F_ST) between temperature-treated pseudomales and non-masculinized females using a comparative genome-scan. Genome-wide associations were identified for the temperature-dependent sex using a genetically all-female population devoid of amh-ΔY.

RESULTS

Twenty-two thousand three hundred ninety-two SNPs were interrogated for the comparison of temperature-treated pseudomales and females, which revealed the largest differentiation on LG23. Outlier F_ST-values (0.35-0.44) were determined for six SNPs in the genomic interval (9,190,077-11,065,693) harbouring the amh gene (9,602,693-9,605,808), exceeding the genome-wide low F_ST of 0.013. Association analysis with a set of 9104 selected SNPs confirmed that the same genomic region on LG23 exerts a significant effect on the temperature-dependent sex.

CONCLUSIONS

This study highlights the role of LG23 in sex determination, harbouring major determinants for temperature-dependent sex reversal in Nile tilapia. Furthermore F_ST outlier detection proves a powerful tool for detection of sex-determining regions in fish genomes.