Fine-scale structural variation of the human genome

Genome-wide detection of structural variation in some sheep breeds using whole-genome long-read sequencing data

Genomic structural variants (SVs) constitute a significant proportion of genetic variation in the genome. The rapid development of long-reads sequencing has facilitated the detection of long-fragment SVs. There is no published study to detect SVs using long-read data from sheep. We applied a long-read mapping approach to detect SVs and characterized a total of 30,771 insertions, deletions, inversions and translocations. We identified 716, 916, 842 and 303 specific SVs in Southdown sheep, Alpine merino sheep, Qilian White Tibetan sheep and Oula sheep, respectively. We annotated these SVs and found that these SV-related genes were primarily enriched in the well-established pathways involved in the regulation of the immune system, growth and development and environmental adaptability. We detected and annotated SVs based on NGS resequencing data to validate the accuracy based on third-generation detection. Moreover, five candidate SVs were verified using the PCR method in 50 sheep. Our study is the first to use a long-reads sequencing approach to construct a novel structural variation map in sheep. We have completed a preliminary exploration of the potential effects of SVs on sheep.

The chromosome-level assembly of the wild diploid alfalfa genome provides insights into the full landscape of genomic variations between cultivated and wild alfalfa

Alfalfa (Medicago sativa L.) is one of the most important forage legumes in the world, including autotetraploid (M. sativa ssp. sativa) and diploid alfalfa (M. sativa ssp. caerulea, progenitor of autotetraploid alfalfa). Here, we reported a high-quality genome of ZW0012 (diploid alfalfa, 769 Mb, contig N50 = 5.5 Mb), which was grouped into the Northern group in population structure analysis, suggesting that our genome assembly filled a major gap among the members of M. sativa complex. During polyploidization, large phenotypic differences occurred between diploids and tetraploids, and the genetic information underlying its massive phenotypic variations remains largely unexplored. Extensive structural variations (SVs) were identified between ZW0012 and XinJiangDaYe (an autotetraploid alfalfa with released genome). We identified 71 ZW0012-specific PAV genes and 1296 XinJiangDaYe-specific PAV genes, mainly involved in defence response, cell growth, and photosynthesis. We have verified the positive roles of MsNCR1 (a XinJiangDaYe-specific PAV gene) in nodulation using an Agrobacterium rhizobia-mediated transgenic method. We also demonstrated that MsSKIP23_1 and MsFBL23_1 (two XinJiangDaYe-specific PAV genes) regulated leaf size by transient overexpression and virus-induced gene silencing analysis. Our study provides a high-quality reference genome of an important diploid alfalfa germplasm and a valuable resource of variation landscape between diploid and autotetraploid, which will facilitate the functional gene discovery and molecular-based breeding for the cultivars in the future.

Mutational landscape of HSP family on human breast cancer

Breast cancer (BRCA) is a prevalent malignancy with the highest incidence among females. BRCA can be categorized into five intrinsic molecular subtypes (LumA, LumB, HER2, Basal, and Normal), each characterized by varying molecular and clinical features determined by the expression of intrinsic genes (PAM50). The Heat Shock Protein (HSP) family is composed of 95 genes evolutionary conservated, they have critical roles in proteostasis in both normal and cancerous processes. Many studies have linked HSP to the development and spread of cancer. They modulate the activity of multiple proteins expressed by oncogenes and anti-oncogenes through a range of interactions. In this study, we evaluate the mutational changes that HSP undergoes in BRCA mainly from the TCGA database. We observe that Copy Number Variations (CNV) are the more frequent events analyzed surpassing the occurrence of point mutations, indels, and translation start site mutations. The Basal subtype showcased the highest count of amplified CNV, including subtype-specific changes, whereas the Luminals tumors accumulated the greatest number of deletion CNV. Meanwhile, the HER2 subtype exhibited a comparatively lower frequency of CNV alterations when compared to the other subtypes. This study integrates CNV and expression data, finding associations between these two variables and the influence of CNV on the deregulation of HSP expression. To enhance the role of HSP as a risk predictor in BRCA, we succeeded in identifying CNV profiles as a prognostic marker. We included Artificial Intelligence to improve the clustering of patients, and we achieved a molecular CNV signature as a significant risk factor independent of known classic markers, including molecular subtypes PAM50. This research enhances the comprehension of HSP DNA alterations in BRCA and its relation with predicting the risk of affected individuals providing insights to develop guide personalized treatment strategies.

A comparative study of structural variant calling in WGS from Alzheimer's disease families

Detecting structural variants (SVs) in whole-genome sequencing poses significant challenges. We present a protocol for variant calling, merging, genotyping, sensitivity analysis, and laboratory validation for generating a high-quality SV call set in whole-genome sequencing from the Alzheimer's Disease Sequencing Project comprising 578 individuals from 111 families. Employing two complementary pipelines, Scalpel and Parliament, for SV/indel calling, we assessed sensitivity through sample replicates (N = 9) with in silico variant spike-ins. We developed a novel metric, D-score, to evaluate caller specificity for deletions. The accuracy of deletions was evaluated by Sanger sequencing. We generated a high-quality call set of 152,301 deletions of diverse sizes. Sanger sequencing validated 114 of 146 detected deletions (78.1%). Scalpel excelled in accuracy for deletions ≤100 bp, whereas Parliament was optimal for deletions >900 bp. Overall, 83.0% and 72.5% of calls by Scalpel and Parliament were validated, respectively, including all 11 deletions called by both Parliament and Scalpel between 101 and 900 bp. Our flexible protocol successfully generated a high-quality deletion call set and a truth set of Sanger sequencing-validated deletions with precise breakpoints spanning 1-17,000 bp.

Single variant, yet "double trouble": TSC and KBG syndrome because of a large de novo inversion

Despite the advances in high-throughput sequencing, many rare disease patients remain undiagnosed. In particular, the patients with well-defined clinical phenotypes and established clinical diagnosis, yet missing or partial genetic diagnosis, may hold a clue to more complex genetic mechanisms of a disease that could be missed by available clinical tests. Here, we report a patient with a clinical diagnosis of Tuberous sclerosis, combined with unusual secondary features, but negative clinical tests including TSC1 and TSC2 Short-read whole-genome sequencing combined with advanced bioinformatics analyses were successful in uncovering a de novo pericentric 87-Mb inversion with breakpoints in TSC2 and ANKRD11, which explains the TSC clinical diagnosis, and confirms a second underlying monogenic disorder, KBG syndrome. Our findings illustrate how complex variants, such as large inversions, may be missed by clinical tests and further highlight the importance of well-defined clinical diagnoses in uncovering complex molecular mechanisms of a disease, such as complex variants and "double trouble" effects.

Genetic interrogation for sequence and copy number variants in systemic lupus erythematosus

Early-onset systemic lupus erythematosus presents with a more severe disease and is associated with a greater genetic burden, especially in patients from Black, Asian or Hispanic ancestries. Next-generation sequencing techniques, notably whole exome sequencing, have been extensively used in genomic interrogation studies to identify causal disease variants that are increasingly implicated in the development of autoimmunity. This Review discusses the known casual variants of polygenic and monogenic systemic lupus erythematosus and its implications under certain genetic disparities while suggesting an age-based sequencing strategy to aid in clinical diagnostics and patient management for improved patient care.

Methods of Detection and Mechanisms of Origin of Complex Structural Genome Variations

Based on classical karyotyping, structural genome variations (SVs) have generally been considered to be either "simple" (with one or two breakpoints) or "complex" (with more than two breakpoints). Studying the breakpoints of SVs at nucleotide resolution revealed additional, subtle structural variations, such that even "simple" SVs turned out to be "complex." Genome-wide sequencing methods, such as fosmid and paired-end mapping, short-read and long-read whole genome sequencing, and single-molecule optical mapping, also indicated that the number of SVs per individual was considerably larger than expected from karyotyping and high-resolution chromosomal array-based studies. Interestingly, SVs were detected in studies of cohorts of individuals without clinical phenotypes. The common denominator of all SVs appears to be a failure to accurately repair DNA double-strand breaks (DSBs) or to halt cell cycle progression if DSBs persist. This review discusses the various DSB response mechanisms during the mitotic cell cycle and during meiosis and their regulation. Emphasis is given to the molecular mechanisms involved in the formation of translocations, deletions, duplications, and inversions during or shortly after meiosis I. Recently, CRISPR-Cas9 studies have provided unexpected insights into the formation of translocations and chromothripsis by both breakage-fusion-bridge and micronucleus-dependent mechanisms.

Viral oncogenes, viruses, and cancer: a third-generation sequencing perspective on viral integration into the human genome

The link between viruses and cancer has intrigued scientists for decades. Certain viruses have been shown to be vital in the development of various cancers by integrating viral DNA into the host genome and activating viral oncogenes. These viruses include the Human Papillomavirus (HPV), Hepatitis B and C Viruses (HBV and HCV), Epstein-Barr Virus (EBV), and Human T-Cell Leukemia Virus (HTLV-1), which are all linked to the development of a myriad of human cancers. Third-generation sequencing technologies have revolutionized our ability to study viral integration events at unprecedented resolution in recent years. They offer long sequencing capabilities along with the ability to map viral integration sites, assess host gene expression, and track clonal evolution in cancer cells. Recently, researchers have been exploring the application of Oxford Nanopore Technologies (ONT) nanopore sequencing and Pacific BioSciences (PacBio) single-molecule real-time (SMRT) sequencing in cancer research. As viral integration is crucial to the development of cancer via viruses, third-generation sequencing would provide a novel approach to studying the relationship interlinking viral oncogenes, viruses, and cancer. This review article explores the molecular mechanisms underlying viral oncogenesis, the role of viruses in cancer development, and the impact of third-generation sequencing on our understanding of viral integration into the human genome.

Incorporating CNV analysis improves the yield of exome sequencing for rare monogenic disorders-an important consideration for resource-constrained settings

Exome sequencing (ES) is a recommended first-tier diagnostic test for many rare monogenic diseases. It allows for the detection of both single-nucleotide variants (SNVs) and copy number variants (CNVs) in coding exonic regions of the genome in a single test, and this dual analysis is a valuable approach, especially in limited resource settings. Single-nucleotide variants are well studied; however, the incorporation of copy number variant analysis tools into variant calling pipelines has not been implemented yet as a routine diagnostic test, and chromosomal microarray is still more widely used to detect copy number variants. Research shows that combined single and copy number variant analysis can lead to a diagnostic yield of up to 58%, increasing the yield with as much as 18% from the single-nucleotide variant only pipeline. Importantly, this is achieved with the consideration of computational costs only, without incurring any additional sequencing costs. This mini review provides an overview of copy number variant analysis from exome data and what the current recommendations are for this type of analysis. We also present an overview on rare monogenic disease research standard practices in resource-limited settings. We present evidence that integrating copy number variant detection tools into a standard exome sequencing analysis pipeline improves diagnostic yield and should be considered a significantly beneficial addition, with relatively low-cost implications. Routine implementation in underrepresented populations and limited resource settings will promote generation and sharing of CNV datasets and provide momentum to build core centers for this niche within genomic medicine.

Systematic evaluation of multiple NGS platforms for structural variants detection

Structural variations (SV) are critical genome changes affecting human diseases. Although many hybridization-based methods exist, evaluating SVs through next-generation sequencing (NGS) data is still necessary for broader research exploration. Here, we comprehensively compared the performance of 16 SV callers and multiple NGS platforms using NA12878 whole genome sequencing (WGS) datasets. The results indicated that several SV callers performed well relatively, such as Manta, GRIDSS, LUMPY, TARDIS, FermiKit, and Wham. Meanwhile, all NGS platforms have a similar performance using a single software. Additionally, we found that the source of undetected SVs was mostly from long reads datasets, therefore, the more appropriate strategy for accurate SV detection will be an integration of long and shorter reads in the future. At present, in the period of NGS as a mainstream method in bioinformatics, our study would provide helpful and comprehensive guidelines for specific categories of SV research.

Inversions on human chromosomes

Human chromosome inversions are types of balanced structural variations, making them difficult to analyze. Thanks to PEM (paired-end sequencing and mapping), there has been tremendous progress in studying inversions. Inversions play an important role as an evolutionary factor, contributing to the formation of gonosomes, speciation of chimpanzees and humans, and inv17q21.3 or inv8p23.1 exhibit the features of natural selection. Both inversions have been related to pathogenic phenotype by directly affecting a gene structure (e.g., inv5p15.1q14.1), regulating gene expression (e.g., inv7q21.3q35) and by predisposing to other secondary arrangements (e.g., inv7q11.23). A polymorphism of human inversions is documented by the InvFEST database (a database that stores information about clinical predictions, validations, frequency of inversions, etc.), but only a small fraction of these inversions is validated, and a detailed analysis is complicated by the frequent location of breakpoints within regions of repetitive sequences.

16p11.2 deletion accelerates subpallial maturation and increases variability in human iPSC-derived ventral telencephalic organoids

Inhibitory interneurons regulate cortical circuit activity, and their dysfunction has been implicated in autism spectrum disorder (ASD). 16p11.2 microdeletions are genetically linked to 1% of ASD cases. However, few studies investigate the effects of this microdeletion on interneuron development. Using ventral telencephalic organoids derived from human induced pluripotent stem cells, we have investigated the effect of this microdeletion on organoid size, progenitor proliferation and organisation into neural rosettes, ganglionic eminence marker expression at early developmental timepoints, and expression of the neuronal marker NEUN at later stages. At early stages, deletion organoids exhibited greater variations in size with concomitant increases in relative neural rosette area and the expression of the ventral telencephalic marker COUPTFII, with increased variability in these properties. Cell cycle analysis revealed an increase in total cell cycle length caused primarily by an elongated G1 phase, the duration of which also varied more than normal. At later stages, deletion organoids increased their NEUN expression. We propose that 16p11.2 microdeletions increase developmental variability and may contribute to ASD aetiology by lengthening the cell cycle of ventral progenitors, promoting premature differentiation into interneurons.

Characterization of the immunoglobulin lambda chain locus from diverse populations reveals extensive genetic variation

Immunoglobulins (IGs), crucial components of the adaptive immune system, are encoded by three genomic loci. However, the complexity of the IG loci severely limits the effective use of short read sequencing, limiting our knowledge of population diversity in these loci. We leveraged existing long read whole-genome sequencing (WGS) data, fosmid technology, and IG targeted single-molecule, real-time (SMRT) long-read sequencing (IG-Cap) to create haplotype-resolved assemblies of the IG Lambda (IGL) locus from 6 ethnically diverse individuals. In addition, we generated 10 diploid assemblies of IGL from a diverse cohort of individuals utilizing IG-Cap. From these 16 individuals, we identified significant allelic diversity, including 36 novel IGLV alleles. In addition, we observed highly elevated single nucleotide variation (SNV) in IGLV genes relative to IGL intergenic and genomic background SNV density. By comparing SNV calls between our high quality assemblies and existing short read datasets from the same individuals, we show a high propensity for false-positives in the short read datasets. Finally, for the first time, we nucleotide-resolved common 5-10 Kb duplications in the IGLC region that contain functional IGLJ and IGLC genes. Together these data represent a significant advancement in our understanding of genetic variation and population diversity in the IGL locus.

Deepening the understanding of CNVs on chromosome 15q11-13 by using hiPSCs: An overview

The human α7 neuronal nicotinic acetylcholine receptor gene (CHRNA7) is widely expressed in the central and peripheral nervous systems. This receptor is implicated in both brain development and adult neurogenesis thanks to its ability to mediate acetylcholine stimulus (Ach). Copy number variations (CNVs) of CHRNA7 gene have been identified in humans and are genetically linked to cognitive impairments associated with multiple disorders, including schizophrenia, bipolar disorder, epilepsy, Alzheimer's disease, and others. Currently, α7 receptor analysis has been commonly performed in animal models due to the impossibility of direct investigation of the living human brain. But the use of model systems has shown that there are very large differences between humans and mice when researchers must study the CNVs and, in particular, the CNV of chromosome 15q13.3 where the CHRNA7 gene is present. In fact, human beings present genomic alterations as well as the presence of genes of recent origin that are not present in other model systems as well as they show a very heterogeneous symptomatology that is associated with both their genetic background and the environment where they live. To date, the induced pluripotent stem cells, obtained from patients carrying CNV in CHRNA7 gene, are a good in vitro model for studying the association of the α7 receptor to human diseases. In this review, we will outline the current state of hiPSCs technology applications in neurological diseases caused by CNVs in CHRNA7 gene. Furthermore, we will discuss some weaknesses that emerge from the overall analysis of the published articles.

Extrachromosomal DNA in Cancer

In cancer, complex genome rearrangements and other structural alterations, including the amplification of oncogenes on circular extrachromosomal DNA (ecDNA) elements, drive the formation and progression of tumors. ecDNA is a particularly challenging structural alteration. By untethering oncogenes from chromosomal constraints, it elevates oncogene copy number, drives intratumoral genetic heterogeneity, promotes rapid tumor evolution, and results in treatment resistance. The profound changes in DNA shape and nuclear architecture generated by ecDNA alter the transcriptional landscape of tumors by catalyzing new types of regulatory interactions that do not occur on chromosomes. The current suite of tools for interrogating cancer genomes is well suited for deciphering sequence but has limited ability to resolve the complex changes in DNA structure and dynamics that ecDNA generates. Here, we review the challenges of resolving ecDNA form and function and discuss the emerging tool kit for deciphering ecDNA architecture and spatial organization, including what has been learned to date about how this dramatic change in shape alters tumor development, progression, and drug resistance.

DNA Methylation Episignatures in Neurodevelopmental Disorders Associated with Large Structural Copy Number Variants: Clinical Implications

Large structural chromosomal deletions and duplications, referred to as copy number variants (CNVs), play a role in the pathogenesis of neurodevelopmental disorders (NDDs) through effects on gene dosage. This review focuses on our current understanding of genomic disorders that arise from large structural chromosome rearrangements in patients with NDDs, as well as difficulties in overlap of clinical presentation and molecular diagnosis. We discuss the implications of epigenetics, specifically DNA methylation (DNAm), in NDDs and genomic disorders, and consider the implications and clinical impact of copy number and genomic DNAm testing in patients with suspected genetic NDDs. We summarize evidence of global methylation episignatures in CNV-associated disorders that can be used in the diagnostic pathway and may provide insights into the molecular pathogenesis of genomic disorders. Finally, we discuss the potential for combining CNV and DNAm assessment into a single diagnostic assay.

Recombination, selection and the evolution of tandem gene arrays

Multigene families—immunity genes or sensory receptors, for instance—are often subject to diversifying selection. Allelic diversity may be favored not only through balancing or frequency-dependent selection at individual loci but also by associating different alleles in multicopy gene families. Using a combination of analytical calculations and simulations, we explored a population genetic model of epistatic selection and unequal recombination, where a trade-off exists between the benefit of allelic diversity and the cost of copy abundance. Starting from the neutral case, where we showed that gene copy number is Gamma distributed at equilibrium, we derived also the mean and shape of the limiting distribution under selection. Considering a more general model, which includes variable population size and population substructure, we explored by simulations mean fitness and some summary statistics of the copy number distribution. We determined the relative effects of selection, recombination, and demographic parameters in maintaining allelic diversity and shaping the mean fitness of a population. One way to control the variance of copy number is by lowering the rate of unequal recombination. Indeed, when encoding recombination by a rate modifier locus, we observe exactly this prediction. Finally, we analyzed the empirical copy number distribution of 3 genes in human and estimated recombination and selection parameters of our model.

Detection of Structural Variations and Fusion Genes in Breast Cancer Samples Using Third-Generation Sequencing

Background: Structural variations (SVs) are common genetic alterations in the human genome that could cause different phenotypes and diseases, including cancer. However, the detection of structural variations using the second-generation sequencing was limited by its short read length, which restrained our understanding of structural variations.

Methods: In this study, we developed a 28-gene panel for long-read sequencing and employed it to Oxford Nanopore Technologies and Pacific Biosciences platforms. We analyzed structural variations in the 28 breast cancer-related genes through long-read genomic and transcriptomic sequencing of tumor, para-tumor, and blood samples in 19 breast cancer patients.

Results: Our results showed that some somatic SVs were recurring among the selected genes, though the majority of them occurred in the non-exonic region. We found evidence supporting the existence of hotspot regions for SVs, which extended our previous understanding that they exist only for single nucleotide variations.

Conclusion: In conclusion, we employed long-read genomic and transcriptomic sequencing to identify SVs from breast cancer patients and proved that this approach holds great potential in clinical application.

Copy number alteration signatures as biomarkers in cancer: a review

Within certain cancers, extensive copy number alterations (CNAs) contribute to a complex and heterogenic genomic profile. This makes it difficult to understand and unravel the distinct molecular dynamics shaping the disease while preventing clinically effective patient stratification. CNA signature analysis represents a novel genomic stratification tool for probing this complexity, offering an intricate framework for deriving CNA patterns at the molecular level. This allows the underlying genomic mechanisms of specific cancers to be revealed, leading to the potential identification of therapeutic targets and prognostic associations. This review outlines the molecular and methodological basis of CNA signatures and focuses on recent advances highlighting their clinical utility, limitations and prospective future as novel diagnostic and prognostic cancer biomarkers.

Frequency and clinical significance of chromosomal inversions prenatally diagnosed by second trimester amniocentesis

To compare the frequency and clinical significance of familial and de novo chromosomal inversions during prenatal diagnosis. This was a retrospective study of inversions diagnosed prenatally in an Asian population by applying conventional GTG-banding to amniocyte cultures. Data from 2005 to 2019 were extracted from a single-center laboratory database. The types, frequencies, and inheritance patterns of multiple inversions were analyzed. Pericentric variant inversions of chromosome 9 or Y were excluded. In total, 56 (0.27%) fetuses with inversions were identified in the 15-year database of 21,120 confirmative diagnostic procedures. Pericentric and paracentric inversions accounted for 62.5% (35/56) and 37.5% of the inversions, respectively. Familial inversions accounted for nearly 90% of cases, and de novo mutation was identified in two pericentric and two paracentric cases. Inversions were most frequently identified on chromosomes 1 and 2 (16.1% of all inversions), followed by chromosomes 6, 7, and 10 (8.9% of all cases). The indications for invasive testing were as follows: advanced maternal age (67.3%), abnormal ultrasound findings (2.1%), abnormal serum aneuploidy screening (20.4%), and other indications (10.2%). The mode of inheritance was available for 67.9% of cases (38/56), with 89.5% of inversions being inherited (34/38). A slight preponderance of inheritance in female fetuses was observed. Three patients with inherited inversions opted for termination (two had severe central nervous system lesions and one had thalassemia major). Gestation continued for 53 fetuses, who exhibited no structural defects at birth or significant developmental problems a year after birth. Our study indicates that approximately 90% of prenatally diagnosed inversions involve familial inheritance, are spreading, and behave like founder effect mutations in this isolated population on an island. This finding can help to alleviate anxiety during prenatal counseling, which further underscores the importance of parental chromosomal analysis, further genetic studies, and appropriate counseling in cases where a nonfamilial inversion is diagnosed.

Effects of Long-Term In Vitro Expansion on Genetic Stability and Tumor Formation Capacity of Stem Cells

Stem cell therapeutics are emerging as novel alternative treatments for various neurodegenerative diseases based on their regenerative potentials. However, stem cell transplantation might have side effects such as tumor formation that limit their clinical applications. Especially, in vitro expansion of stem cells might provoke genetic instability and tumorigenic potential. To address this issue, we analyzed genomic alterations of adult human multipotent neural cells (ahMNCs), a type of human adult neural stem cells, after a long-term in vitro culture process (passage 15) using sensitive analysis techniques including karyotyping, array comparative genomic hybridization (aCGH), and whole exome sequencing (WES). Although karyotyping did not find any major abnormalities in chromosomal number or structure, diverse copy number variations (CNVs) and genetic mutations were detected by aCGH and WES in all five independent ahMNCs. However, the number of CNVs and genetic mutations did not increase and many of them did not persist as in vitro culture progressed. Although most observed CNVs and genetic mutations were not shared by all five ahMNCs, nonsynonymous missense mutations at MUC4 were found in three out of five long-term cultured ahMNC lines. The genetic instability did not confer in vivo tumorigenic potential to ahMNCs. Collectively, these results indicate that, although genetic instability can be induced by long-term in vitro expansion of stem cells, it is not sufficient to fully exert tumor formation capacity of stem cells. Other functional effects of such genetic instability need to be further elucidated.

DeepCNV: a deep learning approach for authenticating copy number variations

Copy number variations (CNVs) are an important class of variations contributing to the pathogenesis of many disease phenotypes. Detecting CNVs from genomic data remains difficult, and the most currently applied methods suffer from an unacceptably high false positive rate. A common practice is to have human experts manually review original CNV calls for filtering false positives before further downstream analysis or experimental validation. Here, we propose DeepCNV, a deep learning-based tool, intended to replace human experts when validating CNV calls, focusing on the calls made by one of the most accurate CNV callers, PennCNV. The sophistication of the deep neural network algorithm is enriched with over 10 000 expert-scored samples that are split into training and testing sets. Variant confidence, especially for CNVs, is a main roadblock impeding the progress of linking CNVs with the disease. We show that DeepCNV adds to the confidence of the CNV calls with an optimal area under the receiver operating characteristic curve of 0.909, exceeding other machine learning methods. The superiority of DeepCNV was also benchmarked and confirmed using an experimental wet-lab validation dataset. We conclude that the improvement obtained by DeepCNV results in significantly fewer false positive results and failures to replicate the CNV association results.

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.

Megabase-scale presence-absence variation with Tripsacum origin was under selection during maize domestication and adaptation

BACKGROUND

Structural variants (SVs) significantly drive genome diversity and environmental adaptation for diverse species. Unlike the prevalent small SVs (< kilobase-scale) in higher eukaryotes, large-size SVs rarely exist in the genome, but they function as one of the key evolutionary forces for speciation and adaptation.

RESULTS

In this study, we discover and characterize several megabase-scale presence-absence variations (PAVs) in the maize genome. Surprisingly, we identify a 3.2 Mb PAV fragment that shows high integrity and is present as complete presence or absence in the natural diversity panel. This PAV is embedded within the nucleolus organizer region (NOR), where the suppressed recombination is found to maintain the PAV against the evolutionary variation. Interestingly, by analyzing the sequence of this PAV, we not only reveal the domestication trace from teosinte to modern maize, but also the footprints of its origin from Tripsacum, shedding light on a previously unknown contribution from Tripsacum to the speciation of Zea species. The functional consequence of the Tripsacum segment migration is also investigated, and environmental fitness conferred by the PAV may explain the whole segment as a selection target during maize domestication and improvement.

CONCLUSIONS

These findings provide a novel perspective that Tripsacum contributes to Zea speciation, and also instantiate a strategy for evolutionary and functional analysis of the "fossil" structure variations during genome evolution and speciation.

Discovery of genomic variation across a generation

Over the past 30 years (the timespan of a generation), advances in genomics technologies have revealed tremendous and unexpected variation in the human genome and have provided increasingly accurate answers to long-standing questions of how much genetic variation exists in human populations and to what degree the DNA complement changes between parents and offspring. Tracking the characteristics of these inherited and spontaneous (or de novo) variations has been the basis of the study of human genetic disease. From genome-wide microarray and next-generation sequencing scans, we now know that each human genome contains over 3 million single nucleotide variants when compared with the ~ 3 billion base pairs in the human reference genome, along with roughly an order of magnitude more DNA—approximately 30 megabase pairs (Mb)—being ‘structurally variable’, mostly in the form of indels and copy number changes. Additional large-scale variations include balanced inversions (average of 18 Mb) and complex, difficult-to-resolve alterations. Collectively, ~1% of an individual’s genome will differ from the human reference sequence. When comparing across a generation, fewer than 100 new genetic variants are typically detected in the euchromatic portion of a child’s genome. Driven by increasingly higher-resolution and higher-throughput sequencing technologies, newer and more accurate databases of genetic variation (for instance, more comprehensive structural variation data and phasing of combinations of variants along chromosomes) of worldwide populations will emerge to underpin the next era of discovery in human molecular genetics.

The neuropathology of autism: A systematic review of post-mortem studies of autism and related disorders

Post-mortem studies allow for the direct investigation of brain tissue in those with autism and related disorders. Several review articles have focused on aspects of post-mortem abnormalities but none has brought together the entire post-mortem literature. Here, we systematically review the evidence from post-mortem studies of autism, and of related disorders that present with autistic features. The literature consists of a small body of studies with small sample sizes, but several remarkably consistent findings are evident. Cortical layering is largely undisturbed, but there are consistent reductions in minicolumn numbers and aberrant myelination. Transcriptomics repeatedly implicate abberant synaptic, metabolic, proliferation, apoptosis and immune pathways. Sufficient replicated evidence is available to implicate non-coding RNA, aberrant epigenetic profiles, GABAergic, glutamatergic and glial dysfunction in autism pathogenesis. Overall, the cerebellum and frontal cortex are most consistently implicated, sometimes revealing distinct region-specific alterations. The literature on related disorders such as Rett syndrome, Fragile X and copy number variations (CNVs) predisposing to autism is particularly small and inconclusive. Larger studies, matched for gender, developmental stage, co-morbidities and drug treatment are required.

Nucleic Acid Tests for Clinical Translation

Nucleic acids, including deoxyribonucleic acid (DNA) and ribonucleic acid (RNA), are natural biopolymers composed of nucleotides that store, transmit, and express genetic information. Overexpressed or underexpressed as well as mutated nucleic acids have been implicated in many diseases. Therefore, nucleic acid tests (NATs) are extremely important. Inspired by intracellular DNA replication and RNA transcription, in vitro NATs have been extensively developed to improve the detection specificity, sensitivity, and simplicity. The principles of NATs can be in general classified into three categories: nucleic acid hybridization, thermal-cycle or isothermal amplification, and signal amplification. Driven by pressing needs in clinical diagnosis and prevention of infectious diseases, NATs have evolved to be a rapidly advancing field. During the past ten years, an explosive increase of research interest in both basic research and clinical translation has been witnessed. In this review, we aim to provide comprehensive coverage of the progress to analyze nucleic acids, use nucleic acids as recognition probes, construct detection devices based on nucleic acids, and utilize nucleic acids in clinical diagnosis and other important fields. We also discuss the new frontiers in the field and the challenges to be addressed.

AID and APOBECs as Multifaceted Intrinsic Virus-Restricting Factors: Emerging Concepts in the Light of COVID-19

The AID (activation-induced cytidine deaminase)/APOBEC (apolipoprotein B mRNA editing enzyme catalytic subunit) family with its multifaceted mode of action emerges as potent intrinsic host antiviral system that acts against a variety of DNA and RNA viruses including coronaviruses. All family members are cytosine-to-uracil deaminases that either have a profound role in driving a strong and specific humoral immune response (AID) or restricting the virus itself by a plethora of mechanisms (APOBECs). In this article, we highlight some of the key aspects apparently linking the AID/APOBECs and SARS-CoV-2. Among those is our discovery that APOBEC4 shows high expression in cell types and anatomical parts targeted by SARS-CoV-2. Additional focus is given by us to the lymphoid structures and AID as the master regulator of germinal center reactions, which result in antibody production by plasma and memory B cells. We propose the dissection of the AID/APOBECs gene signature towards decisive determinants of the patient-specific and/or the patient group-specific antiviral response. Finally, the patient-specific mapping of the AID/APOBEC polymorphisms should be considered in the light of COVID-19.

New Insights into Potocki-Shaffer Syndrome: Report of Two Novel Cases and Literature Review

Potocki-Shaffer syndrome (PSS) is a rare non-recurrent contiguous gene deletion syndrome involving chromosome 11p11.2. Current literature implies a minimal region with haploinsufficiency of three genes, ALX4 (parietal foramina), EXT2 (multiple exostoses), and PHF21A (craniofacial anomalies, and intellectual disability). The rest of the PSS phenotype is still not associated with a specific gene. We report a systematic review of the literature and included two novel cases. Because deletions are highly variable in size, we defined three groups of patients considering the PSS-genes involved. We found 23 full PSS cases (ALX4, EXT2, and PHF21A), 14 cases with EXT2-ALX4, and three with PHF21A only. Among the latter, we describe a novel male child showing developmental delay, café-au-lait spots, liner postnatal overgrowth and West-like epileptic encephalopathy. We suggest PSS cases may have epileptic spasms early in life, and PHF21A is likely to be the causative gene. Given their subtle presentation these may be overlooked and if left untreated could lead to a severe type or deterioration in the developmental plateau. If our hypothesis is correct, a timely therapy may ameliorate PSS phenotype and improve patients’ outcomes. Our analysis also shows PHF21A is a candidate for the overgrowth phenotype.

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

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

Genetics of Obesity in East Asians

Obesity has become a public health problem worldwide. Compared with Europe, people in Asia tend to suffer from type 2 diabetes with a lower body mass index (BMI). Genome-wide association studies (GWASs) have identified over 750 loci associated with obesity. Although the majority of GWAS results were conducted in individuals of European ancestry, a recent GWAS in individuals of Asian ancestry has made a significant contribution to the identification of obesity susceptibility loci. Indeed, owing to the multifactorial character of obesity with a strong environmental component, the revealed loci may have distinct contributions in different ancestral genetic backgrounds and in different environments as presented through diet and exercise among other factors. Uncovering novel, yet unrevealed genes in non-European ancestries may further contribute to explaining the missing heritability for BMI. In this review, we aimed to summarize recent advances in obesity genetics in individuals of Asian ancestry. We therefore compared proposed mechanisms underlying susceptibility loci for obesity associated with individuals of European and Asian ancestries and discussed whether known genetic variants might explain ethnic differences in obesity risk. We further acknowledged that GWAS implemented in individuals of Asian ancestries have not only validated the potential role of previously specified obesity susceptibility loci but also exposed novel ones, which have been missed in the initial genetic studies in individuals of European ancestries. Thus, multi-ethnic studies have a great potential not only to contribute to a better understanding of the complex etiology of human obesity but also potentially of ethnic differences in the prevalence of obesity, which may ultimately pave new avenues in more targeted and personalized obesity treatments.

Potential research ethics violations against an indigenous tribe in Ecuador: a mixed methods approach

Background

Biomedical and ethnographic studies among indigenous people are common practice in health and geographical research. Prior health research misconduct has been documented, particularly when obtaining genetic material. The objective of this study was to crossmatch previously published data with the perceptions of the Waorani peoples about the trading of their genetic material and other biological samples.

Methods

We conducted a mixed methods study design using a tailored 15-item questionnaire in 72 participants and in-depth interviews in 55 participants belonging to 20 Waorani communities about their experiences and perceptions of participating in biomedical research projects. Additionally, we conducted a systematic review of the literature in order to crossmatch the published results of studies stating the approval of an ethics committee and individual consent within their work.

Results

A total of 40 men (60%) and 32 women (40%), with a mean age of 57 ± 15 years agreed to be interviewed for inclusion. Five main categories around the violation of good clinical practices were identified, concerning the obtention of blood samples from a recently contacted Waorani native community within the Amazonian region of Ecuador. These themes are related to the lack of adequate communication between community members and researchers as well as the voluntariness to participate in health research. Additionally, over 40 years, a total of 38 manuscripts related to the use of biological samples in Waorani indigenous people were published. The majority of the studies (68%) did not state within their article obtaining research ethics board approval, and 71% did not report obtaining the informed consent of the participants prior to the execution of the project.

Conclusion

Clinical Research on the Waorani community in the Ecuadorian Amazon basin has been performed on several occasions. Unfortunately, the majority of these projects did not follow the appropriate ethical and professional standards in either reporting the results or fulfilling them. The results of our investigation suggest that biological material, including genetic material, has been used by researchers globally, with some omitting the minimum information required to guarantee transparency and good clinical practices. We highlight the importance of stating ethics within research to avoid breaches in research transparency.

Novel InDels of GHR, GHRH, GHRHR and Their Association with Growth Traits in Seven Chinese Sheep Breeds

The GH growth axis plays an important role in the growth and development of animals and runs through the whole life of animals. Many studies have shown that molecular mutations in key genes of the GH axis will affect the growth and development of animals. The purpose of this study was to explore the distribution characteristics of InDels of GHR, GHRH, and GHRHR in seven Chinese sheep populations, and to further explore the relationship between InDels and sheep growth traits. GHR showed high variation in Chinese sheep, and GHR-53 showed the highest minimum allele frequency (MAF). There was only one InDel mutation site in both GHRH and GHRHR. The genotype frequencies of Hu sheep (HS), Tong sheep (TS), and Lanzhou fat-tail sheep (LFTS) were quite different from other breeds. The association between GHR, GHRH, and GHRHR InDels and body size traits in seven varieties were analyzed. The results showed that there was no significant relationship between GHRH and body size traits in the seven sheep populations. There was a positive association between GHR-21 and hip height of LFSH (p < 0.05). GHR-43 reduced body height and chest depth of Small tail han sheep (STHS) and hip width of TS. GHR-44 significantly affected the body weight of HS, the body height of STHS and the head depth of TS. GHR-53 significantly reduced cannon girth of HS, chest of STHS and forehead width of TS. GHRHR-2 significantly reduced the body weight of LFHS. To sum up, this study revealed the effects of GHR, GHRH, and GHRHR InDels on sheep phenotypic traits, which indicated their potential application prospects in the genetic improvement of mutton sheep.

An integrated peach genome structural variation map uncovers genes associated with fruit traits

BACKGROUND

Genome structural variations (SVs) have been associated with key traits in a wide range of agronomically important species; however, SV profiles of peach and their functional impacts remain largely unexplored.

RESULTS

Here, we present an integrated map of 202,273 SVs from 336 peach genomes. A substantial number of SVs have been selected during peach domestication and improvement, which together affect 2268 genes. Genome-wide association studies of 26 agronomic traits using these SVs identify a number of candidate causal variants. A 9-bp insertion in Prupe.4G186800, which encodes a NAC transcription factor, is shown to be associated with early fruit maturity, and a 487-bp deletion in the promoter of PpMYB10.1 is associated with flesh color around the stone. In addition, a 1.67 Mb inversion is highly associated with fruit shape, and a gene adjacent to the inversion breakpoint, PpOFP1, regulates flat shape formation.

CONCLUSIONS

The integrated peach SV map and the identified candidate genes and variants represent valuable resources for future genomic research and breeding in peach.

Genetic Risk Profiling in Parkinson's Disease and Utilizing Genetics to Gain Insight into Disease-Related Biological Pathways

Parkinson's disease (PD) is a complex disorder underpinned by both environmental and genetic factors. The latter only began to be understood around two decades ago, but since then great inroads have rapidly been made into deconvoluting the genetic component of PD. In particular, recent large-scale projects such as genome-wide association (GWA) studies have provided insight into the genetic risk factors associated with genetically ''complex'' PD (PD that cannot readily be attributed to single deleterious mutations). Here, we discuss the plethora of genetic information provided by PD GWA studies and how this may be utilized to generate polygenic risk scores (PRS), which may be used in the prediction of risk and trajectory of PD. We also comment on how pathway-specific genetic profiling can be used to gain insight into PD-related biological pathways, and how this may be further utilized to nominate causal PD genes and potentially druggable therapeutic targets. Finally, we outline the current limits of our understanding of PD genetics and the potential contribution of variation currently uncaptured in genetic studies, focusing here on uncatalogued structural variants.

An integrated personal and population-based Egyptian genome reference

A small number of de novo assembled human genomes have been reported to date, and few have been complemented with population-based genetic variation, which is particularly important for North Africa, a region underrepresented in current genome-wide references. Here, we combine long- and short-read whole-genome sequencing data with recent assembly approaches into a de novo assembly of an Egyptian genome. The assembly demonstrates well-balanced quality metrics and is complemented with variant phasing via linked reads into haploblocks, which we associate with gene expression changes in blood. To construct an Egyptian genome reference, we identify genome-wide genetic variation within a cohort of 110 Egyptian individuals. We show that differences in allele frequencies and linkage disequilibrium between Egyptians and Europeans may compromise the transferability of European ancestry-based genetic disease risk and polygenic scores, substantiating the need for multi-ethnic genome references. Thus, the Egyptian genome reference will be a valuable resource for precision medicine.

PASA: Identifying More Credible Structural Variants of Hedou12

Although plenty of structural variant detecting approaches for human genomes can be looked up in the literatures, little has been acknowledged on the effectiveness of those structural variant softwares for plant genomes. Moreover, it has been demonstrated frequent occurrences for those structural variant detecting softwares to find too many false structural variants. In this paper, we devote to detect deletions, insertions, and inversions, in total of three kinds of structural variants occurring in Hedou12 genome in contrast to Williams82 genome. To find more potential structural variants, we try to develop new principles to detect discordant and split read map sets supporting structural variants. Aiming to enhance the precision of structural variant detections, we propose two new sequencing characteristic based probability models, which use the sequencing parameters of Hedou12 genome as well as the parameters for Hedou12 paired-end reads to be aligned onto Williams82, to evaluate the probability for a potential structural variant to occur in. To remove the false members from those potential structural variants, we propose a set cover problem model to describe formally on which potential structural variants it should accept to achieve as high as possible a probability summation. This will achieve a solution with more credible structural variants, which can be verified by comparing with DELLY version 0.5.8 and LUMPY version 0.2.2.3. Our algorithm has been verified to be able to find deletions, insertions, and inversions in Hedou12 in contrast to Williams82 DELLY as well as LUMPY fails to find.

Precision oncology: A primer for pediatric surgeons from the APSA cancer committee

Although most children with cancer can be cured of their disease, a subset of patients with adverse tumor types or biological features, and those with relapsed or refractory disease have significantly worse prognosis. Furthermore, current cytotoxic therapy is associated with significant late effects. Precision oncology, using molecular therapeutics targeted against unique genetic features of the patient's tumor, offers the potential to transform the multimodal therapy for these patients. Potentiated by advances in sequencing technology and molecular therapeutic development, and accelerated by large-scale multi-institutional basket trials, the field of pediatric precision oncology has entered the mainstream. These novel therapeutics have important implications for surgical decision making, as well as pre- and postoperative care. This review summarizes the current state of precision medicine in pediatric oncology including the active North American and European precision oncology clinical trials. LEVEL OF EVIDENCE: Treatment study Level V.

Recurrent inversion toggling and great ape genome evolution

Inversions play an important role in disease and evolution but are difficult to characterize because their breakpoints map to large repeats. We increased by sixfold the number (n = 1,069) of previously reported great ape inversions by using single-cell DNA template strand and long-read sequencing. We find that the X chromosome is most enriched (2.5-fold) for inversions, on the basis of its size and duplication content. There is an excess of differentially expressed primate genes near the breakpoints of large (>100 kilobases (kb)) inversions but not smaller events. We show that when great ape lineage-specific duplications emerge, they preferentially (approximately 75%) occur in an inverted orientation compared to that at their ancestral locus. We construct megabase-pair scale haplotypes for individual chromosomes and identify 23 genomic regions that have recurrently toggled between a direct and an inverted state over 15 million years. The direct orientation is most frequently the derived state for human polymorphisms that predispose to recurrent copy number variants associated with neurodevelopmental disease.

Discovery of tandem and interspersed segmental duplications using high-throughput sequencing

MOTIVATION

Several algorithms have been developed that use high-throughput sequencing technology to characterize structural variations (SVs). Most of the existing approaches focus on detecting relatively simple types of SVs such as insertions, deletions and short inversions. In fact, complex SVs are of crucial importance and several have been associated with genomic disorders. To better understand the contribution of complex SVs to human disease, we need new algorithms to accurately discover and genotype such variants. Additionally, due to similar sequencing signatures, inverted duplications or gene conversion events that include inverted segmental duplications are often characterized as simple inversions, likewise, duplications and gene conversions in direct orientation may be called as simple deletions. Therefore, there is still a need for accurate algorithms to fully characterize complex SVs and thus improve calling accuracy of more simple variants.

RESULTS

We developed novel algorithms to accurately characterize tandem, direct and inverted interspersed segmental duplications using short read whole genome sequencing datasets. We integrated these methods to our TARDIS tool, which is now capable of detecting various types of SVs using multiple sequence signatures such as read pair, read depth and split read. We evaluated the prediction performance of our algorithms through several experiments using both simulated and real datasets. In the simulation experiments, using a 30× coverage TARDIS achieved 96% sensitivity with only 4% false discovery rate. For experiments that involve real data, we used two haploid genomes (CHM1 and CHM13) and one human genome (NA12878) from the Illumina Platinum Genomes set. Comparison of our results with orthogonal PacBio call sets from the same genomes revealed higher accuracy for TARDIS than state-of-the-art methods. Furthermore, we showed a surprisingly low false discovery rate of our approach for discovery of tandem, direct and inverted interspersed segmental duplications prediction on CHM1 (<5% for the top 50 predictions).

AVAILABILITY AND IMPLEMENTATION

TARDIS source code is available at https://github.com/BilkentCompGen/tardis, and a corresponding Docker image is available at https://hub.docker.com/r/alkanlab/tardis/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-Wide Detection of CNVs and Association With Body Weight in Sheep Based on 600K SNP Arrays

Copy number variations (CNVs) are important genomic structural variations and can give rise to significant phenotypic diversity. Herein, we used high-density 600K SNP arrays to detect CNVs in two synthetic lines of sheep (DS and SHH) and in Hu sheep (a local Chinese breed). A total of 919 CNV regions (CNVRs) were detected with a total length of 48.17 Mb, accounting for 1.96% of the sheep genome. These CNVRs consisted of 730 gains, 102 losses, and 87 complex CNVRs. These CNVRs were significantly enriched in the segmental duplication (SD) region. A CNVR-based cluster analysis of the three breeds revealed that the DS and SHH breeds share a close genetic relationship. Functional analysis revealed that some genes in these CNVRs were also significantly enriched in the olfactory transduction pathway (oas04740), including members of the OR gene family such as OR6C76, OR4Q2, and OR4K14. Using association analyses and previous gene annotations, we determined that a subset of identified genes was likely to be associated with body weight, including FOXF2, MAPK12, MAP3K11, STRBP, and C14orf132. Together, these results offer valuable information that will guide future efforts to explore the genetic basis for body weight in sheep.

Identification of genomic insertion and flanking sequences of the transgenic drought-tolerant maize line "SbSNAC1-382" using the single-molecule real-time (SMRT) sequencing method

Safety assessment of genetically modified (GM) crops is crucial at the product-development phase before GM crops are placed on the market. Determining characteristics of sequences flanking exogenous insertion sequences is essential for the safety assessment and marketing of transgenic crops. In this study, we used genome walking and whole-genome sequencing (WGS) to identify the flanking sequence characteristics of the SbSNAC1 transgenic drought-tolerant maize line "SbSNAC1-382", but both of the two methods failed. Then, we constructed a genomic fosmid library of the transgenic maize line, which contained 4.18×105 clones with an average insertion fragment of 35 kb, covering 5.85 times the maize genome. Subsequently, three positive clones were screened by pairs of specific primers, and one of the three positive clones was sequenced by using single-molecule real-time (SMRT) sequencing technology. More than 1.95 Gb sequence data (~105× coverage) for the sequenced clone were generated. The junction reads mapped to the boundaries of T-DNA, and the flanking sequences in the transgenic line were identified by comparing all sequencing reads with the maize reference genome and the sequence of the transgenic vector. Furthermore, the putative insertion loci and flanking sequences were confirmed by PCR amplification and Sanger sequencing. The results indicated that two copies of the exogenous T-DNA fragments were inserted at the same genomic site, and the exogenous T-DNA fragments were integrated at the position of Chromosome 5 from 177155650 to 177155696 in the transgenic line 382. In this study, we demonstrated the successful application of the SMRT technology for the characterization of genomic insertion and flanking sequences.

VALOR2: characterization of large-scale structural variants using linked-reads

Most existing methods for structural variant detection focus on discovery and genotyping of deletions, insertions, and mobile elements. Detection of balanced structural variants with no gain or loss of genomic segments, for example, inversions and translocations, is a particularly challenging task. Furthermore, there are very few algorithms to predict the insertion locus of large interspersed segmental duplications and characterize translocations. Here, we propose novel algorithms to characterize large interspersed segmental duplications, inversions, deletions, and translocations using linked-read sequencing data. We redesign our earlier algorithm, VALOR, and implement our new algorithms in a new software package, called VALOR2.

Structural variation in the sequencing era

Identifying structural variation (SV) is essential for genome interpretation but has been historically difficult due to limitations inherent to available genome technologies. Detection methods that use ensemble algorithms and emerging sequencing technologies have enabled the discovery of thousands of SVs, uncovering information about their ubiquity, relationship to disease and possible effects on biological mechanisms. Given the variability in SV type and size, along with unique detection biases of emerging genomic platforms, multiplatform discovery is necessary to resolve the full spectrum of variation. Here, we review modern approaches for investigating SVs and proffer that, moving forwards, studies integrating biological information with detection will be necessary to comprehensively understand the impact of SV in the human genome.

Comparison of Arachis monticola with Diploid and Cultivated Tetraploid Genomes Reveals Asymmetric Subgenome Evolution and Improvement of Peanut

Like many important crops, peanut is a polyploid that underwent polyploidization, evolution, and domestication. The wild allotetraploid peanut species Arachis monticola (A. monticola) is an important and unique link from the wild diploid species to cultivated tetraploid species in the Arachis lineage. However, little is known about A. monticola and its role in the evolution and domestication of this important crop. A fully annotated sequence of ≈2.6 Gb A. monticola genome and comparative genomics of the Arachis species is reported. Genomic reconstruction of 17 wild diploids from AA, BB, EE, KK, and CC groups and 30 tetraploids demonstrates a monophyletic origin of A and B subgenomes in allotetraploid peanuts. The wild and cultivated tetraploids undergo asymmetric subgenome evolution, including homoeologous exchanges, homoeolog expression bias, and structural variation (SV), leading to subgenome functional divergence during peanut domestication. Significantly, SV-associated homoeologs tend to show expression bias and correlation with pod size increase from diploids to wild and cultivated tetraploids. Moreover, genomic analysis of disease resistance genes shows the unique alleles present in the wild peanut can be introduced into breeding programs to improve some resistance traits in the cultivated peanuts. These genomic resources are valuable for studying polyploid genome evolution, domestication, and improvement of peanut production and resistance.

Rare copy number variants in over 100,000 European ancestry subjects reveal multiple disease associations

Copy number variants (CNVs) are suggested to have a widespread impact on the human genome and phenotypes. To understand the role of CNVs across human diseases, we examine the CNV genomic landscape of 100,028 unrelated individuals of European ancestry, using SNP and CGH array datasets. We observe an average CNV burden of ~650 kb, identifying a total of 11,314 deletion, 5625 duplication, and 2746 homozygous deletion CNV regions (CNVRs). In all, 13.7% are unreported, 58.6% overlap with at least one gene, and 32.8% interrupt coding exons. These CNVRs are significantly more likely to overlap OMIM genes (2.94-fold), GWAS loci (1.52-fold), and non-coding RNAs (1.44-fold), compared with random distribution (P < 1 × 10⁻³). We uncover CNV associations with four major disease categories, including autoimmune, cardio-metabolic, oncologic, and neurological/psychiatric diseases, and identify several drug-repurposing opportunities. Our results demonstrate robust frequency definition for large-scale rare variant association studies, identify CNVs associated with major disease categories, and illustrate the pleiotropic impact of CNVs in human disease.

Associations of copy number variations (CNVs) with complex traits are challenging to study because of their low frequency. Here, the authors analyse SNP array and array comparative genomic hybridization data of 100,028 individuals and report their associations with immune-related, cardiometabolic and neuropsychiatric diseases as well as cancer.

High-throughput long paired-end sequencing of a Fosmid library by PacBio

BACKGROUND

Large insert paired-end sequencing technologies are important tools for assembling genomes, delineating associated breakpoints and detecting structural rearrangements. To facilitate the comprehensive detection of inter- and intra-chromosomal structural rearrangements or variants (SVs) and complex genome assembly with long repeats and segmental duplications, we developed a new method based on single-molecule real-time synthesis sequencing technology for generating long paired-end sequences of large insert DNA libraries.

RESULTS

A Fosmid vector, pHZAUFOS3, was developed with the following new features: (1) two 18-bp non-palindromic I-SceI sites flank the cloning site, and another two sites are present in the skeleton of the vector, allowing long DNA inserts (and the long paired-ends in this paper) to be recovered as single fragments and the vector (~ 8 kb) to be fragmented into 2-3 kb fragments by I-SceI digestion and therefore was effectively removed from the long paired-ends (5-10 kb); (2) the chloramphenicol (Cm) resistance gene and replicon (oriV), necessary for colony growth, are located near the two sides of the cloning site, helping to increase the proportion of the paired-end fragments to single-end fragments in the paired-end libraries. Paired-end libraries were constructed by ligating the size-selected, mechanically sheared pooled Fosmid DNA fragments to the Ampicillin (Amp) resistance gene fragment and screening the colonies with Cm and Amp. We tested this method on yeast and Setaria italica Yugu1. Fosmid-size paired-ends with an average length longer than 2 kb for each end were generated. The N50 scaffold lengths of the de novo assemblies of the yeast and S. italica Yugu1 genomes were significantly improved. Five large and five small structural rearrangements or assembly errors spanning tens of bp to tens of kb were identified in S. italica Yugu1 including deletions, inversions, duplications and translocations.

CONCLUSIONS

We developed a new method for long paired-end sequencing of large insert libraries, which can efficiently improve the quality of de novo genome assembly and identify large and small structural rearrangements or assembly errors.