Pseudoautosomal Research Articles

37 Genomic improvement of the horse X chromosome and characterization of the pseudoautosomal boundary

The horse reference genome has been improved with the release of EquCab3 and the first Y chromosome reference. The most complex sequences, such as those in the sex chromosomes, are still unresolved. These include large amplicons, repeats, and the pseudoautosomal boundary (PAB). We initiated a comprehensive study specifically to improve the assembly of the horse sex chromosomes. To refine the assembly of complex regions, we are utilizing 3 new technologies: trio-binning, Hi-C and Bionano optical mapping. Trio-binning uses long read sequences from F1 interspecific hybrids and short reads from parent species. High molecular weight blood DNA was extracted from a female hinny and sequenced on 2 PacBio Sequel cells. Paired-end Illumina reads (150bp) for horse ( Twilight ) and donkey ( Willy ) were obtained from SRA. These sequences as well as the hinny long reads were assembled with trio-binning function of the Canu assembler program. The initial assembly is 2.5 Gb separated into 1,757 contigs with an N50 of 41.5 Mb. We have completed Hi-C sequencing and Bionano optical maps for one thoroughbred stallion ( Bravo ). Both technologies are needed to scaffold the trio-binning Canu assembly. Our initial goal was to use this assembly to better define the PAB in the horse. The PAB demarcates the end of the pseudoautosomal region (PAR) where X-Y recombination stops. Despite the evolutionary and biological importance of the PAB, the region has been characterized at molecular level in only a few species. The PAB of the horse is presently not well defined. Previously, we identified and Sanger sequenced 4 BAC clones - 2 spanning PAB-X and PAB-Y. To identify the PAB of the horse, BAC sequences were aligned to the Y assembly, EquCab3 and a 42 Mb-size contig from trio-binning assembly which corresponds to the short arm of the X. We identified a region on both X and Y where X-Y homology drops from over 97% (PAR) to almost zero, indicative of the PAB. This region corresponds to the location of the XKR3Y gene in the Y but is not well-annotated in the X. We also identified a duplication and an inversion in EquCab3 which was not consistent with the corresponding region in the X-BACs, or the new 42 Mb Xp contig, suggesting a mis-assembly in EquCab3. We believe that these approaches combined will also resolve other complex and repetitive portions in the horse sex chromosomes.

Evolutionary dynamics of the human pseudoautosomal regions.

Recombination between the X and Y human sex chromosomes is limited to the two pseudoautosomal regions (PARs) that present quite distinct evolutionary origins. Despite the crucial importance for male meiosis, genetic diversity patterns and evolutionary dynamics of these regions are poorly understood. In the present study, we analyzed and compared the genetic diversity of the PAR regions using publicly available genomic sequences encompassing both PAR1 and PAR2. Comparisons were performed through allele diversities, linkage disequilibrium status and recombination frequencies within and between X and Y chromosomes. In agreement with previous studies, we confirmed the role of PAR1 as a male-specific recombination hotspot, but also observed similar characteristic patterns of diversity in both regions although male recombination occurs at PAR2 to a much lower extent (at least one recombination event at PAR1 and in ≈1% in normal male meioses at PAR2). Furthermore, we demonstrate that both PARs harbor significantly different allele frequencies between X and Y chromosomes, which could support that recombination is not sufficient to homogenize the pseudoautosomal gene pool or is counterbalanced by other evolutionary forces. Nevertheless, the observed patterns of diversity are not entirely explainable by sexually antagonistic selection. A better understanding of such processes requires new data from intergenerational transmission studies of PARs, which would be decisive on the elucidation of PARs evolution and their role in male-driven heterosomal aneuploidies.

Hypospadias in ring X syndrome

Ring X is a chromosomal anomaly mainly seen in females with turner syndrome and usually present in mosaic form with 45,X cells (45,X/46,X,r(X)) because of their mitotic instability. In males it is an extremely rare finding because large nullisomy for X chromosome material is likely not compatible with survival. Only two cases of male with ring chromosome X were previously reported. We report here a four-year-old male with ring chromosome X characterized using Karyotype, FISH and array CGH and presenting short stature, microcephaly and hypospadias. Molecular investigations showed 923 Kb terminal deletion on the pseudoautosomal region 1 (PAR1) including SHOX gene followed by a duplication of 2.4 Mb. The absence of functional nullisomy because of a second copy of deleted genes was present in chromosome Y PAR1 region may explain the compatibility with survival in our case of male with ring X. Short stature common with the two previously reported cases is likely related to SHOX gene deletion but also to the effect of “ring syndrome”. However, hypospadias was not reported in the previous cases and can be due to the associated duplication outside PAR1 region including in particular PRKX gene coding for a protein involved in urogenital system morphogenesis.

Telomere-to-telomere assembly of the genome of an individual Oikopleura dioica from Okinawa using Nanopore-based sequencing

BackgroundThe larvacean Oikopleura dioica is an abundant tunicate plankton with the smallest (65–70 Mbp) non-parasitic, non-extremophile animal genome identified to date. Currently, there are two genomes available for the Bergen (OdB3) and Osaka (OSKA2016) O. dioica laboratory strains. Both assemblies have full genome coverage and high sequence accuracy. However, a chromosome-scale assembly has not yet been achieved.ResultsHere, we present a chromosome-scale genome assembly (OKI2018_I69) of the Okinawan O. dioica produced using long-read Nanopore and short-read Illumina sequencing data from a single male, combined with Hi-C chromosomal conformation capture data for scaffolding. The OKI2018_I69 assembly has a total length of 64.3 Mbp distributed among 19 scaffolds. 99% of the assembly is contained within five megabase-scale scaffolds. We found telomeres on both ends of the two largest scaffolds, which represent assemblies of two fully contiguous autosomal chromosomes. Each of the other three large scaffolds have telomeres at one end only and we propose that they correspond to sex chromosomes split into a pseudo-autosomal region and X-specific or Y-specific regions. Indeed, these five scaffolds mostly correspond to equivalent linkage groups in OdB3, suggesting overall agreement in chromosomal organization between the two populations. At a more detailed level, the OKI2018_I69 assembly possesses similar genomic features in gene content and repetitive elements reported for OdB3. The Hi-C map suggests few reciprocal interactions between chromosome arms. At the sequence level, multiple genomic features such as GC content and repetitive elements are distributed differently along the short and long arms of the same chromosome.ConclusionsWe show that a hybrid approach of integrating multiple sequencing technologies with chromosome conformation information results in an accurate de novo chromosome-scale assembly of O. dioica’s highly polymorphic genome. This genome assembly opens up the possibility of cross-genome comparison between O. dioica populations, as well as of studies of chromosomal evolution in this lineage.

An 8.22 Mb Assembly and Annotation of the Alpaca (Vicugna pacos) Y Chromosome.

The unique evolutionary dynamics and complex structure make the Y chromosome the most diverse and least understood region in the mammalian genome, despite its undisputable role in sex determination, development, and male fertility. Here we present the first contig-level annotated draft assembly for the alpaca (Vicugna pacos) Y chromosome based on hybrid assembly of short- and long-read sequence data of flow-sorted Y. The latter was also used for cDNA selection providing Y-enriched testis transcriptome for annotation. The final assembly of 8.22 Mb comprised 4.5 Mb of male specific Y (MSY) and 3.7 Mb of the pseudoautosomal region. In MSY, we annotated 15 X-degenerate genes and two novel transcripts, but no transposed sequences. Two MSY genes, HSFY and RBMY, are multicopy. The pseudoautosomal boundary is located between SHROOM2 and HSFY. Comparative analysis shows that the small and cytogenetically distinct alpaca Y shares most of MSY sequences with the larger dromedary and Bactrian camel Y chromosomes. Most of alpaca X-degenerate genes are also shared with other mammalian MSYs, though WWC3Y is Y-specific only in alpaca/camels and the horse. The partial alpaca Y assembly is a starting point for further expansion and will have applications in the study of camelid populations and male biology.

Xp;Yq Unbalanced Translocation with Pseudoautosomal Region Aberrations in a Natural Two-Generation Transmission.

Translocations involving X and Y chromosomes rarely occur in humans and may affect reproductive function. We investigated an Xp:Yq unbalanced translocation with pseudoautosomal region (PAR) aberrations in a natural two-generation transmission. We report the case of an azoospermic male and his fertile mother without any other abnormal clinical phenotypes, except for short stature. Cytogenetic methods, including karyotyping and fluorescence in situ hybridization (FISH), revealed the translocation. Chromosomal microarray comparative genomic hybridization (array-CGH) was used to investigate the regions of Xp partial deletion and Yq partial duplication. Final chromosome karyotypes in the peripheral blood of the infertile male and his mother were 46,Y,der(X)t(X;Y)(p22.33;q11.22) and 46,X,der(X)t(X;Y)(p22.33;q11.22), respectively. Short-stature-homeobox gene deletion was responsible for the short stature in both subjects. PAR aberrations and AZFc duplication may be a direct genetic risk factor for spermatogenesis. This report further supports the use of routine karyotype analysis, FISH-based technology, and array-CGH analysis to identify derivative chromosomes in a complex rearrangement.

Evolutionary stasis of the pseudoautosomal boundary in strepsirrhine primates.

Sex chromosomes are typically comprised of a non-recombining region and a recombining pseudoautosomal region. Accurately quantifying the relative size of these regions is critical for sex-chromosome biology both from a functional and evolutionary perspective. The evolution of the pseudoautosomal boundary (PAB) is well documented in haplorrhines (apes and monkeys) but not in strepsirrhines (lemurs and lorises). Here, we studied the PAB of seven species representing the main strepsirrhine lineages by sequencing a male and a female genome in each species and using sex differences in coverage to identify the PAB. We found that during primate evolution, the PAB has remained unchanged in strepsirrhines whereas several recombination suppression events moved the PAB and shortened the pseudoautosomal region in haplorrhines. Strepsirrhines are well known to have much lower sexual dimorphism than haplorrhines. We suggest that mutations with antagonistic effects between males and females have driven recombination suppression and PAB evolution in haplorrhines.

Meiosis reveals the early steps in the evolution of a neo-XY sex chromosome pair in the African pygmy mouse Mus minutoides.

Sex chromosomes of eutherian mammals are highly different in size and gene content, and share only a small region of homology (pseudoautosomal region, PAR). They are thought to have evolved through an addition-attrition cycle involving the addition of autosomal segments to sex chromosomes and their subsequent differentiation. The events that drive this process are difficult to investigate because sex chromosomes in almost all mammals are at a very advanced stage of differentiation. Here, we have taken advantage of a recent translocation of an autosome to both sex chromosomes in the African pygmy mouse Mus minutoides, which has restored a large segment of homology (neo-PAR). By studying meiotic sex chromosome behavior and identifying fully sex-linked genetic markers in the neo-PAR, we demonstrate that this region shows unequivocal signs of early sex-differentiation. First, synapsis and resolution of DNA damage intermediates are delayed in the neo-PAR during meiosis. Second, recombination is suppressed or largely reduced in a large portion of the neo-PAR. However, the inactivation process that characterizes sex chromosomes during meiosis does not extend to this region. Finally, the sex chromosomes show a dual mechanism of association at metaphase-I that involves the formation of a chiasma in the neo-PAR and the preservation of an ancestral achiasmate mode of association in the non-homologous segments. We show that the study of meiosis is crucial to apprehend the onset of sex chromosome differentiation, as it introduces structural and functional constrains to sex chromosome evolution. Synapsis and DNA repair dynamics are the first processes affected in the incipient differentiation of X and Y chromosomes, and they may be involved in accelerating their evolution. This provides one of the very first reports of early steps in neo-sex chromosome differentiation in mammals, and for the first time a cellular framework for the addition-attrition model of sex chromosome evolution.

Optical Mapping, a Promising Alternative to Gold Standard Cytogenetic Approaches in Acute Lymphoblastic Leukemias: A Blind Comparison on 10 Patients

Acute lymphoblastic leukemias (ALL) are characterized by a large number of genetic abnormalities of prognostic and theranostic interest. Their detection requires the use of several conventional and molecular cytogenetic techniques in order to establish the most exhaustive genetic profile: conventional karyotype, FISH, SNP-Array and RT-MLPA for detection of fusion transcripts. Optical genome mapping technique is based on analysis of ultra-high molecular weight DNA (UHMW: 150 kbp-2.5 Mbp). It allows an automatic, high resolution, genome-wide detection of structural anomalies, including balanced translocations and copy number anomalies, making this technology a promising tool in ALL. However, very few data are available in ALL so far. Our objective was to compare optical mapping (Bionano Genomics) to standard techniques performed routinely (combination of karyotype, FISH, SNP-array and RT-MLPA) in 10 selected B or T-ALL patients carrying at least one major genetic abnormality including t(10;14) TLX1-TRA, t(9;22) BCR-ABL1 (e1a2 and e14a2 transcripts), t(10;11) MLLT10-PICALM, t(11;19) KMT2A-MLLT1, typical hyperdiploid profile, iAMP21, SET-NUP214 fusion, various deletions involving CDKN2A/CDKN2B, IKZF1, PAX5, PTEN, ERG and FBXW7. Abnormalities on SNP-array were retained when allele fraction was > 10% and larger than 1Mbp or > 20 kbp and involving one of the 186 genes defined as potentially relevant in leukemogenesis. A total of 79 anomalies meeting these criteria in SNP-array, and/or found by other standard techniques were identified (45 deletions, 16 whole chromosome gains, 10 duplications or partial gains, and 8 translocations). Bone marrow DMSO pellets from the 10 patients, all treated at the CHU Amiens-Picardie, were used for this study. The UHMW DNA extraction technique provided molecules with an average N50 (≥ 150 kbp) of 286 kbp, and an average coverage of 339x. Direct imaging of the DNA molecules was performed on the Bionano Saphyr system and structural variants were automatically detected using the Rare Variant Pipeline. Abnormalities analysis was blindly performed by two independent operators on the Bionano Access software (v 1.5.2.), one cytogeneticist with no training in optical mapping (operator 1) and the other optical mapping technician without training in cytogenetics (operator 2). Among the 79 relevant anomalies, 72 (91.1%) were found by operator 1 and 71 (89.9%) by operator 2. All (8/8) translocations, 43/45 (95.6%) of deletions, 8/10 (80%) of gains were found by both operators whereas 13/16 (81.25%) and 12/16 (75%) of chromosomal gains were found by operator 1 and 2 respectively. Concerning the 8 discrepancies between optical mapping and standard techniques, 4 of them were filtering issues. Three were low allelic frequency large gains clearly visible in raw data but rejected by software's CNV filters, totally, or partially for a trisomy 3 that caused the discrepancy between the two observers. The fourth was a deletion of VPREB1 gene, erroneously eliminated by operators because included in recurrent deletion zones of IGL locus. The 4 others discrepancies were non detected anomalies. Two of them were copy number variations in Xp22.33 (pseudoautosomal region), a non-covered area with the optical mapping technology, one of which was a P2RY8-CRLF2 fusion. Another sub-clonal gain (trisomy 21) found by karyotype and FISH (quantified at 7% of all nuclei) but not by SNP-array, was missed, even in raw data. The last was a near-tetraploid subclone, found only on karyotype, that might be strictly tetraploid and by definition undetectable by quantitative techniques. Interestingly, new anomalies, not detected by standard techniques, were revealed, including a translocation t(11;14)(p13;q11.2) LMO2-TRA, a translocation with inversion t(7;8)(q34;q24.1) TRB-MYC and a translocation t(14;20)(q32.33;q13.13) IGH-CEBPB. At last, in 3 cases, unidentified complex rearrangements by cytogenetics could be resolved using optical mapping. Thus, optical mapping represents an innovative, time and cost effective alternative to the different cytogenetic techniques currently performed routinely for ALL genetic characterization. It enables identification of complex cytogenetic events, including those currently inaccessible to standard techniques such as chromotripsis or certain complex translocations. Disclosures Lee: Bionano Genomics, San Diego: Current Employment. Ferret:Bionano Genomics: Research Funding.

Locating the Sex Determining Region of Linkage Group 12 of Guppy (Poecilia reticulata).

Despite over 100 years of study, the location of the fully sex-linked region of the guppy (Poecilia reticulata) carrying the male-determining locus, and the regions where the XY pair recombine, remain unclear. Previous population genomics studies to determine these regions used small samples from recently bottlenecked captive populations, which increase the false positive rate of associations between individuals’ sexes and SNPs. Using new data from multiple natural populations, we show that a recently proposed candidate for this species’ male-determining gene is probably not completely sex-linked, leaving the maleness factor still unidentified. Variants in the chromosome 12 region carrying the candidate gene sometimes show linkage disequilibrium with the sex-determining factor, but no consistently male-specific variant has yet been found. Our genetic mapping with molecular markers spread across chromosome 12 confirms that this is the guppy XY pair. We describe two families with recombinants between the X and Y chromosomes, which confirm that the male-determining locus is in the region identified by all previous studies, near the terminal pseudo-autosomal region (PAR), which crosses over at a very high rate in males. We correct the PAR marker order, and assign two unplaced scaffolds to the PAR. We also detect a duplication, with one copy in the male-determining region, explaining signals of sex linkage in a more proximal region.

Twenty Years Later: A Comprehensive Review of the X Chromosome Use in Forensic Genetics.

The unique structure of the X chromosome shaped by evolution has led to the present gender-specific genetic differences, which are not shared by its counterpart, the Y chromosome, and neither by the autosomes. In males, recombination between the X and Y chromosomes is limited to the pseudoautosomal regions, PAR1 and PAR2; therefore, in males, the X chromosome is (almost) entirely transmitted to female offspring. On the other hand, the X chromosome is present in females with two copies that recombine along the whole chromosome during female meiosis and that is transmitted to both female and male descendants. These transmission characteristics, besides the obvious clinical impact (sex chromosome aneuploidies are extremely frequent), make the X chromosome an irreplaceable genetic tool for population genetic-based studies as well as for kinship and forensic investigations. In the early 2000s, the number of publications using X-chromosomal polymorphisms in forensic and population genetic applications increased steadily. However, nearly 20 years later, we observe a conspicuous decrease in the rate of these publications. In light of this observation, the main aim of this article is to provide a comprehensive review of the advances and applications of X-chromosomal markers in population and forensic genetics over the last two decades. The foremost relevant topics are addressed as: (i) developments concerning the number and types of markers available, with special emphasis on short tandem repeat (STR) polymorphisms (STR nomenclatures and practical concerns); (ii) overview of worldwide population (frequency) data; (iii) the use of X-chromosomal markers in (complex) kinship testing and the forensic statistical evaluation of evidence; (iv) segregation and mutation studies; and (v) current weaknesses and future prospects.

Multiscale analysis of SRY-positive 46,XX testicular disorder of sex development: Presentation of nine cases.

46,XX testicular disorder of sex development (46,XX TDSD) is a relatively rare condition characterised by the presence of testicular tissue with 46,XX karyotype. The present study aims to reveal the phenotype to genotype correlation in a series of sex-determining region Y (SRY)-positive 46,XX TDSD cases. We present the clinical findings, hormone profiles and genetic test results of six patients with SRY-positive 46,XX TDSD and give the details and follow-up findings of our three of previously published patients. All patients presented common characteristics such as azoospermia, hypergonadotropic hypogonadism and an SRY gene translocated on the terminal part of the short arm of one of the X chromosomes. Mean±standard deviation (SD) height of the patients was 164.78±8.0cm. Five patients had decreased secondary sexual characteristics, and three patients had gynaecomastia with varying degrees. Five of the seven patients revealed a translocation between protein kinase X (PRKX) and inverted protein kinase Y (PRKY) genes, and the remaining two patients showed a translocation between the pseudoautosomal region 1 (PAR1) of X chromosome and the differential region of Y chromosome. X chromosome inactivation (XCI) analysis results demonstrated random and skewed XCI in 5 cases and 1 case, respectively. In brief, we delineate the phenotypic spectrum of patients with SRY-positive 46,XX TDSD and the underlying mechanisms of Xp;Yp translocations.

Array-CGH revealed gain of Yp11.2 in 49,XXXXY and gain of Xp22.33 in 48,XXYY karyotypes of two rare klinefelter variants.

Klinefelter syndrome (KS) variants often share common features with classical syndrome but some of these variants present with a distinct phenotype. The incidence of sex chromosome tetrasomy and pentasomy are very less and generally diagnosed after prepubertal age. The early diagnosis of complex and unclassified syndromes and it's correlation with genotype is necessary for personalized treatment as well as genetic counselling of the affected families. We describe clinical presentation, and genetic diagnosis of two cases of variant KS. Our first case, a 4 year old male child presented with generalized tonic-clonic seizures (GTCSs), delayed milestones and dysmorphic features while case 2, a-21 years old male who had history of seizures and delayed puberty came to our lab for genetic diagnosis. The chromosomal analysis of case 1 and 2 showed 49,XXXXY and 48,XXYY karyotype respectively. The karyotype results were confirmed with fluorescence in situ hybridization (FISH) and array-CGH analysis. The FISH results were found to be consistent with karyotype but the array-CGH results showed the extra gain of region Yp11.2 in case 1 while the extra gain of region Xp22.33 in case 2. The cases were confirmed as variant KS on the basis of additional sex chromosomes and clinical presentation of deteriorated brain development. The present study suggests that the high doses of sex chromosome linked genes including pseudoautosomal region (PAR) caused the abnormal brain development. The combination of molecular techniques should be utilized for the diagnosis of such complex cases to understand the genotype-phenotype correlation and appropriate genetic counseling.

Theoretical and empirical comparisons of expected and realized relationships for the X-chromosome

BackgroundX-chromosomal loci present different inheritance patterns compared to autosomal loci and must be modeled accordingly. Sexual chromosomes are not systematically considered in whole-genome relationship matrices although rules based on genealogical or marker information have been derived. Loci on the X-chromosome could have a significant contribution to the additive genetic variance, in particular for some traits such as those related to reproduction. Thus, accounting for the X-chromosome relationship matrix might be informative to better understand the architecture of complex traits (e.g., by estimating the variance associated to this chromosome) and to improve their genomic prediction. For such applications, previous studies have shown the benefits of combining information from genotyped and ungenotyped individuals.ResultsIn this paper, we start by presenting rules to compute a genomic relationship matrix (GRM) for the X-chromosome (GX) without making any assumption on dosage compensation, and based on coding of gene content with 0/1 for males and 0/1/2 for females. This coding adjusts naturally to previously derived pedigree-based relationships (S) for the X-chromosome. When needed, we propose to accommodate and estimate dosage compensation and genetic heterogeneity across sexes via multiple trait models. Using a Holstein dairy cattle dataset, including males and females, we then empirically illustrate that realized relationships (GX) matches expectations (S). However, GX presents high deviations from S. GX has also a lower dimensionality compared to the autosomal GRM. In particular, individuals are frequently identical along the entire chromosome. Finally, we confirm that the heritability of gene content for markers on the X-chromosome that are estimated by using S is 1, further demonstrating that S and GX can be combined. For the pseudo-autosomal region, we demonstrate that the expected relationships vary according to position because of the sex-gradient. We end by presenting the rules to construct the 'H matrix’ by combining both relationship matrices.ConclusionsThis work shows theoretically and empirically that a pedigree-based relationship matrix built with rules specifically developed for the X-chromosome (S) matches the realized GRM for the X-chromosome. Therefore, applications that combine expected relationships and genotypes for markers on the X-chromosome should use S and GX.

Using GC Content to Compare Recombination Patterns on the Sex Chromosomes and Autosomes of the Guppy,Poecilia reticulata, and Its Close Outgroup Species

Genetic and physical mapping of the guppy (Poecilia reticulata) have shown that recombination patterns differ greatly between males and females. Crossover events occur evenly across the chromosomes in females, but in male meiosis they are restricted to the tip furthest from the centromere of each chromosome, creating very high recombination rates per megabase, as in pseudoautosomal regions of mammalian sex chromosomes. We used GC content to indirectly infer recombination patterns on guppy chromosomes, based on evidence that recombination is associated with GC-biased gene conversion, so that genome regions with high recombination rates should be detectable by high GC content. We used intron sequences and third positions of codons to make comparisons between sequences that are matched, as far as possible, and are all probably under weak selection. Almost all guppy chromosomes, including the sex chromosome (LG12), have very high GC values near their assembly ends, suggesting high recombination rates due to strong crossover localization in male meiosis. Our test does not suggest that the guppy XY pair has stronger crossover localization than the autosomes, or than the homologous chromosome in the close relative, the platyfish (Xiphophorus maculatus). We therefore conclude that the guppy XY pair has not recently undergone an evolutionary change to a different recombination pattern, or reduced its crossover rate, but that the guppy evolved Y-linkage due to acquiring a male-determining factor that also conferred the male crossover pattern. We also identify the centromere ends of guppy chromosomes, which were not determined in the genome assembly.

Reference genome and transcriptome informed by the sex chromosome complement of the sample increase ability to detect sex differences in gene expression from RNA-Seq data

BackgroundHuman X and Y chromosomes share an evolutionary origin and, as a consequence, sequence similarity. We investigated whether the sequence homology between the X and Y chromosomes affects the alignment of RNA-Seq reads and estimates of differential expression. We tested the effects of using reference genomes and reference transcriptomes informed by the sex chromosome complement of the sample’s genome on the measurements of RNA-Seq abundance and sex differences in expression.ResultsThe default genome includes the entire human reference genome (GRCh38), including the entire sequence of the X and Y chromosomes. We created two sex chromosome complement informed reference genomes. One sex chromosome complement informed reference genome was used for samples that lacked a Y chromosome; for this reference genome version, we hard-masked the entire Y chromosome. For the other sex chromosome complement informed reference genome, to be used for samples with a Y chromosome, we hard-masked only the pseudoautosomal regions of the Y chromosome, because these regions are duplicated identically in the reference genome on the X chromosome. We analyzed the transcript abundance in the whole blood, brain cortex, breast, liver, and thyroid tissues from 20 genetic female (46, XX) and 20 genetic male (46, XY) samples. Each sample was aligned twice: once to the default reference genome and then independently aligned to a reference genome informed by the sex chromosome complement of the sample, repeated using two different read aligners, HISAT and STAR. We then quantified sex differences in gene expression using featureCounts to get the raw count estimates followed by Limma/Voom for normalization and differential expression. We additionally created sex chromosome complement informed transcriptome references for use in pseudo-alignment using Salmon. Transcript abundance was quantified twice for each sample: once to the default target transcripts and then independently to target transcripts informed by the sex chromosome complement of the sample.ConclusionsWe show that regardless of the choice of the read aligner, using an alignment protocol informed by the sex chromosome complement of the sample results in higher expression estimates on the pseudoautosomal regions of the X chromosome in both genetic male and genetic female samples, as well as an increased number of unique genes being called as differentially expressed between the sexes. We additionally show that using a pseudo-alignment approach informed on the sex chromosome complement of the sample eliminates Y-linked expression in female XX samples.

Telomere-to-telomere assembly of a complete human X chromosome

After two decades of improvements, the current human reference genome (GRCh38) is the most accurate and complete vertebrate genome ever produced. However, no single chromosome has been finished end to end, and hundreds of unresolved gaps persist1,2. Here we present a human genome assembly that surpasses the continuity of GRCh382, along with a gapless, telomere-to-telomere assembly of a human chromosome. This was enabled by high-coverage, ultra-long-read nanopore sequencing of the complete hydatidiform mole CHM13 genome, combined with complementary technologies for quality improvement and validation. Focusing our efforts on the human X chromosome3, we reconstructed the centromeric satellite DNA array (approximately 3.1 Mb) and closed the 29 remaining gaps in the current reference, including new sequences from the human pseudoautosomal regions and from cancer-testis ampliconic gene families (CT-X and GAGE). These sequences will be integrated into future human reference genome releases. In addition, the complete chromosome X, combined with the ultra-long nanopore data, allowed us to map methylation patterns across complex tandem repeats and satellite arrays. Our results demonstrate that finishing the entire human genome is now within reach, and the data presented here will facilitate ongoing efforts to complete the other human chromosomes.

Genetic architecture of individual variation in recombination rate on the X chromosome in cattle.

Meiotic recombination is an essential biological process that ensures proper chromosome segregation and creates genetic diversity. Individual variation in global recombination rates has been shown to be heritable in several species, and variants significantly associated with this trait have been identified. Recombination on the sex chromosome has often been ignored in these studies although this trait may be particularly interesting as it may correspond to a biological process distinct from that on autosomes. For instance, recombination in males is restricted to the pseudo-autosomal region (PAR). We herein used a large cattle pedigree with more than 100,000 genotyped animals to improve the genetic map of the X chromosome and to study the genetic architecture of individual variation in recombination rate on the sex chromosome (XRR). The length of the genetic map was 46.4 and 121.2 cM in males and females, respectively, but the recombination rate in the PAR was six times higher in males. The heritability of CO counts on the X chromosome was comparable to that of autosomes in males (0.011) but larger than that of autosomes in females (0.024). XRR was highly correlated (0.76) with global recombination rate (GRR) in females, suggesting that both traits might be governed by shared variants. In agreement, a set of eleven previously identified variants associated with GRR had correlated effects on female XRR (0.86). In males, XRR and GRR appeared to be distinct traits, although more accurate CO counts on the PAR would be valuable to confirm these results.

Cross-Species BAC Mapping Highlights Conservation of Chromosome Synteny across Dragon Lizards (Squamata: Agamidae).

Dragon lizards (Squamata: Agamidae) comprise about 520 species in six subfamilies distributed across Asia, Australasia and Africa. Only five species are known to have sex chromosomes. All of them possess ZZ/ZW sex chromosomes, which are microchromosomes in four species from the subfamily Amphibolurinae, but much larger in Phrynocephalus vlangalii from the subfamily Agaminae. In most previous studies of these sex chromosomes, the focus has been on Australian species from the subfamily Amphibolurinae, but only the sex chromosomes of the Australian central bearded dragon (Pogona vitticeps) are well-characterized cytogenetically. To determine the level of synteny of the sex chromosomes of P. vitticeps across agamid subfamilies, we performed cross-species two-colour FISH using two bacterial artificial chromosome (BAC) clones from the pseudo-autosomal regions of P. vitticeps. We mapped these two BACs across representative species from all six subfamilies as well as two species of chameleons, the sister group to agamids. We found that one of these BAC sequences is conserved in macrochromosomes and the other in microchromosomes across the agamid lineages. However, within the Amphibolurinae, there is evidence of multiple chromosomal rearrangements with one of the BACs mapping to the second-largest chromosome pair and to the microchromosomes in multiple species including the sex chromosomes of P. vitticeps. Intriguingly, no hybridization signal was observed in chameleons for either of these BACs, suggesting a likely agamid origin of these sequences. Our study shows lineage-specific evolution of sequences/syntenic blocks and successive rearrangements and reveals a complex history of sequences leading to their association with important biological processes such as the evolution of sex chromosomes and sex determination.

Inherited duplication of the pseudoautosomal region Xq28 in a subject with Gilles de la Tourette syndrome and intellectual disability: a case report

BackgroundTourette syndrome (TS) is a complex neurodevelopmental disorder (NDD) characterized by multiple chronic involuntary motor and vocal tics with onset during childhood or adolescence. Most TS patients present with additional comorbidities, typically attention deficit hyperactivity disorder (ADHD), obsessive- compulsive disorder (OCD), autism spectrum disorder (ASD) and intellectual disability (ID). Both TS and ID are genetically complex disorders that likely occur as a result of the effects of multiple genes interacting with other environmental factors. In addition to single gene mutations and chromosomal disorders, copy number variations (CNVs) are implicated across many NDDs and ID and contribute to their shared genetic etiology. Screening of CNVs using microarray-based Comparative Genomic Hybridization (aCGH) is now routinely performed in all subjects with NDD and ID.Case presentationWe report a case of a 12-year-old girl diagnosed with Gilles de la Tourette Syndrome associated to behavior disorders and intellectual disability in particular with regard to language. Array-CGH analysis showed a CNV of a subtelomeric region Xq28 (gain of 260 kb) inherited from the healthy father. The duplication contains two genes, VAMP7 and SPRY3 of the PAR2 pseudoautosomal region. FISH analysis revealed that the duplicated segment is located on the short arm of a chromosome 13, resulting in a trisomy of the region. In the proband the expression levels of the genes evaluated in the peripheral blood sample are comparable both those of the mother and to those of female control subjects.ConclusionsAlthough the trisomy of the 260 kb region from Xq28 identified in proband is also shared by the healthy father, it is tantalizing to speculate that, together with genetic risk factors inherited from the mother, it may play a role in the development of a form of Tourette syndrome with intellectual disability. This hypothesis is also supported by the fact that both genes present in the duplicated region (VAMP7 and SPRY3) are expressed in the CNS and are implicated in neurotransmission and neurite growth and branching. In addition, similar CNVs have been identified in individuals whose phenotype is associated with autism spectrum disorders or intellectual disability.