The structural maintenance of chromosomes (SMC) proteins constitute the core of critical complexes involved in structural organization of chromosomes. In yeast, the Smc5/6 complex is known to mediate repair of DNA breaks and replication of repetitive genomic regions, including ribosomal DNA loci and telomeres. In mammalian cells, which have diverse genome structure and scale from yeast, the Smc5/6 complex has also been implicated in DNA damage response, but its further function in unchallenged conditions remains elusive. In this study, we addressed the behavior and function of Smc5/6 during the cell cycle. Chromatin fractionation, immunofluorescence, and live-cell imaging analyses indicated that Smc5/6 associates with chromatin during interphase but largely dissociates from chromosomes when they condense in mitosis. Depletion of Smc5 and Smc6 resulted in aberrant mitotic chromosome phenotypes that were accompanied by the abnormal distribution of topoisomerase IIα (topo IIα) and condensins and by chromosome segregation errors. Importantly, interphase chromatin structure indicated by the premature chromosome condensation assay suggested that Smc5/6 is required for the on-time progression of DNA replication and subsequent binding of topo IIα on replicated chromatids. These results indicate an essential role of the Smc5/6 complex in processing DNA replication, which becomes indispensable for proper sister chromatid assembly in mitosis.

Plasmodium species causing malaria in humans are not monophyletic, sharing common ancestors with nonhuman primate parasites. Plasmodium gonderi is one of the few known Plasmodium species infecting African old-world monkeys that are not found in apes. This study reports a de novo assembled P. gonderi genome with complete chromosomes. The P. gonderi genome shares codon usage, syntenic blocks, and other characteristics with the human parasites Plasmodium ovale s.l. and Plasmodium malariae, also of African origin, and the human parasite Plasmodium vivax and species found in nonhuman primates from Southeast Asia. Using phylogenetically aware methods, newly identified syntenic blocks were found enriched with conserved metabolic genes. Regions outside those blocks harbored genes encoding proteins involved in the vertebrate host-Plasmodium relationship undergoing faster evolution. Such genome architecture may have facilitated colonizing vertebrate hosts. Phylogenomic analyses estimated the common ancestor between P. vivax and an African ape parasite P. vivax-like, within the Asian nonhuman primates parasites clade. Time estimates incorporating P. gonderi placed the P. vivax and P. vivax-like common ancestor in the late Pleistocene, a time of active migration of hominids between Africa and Asia. Thus, phylogenomic and time-tree analyses are consistent with an Asian origin for P. vivax and an introduction of P. vivax-like into Africa. Unlike other studies, time estimates for the clade with Plasmodium falciparum, the most lethal human malaria parasite, coincide with their host species radiation, African hominids. Overall, the newly assembled genome presented here has the quality to support comparative genomic investigations in Plasmodium.

Chromatin assembled with centromere protein A (CENP-A) is the epigenetic mark of centromere identity. Using new reference models, we now identify sites of CENP-A and histone H3.1 binding within the megabase, α-satellite repeat-containing centromeres of 23 human chromosomes. The overwhelming majority (97%) of α-satellite DNA is found to be assembled with histone H3.1-containing nucleosomes with wrapped DNA termini. In both G1 and G2 cell cycle phases, the 2-4% of α-satellite assembled with CENP-A protects DNA lengths centered on 133 bp, consistent with octameric nucleosomes with DNA unwrapping at entry and exit. CENP-A chromatin is shown to contain equimolar amounts of CENP-A and histones H2A, H2B, and H4, with no H3. Solid-state nanopore analyses show it to be nucleosomal in size. Thus, in contrast to models for hemisomes that briefly transition to octameric nucleosomes at specific cell cycle points or heterotypic nucleosomes containing both CENP-A and histone H3, human CENP-A chromatin complexes are octameric nucleosomes with two molecules of CENP-A at all cell cycle phases.

Human DNA present in faecal samples can result in a small number of human reads in gut shotgun metagenomic sequencing data. However, it is presently unclear how much personal information can be reconstructed from such reads, and this has not been quantitatively evaluated. Such a quantitative evaluation is necessary to clarify the ethical concerns related to data sharing and to enable efficient use of human genetic information in stool samples, such as for research and forensics. Here we used genomic approaches to reconstruct personal information from the faecal metagenomes of 343 Japanese individuals with associated human genotype data. Genetic sex could be accurately predicted based on the sequencing depth of sex chromosomes for 97.3% of the samples. Individuals could be re-identified from the matched genotype data based on human reads recovered from the faecal metagenomic data with 93.3% sensitivity using a likelihood score-based method. This method also enabled us to predict the ancestries of 98.3% of the samples. Finally, we performed ultra-deep shotgun metagenomic sequencing of five faecal samples as well as whole-genome sequencing of blood samples. Using genotype-calling approaches, we demonstrated that the genotypes of both common and rare variants could be reconstructed from faecal samples. This included clinically relevant variants. Our approach can be used to quantify personal information contained within gut metagenome data.

Common variable immunodeficiency disorder (CVID) is the most common symptomatic primary immunodeficiency in adults, characterized by B-cell abnormalities and inadequate antibody response. CVID patients have considerable autoimmune comorbidity and we therefore hypothesized that genetic susceptibility to CVID may overlap with autoimmune disorders. Here, in the largest genetic study performed in CVID to date, we compare 778 CVID cases with 10,999 controls across 123,127 single-nucleotide polymorphisms (SNPs) on the Immunochip. We identify the first non-HLA genome-wide significant risk locus at CLEC16A (rs17806056, P=2.0 × 10(-9)) and confirm the previously reported human leukocyte antigen (HLA) associations on chromosome 6p21 (rs1049225, P=4.8 × 10(-16)). Clec16a knockdown (KD) mice showed reduced number of B cells and elevated IgM levels compared with controls, suggesting that CLEC16A may be involved in immune regulatory pathways of relevance to CVID. In conclusion, the CLEC16A associations in CVID represent the first robust evidence of non-HLA associations in this immunodeficiency condition.

Rhabdomyosarcoma is a relatively common tumour of the soft tissue, probably due to regulatory disruption of growth and differentiation of skeletal muscle stem cells. Identification of genes differentially expressed in normal skeletal muscle and in rhabdomyosarcoma may help in understanding mechanisms of tumour development, in discovering diagnostic and prognostic markers and in identifying novel targets for drug therapy.

Humans live in symbiosis with a diverse community of microorganisms, which has evolved to carry out many specific tasks that benefit the host, including protection against invading pathogens. Within the chemical diversity of the gastrointestinal tract, small molecules likely constitute chemical cues for the communication between the microbiota and pathogens. Therefore, we sought to investigate if molecules produced by the human gut microbiota show biological activity against the human pathogen Vibrio cholerae. To probe the effects of the gut metabolome on V. cholerae, we investigated its response to small-molecule extracts from human feces, from a complex bacterial community cultivated in vitro, and from culture supernatants of Enterocloster citroniae, Bacteroides thetaiotaomicron, and Bacteroides vulgatus. Using RNA sequencing, we determined the impact of the human gut metabolome on V. cholerae global gene expression. Among the genes downregulated in the presence of the fecal extract, the most overrepresented functional category was cell motility, which accounted for 39% of repressed genes. Repression of V. cholerae motility by the fecal extract was confirmed phenotypically, and E. citroniae extracts reproduced this phenotype. A complex in vitro microbial community led to increased motility, as did extracts from B. vulgatus, a species present in this community. Accordingly, mucin penetration was also repressed by fecal and E. citroniae extracts, suggesting that the phenotypes observed may have implications for host colonization. Together with previous studies, this work shows that small molecules from the gut metabolome may have a widespread, significant impact on microbe-microbe interactions established in the gut environment.

Sex differences in endothelial cell (EC) biology may reflect intrinsic differences driven by chromosomes or sex steroid exposure and gender differences accumulated over life. We analysed EC gene expression data from boy-girl twins at birth and in non-twin adults to detect sex differences at different stages of life, and show that 14-25% of the EC transcriptome is sex-biased. By combining data from both stages of life, we identified sex differences that are present at birth and maintained throughout life, and those that are acquired over life. Promisingly, we found that genes that present with an acquired sex difference in ECs are more likely to be targets of sex steroids. Annotating both gene sets with data from multiple genome-wide association studies (GWAS) revealed that genes with an intrinsic sex difference in ECs are enriched for coronary artery disease GWAS hits. This study underscores the need for treating sex as a biological variable.

In recent years different types of structural variants (SVs) have been discovered in the human genome and their functional impact has become increasingly clear. Inversions, however, are poorly characterized and more difficult to study, especially those mediated by inverted repeats or segmental duplications. Here, we describe the results of a simple and fast inverse PCR (iPCR) protocol for high-throughput genotyping of a wide variety of inversions using a small amount of DNA. In particular, we analyzed 22 inversions predicted in humans ranging from 5.1 kb to 226 kb and mediated by inverted repeat sequences of 1.6-24 kb. First, we validated 17 of the 22 inversions in a panel of nine HapMap individuals from different populations, and we genotyped them in 68 additional individuals of European origin, with correct genetic transmission in ∼ 12 mother-father-child trios. Global inversion minor allele frequency varied between 1% and 49% and inversion genotypes were consistent with Hardy-Weinberg equilibrium. By analyzing the nucleotide variation and the haplotypes in these regions, we found that only four inversions have linked tag-SNPs and that in many cases there are multiple shared SNPs between standard and inverted chromosomes, suggesting an unexpected high degree of inversion recurrence during human evolution. iPCR was also used to check 16 of these inversions in four chimpanzees and two gorillas, and 10 showed both orientations either within or between species, providing additional support for their multiple origin. Finally, we have identified several inversions that include genes in the inverted or breakpoint regions, and at least one disrupts a potential coding gene. Thus, these results represent a significant advance in our understanding of inversion polymorphism in human populations and challenge the common view of a single origin of inversions, with important implications for inversion analysis in SNP-based studies.

The first meiotic division in human oocytes is highly error-prone and contributes to the uniquely high incidence of aneuploidy observed in human pregnancies. A successful meiosis I (MI) division entails separation of homologous chromosome pairs and co-segregation of sister chromatids. For this to happen, sister kinetochores must form attachments to spindle kinetochore-fibres emanating from the same pole. In mouse and budding yeast, sister kinetochores remain closely associated with each other during MI, enabling them to act as a single unified structure. However, whether this arrangement also applies in human meiosis I oocytes was unclear. In this study, we perform high-resolution imaging of over 1900 kinetochores in human oocytes, to examine the geometry and architecture of the human meiotic kinetochore. We reveal that sister kinetochores in MI are not physically fused, and instead individual kinetochores within a pair are capable of forming independent attachments to spindle k-fibres. Notably, with increasing female age, the separation between kinetochores increases, suggesting a degradation of centromeric cohesion and/or changes in kinetochore architecture. Our data suggest that the differential arrangement of sister kinetochores and dual k-fibre attachments may explain the high proportion of unstable attachments that form in MI and thus indicate why human oocytes are prone to aneuploidy, particularly with increasing maternal age.

Aware of the important benefits of human milk, most U.S. women initiate breastfeeding but difficulties with milk supply lead some to quit earlier than intended. Yet, the contribution of maternal physiology to lactation difficulties remains poorly understood. Human milk fat globules, by enveloping cell contents during their secretion into milk, are a rich source of mammary cell RNA. Here, we pair this non-invasive mRNA source with RNA-sequencing to probe the milk fat layer transcriptome during three stages of lactation: colostral, transitional, and mature milk production. The resulting transcriptomes paint an exquisite portrait of human lactation. The resulting transcriptional profiles cluster not by postpartum day, but by milk Na:K ratio, indicating that women sampled during similar postpartum time frames could be at markedly different stages of gene expression. Each stage of lactation is characterized by a dynamic range (10(5)-fold) in transcript abundances not previously observed with microarray technology. We discovered that transcripts for isoferritins and cathepsins are strikingly abundant during colostrum production, highlighting the potential importance of these proteins for neonatal health. Two transcripts, encoding β-casein (CSN2) and α-lactalbumin (LALBA), make up 45% of the total pool of mRNA in mature lactation. Genes significantly expressed across all stages of lactation are associated with making, modifying, transporting, and packaging milk proteins. Stage-specific transcripts are associated with immune defense during the colostral stage, up-regulation of the machinery needed for milk protein synthesis during the transitional stage, and the production of lipids during mature lactation. We observed strong modulation of key genes involved in lactose synthesis and insulin signaling. In particular, protein tyrosine phosphatase, receptor type, F (PTPRF) may serve as a biomarker linking insulin resistance with insufficient milk supply. This study provides the methodology and reference data set to enable future targeted research on the physiological contributors of sub-optimal lactation in humans.

Studies of the Y chromosome in primates, rodents and carnivores provide compelling evidence that the male specific region of Y (MSY) contains functional genes, many of which have specialized roles in spermatogenesis and male-fertility. Little similarity, however, has been found between the gene content and sequence of MSY in different species. This hinders the discovery of species-specific male fertility genes and limits our understanding about MSY evolution in mammals. Here, a detailed MSY gene catalogue was developed for the horse--an odd-toed ungulate. Using direct cDNA selection from horse testis, and sequence analysis of Y-specific BAC clones, 37 horse MSY genes/transcripts were identified. The genes were mapped to the MSY BAC contig map, characterized for copy number, analyzed for transcriptional profiles by RT-PCR, examined for the presence of ORFs, and compared to other mammalian orthologs. We demonstrate that the horse MSY harbors 20 X-degenerate genes with known orthologs in other eutherian species. The remaining 17 genes are acquired or novel and have so far been identified only in the horse or donkey Y chromosomes. Notably, 3 transcripts were found in the heterochromatic part of the Y. We show that despite substantial differences between the sequence, gene content and organization of horse and other mammalian Y chromosomes, the functions of MSY genes are predominantly related to testis and spermatogenesis. Altogether, 10 multicopy genes with testis-specific expression were identified in the horse MSY, and considered likely candidate genes for stallion fertility. The findings establish an important foundation for the study of Y-linked genetic factors governing fertility in stallions, and improve our knowledge about the evolutionary processes that have shaped Y chromosomes in different mammalian lineages.

Colourless microalgae of the Prototheca genus are the only known plants that have consistently been implicated in a range of clinically relevant opportunistic infections in both animals and humans. The Prototheca algae are emerging pathogens, whose incidence has increased importantly over the past two decades. Prototheca wickerhamii is a major human pathogen, responsible for at least 115 cases worldwide. Although the algae are receiving more attention nowadays, there is still a substantial knowledge gap regarding their biology, and pathogenicity in particular. Here we report, for the first time, the complete nuclear genome, organelle genomes, and transcriptome of the P. wickerhamii type strain ATCC 16529.

Telomeres are specialized nucleoprotein complexes, localized at the physical ends of chromosomes, that contribute to the maintenance of genome stability. One of the features of chronic myeloid leukemia (CML) cells is a reduction in telomere length which may result in increased genomic instability and progression of the disease. Aberrant telomere maintenance in CML is not fully understood and other mechanisms such as the alternative lengthening of telomeres (ALT) are involved. In this work, we employed five BCR-ABL1-positive cell lines, namely K562, KU-812, LAMA-84, MEG-A2, and MOLM-1, commonly used in the laboratories to study the link between mutation, copy number, and expression of telomere maintenance genes with the expression, copy number, and activity of BCR-ABL1. Our results demonstrated that the copy number and expression of BCR-ABL1 are crucial for telomere lengthening. We observed a correlation between BCR-ABL1 expression and telomere length as well as shelterins upregulation. Next-generation sequencing revealed pathogenic variants and copy number alterations in major tumor suppressors, such as TP53 and CDKN2A, but not in telomere-associated genes. Taken together, we showed that BCR-ABL1 kinase expression and activity play a crucial role in the maintenance of telomeres in CML cell lines. Our results may help to validate and properly interpret results obtained by many laboratories employing these in vitro models of CML.

It was demonstrated previously that the three-locus RFLP haplotype, TaqI B-TaqI D-TaqI A (B-D-A), at the DRD2 locus constitutes a powerful genetic marker and probably reflects the most ancient dispersal of anatomically modern humans.

Sex chromosomes are generally derived from a pair of classical type-A chromosomes, and relatively few alternative models have been proposed up to now.1,2 B chromosomes (Bs) are supernumerary and dispensable chromosomes with non-Mendelian inheritance found in many plant and animal species3,4 that have often been considered as selfish genetic elements that behave as genome parasites.5,6 The observation that in some species Bs can be either restricted or predominant in one sex7-14 raised the interesting hypothesis that Bs could play a role in sex determination.15 The characterization of putative B master sex-determining (MSD) genes, however, has not yet been provided to support this hypothesis. Here, in Astyanax mexicanus cavefish originating from Pachón cave, we show that Bs are strongly male predominant. Based on a high-quality genome assembly of a B-carrying male, we characterized the Pachón cavefish B sequence and found that it contains two duplicated loci of the putative MSD gene growth differentiation factor 6b (gdf6b). Supporting its role as an MSD gene, we found that the Pachón cavefish gdf6b gene is expressed specifically in differentiating male gonads, and that its knockout induces male-to-female sex reversal in B-carrying males. This demonstrates that gdf6b is necessary for triggering male sex determination in Pachón cavefish. Altogether these results bring multiple and independent lines of evidence supporting the conclusion that the Pachón cavefish B is a "B-sex" chromosome that contains duplicated copies of the gdf6b gene, which can promote male sex determination in this species.

Mammalian sex-chromosomes originated from a pair of autosomes. A step-wise cessation of recombination is necessary for the proper maintenance of sex-determination and, consequently, generates a four strata structure on the X chromosome. Each stratum shows a specific per-site nucleotide sequence difference (p-distance) between the X and Y chromosomes, depending on the time of recombination arrest. Stratum 4 covers the distal half of the human X chromosome short arm and the p-distance of the stratum is approximately 10%, on average. However, a 100-kb region, which includes KALX and VCX, in the middle of stratum 4 shows a significantly lower p-distance (1-5%), suggesting frequent sequence exchanges or gene conversions between the X and Y chromosomes in humans. To examine the evolutionary mechanism for this low p-distance region, sequences of a corresponding region including KALX/Y from seven species of non-human primates were analyzed.

We report the analysis of 335 microsatellite markers genotyped in 110 multiplex families with autism. All families include at least two "affected" siblings, at least one of whom has autism; the remaining affected sibs carry diagnoses of either Asperger syndrome or pervasive developmental disorder. Affected sib-pair analysis yielded multipoint maximum LOD scores (MLS) that reach the accepted threshold for suggestive linkage on chromosomes 5, X, and 19. Nominal evidence for linkage (point-wise P<.05) was obtained on chromosomes 2, 3, 4, 8, 10, 11, 12, 15, 16, 18, and 20, and secondary loci were found on chromosomes 5 and 19. Analysis of families sharing alleles at the putative X chromosomal linked locus and one or more other putative linked loci produced an MLS of 3.56 for the DXS470-D19S174 marker combination. In an effort to increase power to detect linkage, scan statistics were used to evaluate the significance of peak LOD scores based on statistical evidence at adjacent marker loci. This analysis yielded impressive evidence for linkage to autism and autism-spectrum disorders with significant genomewide P values <.05 for markers on chromosomes 5 and 8 and with suggestive linkage evidence for a marker on chromosome 19.

The generalization of the second Chargaff rule states that counts of any string of nucleotides of length k on a single chromosomal strand equal the counts of its inverse (reverse-complement) k-mer. This Inversion Symmetry (IS) holds for many species, both eukaryotes and prokaryotes, for ranges of k which may vary from 7 to 10 as chromosomal lengths vary from 2Mbp to 200 Mbp. The existence of IS has been demonstrated in the literature, and other pair-wise candidate symmetries (e.g. reverse or complement) have been ruled out.

Balanced chromosome abnormalities (BCAs) occur at a high frequency in healthy and diseased individuals, but cost-efficient strategies to identify BCAs and evaluate whether they contribute to a phenotype have not yet become widespread. Here we apply genome-wide mate-pair library sequencing to characterize structural variation in a patient with unclear neurodevelopmental disease (NDD) and complex de novo BCAs at the karyotype level. Nucleotide-level characterization of the clinically described BCA breakpoints revealed disruption of at least three NDD candidate genes (LINC00299, NUP205, PSMD14) that gave rise to abnormal mRNAs and could be assumed as disease-causing. However, unbiased genome-wide analysis of the sequencing data for cryptic structural variation was key to reveal an additional submicroscopic inversion that truncates the schizophrenia- and bipolar disorder-associated brain transcription factor ZNF804A as an equally likely NDD-driving gene. Deep sequencing of fluorescent-sorted wild-type and derivative chromosomes confirmed the clinically undetected BCA. Moreover, deep sequencing further validated a high accuracy of mate-pair library sequencing to detect structural variants larger than 10 kB, proposing that this approach is powerful for clinical-grade genome-wide structural variant detection. Our study supports previous evidence for a role of ZNF804A in NDD and highlights the need for a more comprehensive assessment of structural variation in karyotypically abnormal individuals and patients with neurocognitive disease to avoid diagnostic deception.