Hominoid chromosomal rearrangements on 17q map to complex regions of segmental duplication

Transcriptome innovations in primates revealed by single-molecule long-read sequencing

Transcriptomic diversity greatly contributes to the fundamentals of disease, lineage-specific biology, and environmental adaptation. However, much of the actual isoform repertoire contributing to shaping primate evolution remains unknown. Here, we combined deep long- and short-read sequencing complemented with mass spectrometry proteomics in a panel of lymphoblastoid cell lines (LCLs) from human, three other great apes, and rhesus macaque, producing the largest full-length isoform catalog in primates to date. Around half of the captured isoforms are not annotated in their reference genomes, significantly expanding the gene models in primates. Furthermore, our comparative analyses unveil hundreds of transcriptomic innovations and isoform usage changes related to immune function and immunological disorders. The confluence of these evolutionary innovations with signals of positive selection and their limited impact in the proteome points to changes in alternative splicing in genes involved in immune response as an important target of recent regulatory divergence in primates.

AC2: An Efficient Protein Sequence Compression Tool Using Artificial Neural Networks and Cache-Hash Models

Recently, the scientific community has witnessed a substantial increase in the generation of protein sequence data, triggering emergent challenges of increasing importance, namely efficient storage and improved data analysis. For both applications, data compression is a straightforward solution. However, in the literature, the number of specific protein sequence compressors is relatively low. Moreover, these specialized compressors marginally improve the compression ratio over the best general-purpose compressors. In this paper, we present AC2, a new lossless data compressor for protein (or amino acid) sequences. AC2 uses a neural network to mix experts with a stacked generalization approach and individual cache-hash memory models to the highest-context orders. Compared to the previous compressor (AC), we show gains of 2–9% and 6–7% in reference-free and reference-based modes, respectively. These gains come at the cost of three times slower computations. AC2 also improves memory usage against AC, with requirements about seven times lower, without being affected by the sequences’ input size. As an analysis application, we use AC2 to measure the similarity between each SARS-CoV-2 protein sequence with each viral protein sequence from the whole UniProt database. The results consistently show higher similarity to the pangolin coronavirus, followed by the bat and human coronaviruses, contributing with critical results to a current controversial subject. AC2 is available for free download under GPLv3 license.

Multicolor-FISH Characterization of a Prenatal Mosaicism for a Chromosomal Rearrangement Undetected by Molecular Cytogenetics

Fetal mosaicism for chromosomal rearrangements remains a challenge to diagnose, even in the era of whole-genome sequencing. We present here a case of fetal mosaicism for a chromosomal rearrangement explored in amniocytes and fetal muscle, consisting of a major cell population (95%) with partial monosomy 4q and a minor population (5%) with additional material replacing the 4qter deleted segment. Molecular techniques (MLPA, array-CGH) failed to assess the origin of this material. Only multicolor-FISH identified the additional segment on chromosome 4 as derived from chromosome 17. Due to the poor prognosis, the couple chose to terminate the pregnancy. Because of low-level mosaicism, chromosomal microarray analysis (CMA), now considered as first-tier prenatal genetic analysis, did not allow the identification of the minor cell line. In case of large CNVs (>5 Mb) detected by CMA, karyotyping may be considered to elucidate the mechanism of the underlying rearrangement and eliminate mosaicism.

Single-cell strand sequencing of a macaque genome reveals multiple nested inversions and breakpoint reuse during primate evolution

Rhesus macaque is an Old World monkey that shared a common ancestor with human ∼25 Myr ago and is an important animal model for human disease studies. A deep understanding of its genetics is therefore required for both biomedical and evolutionary studies. Among structural variants, inversions represent a driving force in speciation and play an important role in disease predisposition. Here we generated a genome-wide map of inversions between human and macaque, combining single-cell strand sequencing with cytogenetics. We identified 375 total inversions between 859 bp and 92 Mbp, increasing by eightfold the number of previously reported inversions. Among these, 19 inversions flanked by segmental duplications overlap with recurrent copy number variants associated with neurocognitive disorders. Evolutionary analyses show that in 17 out of 19 cases, the Hominidae orientation of these disease-associated regions is always derived. This suggests that duplicated sequences likely played a fundamental role in generating inversions in humans and great apes, creating architectures that nowadays predispose these regions to disease-associated genetic instability. Finally, we identified 861 genes mapping at 156 inversions breakpoints, with some showing evidence of differential expression in human and macaque cell lines, thus highlighting candidates that might have contributed to the evolution of species-specific features. This study depicts the most accurate fine-scale map of inversions between human and macaque using a two-pronged integrative approach, such as single-cell strand sequencing and cytogenetics, and represents a valuable resource toward understanding of the biology and evolution of primate species.

An evolutionary driver of interspersed segmental duplications in primates

BACKGROUND

The complex interspersed pattern of segmental duplications in humans is responsible for rearrangements associated with neurodevelopmental disease, including the emergence of novel genes important in human brain evolution. We investigate the evolution of LCR16a, a putative driver of this phenomenon that encodes one of the most rapidly evolving human-ape gene families, nuclear pore interacting protein (NPIP).

RESULTS

Comparative analysis shows that LCR16a has independently expanded in five primate lineages over the last 35 million years of primate evolution. The expansions are associated with independent lineage-specific segmental duplications flanking LCR16a leading to the emergence of large interspersed duplication blocks at non-orthologous chromosomal locations in each primate lineage. The intron-exon structure of the NPIP gene family has changed dramatically throughout primate evolution with different branches showing characteristic gene models yet maintaining an open reading frame. In the African ape lineage, we detect signatures of positive selection that occurred after a transition to more ubiquitous expression among great ape tissues when compared to Old World and New World monkeys. Mouse transgenic experiments from baboon and human genomic loci confirm these expression differences and suggest that the broader ape expression pattern arose due to mutational changes that emerged in cis.

CONCLUSIONS

LCR16a promotes serial interspersed duplications and creates hotspots of genomic instability that appear to be an ancient property of primate genomes. Dramatic changes to NPIP gene structure and altered tissue expression preceded major bouts of positive selection in the African ape lineage, suggestive of a gene undergoing strong adaptive evolution.

High-resolution comparative analysis of great ape genomes

Genetic studies of human evolution require high-quality contiguous ape genome assemblies that are not guided by the human reference. We coupled long-read sequence assembly and full-length complementary DNA sequencing with a multiplatform scaffolding approach to produce ab initio chimpanzee and orangutan genome assemblies. By comparing these with two long-read de novo human genome assemblies and a gorilla genome assembly, we characterized lineage-specific and shared great ape genetic variation ranging from single- to mega-base pair-sized variants. We identified ~17,000 fixed human-specific structural variants identifying genic and putative regulatory changes that have emerged in humans since divergence from nonhuman apes. Interestingly, these variants are enriched near genes that are down-regulated in human compared to chimpanzee cerebral organoids, particularly in cells analogous to radial glial neural progenitors.

Inversion variants in human and primate genomes

For many years, inversions have been proposed to be a direct driving force in speciation since they suppress recombination when heterozygous. Inversions are the most common large-scale differences among humans and great apes. Nevertheless, they represent large events easily distinguishable by classical cytogenetics, whose resolution, however, is limited. Here, we performed a genome-wide comparison between human, great ape, and macaque genomes using the net alignments for the most recent releases of genome assemblies. We identified a total of 156 putative inversions, between 103 kb and 91 Mb, corresponding to 136 human loci. Combining literature, sequence, and experimental analyses, we analyzed 109 of these loci and found 67 regions inverted in one or multiple primates, including 28 newly identified inversions. These events overlap with 81 human genes at their breakpoints, and seven correspond to sites of recurrent rearrangements associated with human disease. This work doubles the number of validated primate inversions larger than 100 kb, beyond what was previously documented. We identified 74 sites of errors, where the sequence has been assembled in the wrong orientation, in the reference genomes analyzed. Our data serve two purposes: First, we generated a map of evolutionary inversions in these genomes representing a resource for interrogating differences among these species at a functional level; second, we provide a list of misassembled regions in these primate genomes, involving over 300 Mb of DNA and 1978 human genes. Accurately annotating these regions in the genome references has immediate applications for evolutionary and biomedical studies on primates.

Characterizing polymorphic inversions in human genomes by single-cell sequencing

Identifying genomic features that differ between individuals and cells can help uncover the functional variants that drive phenotypes and disease susceptibilities. For this, single-cell studies are paramount, as it becomes increasingly clear that the contribution of rare but functional cellular subpopulations is important for disease prognosis, management, and progression. Until now, studying these associations has been challenged by our inability to map structural rearrangements accurately and comprehensively. To overcome this, we coupled single-cell sequencing of DNA template strands (Strand-seq) with custom analysis software to rapidly discover, map, and genotype genomic rearrangements at high resolution. This allowed us to explore the distribution and frequency of inversions in a heterogeneous cell population, identify several polymorphic domains in complex regions of the genome, and locate rare alleles in the reference assembly. We then mapped the entire genomic complement of inversions within two unrelated individuals to characterize their distinct inversion profiles and built a nonredundant global reference of structural rearrangements in the human genome. The work described here provides a powerful new framework to study structural variation and genomic heterogeneity in single-cell samples, whether from individuals for population studies or tissue types for biomarker discovery.

Monkey-based research on human disease: the implications of genetic differences

Assertions that the use of monkeys to investigate human diseases is valid scientifically are frequently based on a reported 90-93% genetic similarity between the species. Critical analyses of the relevance of monkey studies to human biology, however, indicate that this genetic similarity does not result in sufficient physiological similarity for monkeys to constitute good models for research, and that monkey data do not translate well to progress in clinical practice for humans. Salient examples include the failure of new drugs in clinical trials, the highly different infectivity and pathology of SIV/HIV, and poor extrapolation of research on Alzheimer's disease, Parkinson's disease and stroke. The major molecular differences underlying these inter-species phenotypic disparities have been revealed by comparative genomics and molecular biology - there are key differences in all aspects of gene expression and protein function, from chromosome and chromatin structure to post-translational modification. The collective effects of these differences are striking, extensive and widespread, and they show that the superficial similarity between human and monkey genetic sequences is of little benefit for biomedical research. The extrapolation of biomedical data from monkeys to humans is therefore highly unreliable, and the use of monkeys must be considered of questionable value, particularly given the breadth and potential of alternative methods of enquiry that are currently available to scientists.

Chromosomal polymorphism in mammals: an evolutionary perspective

Although chromosome rearrangements (CRs) are central to studies of genome evolution, our understanding of the evolutionary consequences of the early stages of karyotypic differentiation (i.e. polymorphism), especially the non-meiotic impacts, is surprisingly limited. We review the available data on chromosomal polymorphisms in mammals so as to identify taxa that hold promise for developing a more comprehensive understanding of chromosomal change. In doing so, we address several key questions: (i) to what extent are mammalian karyotypes polymorphic, and what types of rearrangements are principally involved? (ii) Are some mammalian lineages more prone to chromosomal polymorphism than others? More specifically, do (karyotypically) polymorphic mammalian species belong to lineages that are also characterized by past, extensive karyotype repatterning? (iii) How long can chromosomal polymorphisms persist in mammals? We discuss the evolutionary implications of these questions and propose several research avenues that may shed light on the role of chromosome change in the diversification of mammalian populations and species.

Comparison of a quantitative Real-Time PCR assay and droplet digital PCR for copy number analysis of the CCL4L genes

The controversy surrounding the findings that copy number variation, of the CCL3 encoding genes, influences HIV-1 infection and disease progression has been in part attributed to the variable results obtained from methods used for copy number evaluation. Like CCL3, the genes encoding the CC chemokine CCL4, also a natural ligand of the CCR5 receptor, are found to occur in population-specific multiple copy number and have been shown to play a protective role against HIV-1. This study evaluated the standard method of quantitative Real-Time PCR (qPCR) and droplet digital PCR (ddPCR) for CCL4L gene copy number determination. The CCL4 encoding genes are CCL4, occurring in two copies per diploid genome (pdg), and the non-allelic CCL4L genes, comprised of CCL4L1 and CCL4L2, which are both found in multiple copies pdg. Copy number of CCL4L, CCL4L1 and CCL4L2 was determined in a cohort of HIV-1-uninfected individuals from the South African Black (n=23) and Caucasian (n=32) population groups using qPCR and ddPCR. A stronger correlation between the number of CCL4L copies and the sum of CCL4L1 and CCL4L2 copies generated by ddPCR (r=0.99, p<0.0001) compared to qPCR (r=0.87, p<0.0001) was observed. Real-Time qPCR exhibited greater inaccuracy at higher copy numbers which is particularly relevant to our cohort of Black individuals who have a higher range of CCL4L copies (3-6) compared to Caucasians (0-4) and a higher population median (4 and 2, respectively). Medians and ranges of CCL4L1 (Black: 2, 0-4, Caucasian: 0, 0-2) and CCL4L2 (Black: 2, 1-5, Caucasian: 2, 0-3) were also higher in the Black population. Droplet digital PCR was shown to be a far superior method to qPCR for assessment of CCL4 gene copy number variation, the accuracy of which is essential for studies of the contribution of variable gene copy number to phenotypic outcomes of host infection and disease course.

Copy number variation in chemokine superfamily: the complex scene of CCL3L-CCL4L genes in health and disease

Genome copy number changes (copy number variations: CNVs) include inherited, de novo and somatically acquired deviations from a diploid state within a particular chromosomal segment. CNVs are frequent in higher eukaryotes and associated with a substantial portion of inherited and acquired risk for various human diseases. CNVs are distributed widely in the genomes of apparently healthy individuals and thus constitute significant amounts of population-based genomic variation. Human CNV loci are enriched for immune genes and one of the most striking examples of CNV in humans involves a genomic region containing the chemokine genes CCL3L and CCL4L. The CCL3L-CCL4L copy number variable region (CNVR) shows extensive architectural complexity, with smaller CNVs within the larger ones and with interindividual variation in breakpoints. Furthermore, the individual genes embedded in this CNVR account for an additional level of genetic and mRNA complexity: CCL4L1 and CCL4L2 have identical exonic sequences but produce a different pattern of mRNAs. CCL3L2 was considered previously as a CCL3L1 pseudogene, but is actually transcribed. Since 2005, CCL3L-CCL4L CNV has been associated extensively with various human immunodeficiency virus-related outcomes, but some recent studies called these associations into question. This controversy may be due in part to the differences in alternative methods for quantifying gene copy number and differentiating the individual genes. This review summarizes and discusses the current knowledge about CCL3L-CCL4L CNV and points out that elucidating their complete phenotypic impact requires dissecting the combinatorial genomic complexity posed by various proportions of distinct CCL3L and CCL4L genes among individuals.

Identification of a recurrent microdeletion at 17q23.1q23.2 flanked by segmental duplications associated with heart defects and limb abnormalities

Segmental duplications, which comprise approximately 5%-10% of the human genome, are known to mediate medically relevant deletions, duplications, and inversions through nonallelic homologous recombination (NAHR) and have been suggested to be hot spots in chromosome evolution and human genomic instability. We report seven individuals with microdeletions at 17q23.1q23.2, identified by microarray-based comparative genomic hybridization (aCGH). Six of the seven deletions are approximately 2.2 Mb in size and flanked by large segmental duplications of >98% sequence identity and in the same orientation. One of the deletions is approximately 2.8 Mb in size and is flanked on the distal side by a segmental duplication, whereas the proximal breakpoint falls between segmental duplications. These characteristics suggest that NAHR mediated six out of seven of these rearrangements. These individuals have common features, including mild to moderate developmental delay (particularly speech delay), microcephaly, postnatal growth retardation, heart defects, and hand, foot, and limb abnormalities. Although all individuals had at least mild dysmorphic facial features, there was no characteristic constellation of features that would elicit clinical suspicion of a specific disorder. The identification of common clinical features suggests that microdeletions at 17q23.1q23.2 constitute a novel syndrome. Furthermore, the inclusion in the minimal deletion region of TBX2 and TBX4, transcription factors belonging to a family of genes implicated in a variety of developmental pathways including those of heart and limb, suggests that these genes may play an important role in the phenotype of this emerging syndrome.

New insights into centromere organization and evolution from the white-cheeked gibbon and marmoset

The evolutionary history of alpha-satellite DNA, the major component of primate centromeres, is hardly defined because of the difficulty in its sequence assembly and its rapid evolution when compared with most genomic sequences. By using several approaches, we have cloned, sequenced, and characterized alpha-satellite sequences from two species representing critical nodes in the primate phylogeny: the white-cheeked gibbon, a lesser ape, and marmoset, a New World monkey. Sequence analyses demonstrate that white-cheeked gibbon and marmoset alpha-satellite sequences are formed by units of approximately 171 and approximately 342 bp, respectively, and they both lack the high-order structure found in humans and great apes. Fluorescent in situ hybridization characterization shows a broad dispersal of alpha-satellite in the white-cheeked gibbon genome including centromeric, telomeric, and chromosomal interstitial localizations. On the other hand, centromeres in marmoset appear organized in highly divergent dimers roughly of 342 bp that show a similarity between monomers much lower than previously reported dimers, thus representing an ancient dimeric structure. All these data shed light on the evolution of the centromeric sequences in Primates. Our results suggest radical differences in the structure, organization, and evolution of alpha-satellite DNA among different primate species, supporting the notion that 1) all the centromeric sequence in Primates evolved by genomic amplification, unequal crossover, and sequence homogenization using a 171 bp monomer as the basic seeding unit and 2) centromeric function is linked to relatively short repeated elements, more than higher-order structure. Moreover, our data indicate that complex higher-order repeat structures are a peculiarity of the hominid lineage, showing the more complex organization in humans.

Breakpoint graphs and ancestral genome reconstructions

Recently completed whole-genome sequencing projects marked the transition from gene-based phylogenetic studies to phylogenomics analysis of entire genomes. We developed an algorithm MGRA for reconstructing ancestral genomes and used it to study the rearrangement history of seven mammalian genomes: human, chimpanzee, macaque, mouse, rat, dog, and opossum. MGRA relies on the notion of the multiple breakpoint graphs to overcome some limitations of the existing approaches to ancestral genome reconstructions. MGRA also generates the rearrangement-based characters guiding the phylogenetic tree reconstruction when the phylogeny is unknown.

Characterization of six human disease-associated inversion polymorphisms

The human genome is a highly dynamic structure that shows a wide range of genetic polymorphic variation. Unlike other types of structural variation, little is known about inversion variants within normal individuals because such events are typically balanced and are difficult to detect and analyze by standard molecular approaches. Using sequence-based, cytogenetic and genotyping approaches, we characterized six large inversion polymorphisms that map to regions associated with genomic disorders with complex segmental duplications mapping at the breakpoints. We developed a metaphase FISH-based assay to genotype inversions and analyzed the chromosomes of 27 individuals from three HapMap populations. In this subset, we find that these inversions are less frequent or absent in Asians when compared with European and Yoruban populations. Analyzing multiple individuals from outgroup species of great apes, we show that most of these large inversion polymorphisms are specific to the human lineage with two exceptions, 17q21.31 and 8p23 inversions, which are found to be similarly polymorphic in other great ape species and where the inverted allele represents the ancestral state. Investigating linkage disequilibrium relationships with genotyped SNPs, we provide evidence that most of these inversions appear to have arisen on at least two different haplotype backgrounds. In these cases, discovery and genotyping methods based on SNPs may be confounded and molecular cytogenetics remains the only method to genotype these inversions.

Combinatorial content of CCL3L and CCL4L gene copy numbers influence HIV-AIDS susceptibility in Ukrainian children

OBJECTIVE

CCL3L and CCL4L genes encode HIV-suppressive chemokines, colocalize on chromosome 17q12 and have copy number variation. Copy number variation of CCL3L associates with HIV-AIDS susceptibility. Here, we determined the influence of the combinatorial content of distinct CCL3L and CCL4L genes on HIV-AIDS susceptibility.

METHODS

By designing gene-specific assays, the association between doses of all CCL3L or CCL4L genes or their individual duplicated components (CCL3La/b and CCL4La/b) with HIV-AIDS susceptibility was determined in 298 perinatally exposed Ukrainian children.

RESULTS

The odds of transmission was increased in children with less than two copies of CCL3L or CCL4L, compared with those with at least two copies, and 10-fold higher when both mother and offspring had less than two CCL3L or CCL4L copies, compared with mother-child pairs with at least two copies. The extent of the pair-wise correlations between CCL3La, CCL3Lb, CCL4La and CCL4Lb copy number varied extensively, with an inverse correlation between CCL4L genes that transcribe a classical chemokine (CCL4La) versus aberrantly-spliced transcripts (CCL4Lb). Children possessing only CCL4Lb progressed four times faster to AIDS than those with only CCL4La. A lower content of CCL3L and CCL4L genes that transcribe classical chemokines was associated with enhanced HIV-AIDS susceptibility.

CONCLUSION

Transmission risk is greatest when mother and offspring both have low CCL3L or CCL4L gene doses. The impact on HIV-AIDS susceptibility of the chemokine gene-rich locus on 17q12 is dependent on the balance between the doses of genes conferring protective (CCL3La and CCL4La) versus detrimental (CCL4Lb) effects. Hence, the combinatorial genomic content of distinct genes within a copy number variable region may determine disease susceptibility.

Evolutionary toggling of the MAPT 17q21.31 inversion region

Using comparative sequencing approaches, we investigated the evolutionary history of the European-enriched 17q21.31 MAPT inversion polymorphism. We present a detailed, BAC-based sequence assembly of the inverted human H2 haplotype and compare it to the sequence structure and genetic variation of the corresponding 1.5-Mb region for the noninverted H1 human haplotype and that of chimpanzee and orangutan. We found that inversion of the MAPT region is similarly polymorphic in other great ape species, and we present evidence that the inversions occurred independently in chimpanzees and humans. In humans, the inversion breakpoints correspond to core duplications with the LRRC37 gene family. Our analysis favors the H2 configuration and sequence haplotype as the likely great ape and human ancestral state, with inversion recurrences during primate evolution. We show that the H2 architecture has evolved more extensive sequence homology, perhaps explaining its tendency to undergo microdeletion associated with mental retardation in European populations.

Copy number variation of CCL3-like genes affects rate of progression to simian-AIDS in Rhesus Macaques (Macaca mulatta)

Variation in genes underlying host immunity can lead to marked differences in susceptibility to HIV infection among humans. Despite heavy reliance on non-human primates as models for HIV/AIDS, little is known about which host factors are shared and which are unique to a given primate lineage. Here, we investigate whether copy number variation (CNV) at CCL3-like genes (CCL3L), a key genetic host factor for HIV/AIDS susceptibility and cell-mediated immune response in humans, is also a determinant of time until onset of simian-AIDS in rhesus macaques. Using a retrospective study of 57 rhesus macaques experimentally infected with SIVmac, we find that CCL3L CNV explains approximately 18% of the variance in time to simian-AIDS (p<0.001) with lower CCL3L copy number associating with more rapid disease course. We also find that CCL3L copy number varies significantly (p<10(-6)) among rhesus subpopulations, with Indian-origin macaques having, on average, half as many CCL3L gene copies as Chinese-origin macaques. Lastly, we confirm that CCL3L shows variable copy number in humans and chimpanzees and report on CCL3L CNV within and among three additional primate species. On the basis of our findings we suggest that (1) the difference in population level copy number may explain previously reported observations of longer post-infection survivorship of Chinese-origin rhesus macaques, (2) stratification by CCL3L copy number in rhesus SIV vaccine trials will increase power and reduce noise due to non-vaccine-related differences in survival, and (3) CCL3L CNV is an ancestral component of the primate immune response and, therefore, copy number variation has not been driven by HIV or SIV per se.