Population Structure, Stratification, and Introgression of Human Structural Variation.

Structural variants contribute substantially to genetic diversity and are important evolutionarily and medically, but they are still understudied. Here we present a comprehensive analysis of structural variation in the Human Genome Diversity panel, a high-coverage dataset of 911 samples from 54 diverse worldwide populations. We identify, in total, 126,018 variants, 78% of which were not identified in previous global sequencing projects. Some reach high frequency and are private to continental groups or even individual populations, including regionally restricted runaway duplications and putatively introgressed variants from archaic hominins. By de novo assembly of 25 genomes using linked-read sequencing, we discover 1,643 breakpoint-resolved unique insertions, in aggregate accounting for 1.9 Mb of sequence absent from the GRCh38 reference. Our results illustrate the limitation of a single human reference and the need for high-quality genomes from diverse populations to fully discover and understand human genetic variation.

Genome-wide de novo L1 Retrotransposition Connects Endonuclease Activity with Replication.

L1 retrotransposon-derived sequences comprise approximately 17% of the human genome. Darwinian selective pressures alter L1 genomic distributions during evolution, confounding the ability to determine initial L1 integration preferences. Here, we generated high-confidence datasets of greater than 88,000 engineered L1 insertions in human cell lines that act as proxies for cells that accommodate retrotransposition in vivo. Comparing these insertions to a null model, in which L1 endonuclease activity is the sole determinant dictating L1 integration preferences, demonstrated that L1 insertions are not significantly enriched in genes, transcribed regions, or open chromatin. By comparison, we provide compelling evidence that the L1 endonuclease disproportionately cleaves predominant lagging strand DNA replication templates, while lagging strand 3'-hydroxyl groups may prime endonuclease-independent L1 retrotransposition in a Fanconi anemia cell line. Thus, acquisition of an endonuclease domain, in conjunction with the ability to integrate into replicating DNA, allowed L1 to become an autonomous, interspersed retrotransposon.

Functional dissection of human mitotic genes using CRISPR-Cas9 tiling screens.

The identity of human protein-coding genes is well known, yet our in-depth knowledge of their molecular functions and domain architecture remains limited by shortcomings in homology-based predictions and experimental approaches focused on whole-gene depletion. To bridge this knowledge gap, we developed a method that leverages CRISPR-Cas9-induced mutations across protein-coding genes for the a priori identification of functional regions at the sequence level. As a test case, we applied this method to 48 human mitotic genes, revealing hundreds of regions required for cell proliferation, including domains that were experimentally characterized, ones that were predicted based on homology, and novel ones. We validated screen outcomes for 15 regions, including amino acids 387-402 of Mad1, which were previously uncharacterized but contribute to Mad1 kinetochore localization and chromosome segregation fidelity. Altogether, we demonstrate that CRISPR-Cas9-based tiling mutagenesis identifies key functional domains in protein-coding genes de novo, which elucidates separation of function mutants and allows functional annotation across the human proteome.

N6-Methyladenine DNA Modification in the Human Genome.

DNA N6-methyladenine (6mA) modification is the most prevalent DNA modification in prokaryotes, but whether it exists in human cells and whether it plays a role in human diseases remain enigmatic. Here, we showed that 6mA is extensively present in the human genome, and we cataloged 881,240 6mA sites accounting for ∼0.051% of the total adenines. [G/C]AGG[C/T] was the most significantly associated motif with 6mA modification. 6mA sites were enriched in the coding regions and mark actively transcribed genes in human cells. DNA 6mA and N6-demethyladenine modification in the human genome were mediated by methyltransferase N6AMT1 and demethylase ALKBH1, respectively. The abundance of 6mA was significantly lower in cancers, accompanied by decreased N6AMT1 and increased ALKBH1 levels, and downregulation of 6mA modification levels promoted tumorigenesis. Collectively, our results demonstrate that DNA 6mA modification is extensively present in human cells and the decrease of genomic DNA 6mA promotes human tumorigenesis.

Are non-protein coding RNAs junk or treasure?: An attempt to explain and reconcile opposing viewpoints of whether the human genome is mostly transcribed into non-functional or functional RNAs.

The human genome project's lasting legacies are the emerging insights into human physiology and disease, and the ascendance of biology as the dominant science of the 21st century. Sequencing revealed that >90% of the human genome is not coding for proteins, as originally thought, but rather is overwhelmingly transcribed into non-protein coding, or non-coding, RNAs (ncRNAs). This discovery initially led to the hypothesis that most genomic DNA is "junk", a term still championed by some geneticists and evolutionary biologists. In contrast, molecular biologists and biochemists studying the vast number of transcripts produced from most of this genome "junk" often surmise that these ncRNAs have biological significance. What gives? This essay contrasts the two opposing, extant viewpoints, aiming to explain their bases, which arise from distinct reference frames of the underlying scientific disciplines. Finally, it aims to reconcile these divergent mindsets in hopes of stimulating synergy between scientific fields.

Analysis of missense variants in the human genome reveals widespread gene-specific clustering and improves prediction of pathogenicity.

We used a machine learning approach to analyze the within-gene distribution of missense variants observed in hereditary conditions and cancer. When applied to 840 genes from the ClinVar database, this approach detected a significant non-random distribution of pathogenic and benign variants in 387 (46%) and 172 (20%) genes, respectively, revealing that variant clustering is widespread across the human exome. This clustering likely occurs as a consequence of mechanisms shaping pathogenicity at the protein level, as illustrated by the overlap of some clusters with known functional domains. We then took advantage of these findings to develop a pathogenicity predictor, MutScore, that integrates qualitative features of DNA substitutions with the new additional information derived from this positional clustering. Using a random forest approach, MutScore was able to identify pathogenic missense mutations with very high accuracy, outperforming existing predictive tools, especially for variants associated with autosomal-dominant disease and cancer. Thus, the within-gene clustering of pathogenic and benign DNA changes is an important and previously underappreciated feature of the human exome, which can be harnessed to improve the prediction of pathogenicity and disambiguation of DNA variants of uncertain significance.

Benchmarking multi-platform sequencing technologies for human genome assembly.

Genome assembly is a computational technique that involves piecing together deoxyribonucleic acid (DNA) fragments generated by sequencing technologies to create a comprehensive and precise representation of the entire genome. Generating a high-quality human reference genome is a crucial prerequisite for comprehending human biology, and it is also vital for downstream genomic variation analysis. Many efforts have been made over the past few decades to create a complete and gapless reference genome for humans by using a diverse range of advanced sequencing technologies. Several available tools are aimed at enhancing the quality of haploid and diploid human genome assemblies, which include contig assembly, polishing of contig errors, scaffolding and variant phasing. Selecting the appropriate tools and technologies remains a daunting task despite several studies have investigated the pros and cons of different assembly strategies. The goal of this paper was to benchmark various strategies for human genome assembly by combining sequencing technologies and tools on two publicly available samples (NA12878 and NA24385) from Genome in a Bottle. We then compared their performances in terms of continuity, accuracy, completeness, variant calling and phasing. We observed that PacBio HiFi long-reads are the optimal choice for generating an assembly with low base errors. On the other hand, we were able to produce the most continuous contigs with Oxford Nanopore long-reads, but they may require further polishing to improve on quality. We recommend using short-reads rather than long-reads themselves to improve the base accuracy of contigs from Oxford Nanopore long-reads. Hi-C is the best choice for chromosome-level scaffolding because it can capture the longest-range DNA connectedness compared to 10× linked-reads and Bionano optical maps. However, a combination of multiple technologies can be used to further improve the quality and completeness of genome assembly. For diploid assembly, hifiasm is the best tool for human diploid genome assembly using PacBio HiFi and Hi-C data. Looking to the future, we expect that further advancements in human diploid assemblers will leverage the power of PacBio HiFi reads and other technologies with long-range DNA connectedness to enable the generation of high-quality, chromosome-level and haplotype-resolved human genome assemblies.

Unveiling human origins of replication using deep learning: accurate prediction and comprehensive analysis.

Accurate identification of replication origins (ORIs) is crucial for a comprehensive investigation into the progression of human cell growth and cancer therapy. Here, we proposed a computational approach Ori-FinderH, which can efficiently and precisely predict the human ORIs of various lengths by combining the Z-curve method with deep learning approach. Compared with existing methods, Ori-FinderH exhibits superior performance, achieving an area under the receiver operating characteristic curve (AUC) of 0.9616 for K562 cell line in 10-fold cross-validation. In addition, we also established a cross-cell-line predictive model, which yielded a further improved AUC of 0.9706. The model was subsequently employed as a fitness function to support genetic algorithm for generating artificial ORIs. Sequence analysis through iORI-Euk revealed that a vast majority of the created sequences, specifically 98% or more, incorporate at least one ORI for three cell lines (Hela, MCF7 and K562). This innovative approach could provide more efficient, accurate and comprehensive information for experimental investigation, thereby further advancing the development of this field.

Towards a Human Cell Atlas: Taking Notes from the Past.

Comprehensively characterizing the cellular composition and organization of tissues has been a long-term scientific challenge that has limited our ability to study fundamental and clinical aspects of human physiology. The Human Cell Atlas (HCA) is a global collaborative effort to create a reference map of all human cells as a basis for both understanding human health and diagnosing, monitoring, and treating disease. Many aspects of the HCA are analogous to the Human Genome Project (HGP), whose completion presents a major milestone in modern biology. To commemorate the HGP's 20-year anniversary of completion, we discuss the launch of the HCA in light of the HGP, and highlight recent progress by the HCA consortium.

The Need for a Human Pangenome Reference Sequence.

The reference human genome sequence is inarguably the most important and widely used resource in the fields of human genetics and genomics. It has transformed the conduct of biomedical sciences and brought invaluable benefits to the understanding and improvement of human health. However, the commonly used reference sequence has profound limitations, because across much of its span, it represents the sequence of just one human haplotype. This single, monoploid reference structure presents a critical barrier to representing the broad genomic diversity in the human population. In this review, we discuss the modernization of the reference human genome sequence to a more complete reference of human genomic diversity, known as a human pangenome.

Exome variant discrepancies due to reference-genome differences.

Despite release of the GRCh38 human reference genome more than seven years ago, GRCh37 remains more widely used by most research and clinical laboratories. To date, no study has quantified the impact of utilizing different reference assemblies for the identification of variants associated with rare and common diseases from large-scale exome-sequencing data. By calling variants on both the GRCh37 and GRCh38 references, we identified single-nucleotide variants (SNVs) and insertion-deletions (indels) in 1,572 exomes from participants with Mendelian diseases and their family members. We found that a total of 1.5% of SNVs and 2.0% of indels were discordant when different references were used. Notably, 76.6% of the discordant variants were clustered within discrete discordant reference patches (DISCREPs) comprising only 0.9% of loci targeted by exome sequencing. These DISCREPs were enriched for genomic elements including segmental duplications, fix patch sequences, and loci known to contain alternate haplotypes. We identified 206 genes significantly enriched for discordant variants, most of which were in DISCREPs and caused by multi-mapped reads on the reference assembly that lacked the variant call. Among these 206 genes, eight are implicated in known Mendelian diseases and 53 are associated with common phenotypes from genome-wide association studies. In addition, variant interpretations could also be influenced by the reference after lifting-over variant loci to another assembly. Overall, we identified genes and genomic loci affected by reference assembly choice, including genes associated with Mendelian disorders and complex human diseases that require careful evaluation in both research and clinical applications.

Fast heritability estimation based on MINQUE and batch training.

Heritability, the proportion of phenotypic variance explained by genome-wide single nucleotide polymorphisms (SNPs) in unrelated individuals, is an important measure of the genetic contribution to human diseases and plays a critical role in studying the genetic architecture of human diseases. Linear mixed model (LMM) has been widely used for SNP heritability estimation, where variance component parameters are commonly estimated by using a restricted maximum likelihood (REML) method. REML is an iterative optimization algorithm, which is computationally intensive when applied to large-scale datasets (e.g. UK Biobank). To facilitate the heritability analysis of large-scale genetic datasets, we develop a fast approach, minimum norm quadratic unbiased estimator (MINQUE) with batch training, to estimate variance components from LMM (LMM.MNQ.BCH). In LMM.MNQ.BCH, the parameters are estimated by MINQUE, which has a closed-form solution for fast computation and has no convergence issue. Batch training has also been adopted in LMM.MNQ.BCH to accelerate the computation for large-scale genetic datasets. Through simulations and real data analysis, we demonstrate that LMM.MNQ.BCH is much faster than two existing approaches, GCTA and BOLT-REML.

Kinesin family member 12-related hepatopathy: A generally indolent disorder with elevated gamma-glutamyl-transferase activity.

Exome sequencing (ES) has identified biallelic kinesin family member 12 (KIF12) mutations as underlying neonatal cholestatic liver disease. We collected information on onset and progression of this entity. Among consecutively referred pediatric patients at our centers, diagnostic ES identified 4 patients with novel, biallelic KIF12 variants using the human GRCh38 reference sequence, as KIF12 remains incompletely annotated in the older reference sequence GRCh37. A review of these and of 21 reported patients with KIF12 variants found that presentation with elevated serum transaminase activity in the context of trivial respiratory infection, without clinical features of liver disease, was more common (n = 18) than manifest cholestatic disease progressing rapidly to liver transplantation (LT; n = 7). Onset of liver disease was at age <1 year in 15 patients; LT was more common in this group. Serum gamma-glutamyl transpeptidase activity (GGT) was elevated in all patients, and total bilirubin was elevated in 15 patients. Liver fibrosis or cirrhosis was present in 14 of 18 patients who were biopsied. The 16 different pathogenic variants and 11 different KIF12 genotypes found were not correlated with age of onset or progression to LT. Identification of biallelic pathogenic KIF12 variants distinguishes KIF12-related disease from other entities with elevated GGT.

The Human Genome Organisation (HUGO) and a vision for Ecogenomics: the Ecological Genome Project.

BACKGROUND: The following outlines ethical reasons for widening the Human Genome Organisation's (HUGO) mandate to include ecological genomics. MAIN: The environment influences an organism's genome through ambient factors in the biosphere (e.g. climate and UV radiation), as well as the agents it comes into contact with, i.e. the epigenetic and mutagenic effects of inanimate chemicals and pollution, and pathogenic organisms. Emerging scientific consensus is that social determinants of health, environmental conditions and genetic factors work together to influence the risk of many complex illnesses. That paradigm can also explain the environmental and ecological determinants of health as factors that underlie the (un)healthy ecosystems on which communities rely. We suggest that The Ecological Genome Project is an aspirational opportunity to explore connections between the human genome and nature. We propose consolidating a view of Ecogenomics to provide a blueprint to respond to the environmental challenges that societies face. This can only be achieved by interdisciplinary engagement between genomics and the broad field of ecology and related practice of conservation. In this respect, the One Health approach is a model for environmental orientated work. The idea of Ecogenomics-a term that has been used to relate to a scientific field of ecological genomics-becomes the conceptual study of genomes within the social and natural environment. CONCLUSION: The HUGO Committee on Ethics, Law and Society (CELS) recommends that an interdisciplinary One Health approach should be adopted in genomic sciences to promote ethical environmentalism. This perspective has been reviewed and endorsed by the HUGO CELS and the HUGO Executive Board.

Transcriptome driven discovery of novel candidate genes for human neurological disorders in the telomer-to-telomer genome assembly era.

BACKGROUND: With the first complete draft of a human genome, the Telomere-to-Telomere Consortium unlocked previously concealed genomic regions for genetic analyses. These regions harbour nearly 2000 potential novel genes with unknown function. In order to uncover candidate genes associated with human neurological pathologies, a comparative transcriptome study using the T2T-CHM13 and the GRCh38 genome assemblies was conducted on previously published datasets for eight distinct human neurological disorders. RESULTS: The analysis of differential expression in RNA sequencing data led to the identification of 336 novel candidate genes linked to human neurological disorders. Additionally, it was revealed that, on average, 3.6% of the differentially expressed genes detected with the GRCh38 assembly may represent potential false positives. Among the noteworthy findings, two novel genes were discovered, one encoding a pore-structured protein and the other a highly ordered ß-strand-rich protein. These genes exhibited upregulation in multiple epilepsy datasets and hold promise as candidate genes potentially modulating the progression of the disease. Furthermore, an analysis of RNA derived from white matter lesions in multiple sclerosis patients indicated significant upregulation of 26 rRNA encoding genes. Additionally, putative pathology related genes were identified for Alzheimer's disease, amyotrophic lateral sclerosis, glioblastoma, glioma, and conditions resulting from the m.3242 A > G mtDNA mutation. CONCLUSION: The results presented here underline the potential of the T2T-CHM13 assembly in facilitating the discovery of candidate genes from transcriptome data in the context of human disorders. Moreover, the results demonstrate the value of remapping sequencing data to a superior genome assembly. Numerous potential pathology related genes, either as causative factors or related elements, have been unveiled, warranting further experimental validation.

Range-limited Heaps' law for functional DNA words in the human genome.

Heaps' or Herdan-Heaps' law is a linguistic law describing the relationship between the vocabulary/dictionary size (type) and word counts (token) to be a power-law function. Its existence in genomes with certain definition of DNA words is unclear partly because the dictionary size in genome could be much smaller than that in a human language. We define a DNA word as a coding region in a genome that codes for a protein domain. Using human chromosomes and chromosome arms as individual samples, we establish the existence of Heaps' law in the human genome within limited range. Our definition of words in a genomic or proteomic context is different from other definitions such as over-represented k-mers which are much shorter in length. Although an approximate power-law distribution of protein domain sizes due to gene duplication and the related Zipf's law is well known, their translation to the Heaps' law in DNA words is not automatic. Several other animal genomes are shown herein also to exhibit range-limited Heaps' law with our definition of DNA words, though with various exponents. When tokens were randomly sampled and sample sizes reach to the maximum level, a deviation from the Heaps' law was observed, but a quadratic regression in log-log type-token plot fits the data perfectly. Investigation of type-token plot and its regression coefficients could provide an alternative narrative of reusage and redundancy of protein domains as well as creation of new protein domains from a linguistic perspective.

Virus-like insertions with sequence signatures similar to those of endogenous nonretroviral RNA viruses in the human genome.

Understanding the genetics and taxonomy of ancient viruses will give us great insights into not only the origin and evolution of viruses but also how viral infections played roles in our evolution. Endogenous viruses are remnants of ancient viral infections and are thought to retain the genetic characteristics of viruses from ancient times. In this study, we used machine learning of endogenous RNA virus sequence signatures to identify viruses in the human genome that have not been detected or are already extinct. Here, we show that the k-mer occurrence of ancient RNA viral sequences remains similar to that of extant RNA viral sequences and can be differentiated from that of other human genome sequences. Furthermore, using this characteristic, we screened RNA viral insertions in the human reference genome and found virus-like insertions with phylogenetic and evolutionary features indicative of an exogenous origin but lacking homology to previously identified sequences. Our analysis indicates that animal genomes still contain unknown virus-derived sequences and provides a glimpse into the diversity of the ancient virosphere.

Convergent Mutations and Single Nucleotide Variants in Mitochondrial Genomes of Modern Humans and Neanderthals.

The genetic contributions of Neanderthals to the modern human genome have been evidenced by the comparison of present-day human genomes with paleogenomes. Neanderthal signatures in extant human genomes are attributed to intercrosses between Neanderthals and archaic anatomically modern humans (AMHs). Although Neanderthal signatures are well documented in the nuclear genome, it has been proposed that there is no contribution of Neanderthal mitochondrial DNA to contemporary human genomes. Here we show that modern human mitochondrial genomes contain 66 potential Neanderthal signatures, or Neanderthal single nucleotide variants (N-SNVs), of which 36 lie in coding regions and 7 result in nonsynonymous changes. Seven N-SNVs are associated with traits such as cycling vomiting syndrome, Alzheimer's disease and Parkinson's disease, and two N-SNVs are associated with intelligence quotient. Based on recombination tests, principal component analysis (PCA) and the complete absence of these N-SNVs in 41 archaic AMH mitogenomes, we conclude that convergent evolution, and not recombination, explains the presence of N-SNVs in present-day human mitogenomes.

Heritable human genome editing: correction, selection and treatment.

Heritable human genome editing (HHGE) to correct a nuclear gene sequence that would result in a serious genetic condition in a future child is presented as 'treatment' in various ethics and policy materials, and as morally preferable to the 'selection' practice of preimplantation genetic testing (PGT), which is subject to the disability critique. However, whether HHGE is 'treatment' for a future child, or another form of 'selection', or whether HHGE instead 'treats' prospective parents, are now central questions in the debate regarding its possible legalisation. This article argues that the idea of 'treatment' for a future child is largely a proxy for 'seriousness of purpose', intended to distinguish HHGE to avoid serious genetic conditions from less obviously justifiable uses; that HHGE is best understood, and morally justified, as a form of 'treatment' for prospective parents who strongly desire an unaffected genetically related child and who have no, or poor, options to achieve this; that HHGE would be morally permissible if consistent with that child's welfare; that legalisation is supportable with reference to the right to respect for private and family life under Article 8 of the European Convention on Human Rights; and that HHGE is morally distinguishable from PGT.

Cultivating DNA Sequencing Technology After the Human Genome Project.

When the Human Genome Project was completed in 2003, automated Sanger DNA sequencing with fluorescent dye labels was the dominant technology. Several nascent alternative methods based on older ideas that had not been fully developed were the focus of technical researchers and companies. Funding agencies recognized the dynamic nature of technology development and that, beyond the Human Genome Project, there were growing opportunities to deploy DNA sequencing in biological research. Consequently, the National Human Genome Research Institute of the National Institutes of Health created a program-widely known as the Advanced Sequencing Technology Program-that stimulated all stages of development of new DNA sequencing methods, from innovation to advanced manufacturing and production testing, with the goal of reducing the cost of sequencing a human genome first to $100,000 and then to $1,000. The events of this period provide a powerful example of how judicious funding of academic and commercial partners can rapidly advance core technology developments that lead to profound advances across the scientific landscape.