Genome maps across 26 human populations reveal population-specific patterns of structural variation

Next-generation data filtering in the genomics era

Genomic data are ubiquitous across disciplines, from agriculture to biodiversity, ecology, evolution and human health. However, these datasets often contain noise or errors and are missing information that can affect the accuracy and reliability of subsequent computational analyses and conclusions. A key step in genomic data analysis is filtering - removing sequencing bases, reads, genetic variants and/or individuals from a dataset - to improve data quality for downstream analyses. Researchers are confronted with a multitude of choices when filtering genomic data; they must choose which filters to apply and select appropriate thresholds. To help usher in the next generation of genomic data filtering, we review and suggest best practices to improve the implementation, reproducibility and reporting standards for filter types and thresholds commonly applied to genomic datasets. We focus mainly on filters for minor allele frequency, missing data per individual or per locus, linkage disequilibrium and Hardy-Weinberg deviations. Using simulated and empirical datasets, we illustrate the large effects of different filtering thresholds on common population genetics statistics, such as Tajima's D value, population differentiation (FST), nucleotide diversity (π) and effective population size (Ne).

Novel and Recurrent Copy Number Variants in ABCA4-Associated Retinopathy

ABCA4 is the most frequently mutated gene leading to inherited retinal disease (IRD) with over 2200 pathogenic variants reported to date. Of these, ~1% are copy number variants (CNVs) involving the deletion or duplication of genomic regions, typically >50 nucleotides in length. An in-depth assessment of the current literature based on the public database LOVD, regarding the presence of known CNVs and structural variants in ABCA4, and additional sequencing analysis of ABCA4 using single-molecule Molecular Inversion Probes (smMIPs) for 148 probands highlighted recurrent and novel CNVs associated with ABCA4-associated retinopathies. An analysis of the coverage depth in the sequencing data led to the identification of eleven deletions (six novel and five recurrent), three duplications (one novel and two recurrent) and one complex CNV. Of particular interest was the identification of a complex defect, i.e., a 15.3 kb duplicated segment encompassing exon 31 through intron 41 that was inserted at the junction of a downstream 2.7 kb deletion encompassing intron 44 through intron 47. In addition, we identified a 7.0 kb tandem duplication of intron 1 in three cases. The identification of CNVs in ABCA4 can provide patients and their families with a genetic diagnosis whilst expanding our understanding of the complexity of diseases caused by ABCA4 variants.

Optical Genome Mapping Reveals Genomic Alterations upon Gene Editing in hiPSCs: Implications for Neural Tissue Differentiation and Brain Organoid Research

Genome editing, notably CRISPR (cluster regularly interspaced short palindromic repeats)/Cas9 (CRISPR-associated protein 9), has revolutionized genetic engineering allowing for precise targeted modifications. This technique's combination with human induced pluripotent stem cells (hiPSCs) is a particularly valuable tool in cerebral organoid (CO) research. In this study, CRISPR/Cas9-generated fluorescently labeled hiPSCs exhibited no significant morphological or growth rate differences compared with unedited controls. However, genomic aberrations during gene editing necessitate efficient genome integrity assessment methods. Optical genome mapping, a high-resolution genome-wide technique, revealed genomic alterations, including chromosomal copy number gain and losses affecting numerous genes. Despite these genomic alterations, hiPSCs retain their pluripotency and capacity to generate COs without major phenotypic changes but one edited cell line showed potential neuroectodermal differentiation impairment. Thus, this study highlights optical genome mapping in assessing genome integrity in CRISPR/Cas9-edited hiPSCs emphasizing the need for comprehensive integration of genomic and morphological analysis to ensure the robustness of hiPSC-based models in cerebral organoid research.

Integrating Optical Genome Mapping and Whole Genome Sequencing in Somatic Structural Variant Detection

Structural variants drive tumorigenesis by disrupting normal gene function through insertions, inversions, translocations, and copy number changes, including deletions and duplications. Detecting structural variants is crucial for revealing their roles in tumor development, clinical outcomes, and personalized therapy. Presently, most studies rely on short-read data from next-generation sequencing that aligns back to a reference genome to determine if and, if so, where a structural variant occurs. However, structural variant discovery by short-read sequencing is challenging, primarily because of the difficulty in mapping regions of repetitive sequences. Optical genome mapping (OGM) is a recent technology used for imaging and assembling long DNA strands to detect structural variations. To capture the structural variant landscape more thoroughly in the human genome, we developed an integrated pipeline that combines Bionano OGM and Illumina whole-genome sequencing and applied it to samples from 29 pediatric B-ALL patients. The addition of OGM allowed us to identify 511 deletions, 506 insertions, 93 duplications/gains, and 145 translocations that were otherwise missed in the short-read data. Moreover, we identified several novel gene fusions, the expression of which was confirmed by RNA sequencing. Our results highlight the benefit of integrating OGM and short-read detection methods to obtain a comprehensive analysis of genetic variation that can aid in clinical diagnosis, provide new therapeutic targets, and improve personalized medicine in cancers driven by structural variation.

Optical Genome Mapping as a Potential Routine Clinical Diagnostic Method

Chromosome analysis (CA) and chromosomal microarray analysis (CMA) have been successfully used to diagnose genetic disorders. However, many conditions remain undiagnosed due to limitations in resolution (CA) and detection of only unbalanced events (CMA). Optical genome mapping (OGM) has the potential to address these limitations by capturing both structural variants (SVs) resulting in copy number changes and balanced rearrangements with high resolution. In this study, we investigated OGM's concordance using 87 SVs previously identified by CA, CMA, or Southern blot. Overall, OGM was 98% concordant with only three discordant cases: (1) uncalled translocation with one breakpoint in a centromere; (2) uncalled duplication with breakpoints in the pseudoautosomal region 1; and (3) uncalled mosaic triplication originating from a marker chromosome. OGM provided diagnosis for three previously unsolved cases: (1) disruption of the SON gene due to a balanced reciprocal translocation; (2) disruption of the NBEA gene due to an inverted insertion; (3) disruption of the TSC2 gene due to a mosaic deletion. We show that OGM is a valid method for the detection of many types of SVs in a single assay and is highly concordant with legacy cytogenomic methods; however, it has limited SV detection capabilities in centromeric and pseudoautosomal regions.

Cribriform-morular Thyroid Carcinoma Arising in a Medulloblastoma Survivor: Two Metachronous Tumors Shared with the Activation of the Wnt Signaling Pathway

Wnt signaling pathway activation is involved in the pathogenesis of a series of malignant tumors and is characterized by the nuclear accumulation of β-catenin protein. The occurrence of two or more Wnt pathway-associated tumors in a single individual is uncommon and generally attributed to inherited cancer syndrome, especially familial adenomatous polyposis (FAP). Herein, we presented a rare case of a child who suffered from the occurrence of Wnt-activated medulloblastoma and cribriform-morular thyroid carcinoma (CMTC) within a 9-year interval. She had no history of FAP and harbored an unexpected somatic mutation of the APC gene in the CMTC tumor. The potential agents involved in the pathogenesis of the two molecular-linked tumors other than FAP were discussed in this report.

The complete diploid reference genome of RPE-1 identifies human phased epigenetic landscapes

Comparative analysis of recent human genome assemblies highlights profound sequence divergence that peaks within polymorphic loci such as centromeres. This raises the question about the adequacy of relying on human reference genomes to accurately analyze sequencing data derived from experimental cell lines. Here, we generated the complete diploid genome assembly for the human retinal epithelial cells (RPE-1), a widely used non-cancer laboratory cell line with a stable karyotype, to use as matched reference for multi-omics sequencing data analysis. Our RPE1v1.0 assembly presents completely phased haplotypes and chromosome-level scaffolds that span centromeres with ultra-high base accuracy (>QV60). We mapped the haplotype-specific genomic variation specific to this cell line including t(Xq;10q), a stable 73.18 Mb duplication of chromosome 10 translocated onto the microdeleted chromosome X telomere t(Xq;10q). Polymorphisms between haplotypes of the same genome reveals genetic and epigenetic variation for all chromosomes, especially at centromeres. The RPE-1 assembly as matched reference genome improves mapping quality of multi-omics reads originating from RPE-1 cells with drastic reduction in alignments mismatches compared to using the most complete human reference to date (CHM13). Leveraging the accuracy achieved using a matched reference, we were able to identify the kinetochore sites at base pair resolution and show unprecedented variation between haplotypes. This work showcases the use of matched reference genomes for multiomics analyses and serves as the foundation for a call to comprehensively assemble experimentally relevant cell lines for widespread application.

Benchmarking long-read genome sequence alignment tools for human genomics applications

Background

The utility of long-read genome sequencing platforms has been shown in many fields including whole genome assembly, metagenomics, and amplicon sequencing. Less clear is the applicability of long reads to reference-guided human genomics, which is the foundation of genomic medicine. Here, we benchmark available platform-agnostic alignment tools on datasets from nanopore and single-molecule real-time platforms to understand their suitability in producing a genome representation.

Results

For this study, we leveraged publicly-available data from sample NA12878 generated on Oxford Nanopore and sample NA24385 on Pacific Biosciences platforms. We employed state of the art sequence alignment tools including GraphMap2, long-read aligner (LRA), Minimap2, CoNvex Gap-cost alignMents for Long Reads (NGMLR), and Winnowmap2. Minimap2 and Winnowmap2 were computationally lightweight enough for use at scale, while GraphMap2 was not. NGMLR took a long time and required many resources, but produced alignments each time. LRA was fast, but only worked on Pacific Biosciences data. Each tool widely disagreed on which reads to leave unaligned, affecting the end genome coverage and the number of discoverable breakpoints. No alignment tool independently resolved all large structural variants (1,001-100,000 base pairs) present in the Database of Genome Variants (DGV) for sample NA12878 or the truthset for NA24385.

Conclusions

These results suggest a combined approach is needed for LRS alignments for human genomics. Specifically, leveraging alignments from three tools will be more effective in generating a complete picture of genomic variability. It should be best practice to use an analysis pipeline that generates alignments with both Minimap2 and Winnowmap2 as they are lightweight and yield different views of the genome. Depending on the question at hand, the data available, and the time constraints, NGMLR and LRA are good options for a third tool. If computational resources and time are not a factor for a given case or experiment, NGMLR will provide another view, and another chance to resolve a case. LRA, while fast, did not work on the nanopore data for our cluster, but PacBio results were promising in that those computations completed faster than Minimap2. Due to its significant burden on computational resources and slow run time, Graphmap2 is not an ideal tool for exploration of a whole human genome generated on a long-read sequencing platform.

Comparative Benchmarking of Optical Genome Mapping and Chromosomal Microarray Reveals High Technological Concordance in CNV Identification and Additional Structural Variant Refinement

The recommended practice for individuals suspected of a genetic etiology for disorders including unexplained developmental delay/intellectual disability (DD/ID), autism spectrum disorders (ASD), and multiple congenital anomalies (MCA) involves a genetic testing workflow including chromosomal microarray (CMA), Fragile-X testing, karyotype analysis, and/or sequencing-based gene panels. Since genomic imbalances are often found to be causative, CMA is recommended as first tier testing for many indications. Optical genome mapping (OGM) is an emerging next generation cytogenomic technique that can detect not only copy number variants (CNVs), triploidy and absence of heterozygosity (AOH) like CMA, but can also define the location of duplications, and detect other structural variants (SVs), including balanced rearrangements and repeat expansions/contractions. This study compares OGM to CMA for clinically reported genomic variants, some of these samples also have structural characterization by fluorescence in situ hybridization (FISH). OGM was performed on IRB approved, de-identified specimens from 55 individuals with genomic abnormalities previously identified by CMA (61 clinically reported abnormalities). SVs identified by OGM were filtered by a control database to remove polymorphic variants and against an established gene list to prioritize clinically relevant findings before comparing with CMA and FISH results. OGM results showed 100% concordance with CMA findings for pathogenic variants and 98% concordant for all pathogenic/likely pathogenic/variants of uncertain significance (VUS), while also providing additional insight into the genomic structure of abnormalities that CMA was unable to provide. OGM demonstrates equivalent performance to CMA for CNV and AOH detection, enhanced by its ability to determine the structure of the genome. This work adds to an increasing body of evidence on the analytical validity and ability to detect clinically relevant abnormalities identified by CMA. Moreover, OGM identifies translocations, structures of duplications and complex CNVs intractable by CMA, yielding additional clinical utility.

Non-Random Enrichment of Single-Nucleotide Polymorphisms Associated with Clopidogrel Resistance within Risk Loci Linked to the Severity of Underlying Cardiovascular Diseases: The Role of Admixture

Cardiovascular disease (CVD) is one of the leading causes of death in Puerto Rico, where clopidogrel is commonly prescribed to prevent ischemic events. Genetic contributors to both a poor clopidogrel response and the severity of CVD have been identified mainly in Europeans. However, the non-random enrichment of single-nucleotide polymorphisms (SNPs) associated with clopidogrel resistance within risk loci linked to underlying CVDs, and the role of admixture, have yet to be tested. This study aimed to assess the possible interaction between genetic biomarkers linked to CVDs and those associated with clopidogrel resistance among admixed Caribbean Hispanics. We identified 50 SNPs significantly associated with CVDs in previous genome-wide association studies (GWASs). These SNPs were combined with another ten SNPs related to clopidogrel resistance in Caribbean Hispanics. We developed Python scripts to determine whether SNPs related to CVDs are in close proximity to those associated with the clopidogrel response. The average and individual local ancestry (LAI) within each locus were inferred, and 60 random SNPs with their corresponding LAIs were generated for enrichment estimation purposes. Our results showed no CVD-linked SNPs in close proximity to those associated with the clopidogrel response among Caribbean Hispanics. Consequently, no genetic loci with a dual predictive role for the risk of CVD severity and clopidogrel resistance were found in this population. Native American ancestry was the most enriched within the risk loci linked to CVDs in this population. The non-random enrichment of disease susceptibility loci with drug-response SNPs is a new frontier in Precision Medicine that needs further attention.

Long-Read DNA Sequencing: Recent Advances and Remaining Challenges

DNA sequencing has revolutionized medicine over recent decades. However, analysis of large structural variation and repetitive DNA, a hallmark of human genomes, has been limited by short-read technology, with read lengths of 100-300 bp. Long-read sequencing (LRS) permits routine sequencing of human DNA fragments tens to hundreds of kilobase pairs in size, using both real-time sequencing by synthesis and nanopore-based direct electronic sequencing. LRS permits analysis of large structural variation and haplotypic phasing in human genomes and has enabled the discovery and characterization of rare pathogenic structural variants and repeat expansions. It has also recently enabled the assembly of a complete, gapless human genome that includes previously intractable regions, such as highly repetitive centromeres and homologous acrocentric short arms. With the addition of protocols for targeted enrichment, direct epigenetic DNA modification detection, and long-range chromatin profiling, LRS promises to launch a new era of understanding of genetic diversity and pathogenic mutations in human populations.

Exome sequencing identified rare recurrent copy number variants and hereditary breast cancer susceptibility

Copy number variants (CNVs) are a major source of genetic variation and can disrupt genes or affect gene dosage. They are known to be causal or underlie predisposition to various diseases. However, the role of CNVs in inherited breast cancer susceptibility has not been thoroughly investigated. To address this, we performed whole-exome sequencing based analysis of rare CNVs in 98 high-risk Northern Finnish breast cancer cases. After filtering, selected candidate alleles were validated and characterized with a combination of orthogonal methods, including PCR-based approaches, optical genome mapping and long-read sequencing. This revealed three recurrent alterations: a 31 kb deletion co-occurring with a retrotransposon insertion (delins) in RAD52, a 13.4 kb deletion in HSD17B14 and a 64 kb partial duplication of RAD51C. Notably, all these genes encode proteins involved in pathways previously identified as essential for breast cancer development. Variants were genotyped in geographically matched cases and controls (altogether 278 hereditary and 1983 unselected breast cancer cases, and 1229 controls). The RAD52 delins and HSD17B14 deletion both showed significant enrichment among cases with indications of hereditary disease susceptibility. RAD52 delins was identified in 7/278 cases (2.5%, P = 0.034, OR = 2.86, 95% CI = 1.10-7.45) and HSD17B14 deletion in 8/278 cases (2.9%, P = 0.014, OR = 3.28, 95% CI = 1.31-8.23), the frequency of both variants in the controls being 11/1229 (0.9%). This suggests a role for RAD52 and HSD17B14 in hereditary breast cancer susceptibility. The RAD51C duplication was very rare, identified only in 2/278 of hereditary cases and 2/1229 controls (P = 0.157, OR = 4.45, 95% CI = 0.62-31.70). The identification of recurrent CNVs in these genes, and especially the relatively high frequency of RAD52 and HSD17B14 alterations in the Finnish population, highlights the importance of studying CNVs alongside single nucleotide variants when searching for genetic factors underlying hereditary disease predisposition.

Genetic variation in the immunoglobulin heavy chain locus shapes the human antibody repertoire

Variation in the antibody response has been linked to differential outcomes in disease, and suboptimal vaccine and therapeutic responsiveness, the determinants of which have not been fully elucidated. Countering models that presume antibodies are generated largely by stochastic processes, we demonstrate that polymorphisms within the immunoglobulin heavy chain locus (IGH) impact the naive and antigen-experienced antibody repertoire, indicating that genetics predisposes individuals to mount qualitatively and quantitatively different antibody responses. We pair recently developed long-read genomic sequencing methods with antibody repertoire profiling to comprehensively resolve IGH genetic variation, including novel structural variants, single nucleotide variants, and genes and alleles. We show that IGH germline variants determine the presence and frequency of antibody genes in the expressed repertoire, including those enriched in functional elements linked to V(D)J recombination, and overlapping disease-associated variants. These results illuminate the power of leveraging IGH genetics to better understand the regulation, function, and dynamics of the antibody response in disease.

Step by step: towards a better understanding of the genetic architecture of Alzheimer's disease

Alzheimer's disease (AD) is considered to have a large genetic component. Our knowledge of this component has progressed over the last 10 years, thanks notably to the advent of genome-wide association studies and the establishment of large consortia that make it possible to analyze hundreds of thousands of cases and controls. The characterization of dozens of chromosomal regions associated with the risk of developing AD and (in some loci) the causal genes responsible for the observed disease signal has confirmed the involvement of major pathophysiological pathways (such as amyloid precursor protein metabolism) and opened up new perspectives (such as the central role of microglia and inflammation). Furthermore, large-scale sequencing projects are starting to reveal the major impact of rare variants - even in genes like APOE - on the AD risk. This increasingly comprehensive knowledge is now being disseminated through translational research; in particular, the development of genetic risk/polygenic risk scores is helping to identify the subpopulations more at risk or less at risk of developing AD. Although it is difficult to assess the efforts still needed to comprehensively characterize the genetic component of AD, several lines of research can be improved or initiated. Ultimately, genetics (in combination with other biomarkers) might help to redefine the boundaries and relationships between various neurodegenerative diseases.

High level of complexity and global diversity of the 3q29 locus revealed by optical mapping and long-read sequencing

BACKGROUND

High sequence identity between segmental duplications (SDs) can facilitate copy number variants (CNVs) via non-allelic homologous recombination (NAHR). These CNVs are one of the fundamental causes of genomic disorders such as the 3q29 deletion syndrome (del3q29S). There are 21 protein-coding genes lost or gained as a result of such recurrent 1.6-Mbp deletions or duplications, respectively, in the 3q29 locus. While NAHR plays a role in CNV occurrence, the factors that increase the risk of NAHR at this particular locus are not well understood.

METHODS

We employed an optical genome mapping technique to characterize the 3q29 locus in 161 unaffected individuals, 16 probands with del3q29S and their parents, and 2 probands with the 3q29 duplication syndrome (dup3q29S). Long-read sequencing-based haplotype resolved de novo assemblies from 44 unaffected individuals, and 1 trio was used for orthogonal validation of haplotypes and deletion breakpoints.

RESULTS

In total, we discovered 34 haplotypes, of which 19 were novel haplotypes. Among these 19 novel haplotypes, 18 were detected in unaffected individuals, while 1 novel haplotype was detected on the parent-of-origin chromosome of a proband with the del3q29S. Phased assemblies from 44 unaffected individuals enabled the orthogonal validation of 20 haplotypes. In 89% (16/18) of the probands, breakpoints were confined to paralogous copies of a 20-kbp segment within the 3q29 SDs. In one del3q29S proband, the breakpoint was confined to a 374-bp region using long-read sequencing. Furthermore, we categorized del3q29S cases into three classes and dup3q29S cases into two classes based on breakpoints. Finally, we found no evidence of inversions in parent-of-origin chromosomes.

CONCLUSIONS

We have generated the most comprehensive haplotype map for the 3q29 locus using unaffected individuals, probands with del3q29S or dup3q29S, and available parents, and also determined the deletion breakpoint to be within a 374-bp region in one proband with del3q29S. These results should provide a better understanding of the underlying genetic architecture that contributes to the etiology of del3q29S and dup3q29S.

A review of the pangenome: how it affects our understanding of genomic variation, selection and breeding in domestic animals?

As large-scale genomic studies have progressed, it has been revealed that a single reference genome pattern cannot represent genetic diversity at the species level. While domestic animals tend to have complex routes of origin and migration, suggesting a possible omission of some population-specific sequences in the current reference genome. Conversely, the pangenome is a collection of all DNA sequences of a species that contains sequences shared by all individuals (core genome) and is also able to display sequence information unique to each individual (variable genome). The progress of pangenome research in humans, plants and domestic animals has proved that the missing genetic components and the identification of large structural variants (SVs) can be explored through pangenomic studies. Many individual specific sequences have been shown to be related to biological adaptability, phenotype and important economic traits. The maturity of technologies and methods such as third-generation sequencing, Telomere-to-telomere genomes, graphic genomes, and reference-free assembly will further promote the development of pangenome. In the future, pangenome combined with long-read data and multi-omics will help to resolve large SVs and their relationship with the main economic traits of interest in domesticated animals, providing better insights into animal domestication, evolution and breeding. In this review, we mainly discuss how pangenome analysis reveals genetic variations in domestic animals (sheep, cattle, pigs, chickens) and their impacts on phenotypes and how this can contribute to the understanding of species diversity. Additionally, we also go through potential issues and the future perspectives of pangenome research in livestock and poultry.

Copy Number Variation and Osteoporosis

PURPOSE OF REVIEW

The purpose of this review is to summarize recent findings on copy number variations and susceptibility to osteoporosis.

RECENT FINDINGS

Osteoporosis is highly influenced by genetic factors, including copy number variations (CNVs). The development and accessibility of whole genome sequencing methods has accelerated the study of CNVs and osteoporosis. Recent findings include mutations in novel genes and validation of previously known pathogenic CNVs in monogenic skeletal diseases. Identification of CNVs in genes previously associated with osteoporosis (e.g. RUNX2, COL1A2, and PLS3) has confirmed their importance in bone remodelling. This process has been associated also with the ETV1-DGKB, AGBL2, ATM, and GPR68 genes, identified by comparative genomic hybridisation microarray studies. Importantly, studies in patients with bone pathologies have associated bone disease with the long non-coding RNA LINC01260 and enhancer sequences residing in the HDAC9 gene. Further functional investigation of genetic loci harbouring CNVs associated with skeletal phenotypes will reveal their role as molecular drivers of osteoporosis.

Pan-genome inversion index reveals evolutionary insights into the subpopulation structure of Asian rice

Understanding and exploiting genetic diversity is a key factor for the productive and stable production of rice. Here, we utilize 73 high-quality genomes that encompass the subpopulation structure of Asian rice (Oryza sativa), plus the genomes of two wild relatives (O. rufipogon and O. punctata), to build a pan-genome inversion index of 1769 non-redundant inversions that span an average of ~29% of the O. sativa cv. Nipponbare reference genome sequence. Using this index, we estimate an inversion rate of ~700 inversions per million years in Asian rice, which is 16 to 50 times higher than previously estimated for plants. Detailed analyses of these inversions show evidence of their effects on gene expression, recombination rate, and linkage disequilibrium. Our study uncovers the prevalence and scale of large inversions (≥100 bp) across the pan-genome of Asian rice and hints at their largely unexplored role in functional biology and crop performance.

Local Ancestry-Informed Candidate Pathway Analysis of Warfarin Stable Dose in Latino Populations

Accuracy of warfarin dose prediction algorithms may be improved by including data from diverse populations in genetic studies of dose variability. Here, we surveyed single nucleotide polymorphisms in vitamin K-related genetic pathways for association with warfarin dose requirements in two admixed Latino populations in standard-principal component adjusted and contemporary-local ancestry adjusted regression models. A total of five variants from vitamin K-related genes/pathways were associated with warfarin dose in both cohorts (P < 0.0125) in standard models. Local ancestry-adjusted analysis unveiled 35 associated variants with absolute effects ranging from β = 9.04 ( ±2.23) to 39.18 ( ±10.89) per ancestral allele in the discovery cohort and β = 6.47 (± 2.02) to 17.82 (± 6.83) in the replication cohort. Importantly, we demonstrate the technical validity of the Tractor model in cohorts with admixed ancestry from three founder populations and bring attention to the technical hurdles obstructing the inclusion of diverse, especially admixed, populations in pharmacogenomic research.

Detection of cryptic balanced chromosomal rearrangements using high-resolution optical genome mapping

BACKGROUND

Chromosomal rearrangements have profound consequences in diverse human genetic diseases. Currently, the detection of balanced chromosomal rearrangements (BCRs) mainly relies on routine cytogenetic G-banded karyotyping. However, cryptic BCRs are hard to detect by karyotyping, and the risk of miscarriage or delivering abnormal offspring with congenital malformations in carrier couples is significantly increased. In the present study, we aimed to investigate the potential of single-molecule optical genome mapping (OGM) in unravelling cryptic chromosomal rearrangements.

METHODS

Eleven couples with normal karyotypes that had abortions/affected offspring with unbalanced rearrangements were enrolled. Ultra-high-molecular-weight DNA was isolated from peripheral blood cells and processed via OGM. The genome assembly was performed followed by variant calling and annotation. Meanwhile, multiple detection strategies, including FISH, long-range-PCR amplicon-based next-generation sequencing and Sanger sequencing were implemented to confirm the results obtained from OGM.

RESULTS

High-resolution OGM successfully detected cryptic reciprocal translocation in all recruited couples, which was consistent with the results of FISH and sequencing. All high-confidence cryptic chromosomal translocations detected by OGM were confirmed by sequencing analysis of rearrangement breakpoints. Moreover, OGM revealed additional complex rearrangement events such as inverted aberrations, further refining potential genetic interpretation.

CONCLUSION

To the best of our knowledge, this is the first study wherein OGM facilitate the rapid and robust detection of cryptic chromosomal reciprocal translocations in clinical practice. With the excellent performance, our findings suggest that OGM is well qualified as an accurate, comprehensive and first-line method for detecting cryptic BCRs in routine clinical testing.

Genomic technologies to improve variation identification in undiagnosed diseases

Human genome variation has increasingly posed challenges and opportunities for patients, medical providers, and an increasing group of stakeholders including advocacy groups, disadvantaged communities, public health experts, and scientists. Here, advances in genomic sequencing and mapping technologies are discussed with particular attention to the increasing ability to detect personal and population genome variation and the potential for accurate integration of variation into health and disease-related care. Genome mapping, one technique used to create genome map scaffolds, has now been combined with long read sequencing. New technologies have led to improved variation detection, including cryptic structural variation and diverse variants with different degrees of disease association. Combined with advances in automated and medical interpretations, variation detection is increasingly being applied in healthcare. These advances promise to make disease diagnostics more rapid, and potentially more accessible, to those with medical needs. Consequentially, the need for medical genetics and genomics experts is increasing. Here, the opportunities and potential challenges for application of genome-scale variation detection in disease are examined. (<300 words).

The Genetics of Intellectual Disability

Intellectual disability (ID) has a prevalence of ~2-3% in the general population, having a large societal impact. The underlying cause of ID is largely of genetic origin; however, identifying this genetic cause has in the past often led to long diagnostic Odysseys. Over the past decades, improvements in genetic diagnostic technologies and strategies have led to these causes being more and more detectable: from cytogenetic analysis in 1959, we moved in the first decade of the 21st century from genomic microarrays with a diagnostic yield of ~20% to next-generation sequencing platforms with a yield of up to 60%. In this review, we discuss these various developments, as well as their associated challenges and implications for the field of ID, which highlight the revolutionizing shift in clinical practice from a phenotype-first into genotype-first approach.

Detection of Transgene Location in the CYP2A13/2B6/2F1-transgenic Mouse Model using Optical Genome Mapping Technology

Most transgenic mouse models are generated through random integration of the transgene. The location of the transgene provides valuable information for assessing potential effects of the transgenesis on the host and for designing genotyping protocols that can amplify across the integration site, but it is challenging to identify. Here, we report the successful utility of optical genome mapping technology to identify the transgene insertion site in a CYP2A13/2B6/2F1-transgenic mouse model, which produces three human cytochrome P450 (P450) enzymes (CYP2A13, CYP2B6, and CYP2F1) that are encoded by neighboring genes on human chromosome 19. These enzymes metabolize many drugs, respiratory toxicants, and chemical carcinogens. Initial efforts to identify candidate insertion sites by whole genome sequencing was unsuccessful, apparently because the transgene is located in a region of the mouse genome that contains highly repetitive sequences. Subsequent utility of the optical genome mapping approach, which compares genome-wide marker distribution between the transgenic mouse genome and a reference mouse (GRCm38) or human (GRCh38) genome, localized the insertion site to mouse chromosome 14, between two marker positions at 4451324 base pair and 4485032 base pair. A transgene-mouse genome junction sequence was further identified through long-polymerase chain reaction amplification and DNA sequencing at GRCm38 Chr.14:4484726. The transgene insertion (∼2.4 megabase pair) contained 5-7 copies of the human transgenes, which replaced a 26.9-33.4 kilobase pair mouse genomic region, including exons 1-4 of Gm3182, a predicted and highly redundant gene. Finally, the sequencing results enabled the design of a new genotyping protocol that can distinguish between hemizygous and homozygous CYP2A13/2B6/2F1-transgenic mice. SIGNIFICANCE STATEMENT: This study characterizes the genomic structure of, and provides a new genotyping method for, a transgenic mouse model that expresses three human P450 enzymes, CYP2A13, CYP2B6, and CYP2F1, that are important in xenobiotic metabolism and toxicity. The demonstrated success in applying the optical genome mapping technology for identification of transgene insertion sites should encourage others to do the same for other transgenic models generated through random integration, including most of the currently available human P450 transgenic mouse models.

CONGA: Copy number variation genotyping in ancient genomes and low-coverage sequencing data

To date, ancient genome analyses have been largely confined to the study of single nucleotide polymorphisms (SNPs). Copy number variants (CNVs) are a major contributor of disease and of evolutionary adaptation, but identifying CNVs in ancient shotgun-sequenced genomes is hampered by typical low genome coverage (<1×) and short fragments (<80 bps), precluding standard CNV detection software to be effectively applied to ancient genomes. Here we present CONGA, tailored for genotyping CNVs at low coverage. Simulations and down-sampling experiments suggest that CONGA can genotype deletions >1 kbps with F-scores >0.75 at ≥1×, and distinguish between heterozygous and homozygous states. We used CONGA to genotype 10,002 outgroup-ascertained deletions across a heterogenous set of 71 ancient human genomes spanning the last 50,000 years, produced using variable experimental protocols. A fraction of these (21/71) display divergent deletion profiles unrelated to their population origin, but attributable to technical factors such as coverage and read length. The majority of the sample (50/71), despite originating from nine different laboratories and having coverages ranging from 0.44×-26× (median 4×) and average read lengths 52-121 bps (median 69), exhibit coherent deletion frequencies. Across these 50 genomes, inter-individual genetic diversity measured using SNPs and CONGA-genotyped deletions are highly correlated. CONGA-genotyped deletions also display purifying selection signatures, as expected. CONGA thus paves the way for systematic CNV analyses in ancient genomes, despite the technical challenges posed by low and variable genome coverage.

Recent advances and current challenges in population genomics of structural variation in animals and plants

The field of population genomics has seen a surge of studies on genomic structural variation over the past two decades. These studies witnessed that structural variation is taxonomically ubiquitous and represent a dominant form of genetic variation within species. Recent advances in technology, especially the development of long-read sequencing platforms, have enabled the discovery of structural variants (SVs) in previously inaccessible genomic regions which unlocked additional structural variation for population studies and revealed that more SVs contribute to evolution than previously perceived. An increasing number of studies suggest that SVs of all types and sizes may have a large effect on phenotype and consequently major impact on rapid adaptation, population divergence, and speciation. However, the functional effect of the vast majority of SVs is unknown and the field generally lacks evidence on the phenotypic consequences of most SVs that are suggested to have adaptive potential. Non-human genomes are heavily under-represented in population-scale studies of SVs. We argue that more research on other species is needed to objectively estimate the contribution of SVs to evolution. We discuss technical challenges associated with SV detection and outline the most recent advances towards more representative reference genomes, which opens a new era in population-scale studies of structural variation.

The 22q11.2 Low Copy Repeats

LCR22s are among the most complex loci in the human genome and are susceptible to nonallelic homologous recombination. This can lead to a variety of genomic disorders, including deletions, duplications, and translocations, of which the 22q11.2 deletion syndrome is the most common in humans. Interrogating these phenomena is difficult due to the high complexity of the LCR22s and the inaccurate representation of the LCRs across different reference genomes. Optical mapping techniques, which provide long-range chromosomal maps, could be used to unravel the complex duplicon structure. These techniques have already uncovered the hypervariability of the LCR22-A haplotype in the human population. Although optical LCR22 mapping is a major step forward, long-read sequencing approaches will be essential to reach nucleotide resolution of the LCR22s and map the crossover sites. Accurate maps and sequences are needed to pinpoint potential predisposing alleles and, most importantly, allow for genotype-phenotype studies exploring the role of the LCR22s in health and disease. In addition, this research might provide a paradigm for the study of other rare genomic disorders.

Complex genomic rearrangements: an underestimated cause of rare diseases

Complex genomic rearrangements (CGRs) are known contributors to disease but are often missed during routine genetic screening. Identifying CGRs requires (i) identifying copy number variants (CNVs) concurrently with inversions, (ii) phasing multiple breakpoint junctions incis, as well as (iii) detecting and resolving structural variants (SVs) within repeats. We demonstrate how combining cytogenetics and new sequencing methodologies is being successfully applied to gain insights into the genomic architecture of CGRs. In addition, we review CGR patterns and molecular features revealed by studying constitutional genomic disorders. These data offer invaluable lessons to individuals interested in investigating CGRs, evaluating their clinical relevance and frequency, as well as assessing their impact(s) on rare genetic diseases.

Optical mapping compendium of structural variants across global cattle breeds

Structural variants (SV) have been linked to important bovine disease phenotypes, but due to the difficulty of their accurate detection with standard sequencing approaches, their role in shaping important traits across cattle breeds is largely unexplored. Optical mapping is an alternative approach for mapping SVs that has been shown to have higher sensitivity than DNA sequencing approaches. The aim of this project was to use optical mapping to develop a high-quality database of structural variation across cattle breeds from different geographical regions, to enable further study of SVs in cattle. To do this we generated 100X Bionano optical mapping data for 18 cattle of nine different ancestries, three continents and both cattle sub-species. In total we identified 13,457 SVs, of which 1,200 putatively overlap coding regions. This resource provides a high-quality set of optical mapping-based SV calls that can be used across studies, from validating DNA sequencing-based SV calls to prioritising candidate functional variants in genetic association studies and expanding our understanding of the role of SVs in cattle evolution.

Measurement(s)	Optical Mapping
Technology Type(s)	Optical Mapping
Factor Type(s)	Structural variants
Sample Characteristic - Organism	Bos taurus
Sample Characteristic - Location	United Kingdom • Kenya • Zambia • Uganda • Brazil • Nigeria

Identification of a familial complex chromosomal rearrangement by optical genome mapping

BACKGROUND

Complex chromosomal rearrangements (CCRs) are rare chromosomal structural variations, containing a variety of rearrangements such as translocation, inversion and/or insertion. With the development of cytogenetic and molecular genetic techniques, some chromosomal rearrangements that were initially considered to be simple reciprocal translocations in the past might eventually involve more complex chromosomal rearrangements.

CASE PRESENTATION

In this case, a pregnant woman, who had a spontaneous abortion last year, had abnormal prenatal test results again in the second pregnancy. Applying a combination of genetic methods including karyotype analysis, chromosomal microarray analysis, fluorescence in situ hybridization and optical genome mapping confirmed that the pregnant woman was a carrier of a CCR involving three chromosomes and four breakpoints, and the CCR was paternal-origin. Her first and second pregnancy abnormalities were caused by chromosomal microdeletions and microduplications due to the malsegregations of the derivative chromosomes.

CONCLUSIONS

We presented a rare familial CCR involving three chromosomes and four breakpoints. This study provided precise and detailed information for the subsequent reproductive decision-making and genetic counselling of the patient.

Confined migration induces heterochromatin formation and alters chromatin accessibility

During migration, cells often squeeze through small constrictions, requiring extensive deformation. We hypothesized that nuclear deformation associated with such confined migration could alter chromatin organization and function. By studying cells migrating through microfluidic devices that mimic interstitial spaces in vivo, we found that confined migration results in increased H3K9me3 and H3K27me3 heterochromatin marks that persist for days. This "confined migration-induced heterochromatin" (CMiH) was distinct from heterochromatin formation during migration initiation. Confined migration decreased chromatin accessibility at intergenic regions near centromeres and telomeres, suggesting heterochromatin spreading from existing sites. Consistent with the overall decrease in accessibility, global transcription was decreased during confined migration. Intriguingly, we also identified increased accessibility at promoter regions of genes linked to chromatin silencing, tumor invasion, and DNA damage response. Inhibiting CMiH reduced migration speed, suggesting that CMiH promotes confined migration. Together, our findings indicate that confined migration induces chromatin changes that regulate cell migration and other functions.

Epigenetic tumor heterogeneity in the era of single-cell profiling with nanopore sequencing

Nanopore sequencing has brought the technology to the next generation in the science of sequencing. This is achieved through research advancing on: pore efficiency, creating mechanisms to control DNA translocation, enhancing signal-to-noise ratio, and expanding to long-read ranges. Heterogeneity regarding epigenetics would be broad as mutations in the epigenome are sensitive to cause new challenges in cancer research. Epigenetic enzymes which catalyze DNA methylation and histone modification are dysregulated in cancer cells and cause numerous heterogeneous clones to evolve. Detection of this heterogeneity in these clones plays an indispensable role in the treatment of various cancer types. With single-cell profiling, the nanopore sequencing technology could provide a simple sequence at long reads and is expected to be used soon at the bedside or doctor's office. Here, we review the advancements of nanopore sequencing and its use in the detection of epigenetic heterogeneity in cancer.

Oxford Nanopore and Bionano Genomics technologies evaluation for plant structural variation detection

Background

Structural Variations (SVs) are genomic rearrangements derived from duplication, deletion, insertion, inversion, and translocation events. In the past, SVs detection was limited to cytological approaches, then to Next-Generation Sequencing (NGS) short reads and partitioned assemblies. Nowadays, technologies such as DNA long read sequencing and optical mapping have revolutionized the understanding of SVs in genomes, due to the enhancement of the power of SVs detection.

This study aims to investigate performance of two techniques, 1) long-read sequencing obtained with the MinION device (Oxford Nanopore Technologies) and 2) optical mapping obtained with Saphyr device (Bionano Genomics) to detect and characterize SVs in the genomes of the two ecotypes of Arabidopsis thaliana, Columbia-0 (Col-0) and Landsberg erecta 1 (Ler-1).

Results

We described the SVs detected from the alignment of the best ONT assembly and DLE-1 optical maps of A. thaliana Ler-1 against the public reference genome Col-0 TAIR10.1. After filtering (SV > 1 kb), 1184 and 591 Ler-1 SVs were retained from ONT and Bionano technologies respectively. A total of 948 Ler-1 ONT SVs (80.1%) corresponded to 563 Bionano SVs (95.3%) leading to 563 common locations. The specific locations were scrutinized to assess improvement in SV detection by either technology. The ONT SVs were mostly detected near TE and gene features, and resistance genes seemed particularly impacted.

Conclusions

Structural variations linked to ONT sequencing error were removed and false positives limited, with high quality Bionano SVs being conserved. When compared with the Col-0 TAIR10.1 reference genome, most of the detected SVs discovered by both technologies were found in the same locations. ONT assembly sequence leads to more specific SVs than Bionano one, the latter being more efficient to characterize large SVs. Even if both technologies are complementary approaches, ONT data appears to be more adapted to large scale populations studies, while Bionano performs better in improving assembly and describing specificity of a genome compared to a reference.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08499-4.

Multiplex structural variant detection by whole-genome mapping and nanopore sequencing

Identification of structural variants (SVs) breakpoints is important in studying mutations, mutagenic causes, and functional impacts. Next-generation sequencing and whole-genome optical mapping are extensively used in SV discovery and characterization. However, multiple platforms and computational approaches are needed for comprehensive analysis, making it resource-intensive and expensive. Here, we propose a strategy combining optical mapping and cas9-assisted targeted nanopore sequencing to analyze SVs. Optical mapping can economically and quickly detect SVs across a whole genome but does not provide sequence-level information or precisely resolve breakpoints. Furthermore, since only a subset of all SVs is known to affect biology, we attempted to type a subset of all SVs using targeted nanopore sequencing. Using our approach, we resolved the breakpoints of five deletions, five insertions, and an inversion, in a single experiment.

Population-scale long-read sequencing uncovers transposable elements associated with gene expression variation and adaptive signatures in Drosophila

High quality reference genomes are crucial to understanding genome function, structure and evolution. The availability of reference genomes has allowed us to start inferring the role of genetic variation in biology, disease, and biodiversity conservation. However, analyses across organisms demonstrate that a single reference genome is not enough to capture the global genetic diversity present in populations. In this work, we generate 32 high-quality reference genomes for the well-known model species D. melanogaster and focus on the identification and analysis of transposable element variation as they are the most common type of structural variant. We show that integrating the genetic variation across natural populations from five climatic regions increases the number of detected insertions by 58%. Moreover, 26% to 57% of the insertions identified using long-reads were missed by short-reads methods. We also identify hundreds of transposable elements associated with gene expression variation and new TE variants likely to contribute to adaptive evolution in this species. Our results highlight the importance of incorporating the genetic variation present in natural populations to genomic studies, which is essential if we are to understand how genomes function and evolve.

Even in well-studied species, there is still substantial natural genetic variation that has not been characterized. Here, the authors use long read sequencing to discover transposable elements in the Drosophila genome not detected by short read sequencing, and link them to gene expression.

National Genome Initiatives in Europe and the United Kingdom in the Era of Whole-Genome Sequencing: A Comprehensive Review

Identification of genomic variability in population plays an important role in the clinical diagnostics of human genetic diseases. Thanks to rapid technological development in the field of massive parallel sequencing technologies, also known as next-generation sequencing (NGS), complex genomic analyses are now easier and cheaper than ever before, which consequently leads to more effective utilization of these techniques in clinical practice. However, interpretation of data from NGS is still challenging due to several issues caused by natural variability of DNA sequences in human populations. Therefore, development and realization of projects focused on description of genetic variability of local population (often called "national or digital genome") with a NGS technique is one of the best approaches to address this problem. The next step of the process is to share such data via publicly available databases. Such databases are important for the interpretation of variants with unknown significance or (likely) pathogenic variants in rare diseases or cancer or generally for identification of pathological variants in a patient's genome. In this paper, we have compiled an overview of published results of local genome sequencing projects from United Kingdom and Europe together with future plans and perspectives for newly announced ones.

Lower Exome Sequencing Coverage of Ancestrally African Patients in The Cancer Genome Atlas

BACKGROUND

In the United States, cancer disproportionately impacts Black and African American individuals. Identifying genetic factors underlying cancer disparities has been an important research focus and requires data that are equitable in both quantity and quality across racial groups. It is widely recognized that DNA databases quantitatively underrepresent minorities. However, the differences in data quality between racial groups have not been well studied.

METHODS

We compared the qualities of germline and tumor exomes between ancestrally African and European patients in The Cancer Genome Atlas of 7 cancers with at least 50 self-reported Black patients in the context of sequencing depth, tumor purity, and qualities of germline variants and somatic mutations.

RESULTS

Germline and tumor exomes from ancestrally African patients were sequenced at statistically significantly lower depth in 6 out of the 7 cancers. For 3 cancers, most ancestrally European exomes were sequenced in early sample batches at higher depth, whereas ancestrally African exomes were concentrated in later batches and sequenced at much lower depth. For the other 3 cancers, the reasons of lower sequencing coverage of ancestrally African exomes remain unknown. Furthermore, even when the sequencing depths were comparable, African exomes had disproportionally higher percentages of positions with insufficient coverage, likely because of the known European bias in the human reference genome that impacted exome capture kit design.

CONCLUSIONS

Overall and positional lower sequencing depths of ancestrally African exomes in The Cancer Genome Atlas led to underdetection and lower quality of variants, highlighting the need to consider epidemiological factors for future genomics studies.

Towards accurate and reliable resolution of structural variants for clinical diagnosis

Structural variants (SVs) are a major source of human genetic diversity and have been associated with different diseases and phenotypes. The detection of SVs is difficult, and a diverse range of detection methods and data analysis protocols has been developed. This difficulty and diversity make the detection of SVs for clinical applications challenging and requires a framework to ensure accuracy and reproducibility. Here, we discuss current developments in the diagnosis of SVs and propose a roadmap for the accurate and reproducible detection of SVs that includes case studies provided from the FDA-led SEquencing Quality Control Phase II (SEQC-II) and other consortium efforts.

Optical genome mapping identifies rare structural variations as predisposition factors associated with severe COVID-19

Impressive global efforts have identified both rare and common gene variants associated with severe COVID-19 using sequencing technologies. However, these studies lack the sensitivity to accurately detect several classes of variants, especially large structural variants (SVs), which account for a substantial proportion of genetic diversity including clinically relevant variation. We performed optical genome mapping on 52 severely ill COVID-19 patients to identify rare/unique SVs as decisive predisposition factors associated with COVID-19. We identified 7 SVs involving genes implicated in two key host-viral interaction pathways: innate immunity and inflammatory response, and viral replication and spread in nine patients, of which SVs in STK26 and DPP4 genes are the most intriguing candidates. This study is the first to systematically assess the potential role of SVs in the pathogenesis of COVID-19 severity and highlights the need to evaluate SVs along with sequencing variants to comprehensively associate genomic information with interindividual variability in COVID-19 phenotypes.

Finding Overlapping Rmaps via Clustering

Optical mapping has been largely automated, and first produces single molecule restriction maps, called Rmaps, which are assembled to generate genome wide optical maps. Since the location and orientation of each Rmap is unknown, the first problem in the analysis of this data is finding related Rmaps, i.e., pairs of Rmaps that share the same orientation and have significant overlap in their genomic location. Although heuristics for identifying related Rmaps exist, they all require quantization of the data which leads to a loss in the precision. In this paper, we propose a Gaussian mixture modelling clustering based method, which we refer to as O, that finds overlapping Rmaps without quantization. Using both simulated and real datasets, we show that OMclust substantially improves the precision (from 48.3% to 73.3%) over the state-of-the art methods while also reducing CPU time and memory consumption. Further, we integrated OMclust into the error correction methods (Elmeri and Comet) to demonstrate the increase in the performance of these methods. When OMclust was combined with Comet to error correct Rmap data generated from human DNA, it was able to error correct close to 3x more Ramps, and reduced the CPU time by more than 35x.

Emerging Insights Into Chronic Renal Disease Pathogenesis in Hypertension From Human and Animal Genomic Studies

The pathogenic links between elevated blood pressure and chronic kidney disease remain obscure. This article examines progress in population genetics and in animal models of hypertension and chronic kidney disease. It also provides a critique of the application of genome-wide association studies to understanding the heritability of renal function. Emerging themes identified indicate that heritable risk of chronic kidney disease in hypertension can arise from genetic variation in (1) glomerular and tubular protein handling mechanisms; (2) autoregulatory capacity of the renal vasculature; and (3) innate and adaptive immune mechanisms. Increased prevalence of hypertension-associated chronic kidney disease that occurs with aging may reflect amplification of heritable risks by normal aging processes affecting immunity and autoregulation.

Long-read technologies identify a hidden inverted duplication in a family with choroideremia

The lack of molecular diagnoses in rare genetic diseases can be explained by limitations of current standard genomic technologies. Upcoming long-read techniques have complementary strengths to overcome these limitations, with a particular strength in identifying structural variants. By using optical genome mapping and long-read sequencing, we aimed to identify the pathogenic variant in a large family with X-linked choroideremia. In this family, aberrant splicing of exon 12 of the choroideremia gene CHM was detected in 2003, but the underlying genomic defect remained elusive. Optical genome mapping and long-read sequencing approaches now revealed an intragenic 1,752 bp inverted duplication including exon 12 and surrounding regions, located downstream of the wild-type copy of exon 12. Both breakpoint junctions were confirmed with Sanger sequencing and segregate with the X-linked inheritance in the family. The breakpoint junctions displayed sequence microhomology suggestive for an erroneous replication mechanism as the origin of the structural variant. The inverted duplication is predicted to result in a hairpin formation of the pre-mRNA with the wild-type exon 12, leading to exon skipping in the mature mRNA. The identified inverted duplication is deemed the hidden pathogenic cause of disease in this family. Our study shows that optical genome mapping and long-read sequencing have significant potential for the identification of (hidden) structural variants in rare genetic diseases.

Application of full-genome analysis to diagnose rare monogenic disorders

Current genetic tests for rare diseases provide a diagnosis in only a modest proportion of cases. The Full-Genome Analysis method, FGA, combines long-range assembly and whole-genome sequencing to detect small variants, structural variants with breakpoint resolution, and phasing. We built a variant prioritization pipeline and tested FGA’s utility for diagnosis of rare diseases in a clinical setting. FGA identified structural variants and small variants with an overall diagnostic yield of 40% (20 of 50 cases) and 35% in exome-negative cases (8 of 23 cases), 4 of these were structural variants. FGA detected and mapped structural variants that are missed by short reads, including non-coding duplication, and phased variants across long distances of more than 180 kb. With the prioritization algorithm, longer DNA technologies could replace multiple tests for monogenic disorders and expand the range of variants detected. Our study suggests that genomes produced from technologies like FGA can improve variant detection and provide higher resolution genome maps for future application.

Characterization of full-length LINE-1 insertions in 154 genomes

Long interspersed nuclear elements (LINEs) are retrotransposons that contribute to genetic variation in the human genome. LINE-1 elements in larger-scale studies are challenging to identify using sequencing technologies due to cost and scalability. We developed an approach using optical mapping for detection of full-length LINE-1 insertions and 10× sequencing for confirmation. We found 51 true positive full-length LINE-1 insertions, of which 4 are novel insertions, in NA12878. Repeating our analysis on a larger sample set representing 26 populations, we identified 329 full-length LINE-1 elements, of which 123 are novel. 24.8% of these 329 LINE-1 insertions were shared amongst all 5 superpopulations (AFR, AMR, EUR, EAS, SAS). The African superpopulation has a higher percentage of population-specific LINE-1 insertions than any other superpopulation. These data indicate that our approach can provide high-speed, cost-effective, and increased accuracy for LINE-1 detection. These data also provide an insight into variations of LINE-1 elements between different populations.

Multiplatform discovery and regulatory function analysis of structural variations in non-small cell lung carcinoma

Non-small cell lung carcinoma (NSCLC), the most common form of lung cancer, is the leading cause of cancer-related death worldwide. We perform whole-genome sequencing (WGS) on samples from 43 primary patients with NSCLC and matched normal samples and analyze their matched open chromatin data and transcriptome data. Our results indicate that next-generation sequencing (NGS) and the Bionano Genomics (BNG) platform should be viewed as complementary technologies in terms of structural variations detection. By creating a framework integrating these two platforms, we detect high-technical-confidence somatic structural variations (SVs) in NSCLC cases, which could aid in the efficient investigation of new candidate oncogenes, such as TRIO and SESTD1. Our findings highlight the impact of somatic SVs on NSCLC oncogenesis and lay a foundation for exploring associations among somatic SVs, gene expression, and regulatory networks in patients with NSCLC.

Genomic regions associated with microdeletion/microduplication syndromes exhibit extreme diversity of structural variation

Segmental duplications (SDs) are a class of long, repetitive DNA elements whose paralogs share a high level of sequence similarity with each other. SDs mediate chromosomal rearrangements that lead to structural variation in the general population as well as genomic disorders associated with multiple congenital anomalies, including the 7q11.23 (Williams-Beuren Syndrome, WBS), 15q13.3, and 16p12.2 microdeletion syndromes. Population-level characterization of SDs has generally been lacking because most techniques used for analyzing these complex regions are both labor and cost intensive. In this study, we have used a high-throughput technique to genotype complex structural variation with a single molecule, long-range optical mapping approach. We characterized SDs and identified novel structural variants (SVs) at 7q11.23, 15q13.3, and 16p12.2 using optical mapping data from 154 phenotypically normal individuals from 26 populations comprising five super-populations. We detected several novel SVs for each locus, some of which had significantly different prevalence between populations. Additionally, we localized the microdeletion breakpoints to specific paralogous duplicons located within complex SDs in two patients with WBS, one patient with 15q13.3, and one patient with 16p12.2 microdeletion syndromes. The population-level data presented here highlights the extreme diversity of large and complex SVs within SD-containing regions. The approach we outline will greatly facilitate the investigation of the role of inter-SD structural variation as a driver of chromosomal rearrangements and genomic disorders.

Towards improved genetic diagnosis of human differences of sex development

Despite being collectively among the most frequent congenital developmental conditions worldwide, differences of sex development (DSD) lack recognition and research funding. As a result, what constitutes optimal management remains uncertain. Identification of the individual conditions under the DSD umbrella is challenging and molecular genetic diagnosis is frequently not achieved, which has psychosocial and health-related repercussions for patients and their families. New genomic approaches have the potential to resolve this impasse through better detection of protein-coding variants and ascertainment of under-recognized aetiology, such as mosaic, structural, non-coding or epigenetic variants. Ultimately, it is hoped that better outcomes data, improved understanding of the molecular causes and greater public awareness will bring an end to the stigma often associated with DSD.

Towards population-scale long-read sequencing

Long-read sequencing technologies have now reached a level of accuracy and yield that allows their application to variant detection at a scale of tens to thousands of samples. Concomitant with the development of new computational tools, the first population-scale studies involving long-read sequencing have emerged over the past 2 years and, given the continuous advancement of the field, many more are likely to follow. In this Review, we survey recent developments in population-scale long-read sequencing, highlight potential challenges of a scaled-up approach and provide guidance regarding experimental design. We provide an overview of current long-read sequencing platforms, variant calling methodologies and approaches for de novo assemblies and reference-based mapping approaches. Furthermore, we summarize strategies for variant validation, genotyping and predicting functional impact and emphasize challenges remaining in achieving long-read sequencing at a population scale.

Expanding the conservation genomics toolbox: Incorporating structural variants to enhance genomic studies for species of conservation concern

Structural variants (SVs) are large rearrangements (>50 bp) within the genome that impact gene function and the content and structure of chromosomes. As a result, SVs are a significant source of functional genomic variation, that is, variation at genomic regions underpinning phenotype differences, that can have large effects on individual and population fitness. While there are increasing opportunities to investigate functional genomic variation in threatened species via single nucleotide polymorphism (SNP) data sets, SVs remain understudied despite their potential influence on fitness traits of conservation interest. In this future-focused Opinion, we contend that characterizing SVs offers the conservation genomics community an exciting opportunity to complement SNP-based approaches to enhance species recovery. We also leverage the existing literature-predominantly in human health, agriculture and ecoevolutionary biology-to identify approaches for readily characterizing SVs and consider how integrating these into the conservation genomics toolbox may transform the way we manage some of the world's most threatened species.