Copy number variation of horse Y chromosome genes in normal equine populations and in horses with abnormal sex development and subfertility: relationship of copy number variations with Y haplogroups.

Structural rearrangements like copy number variations in the male-specific Y chromosome have been associated with male fertility phenotypes in human and mouse but have been sparsely studied in other mammalian species. Here, we designed digital droplet PCR assays for 7 horse male-specific Y chromosome multicopy genes and SRY and evaluated their absolute copy numbers in 209 normal male horses of 22 breeds, 73 XY horses with disorders of sex development and/or infertility, 5 Przewalski's horses and 2 kulans. This established baseline copy number for these genes in horses. The TSPY gene showed the highest copy number and was the most copy number variable between individuals and breeds. SRY was a single-copy gene in most horses but had 2-3 copies in some indigenous breeds. Since SRY is flanked by 2 copies of RBMY, their copy number variations were interrelated and may lead to SRY-negative XY disorders of sex development. The Przewalski's horse and kulan had 1 copy of SRY and RBMY. TSPY and ETSTY2 showed significant copy number variations between cryptorchid and normal males (P < 0.05). No significant copy number variations were observed in subfertile/infertile males. Notably, copy number of TSPY and ETSTY5 differed between successive male generations and between cloned horses, indicating germline and somatic mechanisms for copy number variations. We observed no correlation between male-specific Y chromosome gene copy number variations and male-specific Y chromosome haplotypes. We conclude that the ampliconic male-specific Y chromosome reference assembly has deficiencies and further studies with an improved male-specific Y chromosome assembly are needed to determine selective constraints over horse male-specific Y chromosome gene copy number and their relation to stallion reproduction and male biology.

A transposon expression burst accompanies the activation of Y-chromosome fertility genes during Drosophila spermatogenesis.

Transposable elements (TEs) must replicate in germline cells to pass novel insertions to offspring. In Drosophila melanogaster ovaries, TEs can exploit specific developmental windows of opportunity to evade host silencing and increase their copy numbers. However, TE activity and host silencing in the distinct cell types of Drosophila testis are not well understood. Here, we reanalyze publicly available single-cell RNA-seq datasets to quantify TE expression in the distinct cell types of the Drosophila testis. We develop a method for identification of TE and host gene expression modules and find that a distinct population of early spermatocytes expresses a large number of TEs at much higher levels than other germline and somatic components of the testes. This burst of TE expression coincides with the activation of Y chromosome fertility factors and spermatocyte-specific transcriptional regulators, as well as downregulation of many components of the piRNA pathway. The TEs expressed by this cell population are specifically enriched on the Y chromosome and depleted on the X chromosome, relative to other active TEs. These data suggest that some TEs may achieve high insertional activity in males by exploiting a window of opportunity for mobilization created by the activation of spermatocyte-specific and Y chromosome-specific transcriptional programs.

Dynamic sex chromosome expression in Drosophila male germ cells.

Given their copy number differences and unique modes of inheritance, the evolved gene content and expression of sex chromosomes is unusual. In many organisms the X and Y chromosomes are inactivated in spermatocytes, possibly as a defense mechanism against insertions into unpaired chromatin. In addition to current sex chromosomes, Drosophila has a small gene-poor X-chromosome relic (4th) that re-acquired autosomal status. Here we use single cell RNA-Seq on fly larvae to demonstrate that the single X and pair of 4th chromosomes are specifically inactivated in primary spermatocytes, based on measuring all genes or a set of broadly expressed genes in testis we identified. In contrast, genes on the single Y chromosome become maximally active in primary spermatocytes. Reduced X transcript levels are due to failed activation of RNA-Polymerase-II by phosphorylation of Serine 2 and 5.

Bayesian Inference in Y-Linked Two-Sex Branching Processes with Mutations: ABC Approach.

A Y-linked two-sex branching process with mutations and blind choice of males is a suitable model for analyzing the evolution of the number of carriers of a Y-linked allele and its mutations. Such a model considers a two-sex monogamous population in which each female chooses her partner from among the male population without caring about his type (i.e., the allele he carries). In this work, we deal with the problem of estimating the main parameters of these models by developing Bayesian inference in a parametric framework. First, we consider as a sample scheme the observation of the total number of females and males up to some generation as well as the number of males of each genotype in the last generation. Subsequently, we introduce the information on the mutated males in only the last generation, obtaining in this way a second sample scheme. For both samples, we apply the Approximate Bayesian Computation (ABC) method to approximate the posterior distributions of the main parameters of the model. The accuracy of the procedure based on these samples is illustrated and discussed by way of simulated examples.

X- and Y-Linked Chromatin-Modifying Genes as Regulators of Sex-Specific Cancer Incidence and Prognosis.

Biological sex profoundly conditions organismal development and physiology, imposing wide-ranging effects on cell signaling, metabolism, and immune response. These effects arise from sex-specified differences in hormonal exposure, and from intrinsic genetic and epigenetic differences associated with the presence of an XX versus XY chromosomal complement. In addition, biological sex is now recognized to be a determinant of the incidence, presentation, and therapeutic response of multiple forms of cancer, including cancers not specifically associated with male or female anatomy. Although multiple factors contribute to sex-based differences in cancer, a growing body of research emphasizes a role for differential activity of X- and Y-linked tumor-suppressor genes in males and females. Among these, the X-linked KDM6A/UTX and KDM5C/JARID1C/SMCX, and their Y-linked paralogs UTY/KDM6C and KDM5D/JARID1D/SMCY encode lysine demethylases. These epigenetic modulators profoundly influence gene expression, based on enzymatic activity in demethylating H3K27me3 and H3K4me3, and nonenzymatic scaffolding roles for large complexes that open and close chromatin for transcription. In a growing number of cases, mutations affecting these proteins have been recognized to strongly influence cancer risk, prognosis, and response to specific therapies. However, sex-specific patterns of mutation, expression, and activity of these genes, coupled with tissue-specific requirement for their function as tumor suppressors, together exemplify the complex relationship between sex and cancer vulnerabilities. In this review, we summarize and discuss the current state of the literature on the roles of these proteins in contributing to sex bias in cancer, and the status of clinical agents relevant to their function.

Sex-Associated Gene Expression Alterations Correlate With Esophageal Cancer Survival.

OBJECTIVES: Esophageal cancer (EC) is a significant cause of cancer death with 5-year survival of 10%-15% and males more frequently affected. Genetic evaluation for loci highlighting risk has been performed, but survival data are limited. The Cancer Genome Atlas (TCGA) data sets allow for potential prognostic marker assessment in large patient cohorts. The study aimed to use the TCGA EC data set to assess whether survival varies by sex and explore genetic alterations that may explain variation observed. METHODS: TCGA clinical/RNA-seq data sets (n = 185, 158 males/27 females) were downloaded from the cancer genome browser. Data analysis/figure preparation was performed in R and GraphPad Prism 7. Survival analysis was performed using the survival package. Text mining of PubMed was performed using the tm, RISmed, and wordcloud packages. Pathway analysis was performed using the Reactome database. RESULTS: In EC, male sex/high tumor grade reduced overall survival (hazard ratio = 2.27 [0.99-5.24] for M vs F and 2.49 [0.89-6.92] for low vs high grade, respectively) and recurrence-free survival (hazard ratio = 4.09 [0.98-17.03] for M vs F and 3.36 [0.81-14.01] for low vs high grade, respectively). To investigate the genetic basis for sex-based survival differences in EC, corresponding gene expression data were analyzed. Sixty-nine genes were dysregulated at the P < 0.01 level by the Wilcox test, 33% were X-chromosome genes, and 7% were Y-chromosome genes. DISCUSSION: Female sex potentially confers an EC survival advantage. Importantly, we demonstrate a genetic/epigenetic basis for these survival differences that are independent of lifestyle-associated risk factors overrepresented in males. Further research may lead to novel concepts in treating/measuring EC aggressiveness by sex.

Dosage regulation, and variation in gene expression and copy number of human Y chromosome ampliconic genes.

The Y chromosome harbors nine multi-copy ampliconic gene families expressed exclusively in testis. The gene copies within each family are >99% identical to each other, which poses a major challenge in evaluating their copy number. Recent studies demonstrated high variation in Y ampliconic gene copy number among humans. However, how this variation affects expression levels in human testis remains understudied. Here we developed a novel computational tool Ampliconic Copy Number Estimator (AmpliCoNE) that utilizes read sequencing depth information to estimate Y ampliconic gene copy number per family. We applied this tool to whole-genome sequencing data of 149 men with matched testis expression data whose samples are part of the Genotype-Tissue Expression (GTEx) project. We found that the Y ampliconic gene families with low copy number in humans were deleted or pseudogenized in non-human great apes, suggesting relaxation of functional constraints. Among the Y ampliconic gene families, higher copy number leads to higher expression. Within the Y ampliconic gene families, copy number does not influence gene expression, rather a high tolerance for variation in gene expression was observed in testis of presumably healthy men. No differences in gene expression levels were found among major Y haplogroups. Age positively correlated with expression levels of the HSFY and PRY gene families in the African subhaplogroup E1b, but not in the European subhaplogroups R1b and I1. We also found that expression of five Y ampliconic gene families is coordinated with that of their non-Y (i.e. X or autosomal) homologs. Indeed, five ampliconic gene families had consistently lower expression levels when compared to their non-Y homologs suggesting dosage regulation, while the HSFY family had higher expression levels than its X homolog and thus lacked dosage regulation.

Copy number variation of ZNF280AY across 21 cattle breeds and its association with the reproductive traits of Holstein and Simmental bulls.

The mammalian Y chromosome gene families in the ampliconic region are expressed predominantly or exclusively in the testis, and their copy number variations (CNV) are significantly associated with male reproductive traits, suggesting they have important roles in spermatogenesis and testicular development. ZNF280AY (zinc finger protein 280A, Y-linked) is a member of the zinc finger protein family and has been identified as a bovid-specific Y-chromosome gene. The current study applied a reliable quantitative real-time PCR method to estimate the CNV of ZNF280AY in 715 bulls across 21 cattle breeds and to further investigate the association of the CNV of ZNF280AY with bull reproductive traits and ZNF280AY mRNA expression levels in adult testis. The results revealed that the median copy number of ZNF280AY was 47, and the copy number varied from 11 to 154, showing significant CNV between and within the investigated cattle breeds. In addition, all 715 bulls were classified into Y1, Y2, and Y3 lineage groups based on a rapid genotyping method described previously. Pairwise comparisons indicated that bulls belonging to the Y1 lineage had a significantly lower median copy number (40) than bulls belonging to the Y2 (52) and Y3 lineages (57). Association analysis revealed that the CNV of ZNF280AY was correlated negatively with the percentage of normal sperm and sperm concentration in Holstein bulls, whereas no significant correlation was observed with ejaculation volume, total sperm count, sperm motility, postthaw motility (PTM), and scrotal circumference in Holstein and Simmental bulls. Furthermore, no correlation was observed between ZNF280AY copy number and ZNF280AY mRNA expression levels in the testis. The current study suggests that the CNV of the ZNF280AY gene family is associated with male reproductive traits and may serve as a valuable marker for early bull fertility selection in Holstein breeding programs.

Y chromosome gene copy number and lack of autism phenotype in a male with an isodicentric Y chromosome and absent NLGN4Y expression.

We describe a unique male with a dicentric Y chromosome whose phenotype was compared to that of males with 47,XYY (XYY). The male Y-chromosome aneuploidy XYY is associated with physical, behavioral/cognitive phenotypes, and autism spectrum disorders. We hypothesize that increased risk for these phenotypes is caused by increased copy number/overexpression of Y-encoded genes. Specifically, an extra copy of the neuroligin gene NLGN4Y might elevate the risk of autism in boys with XYY. We present a unique male with the karyotype 46,X,idic(Y)(q11.22), which includes duplication of the Y short arm and proximal long arm and deletion of the distal long arm, evaluated his physical, behavioral/cognitive, and neuroimaging/magnetoencephalography (MEG) phenotypes, and measured blood RNA expression of Y genes. The proband had tall stature and cognitive function within the typical range, without autism features. His blood RNA showed twofold increase in expression of Yp genes versus XY controls, and absent expression of deleted Yq genes, including NLGN4Y. The M100 latencies were similar to findings in typically developing males. In summary, the proband had overexpression of a subset of Yp genes, absent NLGN4Y expression, without ASD findings or XYY-MEG latency findings. These results are consistent with a role for NLGN4Y overexpression in the etiology of behavioral phenotypes associated with XYY. Further investigation of NLGN4Y as an ASD risk gene in XYY is warranted. The genotype and phenotype(s) of this subject may also provide insight into how Y chromosome genes contribute to normal male development and the male predominance in ASD.

Low prevalence of male microchimerism in women with Mayer-Rokitansky-Küster-Hauser syndrome.

STUDY QUESTION: Is there an increased prevalence of male microchimerism in women with Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome, as evidence of fetal exposure to blood and anti-Müllerian hormone (AMH) from a (vanished) male co-twin resulting in regression of the Müllerian duct derivatives? SUMMARY ANSWER: Predominant absence of male microchimerism in adult women with MRKH syndrome does not support our hypothesis that intrauterine blood exchange with a (vanished) male co-twin is the pathophysiological mechanism. WHAT IS KNOWN ALREADY: The etiology of MRKH is unclear. Research on the phenotype analogous condition in cattle (freemartinism) has yielded the hypothesis that Müllerian duct development is inhibited by exposure to AMH in utero. In cattle, the male co-twin has been identified as the source for AMH, which is transferred via placental blood exchange. In human twins, a similar exchange of cellular material has been documented by detection of chimerism, but it is unknown whether this has clinical consequences. STUDY DESIGN, SIZE, DURATION: An observational case-control study was performed to compare the presence of male microchimerism in women with MRKH syndrome and control women. Through recruitment via the Dutch patients' association of women with MRKH (comprising 300 members who were informed by email or regular mail), we enrolled 96 patients between January 2017 and July 2017. The control group consisted of 100 women who reported never having been pregnant. PARTICIPANTS/MATERIALS, SETTING, METHODS: After written informed consent, peripheral blood samples were obtained by venipuncture, and genomic DNA was extracted. Male microchimerism was detected by Y-chromosome-specific real-time quantitative PCR, with use of DYS14 marker. Possible other sources for microchimerism, for example older brothers, were evaluated using questionnaire data. MAIN RESULTS AND THE ROLE OF CHANCE: The final analysis included 194 women: 95 women with MRKH syndrome with a mean age of 40.9 years and 99 control women with a mean age of 30.2 years. In total, 54 women (56.8%) were identified as having typical MRKH syndrome, and 41 women (43.2%) were identified as having atypical MRKH syndrome (when extra-genital malformations were present). The prevalence of male microchimerism was significantly higher in the control group than in the MRKH group (17.2% versus 5.3%, P = 0.009). After correcting for age, women in the control group were 5.8 times more likely to have male microchimerism (odds ratio 5.84 (CI 1.59-21.47), P = 0.008). The mean concentration of male microchimerism in the positive samples was 56.0 male genome equivalent per 1 000 000 cells. The prevalence of male microchimerism was similar in women with typical MRKH syndrome and atypical MRKH syndrome (5.6% versus 4.9%, P = 0.884). There were no differences between women with or without microchimerism in occurrence of alternative sources of XY cells, such as older brothers, previous blood transfusion, or history of sexual intercourse. LIMITATIONS, REASON FOR CAUTION: We are not able to draw definitive conclusions regarding the occurrence of AMH exchange during embryologic development in women with MRKH syndrome. Our subject population includes all adult women and therefore is reliant on long-term prevalence of microchimerism. Moreover, we have only tested blood, and, theoretically, the cells may have grafted anywhere in the body during development. It must also be considered that the exchange of AMH may occur without the transfusion of XY cells and therefore cannot be discovered by chimerism detection. WIDER IMPLICATIONS OF THE FINDINGS: This is the first study to test the theory that freemartinism causes the MRKH syndrome in humans. The study aimed to test the presence of male microchimerism in women with MRKH syndrome as a reflection of early fetal exposure to blood and AMH from a male (vanished) co-twin. We found that male microchimerism was only present in 5.3% of the women with MRKH syndrome, a significantly lower percentage than in the control group (17.2%). Our results do not provide evidence for an increased male microchimerism in adult women with MRKH as a product of intrauterine blood exchange. However, the significant difference in favor of the control group is of interest to the ongoing discussion on microchimeric cell transfer and the possible sources of XY cells. STUDY FUNDING/COMPETING INTEREST(S): None. TRIAL REGISTRATION NUMBER: Dutch trial register, NTR5961.

Sex-specific polymorphism of MET1 and ARR17 genes in Populus × sibirica.

The genus Populus is an effective model in tree genetics. This genus includes dioecious species and, recently, whole genome resequencing of P. trichocarpa and P. balsamifera enabled the identification of sex-linked regions and sex-associated single nucleotide polymorphisms (SNPs). These results created new opportunities to study sex determination in poplars. In the present work, we performed deep sequencing of genes encoding METHYLTRANSFERASE1 (MET1) and homolog of ARABIDOPSIS RESPONSE REGULATOR 17 (ARR17), which are localized in a sex-linked region of Populus genome and contain a number of sex-associated SNPs. Amplicon libraries for 38 samples of P. × sibirica (19 males and 19 females) were sequenced on MiSeq Illumina (300 nt paired-end reads) and approximately 4000× coverage was obtained for each sample. In total, from 80 to 179 SNPs were detected in poplar individuals for MET1, and from 16 to 49 SNPs were detected for ARR17. We identified 17 sex-specific SNPs (11 in MET1 and 6 in ARR17) - they were present in all males but absent in all females. For identified sex-specific SNP sites, females were homozygous, while males were heterozygous. Moreover, colocation of sex-specific SNPs confirming the XY sex-determination system of poplars was revealed: in one allelic variant, males had the same nucleotides as females, while in the other, sex-specific SNPs were present. Based on the data obtained, we developed and successfully applied a high-resolution melting-based approach for sex identification in poplars. The developed molecular markers are useful for distinguishing between male and female poplars in scientific research and can also be applied to select male-only genotypes for use in city landscaping and production of paper, pulp, and biofuel.

Rethinking sex determination of non-gonadal tissues.

Evolution of genetic mechanisms of sex determination led to two processes causing sex differences in somatic phenotypes: gonadal differentiation and sex chromosome dosage inequality. In species with heteromorphic sex chromosomes, the sex of the individual is established at the time of formation of the zygote, leading to inherent sex differences in expression of sex chromosome genes beginning as soon as the embryonic transcriptome is activated. The inequality of sex chromosome gene expression causes sexual differentiation of the gonads and of non-gonadal tissues. The difference in gonad type in turn causes lifelong differences in gonadal hormones, which interact with unequal effects of X and Y genes acting within cells. Separating the effects of gonadal hormones and sex chromosomes has been possible using mouse models in which gonadal determination is separated from the sex chromosomes, allowing comparison of XX and XY mice with the same type of gonad. Sex differences caused by gonadal hormones and sex chromosomes affect basic physiology and disease mechanisms in most or all tissues.

Genetic analysis of 29 Y-STR loci in Han population from Dongfang, Southern China.

In the present study, blood samples of 984 unrelated Han individuals were collected from Dongfang, Southern China, after informed consent. A total of 29 Y-chromosomal short tandem repeat (Y-STR) were analyzed, including DYF387S1, DYS19, DYS385ab, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS444, DYS447, DYS448, DYS449, DYS456, DYS458, DYS460, DYS481, DYS508, DYS518, DYS533, DYS576, DYS635, DYS643 and GATAH4. A total of 749 different haplotypes were found among 984 individuals, of which 645 were unique. The haplotype diversity was 0.9988 and the discrimination capacity was 0.7612, while the match probability was 0.0025. The smallest genetic distance (RST = 0.0155) was found between the Dongfang Han population and Guizhou Han population, while the largest genetic distance (RST = 0.1284) was observed with Gansu Tibetan.

Codon usage bias and its influencing factors for Y-linked genes in human.

The non-uniform usage of synonymous codons during translation of a protein is the codon usage bias and is mainly influenced by natural selection and mutation pressure. We have used bioinformatic approaches to analyze codon usage bias of human Y-linked genes. Effective number of codon (ENC) suggested that the overall extent of codon usage bias of genes was low. The relative synonymous codon usage (RSCU) analysis revealed that AGA and CTG codons were over-represented in Y-linked genes. Compositional constraint under mutation pressure influenced the codon usage pattern as revealed by the correspondence analysis (COA). Parity plot suggests that both natural selection and mutation pressure might have influenced the codon usage bias of Y-linked genes.

DNA analysis of root canal-filled teeth.

Teeth are markedly useful as samples for DNA analysis; however, intact teeth are not always available. This study examined the possibility of identifying autosomal and Y-chromosome short tandem repeat (STR) types in samples from 34 teeth (15 intact and 19 root canal filled) that had been preserved for 10-33years after dental extraction. The aim was to explore the feasibility of individual identification by DNA analysis of samples obtained from highly decomposed and skeletonized corpses. Only one out of 24 autosomal STR loci was not identified in two of the 15 intact teeth, whereas all 23 loci of the Y chromosome STR were detected. One or two autosomal STR loci remained unidentified in eight of the 19 root-filled teeth, and four or five of the 23 Y STR loci were undetected in three cases. However, the types were identified in about 20 loci in all samples, and the composition of the root canal filling material did not appear to interfere with the PCR. This study demonstrates that the storage period of the teeth had no influence on our results indicating that root canal filled teeth can be used for DNA analysis.

Role of testosterone and Y chromosome genes for the masculinization of the human brain.

Women with complete androgen insensitivity syndrome (CAIS) have a male (46,XY) karyotype but no functional androgen receptors. Their condition, therefore, offers a unique model for studying testosterone effects on cerebral sex dimorphism. We present MRI data from 16 women with CAIS and 32 male (46,XY) and 32 female (46,XX) controls. METHODS: FreeSurfer software was employed to measure cortical thickness and subcortical structural volumes. Axonal connections, indexed by fractional anisotropy, (FA) were measured with diffusion tensor imaging, and functional connectivity with resting state fMRI. RESULTS: Compared to men, CAIS women displayed a "female" pattern by having thicker parietal and occipital cortices, lower FA values in the right corticospinal, superior and inferior longitudinal tracts, and corpus callosum. Their functional connectivity from the amygdala to the medial prefrontal cortex, was stronger and amygdala-connections to the motor cortex weaker than in control men. CAIS and control women also showed stronger posterior cingulate and precuneus connections in the default mode network. Thickness of the motor cortex, the caudate volume, and the FA in the callosal body followed, however, a "male" pattern. CONCLUSION: Altogether, these data suggest that testosterone modulates the microstructure of somatosensory and visual cortices and their axonal connections to the frontal cortex. Testosterone also influenced functional connections from the amygdala, whereas the motor cortex could, in agreement with our previous reports, be moderated by processes linked to X-chromosome gene dosage. These data raise the question about other genetic factors masculinizing the human brain than the SRY gene and testosterone. Hum Brain Mapp 38:1801-1814, 2017. © 2017 Wiley Periodicals, Inc.

Extremely low nucleotide diversity in the X-linked region of papaya caused by a strong selective sweep.

BACKGROUND: The papaya Y-linked region showed clear population structure, resulting in the detection of the ancestral male population that domesticated hermaphrodite papayas were selected from. The same populations were used to study nucleotide diversity and population structure in the X-linked region. RESULTS: Diversity is very low for all genes in the X-linked region in the wild dioecious population, with nucleotide diversity π syn = 0.00017, tenfold lower than the autosomal region (π syn = 0.0017) and 12-fold lower than the Y-linked region (π syn = 0.0021). Analysis of the X-linked sequences shows an undivided population, suggesting a geographically wide diversity-reducing event, whereas two subpopulations were observed in the autosomes separating gynodioecy and dioecy and three subpopulations in the Y-linked region separating three male populations. The extremely low diversity in the papaya X-linked region was probably caused by a recent, strong selective sweep before domestication, involving either the spread of a recessive mutation in an X-linked gene that is beneficial to males or a partially dominant mutation that benefitted females or both sexes. Nucleotide diversity in the domesticated X samples is about half that in the wild Xs, probably due to the bottleneck when hermaphrodites were selected during domestication. CONCLUSIONS: The extreme low nucleotide diversity in the papaya X-linked region is much greater than observed in humans, great apes, and the neo-X chromosome of Drosophila miranda, which show the expected pattern of Y-linked genes < X-linked genes < autosomal genes; papaya shows an unprecedented pattern of X-linked genes < autosomal genes < Y-linked genes.

Evaluation of Y chromosomal SNP haplogrouping in the HID-Ion AmpliSeq™ Identity Panel.

The Y chromosomal haplogroup determined from single nucleotide polymorphism (SNP) combinations is a valuable genetic marker to study ancestral male lineage and ethical distribution. Next-generation sequencing has been developed for widely diverse genetics fields. For this study, we demonstrate 34 Y-SNP typing employing the Ion PGM™ system to perform haplogrouping. DNA libraries were constructed using the HID-Ion AmpliSeq™ Identity Panel. Emulsion PCR was performed, then DNA sequences were analyzed on the Ion 314 and 316 Chip Kit v2. Some difficulties became apparent during the analytic processes. No-call was reported at rs2032599 and M479 in six samples, in which the least coverage was observed at M479. A minor misreading occurred at rs2032631 and M479. A real time PCR experiment using other pairs of oligonucleotide primers showed that these events might result from the flanking sequence. Finally, Y haplogroup was determined completely for 81 unrelated males including Japanese (n=59) and Malay (n=22) subjects. The allelic divergence differed between the two populations. In comparison with the conventional Sanger method, next-generation sequencing provides a comprehensive SNP analysis with convenient procedures, but further system improvement is necessary.

Population genetic study for 24 STR loci and Y indel (GlobalFiler™ PCR Amplification kit and PowerPlex® Fusion system) in 1000 Korean individuals.

Allele frequencies for 23 autosomal short tandem repeat loci (D3S1358, vWA, D16S539, CSF1PO, TPOX, D8S1179, D21S11, D18S51, TH01, FGA, D5S818, D13S317, D7S820, D2S441, D19S433, D22S1045, D10S1248, D1S1656, D12S391, D2S1338, SE33, Penta D, Penta E), 1 Y-chromosome short tandem repeat locus (DYS391) and Y indel were obtained from 1000 unrelated individuals of the Korean population.

Ancestral Y-linked genes were maintained by translocation to the X and Y chromosomes fused to an autosomal pair in the Okinawa spiny rat Tokudaia muenninki.

Two species of the genus Tokudaia lack the Y chromosome and SRY, but several Y-linked genes have been rescued by translocation or transposition to other chromosomes. Tokudaia muenninki is the only species in the genus that maintains the Y owing to sex chromosome-autosome fusions. According to previous studies, many SRY pseudocopies and other Y-linked genes have evolved by excess duplication in this species. Using RNA-seq and RT-PCR, we found that ZFY, EIF2S3Y, TSPY, UTY, DDX3Y, USP9Y, and RBMY, but not UBA1Y, had high deduced amino acid sequence similarity and similar expression patterns with other rodents, suggesting that these genes were functional. Based on FISH and quantitative real-time PCR, all of the genes except for UTY and DDX3Y were amplified on the X and Y chromosomes with approximately 10-66 copies in the male genome. In a comparative analysis of the 372.4-kb BAC sequence and Y-linked gene transcripts from T. muenninki with the mouse Y genomic sequence, we observed that multiple-copy genes in the ancestral Y genome were nonfunctional, indicating that the gene functions were assumed by amplified copies. We also found a LTR sequence at the distal end of a SRY duplication unit, suggesting that unequal sister chromatid exchange mediated by retrotransposable elements could have been involved in SRY amplification. Our results revealed that the Y-linked genes were rescued from degeneration via translocations to other sex chromosomal regions and amplification events in T. muenninki.