Y-chromosome Sequences Research Articles

Y Chromosome Loss and Implications for Oncology.

The Y chromosome has recognized functions in promoting male sex determination and regulating aspects of fertility. However, recent work has demonstrated important roles for the Y chromosome and Y-encoded genes in multiple domains of male health, including cancer. It is well established that males experience shorter lifespans than females, and this sex bias on overall mortality is accentuated in populations with longer life expectancy, in part related to elevated rates of cancer. The majority of human malignancies exhibit a sex bias with elevated frequencies in males. For many of these cancer types, the disparity has not been explained by environmental risk factors such as tobacco use. Notably, loss of the Y chromosome (LOY) detected in blood cells, termed mosaic LOY, is a common event that is related to advancing age and is associated with a shortened lifespan. Mosaic LOY is linked to increased incidence and mortality across a range of malignancies. Furthermore, tumors arising in different anatomic sites exhibit different frequencies of partial or complete Y chromosome loss. Causal oncogenic or tumor-suppressive roles have been documented for several Y-encoded genes, such as lysine-specific demethylase 5 D, that exert pleiotropic effects on cellular functions by virtue of genome-wide regulation of gene activity. In this review, we discuss aspects of the Y chromosome relevant to oncology. The recent completion of the entire human Y-chromosome sequence provides a reference map of Y-encoded genes and regulatory elements to enable causal molecular studies that may explain and exploit the marked disparity in male cancer risk and mortality.

Paternal genetic structure analysis of the modern Han populations based on Y-SNP and Y-STR.

The Han populations represent the largest ethnic group in China. Previous studies have primarily focused on investigating their genetic origins, migration and integration, as well as paternal genetic relationships within specific regional Han populations. However, a comprehensive analysis of the global paternal genetic structure of Han populations is lacking. In this study, we performed Y-chromosome sequencing on 362 unrelated male samples from Chinese Han individuals collected from Qinghai, Sichuan and Liaoning provinces. We then integrated relevant data from reported studies. Our final dataset comprised 1830 samples from 16 Han populations across 15 provinces in China, encompassing information on 89 Y-SNPs and 16 Y-STRs. Statistical analyses were conducted to assess Y-STR haplotype diversity (HD) and Y-SNP haplogroup frequencies. Additionally, we employed principal component analysis (PCA), phylogenetic tree and haplotype network to explore genetic differentiation within Han populations and the genetic relationships between Han populations and ethnic minorities surrounding them. Our results demonstrated that the O-M175 haplogroup represents the predominant paternal lineage in Han populations, with frequencies ranging from 60.53% (Qinghai Han) to 92.7% (Guangdong Han). Moreover, the subclades downstream of O-M175 showed distinct regional variations in their distribution patterns. The O2-M122 haplogroup was prevalent in all Han populations and demonstrated a gradual decline in frequency from north to south. Conversely, the distribution frequency of the O1b-M268 haplogroup decreased from south to north, particularly showed significant presence among Han populations in the Lingnan region. Haplogroup O1a-M119 distributed more frequently in the central Han populations. Our findings revealed that Chinese Han populations can be categorized into three subgroups: northern, central, and southern. Notably, there were significant differences among Han in Qinghai and other regions. Regarding the genetic relationships between Han populations and surrounding ethnic minorities, we observed a closer genetic affinity between different Han populations, but northern Han demonstrated a stronger relationship with the Hui ethnic group, while southern Han exhibited a closer connection with the Gelao and Li ethnic groups. In summary, this study presented a systematic analysis of haplogroup distribution, genetic substructure of Han populations and genetic relationships between Han populations and surrounding ethnic minorities based on 89 Y-SNPs and 16 Y-STRs systematically. Our research supplemented valuable insights into population genetics and forensic genetics, and provided data support for the forensic application of Y chromosome. The integration of Y-SNP haplogroups with Y-STR haplotypes offers enhanced understanding of the genetic substructure within Han populations, which holds significance for both population genetics research and forensic science applications.

Multiple founding paternal lineages inferred from the newly-developed 639-plex Y-SNP panel suggested the complex admixture and migration history of Chinese people

BackgroundNon-recombining regions of the Y-chromosome recorded the evolutionary traces of male human populations and are inherited haplotype-dependently and male-specifically. Recent whole Y-chromosome sequencing studies have identified previously unrecognized population divergence, expansion and admixture processes, which promotes a better understanding and application of the observed patterns of Y-chromosome genetic diversity.ResultsHere, we developed one highest-resolution Y-chromosome single nucleotide polymorphism (Y-SNP) panel targeted for uniparental genealogy reconstruction and paternal biogeographical ancestry inference, which included 639 phylogenetically informative SNPs. We genotyped these loci in 1033 Chinese male individuals from 33 ethnolinguistically diverse populations and identified 256 terminal Y-chromosomal lineages with frequency ranging from 0.0010 (singleton) to 0.0687. We identified six dominant common founding lineages associated with different ethnolinguistic backgrounds, which included O2a2b1a1a1a1a1a1a1-M6539, O2a1b1a1a1a1a1a1-F17, O2a2b1a1a1a1a1b1a1b-MF15397, O2a2b2a1b1-A16609, O1b1a1a1a1b2a1a1-F2517, and O2a2b1a1a1a1a1a1-F155. The AMOVA and nucleotide diversity estimates revealed considerable differences and high genetic diversity among ethnolinguistically different populations. We constructed one representative phylogenetic tree among 33 studied populations based on the haplogroup frequency spectrum and sequence variations. Clustering patterns in principal component analysis and multidimensional scaling results showed a genetic differentiation between Tai-Kadai-speaking Li, Mongolic-speaking Mongolian, and other Sinitic-speaking Han Chinese populations. Phylogenetic topology inferred from the BEAST and Network relationships reconstructed from the popART further showed the founding lineages from culturally/linguistically diverse populations, such as C2a/C2b was dominant in Mongolian people and O1a/O1b was dominant in island Li people. We also identified many lineages shared by more than two ethnolinguistically different populations with a high proportion, suggesting their extensive admixture and migration history.ConclusionsOur findings indicated that our developed high-resolution Y-SNP panel included major dominant Y-lineages of Chinese populations from different ethnic groups and geographical regions, which can be used as the primary and powerful tool for forensic practice. We should emphasize the necessity and importance of whole sequencing of more ethnolinguistically different populations, which can help identify more unrecognized population-specific variations for the promotion of Y-chromosome-based forensic applications.

PSI-4 Investigation of a Naturally Occurring XX Male Bull

Abstract The Y chromosome–linked gene sex-determining region (SRY) is believed to be the master initiator of male sex determination in almost all eutherian mammals, functioning to upregulate expression of its direct target gene SRY-related HMG box–containing gene 9 (SOX9). XX maleness in humans is a rare syndrome with a frequency of 1 in 20,000–25,000 males. Approximately 90% of these cases carry some amount of the Y chromosome sequences due to an illegitimate recombination between X and Y chromosomes, where most of these cases have normal male genitalia. A more severe manifestation, the remaining 10% do not have any Y-chromosome sequences, and typically have ambiguous genitalia. In this study, we report a case of an infertile bull. Physical examination of this purebred Holstein bull showed a phenotypically androgynous animal with a normal penis, cryptorchidism (2 diminutive testicles, only one descended), and with a skeletal conformation and stance usually associated with female individuals. Lymphocytes from heparin sodium anticoagulated blood were cultured, harvested, and prepared according to conventional cell culture methods revealed a normal 60, XX karyotype with no numerical or structural chromosomal aberrations. PCR analyses were negative for the SRY gene, and positive for androgen receptor gene on the X chromosome. No sperm were found in seminal fluid collected in a routine breeding soundness exam at 17 months of age. The bull was slaughtered at 18 months of age, and its reproductive tract was collected. The testicle weighed 18.9 g (descended), and 42.4 g (abdominal). Neither testicle/epididymal unit was appropriate size or fully developed. Two very enlarged vesicular glands (43.36 g and 52.52 g, normally each ~ 20 g) were observed alongside what appeared to be uterine tissue. Histology showed no tubular grooves in the ampulla and an epididymis devoid of sperm. The development of a male phenotype in SRY-negative XX individuals suggests that in addition to SRY, there are autosomal or X-linked genes involved in the sex-determining pathway. In humans, a dose change of the SOX9 gene, loss of function mutations in the WNT4, FoxL2, and RSPO1 testicular suppressing genes, and genomic rearrangements within the regulatory region of the SRY-box transcription factor 3 (SOX3) gene have all been associated with XX male sex reversal. To investigate whether a mutation or genomic rearrangement was responsible for the male phenotype in the SRY-negative 60, XX bull, DNA was collected from the bull and his sire and dam for long-read trio sequencing. These data may provide evidence for either a loss of function mutation in a gene normally inhibiting testes formation in XX individuals, or a gain of function mutation in a gene downstream of SRY in the sex determining pathway.

Genetic insights into the social organization of Neanderthals

Genomic analyses of Neanderthals have previously provided insights into their population history and relationship to modern humans1–8, but the social organization of Neanderthal communities remains poorly understood. Here we present genetic data for 13 Neanderthals from two Middle Palaeolithic sites in the Altai Mountains of southern Siberia: 11 from Chagyrskaya Cave9,10 and 2 from Okladnikov Cave11—making this one of the largest genetic studies of a Neanderthal population to date. We used hybridization capture to obtain genome-wide nuclear data, as well as mitochondrial and Y-chromosome sequences. Some Chagyrskaya individuals were closely related, including a father–daughter pair and a pair of second-degree relatives, indicating that at least some of the individuals lived at the same time. Up to one-third of these individuals’ genomes had long segments of homozygosity, suggesting that the Chagyrskaya Neanderthals were part of a small community. In addition, the Y-chromosome diversity is an order of magnitude lower than the mitochondrial diversity, a pattern that we found is best explained by female migration between communities. Thus, the genetic data presented here provide a detailed documentation of the social organization of an isolated Neanderthal community at the easternmost extent of their known range.

CSYseq: The first Y-chromosome sequencing tool typing a large number of Y-SNPs and Y-STRs to unravel worldwide human population genetics.

Male-specific Y-chromosome (chrY) polymorphisms are interesting components of the DNA for population genetics. While single nucleotide polymorphisms (Y-SNPs) indicate distant evolutionary ancestry, short tandem repeats (Y-STRs) are able to identify close familial kinships. Detailed chrY analysis provides thus both biogeographical background information as paternal lineage identification. The rapid advancement of high-throughput massive parallel sequencing (MPS) technology in the past decade has revolutionized genetic research. Using MPS, single-base information of both Y-SNPs as Y-STRs can be analyzed in a single assay typing multiple samples at once. In this study, we present the first extensive chrY-specific targeted resequencing panel, the ‘CSYseq’, which simultaneously identifies slow mutating Y-SNPs as evolution markers and rapid mutating Y-STRs as patrilineage markers. The panel was validated by paired-end sequencing of 130 males, distributed over 65 deep-rooted pedigrees covering 1,279 generations. The CSYseq successfully targets 15,611 Y-SNPs including 9,014 phylogenetic informative Y-SNPs to identify 1,443 human evolutionary Y-subhaplogroup lineages worldwide. In addition, the CSYseq properly targets 202 Y-STRs, including 81 slow, 68 moderate, 27 fast and 26 rapid mutating Y-STRs to individualize close paternal relatives. The targeted chrY markers cover a high average number of reads (Y-SNP = 717, Y-STR = 150), easy interpretation, powerful discrimination capacity and chrY specificity. The CSYseq is interesting for research on different time scales: to identify evolutionary ancestry, to find distant family and to discriminate closely related males. Therefore, this panel serves as a unique tool valuable for a wide range of genetic-genealogical applications in interdisciplinary research within evolutionary, population, molecular, medical and forensic genetics.

Optimization of diagnostic strategy for non-invasive cell-free foetal RHD determination from maternal plasma.

The aim of the study was to optimize routine non-invasiveprenatal detection of fetal RHD gene from plasma of RhD-negative pregnant women (the median of gestational age was 25weeks, range 10-38) to detect RhD materno-fetal incompatibility and to avoid the redundant immunoprophylaxis. Initially only one exon of RHD gene (exon 10) was investigated in 281 plasma samples (144 verified after delivery), in the second phase three RHD exons (5, 7, 10) were analyzed in 246 samples of plasma and maternal genomic DNA (204 verified) by real-time PCR method. Detection of Y-chromosomal sequence DYS-14 and five X-chromosomal insertion/deletion polymorphisms was used to confirm the fetal cfDNA detectability in plasma. Specific polymorphisms in RHD gene were detected by sequence-specific primer PCR in nine samples. When only the RHD exon 10 was tested, 2·8% of verified samples were false positive and 3·5% false negative. With three RHD exons (5, 7, 10) and maternal genomic DNA testing, only one case was false negative (0·5%). Nine samples were inconclusive due to RHD-positive results in maternal genomic DNA. These samples were analyzed for specific mutations in RHD gene. Combination of both methods for fetal cfDNA verification succeeded in 75% of tested group. Implementation of analysis of three RHD exons and maternal genomic DNA to routine practice lowers dramatically the ratio of false positive and negative results. This method enables more accurate determination of fetalRHD status with the reduction of unnecessary medical care and RhD immunoprophylaxis.

Shared paternal ancestry of Han, Tai-Kadai-speaking, and Austronesian-speaking populations as revealed by the high resolution phylogeny of O1a-M119 and distribution of its sub-lineages within China.

The aim of this research was to explore the origin, diversification, and demographic history of O1a-M119 over the past 10,000 years, as well as its role during the formation of East Asian and Southeast Asian populations, particularly the Han, Tai-Kadai-speaking, and Austronesian-speaking populations. Y-chromosome sequences (n = 141) of the O1a-M119 lineage, including 17 newly generated in this study, were used to reconstruct a revised phylogenetic tree with age estimates, and identify sub-lineages. The geographic distribution of 12 O1a-M119 sub-lineages was summarized, based on 7325 O1a-M119 individuals identified among 60,009 Chinese males. A revised phylogenetic tree, age estimation, and distribution maps indicated continuous expansion of haplogroup O1a-M119 over the past 10,000 years, and differences in demographic history across geographic regions. We propose several sub-lineages of O1a-M119 as founding paternal lineages of Han, Tai-Kadai-speaking, and Austronesian-speaking populations. The sharing of several young O1a-M119 sub-lineages with expansion times less than 6000 years between these three population groups supports a partial common ancestry for them in the Neolithic Age; however, the paternal genetic divergence pattern is much more complex than previous hypotheses based on ethnology, archeology, and linguistics. Our analyses contribute to a better understanding of the demographic history of O1a-M119 sub-lineages over the past 10,000 years during the emergence of Han, Austronesians, Tai-Kadai-speaking populations. The data described in this study will assist in understanding of the history of Han, Tai-Kadai-speaking, and Austronesian-speaking populations from ethnology, archeology, and linguistic perspectives in the future.

Detection of cell-free foetal DNA fraction in female-foetus bearing pregnancies using X-chromosomal insertion/deletion polymorphisms examined by digital droplet PCR

In families with X-linked recessive diseases, foetal sex is determined prenatally by detection of Y-chromosomal sequences in cell-free foetal DNA (cffDNA) in maternal plasma. The same procedure is used to confirm the cffDNA presence during non-invasive prenatal RhD incompatibility testing but there are no generally accepted markers for the detection of cffDNA fraction in female-foetus bearing pregnancies. We present a methodology allowing the detection of paternal X-chromosomal alleles on maternal background and the confirmation of female sex of the foetus by positive amplification signals. Using digital droplet PCR (ddPCR) we examined X-chromosomal INDEL (insertion/deletion) polymorphisms: rs2307932, rs16397, rs16637, rs3048996, rs16680 in buccal swabs of 50 females to obtain the population data. For all INDELs, we determined the limits of detection for each ddPCR assay. We examined the cffDNA from 63 pregnant women bearing Y-chromosome negative foetuses. The analysis with this set of INDELs led to informative results in 66.67% of examined female-foetus bearing pregnancies. Although the population data predicted higher informativity (74%) we provided the proof of principle of this methodology. We successfully applied this methodology in prenatal diagnostics in a family with Wiscott–Aldrich syndrome and in pregnancies tested for the risk of RhD incompatibility.

Mutation Rate Variability across Human Y-Chromosome Haplogroups.

A common assumption in dating patrilineal events using Y-chromosome sequencing data is that the Y-chromosome mutation rate is invariant across haplogroups. Previous studies revealed interhaplogroup heterogeneity in phylogenetic branch length. Whether this heterogeneity is caused by interhaplogroup mutation rate variation or nongenetic confounders remains unknown. Here, we analyzed whole-genome sequences from cultured cells derived from >1,700 males. We confirmed the presence of branch length heterogeneity. We demonstrate that sex-chromosome mutations that appear within cell lines, which likely occurred somatically or in vitro (and are thus not influenced by nongenetic confounders) are informative for germline mutational processes. Using within-cell-line mutations, we computed a relative Y-chromosome somatic mutation rate, and uncovered substantial variation (up to 83.3%) in this proxy for germline mutation rate among haplogroups. This rate positively correlates with phylogenetic branch length, indicating that interhaplogroup mutation rate variation is a likely cause of branch length heterogeneity.

Massively parallel sequencing of sex-chromosomal STRs in Saudi Arabia reveals patrilineage-associated sequence variants

Massively parallel sequencing (MPS) of forensic STRs has the potential to reveal additional allele diversity compared to conventional capillary electrophoresis (CE) typing strategies, but population studies are currently relatively few in number. The Verogen ForenSeq™ DNA Signature Prep Kit includes both Y-STRs and X-STRs among its targeted loci, and here we report the sequences of these loci, analysed using Verogen’s ForenSeq™ Universal Analysis Software (UAS) v1.3 and STRait Razor v3.0, in a representative sample of 89 Saudi Arabian males. We identified 56 length variants (equivalent to CE alleles) and 75 repeat sequence sub-variants across the six X-STRs analysed; equivalent figures for the set of 24 Y-STRs were 147 and 192 respectively. We also observed two flanking sequence variants for the X-, and six for the Y-STRs. Recovery of sequence data and concordance with CE data (where available) across the tested loci was good, though rare flanking variation affected interpretation and allele calling at DYF387S1 and DXS7132. Examination of flanking sequences of the Y-STRs revealed five SNPs (L255, M4790, BY7692, Z16708 and S17543) previously shown to define specific haplogroups by Y-chromosome sequencing. These define Y-haplogroups in 62 % of our sample, a proportion that increases to 91 % when haplogroup-associated repeat-sequence motifs are also considered. A population-level comparison of the Saudi Arabian X-STRs with a global sample showed our dataset to be part of a large cluster of populations of West Eurasian and Middle Eastern origin.

Paternal gene pool of Malays in Southeast Asia and its applications for the early expansion of Austronesians.

The origin and differentiation of Austronesian populations and their languages have long fascinated linguists, archeologists, and geneticists. However, the founding process of Austronesians and when they separated from their close relatives, such as the Daic and Austro-Asiatic populations in the mainland of Asia, remain unclear. In this study, we explored the paternal origin of Malays in Southeast Asia and the early differentiation of Austronesians. We generated whole Y-chromosome sequences of 50 Malays and co-analyzed 200 sequences from other Austronesians and related populations. We generated a revised phylogenetic tree with time estimation. We identified six founding paternal lineages among the studied Malays samples. These founding lineages showed a surprisingly coincident expansion age at 5000 to 6000 years ago. We also found numerous mostly close related samples of the founding lineages of Malays among populations from Mainland of Asia. Our analyses provided a refined phylogenetic resolution for the dominant paternal lineages of Austronesians found by previous studies. We suggested that the co-expansion of numerous founding paternal lineages corresponds to the initial differentiation of the most recent common ancestor of modern Austronesians. The splitting time and divergence pattern in perspective of paternal Y-chromosome evidence are highly consistent with the previous theories of ethnologists, linguists, and archeologists.

Paternal origin of Tungusic-speaking populations: Insights from the updated phylogenetic tree of Y-chromosome haplogroup C2a-M86.

Haplogroup C2a-M48 is the predominant paternal lineage of Tungusic-speaking populations, one of the largest population groups in Siberia. Up until now, the origins and dispersal of Tungusic-speaking populations have remained unclear. In this study, the demographic history of Tungusic-speaking populations was explored using the phylogenetic analysis of haplogroup C2a-M86, the major subbranch of C2a-M48. In total, 18 newly generated Y chromosome sequences from C2a-M48 males and 20 previously available Y-chromosome sequences from this haplogroup were analyzed. A highly revised phylogenetic tree of haplogroup C2a-M86 with age estimates was reconstructed. Frequencies of this lineage in the literature were collected and a comprehensive analysis of this lineage in 13 022 individuals from 245 populations in Eurasia was performed. The distribution map of C2a-M48 indicated the most probable area of origin and diffusion route of this paternal lineage in North Eurasia. Most C2a-M86 samples from Tungusic-speaking populations belonged to the sublineage C2a-F5484, which emerged about 3300 years ago. We identified six unique sublineages corresponding to the Manchu, Evenks, Evens, Oroqen, and Daurpopulations; these sublineages diverged gradually over the past 1900 years. Notably, we observed a clear north-south dichotomous structure for sublineages derived from C2a-F5484, consistent with the internal north-south divergence of Tungusic languages and ethnic groups. We identified the important founding paternal haplogroup, C2a-F5484, for Tungusic-speaking populations as well as numerous unique subgroups of this haplogroup. We propose that the timeframe for the divergence of C2a-F5484 corresponds with the early differentiation of ancestral Tungusic-speaking populations.

Ancient Patrilineal Lines and Relatively High ECAY Diversity Preserved in Indigenous Horses Revealed With Novel Y-Chromosome Markers.

Extremely low nucleotide diversity of modern horse Y-chromosome has been reported, and only poor phylogenetic resolution could be resulted from limited Y-chromosome markers. In this study, three types of horse Y-chromosome markers, including Single-nucleotide polymorphisms (SNPs), copy number variants (CNVs), and allele-specific CNVs, were developed by screening more than 300 male horses from 23 indigenous Chinese horse populations and 4 imported horse breeds. Fourteen segregating sites including a novel SNP in the AMELY gene were found in approximately 53 kb of male-specific Y-chromosome sequences. CNVs were detected at 11 of 14 sites, while allele-specific CNVs at 6 polymorphic sites in repeated fragments were also determined. The phylogenetic analyses with the SNPs identified in this study and previously published 51 SNPs obtained mainly from European horses showed that indigenous Chinese horses exhibit much deeper divergence than European and Middle Eastern horses, while individuals of Chinese horses with the C allele of the AMELY gene constituted the most ancient group. Via SNPs, CNVs, and allele-specific CNVs, much higher diversity of paternal lines can be detected than those identified with merely SNPs. Our results indicated that there are ancient paternal horse lines preserved in southwestern China, which sheds new light on the domestication and immigration of horses, and suggest that the priorities of the conservation should be given to the ancient and rare paternal lines. These three marker types provided finer phylogenetic resolution of horse patrilineal lines, and the strategies used in the present study also provide valuable reference for the genetic studies of other mammalian patrilineages.

Phylogenetic analysis of the Y-chromosome haplogroup C2b-F1067, a dominant paternal lineage in Eastern Eurasia.

Human Y-chromosome haplogroup C2b-F1067 is one of the dominant paternal lineages of populations in Eastern Eurasia. In order to explore the origin, diversification, and expansion of this haplogroup, we generated 206 new Y-chromosome sequences from C2b-F1067 males and coanalyzed 220 Y-chromosome sequences of this haplogroup. BEAST software was used to reconstruct a revised phylogenetic tree of haplogroup C2b-F1067 with age estimates. The revised phylogeny of C2b-F1067 included 155 sublineages, 1986 non-private variants, and >6000 private variants. The age estimation suggested that the initial splitting of C2b-F1067 happened at about 32.8 thousand years ago (kya) and the major sublineages of this haplgroup experienced continuous expansion in the most recent 10,000 years. We identified numerous sublineages that were nearly specific for Korean, Mongolian, Chinese, and other ethnic minorities in China. In particular, we evaluated the candidate-specific lineage for the Dayan Khan family and the Confucius family, the descendants of the ruling family of the Chinese Shang dynasty. These findings suggest that ancient populations with varied C2b-F1067 sublineages played an important role during the formation of most modern populations in Eastern Eurasia, and thus eventually became the founding paternal lineages of these populations.

De novo assembly of the olive fruit fly (Bactrocera oleae) genome with linked-reads and long-read technologies minimizes gaps and provides exceptional Y chromosome assembly

BackgroundThe olive fruit fly, Bactrocera oleae, is the most important pest in the olive fruit agribusiness industry. This is because female flies lay their eggs in the unripe fruits and upon hatching the larvae feed on the fruits thus destroying them. The lack of a high-quality genome and other genomic and transcriptomic data has hindered progress in understanding the fly’s biology and proposing alternative control methods to pesticide use.ResultsGenomic DNA was sequenced from male and female Demokritos strain flies, maintained in the laboratory for over 45 years. We used short-, mate-pair-, and long-read sequencing technologies to generate a combined male-female genome assembly (GenBank accession GCA_001188975.2). Genomic DNA sequencing from male insects using 10x Genomics linked-reads technology followed by mate-pair and long-read scaffolding and gap-closing generated a highly contiguous 489 Mb genome with a scaffold N50 of 4.69 Mb and L50 of 30 scaffolds (GenBank accession GCA_001188975.4). RNA-seq data generated from 12 tissues and/or developmental stages allowed for genome annotation. Short reads from both males and females and the chromosome quotient method enabled identification of Y-chromosome scaffolds which were extensively validated by PCR.ConclusionsThe high-quality genome generated represents a critical tool in olive fruit fly research. We provide an extensive RNA-seq data set, and genome annotation, critical towards gaining an insight into the biology of the olive fruit fly. In addition, elucidation of Y-chromosome sequences will advance our understanding of the Y-chromosome’s organization, function and evolution and is poised to provide avenues for sterile insect technique approaches.

Analysis of whole Y-chromosome sequences reveals the Japanese population history in the Jomon period

The Jomon and the Yayoi are considered to be the two major ancestral populations of the modern mainland Japanese. The Jomon people, who inhabited mainland Japan, admixed with Yayoi immigrants from the Asian continent. To investigate the population history in the Jomon period (14,500–2,300 years before present [YBP]), we analyzed whole Y-chromosome sequences of 345 Japanese males living in mainland Japan. A phylogenetic analysis of East Asian Y chromosomes identified a major clade (35.4% of mainland Japanese) consisting of only Japanese Y chromosomes, which seem to have originated from indigenous Jomon people. A Monte Carlo simulation indicated that ~70% of Jomon males had Y chromosomes in this clade. The Bayesian skyline plots of 122 Japanese Y chromosomes in the clade detected a marked decrease followed by a subsequent increase in the male population size from around the end of the Jomon period to the beginning of the Yayoi period (2,300 YBP). The colder climate in the Late to Final Jomon period may have resulted in critical shortages of food for the Jomon people, who were hunter-gatherers, and the rice farming introduced by Yayoi immigrants may have helped the population size of the Jomon people to recover.

Relating Clans Ao and Aisin Gioro from northeast China by whole Y-chromosome sequencing.

The Y-chromosome haplogroup C2b1a3a2-F8951 is the paternal lineage of the Aisin Gioro clan, the most important brother branch of the famous Mongolic-speaking population characteristic haplogroup C2*-Star Cluster (C2b1a3a1-F3796). However, investigations on its internal phylogeny are still limited. In this study, we used whole Y-chromosome sequencing to update its phylogenetic tree. In the revised tree, C2b1a3a2-F8951 and C2*-Star Cluster differentiated 3852 years ago (95% CI = 3295-4497). Approximately 3558 years ago (95% CI = 3013-4144), C2b1a3a2-F8951 was divided into two main subclades, C2b1a3a2a-F14753 and C2b1a3a2b-F5483. Currently, samples of C2b1a3a2-F8951 were mainly from the House of Aisin Gioro clan, the Ao family from Daur and some individuals mainly from northeast China. Although other haplogroups are also found in the Ao family, including C2b1a2-M48, C2b1a3a1-F3796, C2a1b-F845, and N1c-M178, the haplogroup C2b1a3a2-F8951 is still the most distinct genetic component. For haplogroup C2b1a3a2-F8951, the time of the most recent common ancestor of the House of Aisin Gioro clan and the Ao family were both very late, just a few hundred years ago. Some family-specific Y-SNPs of the House of Aisin Gioro and the Ao family were also discovered. This revision evidently improved the resolving power of Y-chromosome phylogeny in northeast Asia, deepening our understanding of the origin of these two families, even the Mongolic-speaking population.

The Red Fox Y-Chromosome in Comparative Context.

While the number of mammalian genome assemblies has proliferated, Y-chromosome assemblies have lagged behind. This discrepancy is caused by biological features of the Y-chromosome, such as its high repeat content, that present challenges to assembly with short-read, next-generation sequencing technologies. Partial Y-chromosome assemblies have been developed for the cat (Felis catus), dog (Canis lupus familiaris), and grey wolf (Canis lupus lupus), providing the opportunity to examine the red fox (Vulpes vulpes) Y-chromosome in the context of closely related species. Here we present a data-driven approach to identifying Y-chromosome sequence among the scaffolds that comprise the short-read assembled red fox genome. First, scaffolds containing genes found on the Y-chromosomes of cats, dogs, and wolves were identified. Next, analysis of the resequenced genomes of 15 male and 15 female foxes revealed scaffolds containing male-specific k-mers and patterns of inter-sex copy number variation consistent with the heterogametic chromosome. Analyzing variation across these two metrics revealed 171 scaffolds containing 3.37 Mbp of putative Y-chromosome sequence. The gene content of these scaffolds is consistent overall with that of the Y-chromosome in other carnivore species, though the red fox Y-chromosome carries more copies of BCORY2 and UBE1Y than has been reported in related species and fewer copies of SRY than in other canids. The assignment of these scaffolds to the Y-chromosome serves to further characterize the content of the red fox draft genome while providing resources for future analyses of canid Y-chromosome evolution.

Y-chromosomal connection between Hungarians and geographically distant populations of the Ural Mountain region and West Siberia

Hungarians who live in Central Europe today are one of the westernmost Uralic speakers. Despite of the proposed Volga-Ural/West Siberian roots of the Hungarian language, the present-day Hungarian gene pool is highly similar to that of the surrounding Indo-European speaking populations. However, a limited portion of specific Y-chromosomal lineages from haplogroup N, sometimes associated with the spread of Uralic languages, link modern Hungarians with populations living close to the Ural Mountain range on the border of Europe and Asia. Here we investigate the paternal genetic connection between these spatially separated populations. We reconstruct the phylogeny of N3a4-Z1936 clade by using 33 high-coverage Y-chromosomal sequences and estimate the coalescent times of its sub-clades. We genotype close to 5000 samples from 46 Eurasian populations to show the presence of N3a4-B539 lineages among Hungarians and in the populations from Ural Mountain region, including Ob-Ugric-speakers from West Siberia who are geographically distant but linguistically closest to Hungarians. This sub-clade splits from its sister-branch N3a4-B535, frequent today among Northeast European Uralic speakers, 4000–5000 ya, which is in the time-frame of the proposed divergence of Ugric languages.