Genomic rearrangements and sporadic disease

Genome-wide association testing beyond SNPs

Decades of genetic association testing in human cohorts have provided important insights into the genetic architecture and biological underpinnings of complex traits and diseases. However, for certain traits, genome-wide association studies (GWAS) for common SNPs are approaching signal saturation, which underscores the need to explore other types of genetic variation to understand the genetic basis of traits and diseases. Copy number variation (CNV) is an important source of heritability that is well known to functionally affect human traits. Recent technological and computational advances enable the large-scale, genome-wide evaluation of CNVs, with implications for downstream applications such as polygenic risk scoring and drug target identification. Here, we review the current state of CNV-GWAS, discuss current limitations in resource infrastructure that need to be overcome to enable the wider uptake of CNV-GWAS results, highlight emerging opportunities and suggest guidelines and standards for future GWAS for genetic variation beyond SNPs at scale.

FANCD2 genome binding is nonrandom and is enriched at large transcriptionally active neural genes prone to copy number variation

Fanconi anemia (FA) is a rare genetic disease characterized by congenital abnormalities and increased risk for bone marrow failure and cancer. Central nervous system defects, including acute and irreversible loss of neurological function and white matter lesions with calcifications, have become increasingly recognized among FA patients, and are collectively referred to as Fanconi Anemia Neurological Syndrome or FANS. The molecular etiology of FANS is poorly understood. In this study, we have used a functional integrative genomics approach to further define the function of the FANCD2 protein and FA pathway. Combined analysis of new and existing FANCD2 ChIP-seq datasets demonstrates that FANCD2 binds nonrandomly throughout the genome with binding enriched at transcription start sites and in broad regions spanning protein-coding gene bodies. FANCD2 demonstrates a strong preference for large neural genes involved in neuronal differentiation, synapse function, and cell adhesion, with many of these genes implicated in neurodevelopmental and neuropsychiatric disorders. Furthermore, FANCD2 binds to regions of the genome that replicate late, undergo mitotic DNA synthesis (MiDAS) under conditions of replication stress, and are hotspots for copy number variation. Our analysis describes an important targeted role for FANCD2 and the FA pathway in the maintenance of large neural gene stability.

Copy Number Variants in 30 Saudi Pediatric Patients with Neurodevelopmental Disorders: From Unknown Significance to Diagnosis

Background

Structural variants (SVs), such as copy number variants (CNVs), insertions, deletions, inversions, and translocations, contribute significantly to genetic diversity and disease etiology. CNVs, which involve the duplication or deletion of DNA segments, are particularly impactful on genes crucial for biological functions and disease processes.

Objective

To reassess unclassified SVs that may be underlying unresolved neurodevelopmental disorders among Saudi patients.

Methodology

In this retrospective study conducted at King Saud Medical City, Riyadh, Saudi Arabia, 30 probands with neurodevelopmental disorders and congenital malformations were examined using next-generation sequencing methods-exome sequencing, gene panels, or SNP arrays (the Illumina platform). Reclassification was aided by online tools such as VarSome and ClinVar, with pathogenicity assessments using the ClinGen CNV Pathogenicity Calculator based on American College of Medical Genetics and Genomics criteria for CNV loss and gain, and dosage sensitivity.

Results

A total of 31 CNVs were analyzed, of which 2 were reclassified: one as benign and the other as pathogenic. The pathogenic CNV, [3p13p12.3 (70411134_75249376) x1], included a deletion of the FOXP1 gene and was associated with an intellectual developmental disorder, language impairment, possible autistic features, psychomotor impairment, developmental regression, and epilepsy.

Conclusion

This study underscores the importance of continuously documenting and revisiting unclassified CNVs in accessible databases to enhance the diagnosis and understanding of complex genotype-phenotype relationships. Reclassifying these CNVs not only accelerates diagnostic processes but also enriches our insight into their significant roles in health and disease.

Genome integrity as a potential index of longevity in Ashkenazi Centenarian's families

The aging process, or senescence, is characterized by age-specific decline in physical and physiological function, and increased frailty and genomic changes, including mutation accumulation. However, the mechanisms through which changes in genomic architecture influence human longevity have remained obscure. Copy number variants (CNVs), an abundant class of genomic variants, offer unique opportunities for understanding age-related genomic changes. Here we report the spectrum of CNVs in a cohort of 670 Ashkenazi Jewish centenarians, their progeny, and unrelated controls. The average ages of these groups were 97.4 ± 2.8, 69.2 ± 9.2, and 66.5 ± 7.0 respectively. For the first time, we compared different size classes of CNVs, from 1 kB to 100 MB in size. Using a high-resolution custom Affymetrix array, targeting 44,639 genomic regions, we identified a total of 12,166, 22,188, and 10,285 CNVs in centenarians, their progeny, and control groups, respectively. Interestingly, the offspring group showed the highest number of unique CNVs, followed by control and centenarians. While both gains and losses were found in all three groups, centenarians showed a significantly higher average number of both total gains and losses relative to their controls (p < 0.0327, 0.0182, respectively). Moreover, centenarians showed a lower total length of genomic material lost, suggesting that they may maintain superior genomic integrity over time. We also observe a significance fold increase of CNVs among the offspring, implying greater genomic integrity and a putative mechanism for longevity preservation. Genomic regions that experienced loss or gains appear to be distributed across many sites in the genome and contain genes involved in DNA transcription, cellular transport, developmental pathways, and metabolic functions. Our findings suggest that the exceptional longevity observed in centenarians may be attributed to the prolonged maintenance of functionally important genes. These genes are intrinsic to specific genomic regions as well as to the overall integrity of the genomic architecture. Additionally, a strong association between longer CNVs and differential gene expression observed in this study supports the notion that genomic integrity could positively influence longevity.

Modelling DNA damage-repair and beyond

The paper presents a review of mechanistic modelling studies of DNA damage and DNA repair, and consequences to follow in mammalian cell nucleus. We hypothesize DNA deletions are consequences of repair of double strand breaks leading to the modifications of genome that play crucial role in long term development of genetic inheritance and diseases. The aim of the paper is to review formation mechanisms underlying naturally occurring DNA deletions in the human genome and their potential relevance for bridging the gap between induced DNA double strand breaks and deletions in damaged human genome from endogenous and exogenous events. The model of the cell nucleus presented enables simulation of DNA damage at molecular level identifying the spectrum of damage induced in all chromosomal territories and loops. Our mechanistic modelling of DNA repair for double stand breaks (DSB), single strand breaks (SSB) and base damage (BD), shows the complexity of DNA damage is responsible for the longer repair times and the reason for the biphasic feature of mammalian cells repair curves. In the absence of experimentally determined data, the mechanistic model of repair predicts the in vivo rate constants for the proteins involved in the repair of DSB, SSB, and of BD.

Comprehensive Analysis of Clinically Relevant Copy Number Alterations (CNAs) Using a 523-Gene Next-Generation Sequencing Panel and NxClinical Software in Solid Tumors

Copy number alterations (CNAs) are significant in tumor initiation and progression. Identifying these aberrations is crucial for targeted therapies and personalized cancer diagnostics. Next-generation sequencing (NGS) methods present advantages in scalability and cost-effectiveness, surpassing limitations associated with reference assemblies and probe capacities in traditional laboratory approaches. This retrospective study evaluated CNAs in 50 FFPE tumor samples (breast cancer, ovarian carcinoma, pancreatic cancer, melanoma, and prostate carcinoma) using Illumina's TruSight Oncology 500 (TSO500) and the Affymetrix Oncoscan Molecular Inversion Probe (OS-MIP) (ThermoFisher Scientific, Waltham, MA, USA). NGS analysis with the NxClinical 6.2 software demonstrated a high sensitivity and specificity (100%) for CNA detection, with a complete concordance rate as compared to the OS-MIP. All 54 known CNAs were identified by NGS, with gains being the most prevalent (63%). Notable CNAs were observed in MYC (18%), TP53 (12%), BRAF (8%), PIK3CA, EGFR, and FGFR1 (6%) genes. The diagnostic parameters exhibited high accuracy, including a positive predictive value, negative predictive value, and overall diagnostic accuracy. This study underscores NxClinical as a reliable software for identifying clinically relevant gene alterations using NGS TSO500, offering valuable insights for personalized cancer treatment strategies based on CNA analysis.

Large-scale genomic rearrangements boost SCRaMbLE in Saccharomyces cerevisiae

Synthetic Chromosome Rearrangement and Modification by LoxP-mediated Evolution (SCRaMbLE) is a promising tool to study genomic rearrangements. However, the potential of SCRaMbLE to study genomic rearrangements is currently hindered, because a strain containing all 16 synthetic chromosomes is not yet available. Here, we construct SparLox83R, a yeast strain containing 83 loxPsym sites distributed across all 16 chromosomes. SCRaMbLE of SparLox83R produces versatile genome-wide genomic rearrangements, including inter-chromosomal events. Moreover, when combined with synthetic chromosomes, SCRaMbLE of hetero-diploids with SparLox83R leads to increased diversity of genomic rearrangements and relatively faster evolution of traits compared to hetero-diploids only with wild-type chromosomes. Analysis of the SCRaMbLEd strain with increased tolerance to nocodazole demonstrates that genomic rearrangements can perturb the transcriptome and 3D genome structure and consequently impact phenotypes. In summary, a genome with sparsely distributed loxPsym sites can serve as a powerful tool for studying the consequence of genomic rearrangements and accelerating strain engineering in Saccharomyces cerevisiae.

Identification of a spontaneously arising variant affecting thermotaxis behavior in a recombinant inbred Caenorhabditis elegans line

Analyses of the contributions of genetic variants in wild strains to phenotypic differences have led to a more complete description of the pathways underlying cellular functions. Causal loci are typically identified via interbreeding of strains with distinct phenotypes in order to establish recombinant inbred lines (RILs). Since the generation of RILs requires growth for multiple generations, their genomes may contain not only different combinations of parental alleles but also genetic changes that arose de novo during the establishment of these lines. Here, we report that in the course of generating RILs between Caenorhabditis elegans strains that exhibit distinct thermotaxis behavioral phenotypes, we identified spontaneously arising variants in the ttx-1 locus. ttx-1 encodes the terminal selector factor for the AFD thermosensory neurons, and loss-of-function mutations in ttx-1 abolish thermotaxis behaviors. The identified genetic changes in ttx-1 in the RIL are predicted to decrease ttx-1 function in part via specifically affecting a subset of AFD-expressed ttx-1 isoforms. Introduction of the relevant missense mutation in the laboratory C. elegans strain via gene editing recapitulates the thermotaxis behavioral defects of the RIL. Our results suggest that spontaneously occurring genomic changes in RILs may complicate identification of loci contributing to phenotypic variation, but that these mutations may nevertheless lead to the identification of important causal molecules and mechanisms.

Effects of Parkin on the Mitochondrial Genome in the Heart and Brain of Mitochondrial Mutator Mice

Mitochondrial dysfunction has been implicated in neurodegenerative diseases like Parkinson's disease (PD). This study investigates the role of Parkin, a protein involved in mitochondrial quality control, and strongly linked to PD, in the context of mitochondrial DNA (mtDNA) mutations. Mitochondrial mutator mice (PolgD257A/D257A ) (Polg) are used and bred with Parkin knockout (PKO) mice or mice with disinhibited Parkin (W402A). In the brain, mtDNA mutations are analyzed in synaptosomes, presynaptic neuronal terminals, which are far from neuronal soma, which likely renders mitochondria there more vulnerable compared with brain homogenate. Surprisingly, PKO results in reduced mtDNA mutations in the brain but increased control region multimer (CRM) in synaptosomes. In the heart, both PKO and W402A lead to increased mutations, with W402A showing more mutations in the heart than PKO. Computational analysis reveals many of these mutations are deleterious. These findings suggest that Parkin plays a tissue-dependent role in regulating mtDNA damage response, with differential effects in the brain and heart. Understanding the specific role of Parkin in different tissues may provide insights into the underlying mechanisms of PD and potential therapeutic strategies. Further investigation into these pathways can enhance the understanding of neurodegenerative diseases associated with mitochondrial dysfunction.

SNV/indel hypermutator phenotype in biallelic RAD51C variant: Fanconi anemia

We previously reported a fetus with Fanconi anemia (FA), complementation group O due to compound heterozygous variants involving RAD51C. Interestingly, the trio exome sequencing analysis also detected eight apparent de novo mosaic variants with variant allele fraction (VAF) ranging between 11.5 and 37%. Here, using whole genome sequencing and a 'home-brew' variant filtering pipeline and DeepMosaic module, we investigated the number and signature of de novo heterozygous and mosaic variants and the hypothesis of a rare phenomenon of hypermutation. Eight-hundred-thirty apparent de novo SNVs and 21 de novo indels had VAFs below 37.41% and were considered postzygotic somatic mosaic variants. The VAFs showed a bimodal distribution, with one component having an average VAF of 25% (range: 18.7-37.41%) (n = 446), representing potential postzygotic first mitotic events, and the other component with an average VAF of 12.5% (range 9.55-18.69%) (n = 384), describing potential second mitotic events. No increased rate of CNV formation was observed. The mutational pattern analysis for somatic single base substitution showed SBS40, SBS5, and SBS3 as the top recognized signatures. SBS3 is a known signature associated with homologous recombination-based DNA damage repair error. Our data demonstrate that biallelic RAD51C variants show evidence for defective genomic DNA damage repair and thereby result in a hypermutator phenotype with the accumulation of postzygotic de novo mutations, at least in the prenatal period. This 'genome hypermutator phenomenon' might contribute to the observed hematological manifestations and the predisposition to tumors in patients with FA. We propose that other FA groups should be investigated for genome-wide de novo variants.

Genome-Wide Detection and Analysis of Copy Number Variation in Anhui Indigenous and Western Commercial Pig Breeds Using Porcine 80K SNP BeadChip

Copy number variation (CNV) is an important class of genetic variations widely associated with the porcine genome, but little is known about the characteristics of CNVs in foreign and indigenous pig breeds. We performed a genome-wide comparison of CNVs between Anhui indigenous pig (AHIP) and Western commercial pig (WECP) breeds based on data from the Porcine 80K SNP BeadChip. After analysis using the PennCNV software, we detected 3863 and 7546 CNVs in the AHIP and WECP populations, respectively. We obtained 225 (loss: 178, gain: 47) and 379 (loss: 293, gain: 86) copy number variation regions (CNVRs) randomly distributed across the autosomes of the AHIP and WECP populations, accounting for 10.90% and 22.57% of the porcine autosomal genome, respectively. Functional enrichment analysis of genes in the CNVRs identified genes related to immunity (FOXJ1, FOXK2, MBL2, TNFRSF4, SIRT1, NCF1) and meat quality (DGAT1, NT5E) in the WECP population; these genes were a loss event in the WECP population. This study provides important information on CNV differences between foreign and indigenous pig breeds, making it possible to provide a reference for future improvement of these breeds and their production performance.

Exome/Genome Sequencing in Undiagnosed Syndromes

Exome sequencing (ES) and genome sequencing (GS) have radically transformed the diagnostic approach to undiagnosed rare/ultrarare Mendelian diseases. Next-generation sequencing (NGS), the technology integral for ES, GS, and most large (100+) gene panels, has enabled previously unimaginable diagnoses, changes in medical management, new treatments, and accurate reproductive risk assessments for patients, as well as new disease gene discoveries. Yet, challenges remain, as most individuals remain undiagnosed with current NGS. Improved NGS technology has resulted in long-read sequencing, which may resolve diagnoses in some patients who do not obtain a diagnosis with current short-read ES and GS, but its effectiveness is unclear, and it is expensive. Other challenges that persist include the resolution of variants of uncertain significance, the urgent need for patients with ultrarare disorders to have access to therapeutics, the need for equity in patient access to NGS-based testing, and the study of ethical concerns. However, the outlook for undiagnosed disease resolution is bright, due to continual advancements in the field.

Prenatally diagnosed 16p11.2 copy number variations by SNP Array: A retrospective case series

OBJECTIVE

The 16p11.2 copy number variations (CNVs) are increasingly recognized as one of the most frequent genomic disorders, with a broad spectrum of phenotypes. The fetal phenotype associated with 16p11.2 CNVs is poorly described. The current study presents prenatal series of 16p11.2 CNVs and provides a better understanding of this submicroscopic imbalance in prenatal diagnosis.

METHOD

Retrospective case series were extracted from a single tertiary referral center performing prenatal single nucleotide polymorphism (SNP) array from April 2017 to December 2021. The maternal demographics, indication for amniocentesis, ultrasound findings, SNP array results, inheritance of the CNVs, and pregnancy outcomes were studied.

RESULTS

We indentified 30 fetuses carrying 16p11.2 CNVs, representing 0.35% (30/8578) of prenatal SNP array results. The series included 17 fetuses with a proximal deletion, 7 with a distal deletion, 4 with a proximal duplication, and 2 with a distal duplication. Prenatal ultrasound anomalies were reported in 80% of these cases. The most common presentation was vertebralanomalies (9/30). Other features noted in more than one fetus were increased nuchal translucency/nuchal fold (NT/NF) (5/30), absent/hypoplastic nasal bone (3/30), polyhydramnios (3/30), ventricular septal defect (VSD) (2/30), unilateral mild ventriculomegaly (2/30), fetal growth restriction (FGR) (2/30), right aortic arch (2/30). All the 9 vertebralanomalies were present in fetuses harboring proximal deletion (9/17). Familial transmission was confirmed in 44% of cases (11/25) and termination of pregnancy was requested in 62.1% (18/29) of cases.

CONCLUSION

The 16p11.2 CNVs can have variable prenatal phenotypes and these CNVs are frequently inherited from parents with a milder or normal phenotype. Our results underline that vertebral deformities were frequent in cases of 16p11.2 proximal deletion, and further demonstrate the incomplete penetrance of the CNVs.

The holistic approach to the CHRNA7 gene, hsa-miR-3158-5p, and 15q13.3 hotspot CNVs in migraineurs

Migraine is a neurological disease characterized by severe headache attacks. Combinations of different genetic variations such as copy number variation (CNV) in a gene and microRNA (miRNA) expression can provide a holistic approach to the disease as a pathophysiological, diagnostic, and therapeutic target. CNVs, the Cholinergic Receptor Nicotinic Alpha 7 Subunit (CHRNA7) gene, and expression of gene-targeting miRNAs (hsa-miR-548e-5p and hsa-miR-3158-5p) in migraineurs (n = 102; with aura, n = 43; without aura, n = 59) and non-migraines (n = 120) aged 15-60 years, comparative, case-control study was conducted. Genetic markers were compared with biochemical parameters (BMI, WBC, Urea, GFR, ESR, CRP, HBG). All analyzes were performed by quantitative Real-Time PCR (q-PCR) and fold change was calculated with the 2-ΔΔCT method. The diagnostic power of the CHRNA7 gene, CNV, and miRNAs were analyzed with the receiver operating curve (ROC). CHRNA7 gene and hsa-miR-3158-5p are down-regulated in migraineurs and the gene is controlled by this miRNA via CNVs (p < .05). Both deletion and duplication were detected in patients with migraine for CVN numbers (p = .05). The number of CNV deletions was higher than duplications. When CHRNA7-CNV-hsa-miR-3158-5p was modeled together in the ROC analysis, the area under the curve (AUC) was 0.805, and the diagnostic power was "good". In migraineurs, the CHRNA7 gene can be controlled by hsa-miR-3158-5p via CNVs to modulate the mechanism of pain. These three genetic markers have diagnostic potential and may be used in antimigraine treatments.

Sequencing individual genomes with recurrent genomic disorder deletions: an approach to characterize genes for autosomal recessive rare disease traits

BACKGROUND

In medical genetics, discovery and characterization of disease trait contributory genes and alleles depends on genetic reasoning, study design, and patient ascertainment; we suggest a segmental haploid genetics approach to enhance gene discovery and molecular diagnostics.

METHODS

We constructed a genome-wide map for nonallelic homologous recombination (NAHR)-mediated recurrent genomic deletions and used this map to estimate population frequencies of NAHR deletions based on large-scale population cohorts and region-specific studies. We calculated recessive disease carrier burden using high-quality pathogenic or likely pathogenic variants from ClinVar and gnomAD. We developed a NIRD (NAHR deletion Impact to Recessive Disease) score for recessive disorders by quantifying the contribution of NAHR deletion to the overall allele load that enumerated all pairwise combinations of disease-causing alleles; we used a Punnett square approach based on an assumption of random mating. Literature mining was conducted to identify all reported patients with defects in a gene with a high NIRD score; meta-analysis was performed on these patients to estimate the representation of NAHR deletions in recessive traits from contemporary human genomics studies. Retrospective analyses of extant clinical exome sequencing (cES) were performed for novel rare recessive disease trait gene and allele discovery from individuals with NAHR deletions.

RESULTS

We present novel genomic insights regarding the genome-wide impact of NAHR recurrent segmental variants on recessive disease burden; we demonstrate the utility of NAHR recurrent deletions to enhance discovery in the challenging context of autosomal recessive (AR) traits and biallelic variation. Computational results demonstrate new mutations mediated by NAHR, involving recurrent deletions at 30 genomic regions, likely drive recessive disease burden for over 74% of loci within these segmental deletions or at least 2% of loci genome-wide. Meta-analyses on 170 literature-reported patients implicate that NAHR deletions are depleted from the ascertained pool of AR trait alleles. Exome reanalysis of personal genomes from subjects harboring recurrent deletions uncovered new disease-contributing variants in genes including COX10, ERCC6, PRRT2, and OTUD7A.

CONCLUSIONS

Our results demonstrate that genomic sequencing of personal genomes with NAHR deletions could dramatically improve allele and gene discovery and enhance clinical molecular diagnosis. Moreover, results suggest NAHR events could potentially enable human haploid genetic screens as an approach to experimental inquiry into disease biology.

Analysis of genomic alterations in cancer associated human pancreatic stellate cells

Pancreatic stellate cells (PSCs) constitute important cells of the pancreatic microenvironment and their close interaction with cancer cells is important in pancreatic cancer. It is currently not known whether PSCs accumulate genetic alterations that contribute to tumor biology. Our aim was to analyze genetic alterations in cancer associated PSCs. PSC DNA was matched to DNA isolated from pancreatic cancer patients’ blood (n = 5) and analyzed by Next-Generation Sequencing (NGS). Bioinformatic analysis was performed using the GATK software and pathogenicity prediction scores. Sanger sequencing was carried out to verify specific genetic alterations in a larger panel of PSCs (n = 50). NGS and GATK analysis identified on average 26 single nucleotide variants in PSC DNA as compared to the matched blood DNA that could be visualized with the Integrative Genomics Viewer. The absence of PDAC driver mutations (KRAS, p53, p16/INK4a, SMAD4) confirmed that PSC isolations were not contaminated with cancer cells. After filtering the variants, using different pathogenicity scores, ten genes were identified (SERPINB2, CNTNAP4, DENND4B, DPP4, FGFBP2, MIGA2, POLE, SNRNP40, TOP2B, and ZDHHC18) in single samples and confirmed by Sanger sequencing. As a proof of concept, functional analysis using control and SERPINB2 knock-out fibroblasts revealed functional effects on growth, migration, and collagen contraction. In conclusion, PSC DNA exhibit a substantial amount of single nucleotide variants that might have functional effects potentially contributing to tumor aggressiveness.

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.

Characterization of non-specific lipid transfer protein (nsLtp) gene families in the Brassica napus pangenome reveals abundance variation

BACKGROUND

Brassica napus is an important agricultural species, improving stress resistance was one of the main breeding goals at present. Non-specific lipid transfer proteins (nsLTPs) are small, basic proteins which are involved in some biotic or abiotic stress responses. B. napus is susceptible to a variety of fungal diseases, so identify the BnLTPs and their expression in disease responses is very important. The common reference genome of B. napus does not contain all B. napus genes because of gene presence/absence variations between individuals. Therefore, it was necessary to search for candidate BnLTP genes in the B. napus pangenome.

RESULTS

In the present study, the BnLTP genes were identified throughout the pangenome, and different BnLTP genes were presented among varieties. Totally, 246 BnLTP genes were identified and could be divided into five types (1, 2, C, D, and G). The classification, phylogenetic reconstruction, chromosome distribution, functional annotation, and gene expression were analyzed. We also identified potential cis-elements that respond to biotic and abiotic stresses in the 2 kb upstream regions of all BnLTP genes. RNA sequencing analysis showed that the BnLTP genes were involved in the response to Sclerotinia sclerotiorum infection. We identified 32 BnLTPs linked to blackleg resistance quantitative trait locus (QTL).

CONCLUSION

The identification and analysis of LTP genes in the B. napus pangenome could help to elucidate the function of BnLTP family members and provide new information for future molecular breeding in B. napus.

Clan genomics: From OMIM phenotypic traits to genes and biology

Clinical characterization of a patient phenotype has been the quintessential approach for elucidating a differential diagnosis and a hypothesis to explore a potential clinical diagnosis. This has resulted in a language of medicine and a semantic ontology, with both specialty- and subspecialty-specific lexicons, that can be challenging to translate and interpret. There is no 'Rosetta Stone' of clinical medicine such as the genetic code that can assist translation and interpretation of the language of genetics. Nevertheless, the information content embodied within a clinical diagnosis can guide management, therapeutic intervention, and potentially prognostic outlook of disease enabling anticipatory guidance for patients and families. Clinical genomics is now established firmly in medical practice. The granularity and informative content of a personal genome is immense. Yet, we are limited in our utility of much of that personal genome information by the lack of functional characterization of the overwhelming majority of computationally annotated genes in the haploid human reference genome sequence. Whereas DNA and the genetic code have provided a 'Rosetta Stone' to translate genetic variant information, clinical medicine, and clinical genomics provide the context to understand human biology and disease. A path forward will integrate deep phenotyping, such as available in a clinical synopsis in the Online Mendelian Inheritance in Man (OMIM) entries, with personal genome analyses.

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.

Adverse biological effects of ingested polystyrene microplastics using Drosophila melanogaster as a model in vivo organism

The abundant presence and extensive use of polystyrene microplastics (PSMPs) has recently become a serious environmental concern, as impact of exposure to these substances on human health remains unknown. While in vitro studies yield data on adverse effect of PSMPs, in vivo approaches are more relevant for risk assessment. Drosophila melanogaster is one of the most genetically and experimentally accessible model organisms used in biology as an in vivo model. D. melanogaster was selected as a representative in vivo model organism to examine the genotoxic potential of PSMPs at 5 concentrations of three different sizes namely 4, 10, or 20 µm. In particular, the wing somatic mutation and recombination test (SMART), a scalable, time-efficient in vivo assay developed to study genotoxicity of various compounds in a rapid manner at low costs was used. The third-instar Drosophila larvae were exposed to PSMPs through food at 5 concentrations ranging from 0.01-10 mM. Viability (lethality), larval length, morphological deformations, locomotor activity (climbing behavior), and genotoxic effects were the end-points measured. Exposure to PSMPs at 4, 10, or 20 µm produced significant morphological defects, impaired climbing behavior, and genotoxicity as evidenced by the SMART test demonstrating induction of somatic recombination. Significant increases were observed in the frequency of total spots, suggesting that PSMPs might induce genotoxic activity predominantly via initiation of somatic DNA recombination in a concentration-dependent manner.

Identification of Copy Number Variants in a Southern Chinese Cohort of Patients with Congenital Scoliosis

Congenital scoliosis (CS) is a lateral curvature of the spine resulting from congenital vertebral malformations (CVMs) and affects 0.5–1/1000 live births. The copy number variant (CNV) at chromosome 16p11.2 has been implicated in CVMs and recent studies identified a compound heterozygosity of 16p11.2 microdeletion and TBX6 variant/haplotype causing CS in multiple cohorts, which explains about 5–10% of the affected cases. Here, we studied the genetic etiology of CS by analyzing CNVs in a cohort of 67 patients with congenital hemivertebrae and 125 family controls. We employed both candidate gene and family-based approaches to filter CNVs called from whole exome sequencing data. This identified 12 CNVs in four scoliosis-associated genes (TBX6, NOTCH2, DSCAM, and SNTG1) as well as eight recessive and 64 novel rare CNVs in 15 additional genes. Some candidates, such as DHX40, NBPF20, RASA2, and MYSM1, have been found to be associated with syndromes with scoliosis or implicated in bone/spine development. In particular, the MYSM1 mutant mouse showed spinal deformities. Our findings suggest that, in addition to the 16p11.2 microdeletion, other CNVs are potentially important in predisposing to CS.

Genome diversity and instability in human germ cells and preimplantation embryos

Genome diversity is essential for evolution and is of fundamental importance to human health. Generating genome diversity requires phases of DNA damage and repair that can cause genome instability. Humans have a high incidence of de novo congenital disorders compared to other organisms. Recent access to eggs, sperm and preimplantation embryos is revealing unprecedented rates of genome instability that may result in infertility and de novo mutations that cause genomic imbalance in at least 70% of conceptions. The error type and incidence of de novo mutations differ during developmental stages and are influenced by differences in male and female meiosis. In females, DNA repair is a critical factor that determines fertility and reproductive lifespan. In males, aberrant meiotic recombination causes infertility, embryonic failure and pregnancy loss. Evidence suggest germ cells are remarkably diverse in the type of genome instability that they display and the DNA damage responses they deploy. Additionally, the initial embryonic cell cycles are characterized by a high degree of genome instability that cause congenital disorders and may limit the use of CRISPR-Cas9 for heritable genome editing.

A novel homozygous SLC13A5 whole-gene deletion generated by Alu/Alu-mediated rearrangement in an Iraqi family with epileptic encephalopathy

Biallelic loss-of-function (LoF) of SLC13A5 (solute carrier family 13, member 5) induced deficiency in sodium/citrate transporter (NaCT) causes autosomal recessive developmental epileptic encephalopathy 25 with hypoplastic amelogenesis imperfecta (DEE25; MIM #615905). Many pathogenic SLC13A5 single nucleotide variants (SNVs) and small indels have been described; however, no cases with copy number variants (CNVs) have been sufficiently investigated. We describe a consanguineous Iraqi family harboring an 88.5 kb homozygous deletion including SLC13A5 in Chr17p13.1. The three affected male siblings exhibit neonatal-onset epilepsy with fever-sensitivity, recurrent status epilepticus, global developmental delay/intellectual disability (GDD/ID), and other variable neurological findings as shared phenotypical features of DEE25. Two of the three affected subjects exhibit hypoplastic amelogenesis imperfecta (AI), while the proband shows no evidence of dental abnormalities or AI at 2 years of age with apparently unaffected primary dentition. Characterization of the genomic architecture at this locus revealed evidence for genomic instability generated by an Alu/Alu-mediated rearrangement; confirmed by break-point junction Sanger sequencing. This multiplex family from a distinct population elucidates the phenotypic consequence of complete LoF of SLC13A5 and illustrates the importance of read-depth-based CNV detection in comprehensive exome sequencing analysis to solve cases that otherwise remain molecularly unsolved.

The Copy Number Variation of OsMTD1 Regulates Rice Plant Architecture

Copy number variation (CNV) may have phenotypic effects by altering the expression level of the gene(s) or regulatory element(s) contained. It is believed that CNVs play pivotal roles in controlling plant architecture and other traits in plant. However, the effects of CNV contributing to special traits remain largely unknown. Here we report a CNV involved in rice architecture by modulating tiller number and leaf angle. In the genome of Oryza sativa ssp. japonica cv. Nipponbare, we found a locus Loc_Os08g34249 is derived from a 13,002-bp tandem duplication in the nearby region of OsMTD1, a gene regulating tillering in rice. Further survey of 230 rice cultivars showed that the duplication occurred in only 13 japonica rice cultivars. Phenotypic investigation indicated that this CNV region may contribute to tiller number. Moreover, we revealed that OsMTD1 not only influences rice tiller number and leaf angle, but also represses pri-miR156f transcription in the CNV region. Intriguingly, this CNV performs function through both the dosage and position effects on OsMTD1 and pri-miR156f. Thus, our work identified a CNV and revealed a molecular regulatory basis for its effects on plant architecture, implying this CNV may possess importance and application potential in molecular breeding in rice.

Gene diversity explains variation in biological features of insect killing fungus, Beauveria bassiana

Beauveria bassiana is a species complex whose isolates show considerable natural genetic variability. However, little is known about how this genetic diversity affects the fungus performance. Herein, we characterized the diversity of genes involved in various mechanisms of the infective cycle of 42 isolates that have different growth rates, thermotolerance and virulence. The analysed genes showed general genetic diversity measured as non-synonymous changes (NSC) and copy number variation (CNV), with most of them being subjected to positive episodic diversifying selection. Correlation analyses between NSC or CNV and the isolate virulence, thermotolerance and growth rate revealed that various genes shaped the biological features of the fungus. Lectin-like, mucin signalling, Biotrophy associated and chitinase genes NSCs correlated with the three biological features of B. bassiana. In addition, other genes (i.e. DNA photolyase and cyclophilin B) that had relatively conserved sequences, had variable CNs across the isolates which were correlated with the variability of either virulence or thermotolerance of B. bassiana isolates. The data obtained is important for a better understanding of population structure, ecological and potential impact when isolates are used as mycoinsecticides and can justify industrialization of new isolates.

Analysis of copy number variation in dogs implicates genomic structural variation in the development of anterior cruciate ligament rupture

Anterior cruciate ligament (ACL) rupture is an important condition of the human knee. Second ruptures are common and societal costs are substantial. Canine cranial cruciate ligament (CCL) rupture closely models the human disease. CCL rupture is common in the Labrador Retriever (5.79% prevalence), ~100-fold more prevalent than in humans. Labrador Retriever CCL rupture is a polygenic complex disease, based on genome-wide association study (GWAS) of single nucleotide polymorphism (SNP) markers. Dissection of genetic variation in complex traits can be enhanced by studying structural variation, including copy number variants (CNVs). Dogs are an ideal model for CNV research because of reduced genetic variability within breeds and extensive phenotypic diversity across breeds. We studied the genetic etiology of CCL rupture by association analysis of CNV regions (CNVRs) using 110 case and 164 control Labrador Retrievers. CNVs were called from SNPs using three different programs (PennCNV, CNVPartition, and QuantiSNP). After quality control, CNV calls were combined to create CNVRs using ParseCNV and an association analysis was performed. We found no strong effect CNVRs but found 46 small effect (max(T) permutation P<0.05) CCL rupture associated CNVRs in 22 autosomes; 25 were deletions and 21 were duplications. Of the 46 CCL rupture associated CNVRs, we identified 39 unique regions. Thirty four were identified by a single calling algorithm, 3 were identified by two calling algorithms, and 2 were identified by all three algorithms. For 42 of the associated CNVRs, frequency in the population was <10% while 4 occurred at a frequency in the population ranging from 10–25%. Average CNVR length was 198,872bp and CNVRs covered 0.11 to 0.15% of the genome. All CNVRs were associated with case status. CNVRs did not overlap previous canine CCL rupture risk loci identified by GWAS. Associated CNVRs contained 152 annotated genes; 12 CNVRs did not have genes mapped to CanFam3.1. Using pathway analysis, a cluster of 19 homeobox domain transcript regulator genes was associated with CCL rupture (P = 6.6E-13). This gene cluster influences cranial-caudal body pattern formation during embryonic limb development. Clustered genes were found in 3 CNVRs on chromosome 14 (HoxA), 28 (NKX6-2), and 36 (HoxD). When analysis was limited to deletion CNVRs, the association was strengthened (P = 8.7E-16). This study suggests a component of the polygenic risk of CCL rupture in Labrador Retrievers is associated with small effect CNVs and may include aspects of stifle morphology regulated by homeobox domain transcript regulator genes.

Genetic mutations contributing to non-obstructive azoospermia

Non-obstructive azoospermia is a distinct diagnosis within male infertility in which no sperm is found in the ejaculate as a result of spermatogenesis failure. Because of the increased prevalence of genetic abnormalities in men with non-obstructive azoospermia, male infertility guidelines recommend screening for karyotype abnormalities and Y chromosome microdeletions in this population. Numerous karyotype abnormalities may be present resulting in impaired spermatogenesis, including: Klinefelter syndrome, translocations, and deletions. Y chromosome microdeletions of the AZFa, AZFb, AZFc subregions all can also result in non-obstructive azoospermia with the possibility of sperm being present if only the AZFc subregion is deleted. While these are the two genetic tests recommended by the guidelines, nearly 50%-80% of non-obstructive azoospermia has no identifiable cause and is deemed idiopathic. Several other genetic defects can lead to non-obstructive azoospermia including Kallmann syndrome, mild androgen insensitivity syndrome, and TEX11. While many additional candidate genes have been proposed, many have yet to be verified or are so infrequent in the population that screening is cost-ineffective. Much research is still required in the genetics of non-obstructive azoospermia and will require multi-institutional initiatives to better understand the genetics of condition.

Probability distribution of copy number alterations along the genome: an algorithm to distinguish different tumour profiles

Copy number alterations (CNAs) comprise deletions or amplifications of fragments of genomic material that are particularly common in cancer and play a major contribution in its development and progression. High resolution microarray-based genome-wide technologies have been widely used to detect CNAs, generating complex datasets that require further steps to allow for the determination of meaningful results. In this work, we propose a methodology to determine common regions of CNAs from these datasets, that in turn are used to infer the probability distribution of disease profiles in the population. This methodology was validated using simulated data and assessed using real data from Head and Neck Squamous Cell Carcinoma and Lung Adenocarcinoma, from the TCGA platform. Probability distribution profiles were produced allowing for the distinction between different phenotypic groups established within that cohort. This method may be used to distinguish between groups in the diseased population, within well-established degrees of confidence. The application of such methods may be of greater value in the clinical context both as a diagnostic or prognostic tool and, even as a useful way for helping to establish the most adequate treatment and care plans.

Retrospective clinical trial experimentally validates glioblastoma genome-wide pattern of DNA copy-number alterations predictor of survival

Modeling of genomic profiles from the Cancer Genome Atlas (TCGA) by using recently developed mathematical frameworks has associated a genome-wide pattern of DNA copy-number alterations with a shorter, roughly one-year, median survival time in glioblastoma (GBM) patients. Here, to experimentally test this relationship, we whole-genome sequenced DNA from tumor samples of patients. We show that the patients represent the U.S. adult GBM population in terms of most normal and disease phenotypes. Intratumor heterogeneity affects ≈ 11 % and profiling technology and reference human genome specifics affect <1% of the classifications of the tumors by the pattern, where experimental batch effects normally reduce the reproducibility, i.e., precision, of classifications based upon between one to a few hundred genomic loci by >30%. With a 2.25-year Kaplan-Meier median survival difference, a 3.5 univariate Cox hazard ratio, and a 0.78 concordance index, i.e., accuracy, the pattern predicts survival better than and independent of age at diagnosis, which has been the best indicator since 1950. The prognostic classification by the pattern may, therefore, help to manage GBM pseudoprogression. The diagnostic classification may help drugs progress to regulatory approval. The therapeutic predictions, of previously unrecognized targets that are correlated with survival, may lead to new drugs. Other methods missed this relationship in the roughly 3B-nucleotide genomes of the small, order of magnitude of 100, patient cohorts, e.g., from TCGA. Previous attempts to associate GBM genotypes with patient phenotypes were unsuccessful. This is a proof of principle that the frameworks are uniquely suitable for discovering clinically actionable genotype-phenotype relationships.

Cellular and genomic approaches for exploring structural chromosomal rearrangements

Human chromosomes are arranged in a linear and conserved sequence order that undergoes further spatial folding within the three-dimensional space of the nucleus. Although structural variations in this organization are an important source of natural genetic diversity, cytogenetic aberrations can also underlie a number of human diseases and disorders. Approaches for studying chromosome structure began half a century ago with karyotyping of Giemsa-banded chromosomes and has now evolved to encompass high-resolution fluorescence microscopy, reporter-based assays, and next-generation DNA sequencing technologies. Here, we provide a general overview of experimental methods at different resolution and sensitivity scales and discuss how they can be complemented to provide synergistic insight into the study of human chromosome structural rearrangements. These approaches range from kilobase-level resolution DNA fluorescence in situ hybridization (FISH)-based imaging approaches of individual cells to genome-wide sequencing strategies that can capture nucleotide-level information from diverse sample types. Technological advances coupled to the combinatorial use of multiple methods have resulted in the discovery of new rearrangement classes along with mechanistic insights into the processes that drive structural alterations in the human genome.

A spontaneous complex structural variant in rcan-1 increases exploratory behavior and laboratory fitness of Caenorhabditis elegans

Over long evolutionary timescales, major changes to the copy number, function, and genomic organization of genes occur, however, our understanding of the individual mutational events responsible for these changes is lacking. In this report, we study the genetic basis of adaptation of two strains of C. elegans to laboratory food sources using competition experiments on a panel of 89 recombinant inbred lines (RIL). Unexpectedly, we identified a single RIL with higher relative fitness than either of the parental strains. This strain also displayed a novel behavioral phenotype, resulting in higher propensity to explore bacterial lawns. Using bulk-segregant analysis and short-read resequencing of this RIL, we mapped the change in exploration behavior to a spontaneous, complex rearrangement of the rcan-1 gene that occurred during construction of the RIL panel. We resolved this rearrangement into five unique tandem inversion/duplications using Oxford Nanopore long-read sequencing. rcan-1 encodes an ortholog to human RCAN1/DSCR1 calcipressin gene, which has been implicated as a causal gene for Down syndrome. The genomic rearrangement in rcan-1 creates two complete and two truncated versions of the rcan-1 coding region, with a variety of modified 5’ and 3’ non-coding regions. While most copy-number variations (CNVs) are thought to act by increasing expression of duplicated genes, these changes to rcan-1 ultimately result in the reduction of its whole-body expression due to changes in the upstream regions. By backcrossing this rearrangement into a common genetic background to create a near isogenic line (NIL), we demonstrate that both the competitive advantage and exploration behavioral changes are linked to this complex genetic variant. This NIL strain does not phenocopy a strain containing an rcan-1 loss-of-function allele, which suggests that the residual expression of rcan-1 is necessary for its fitness effects. Our results demonstrate how colonization of new environments, such as those encountered in the laboratory, can create evolutionary pressure to modify gene function. This evolutionary mismatch can be resolved by an unexpectedly complex genetic change that simultaneously duplicates and diversifies a gene into two uniquely regulated genes. Our work shows how complex rearrangements can act to modify gene expression in ways besides increased gene dosage.

Evolution acts on genetic variants that modify phenotypes that increase the likelihood of staying alive and passing on these genetic changes to subsequent generations (i.e. fitness). There is general interest in understanding the types of genetic variants that can increase fitness in specific environments. One route that fitness can be increased is through changes in behavior, such as finding new food sources. Here, we identify a spontaneous genetic change that increases exploration behavior and fitness of animals in laboratory environments. Interestingly, this genetic change is not a simple genetic change that deletes or changes the sequence of a protein product, but rather a complex structural variant that simultaneously duplicates the rcan-1 gene and also modifies its expression in a number of tissues. Our work demonstrates how a complex structural change can duplicate a gene, modify the DNA control regions that determine its cellular sites of action, and confer a fitness advantage that could lead to its spread in a population.

Evaluation of the functional effects of genetic variants‒missense and nonsense SNPs, indels and copy number variations‒in the gene encoding human deoxyribonuclease I potentially implicated in autoimmunity

Genetic variants, such as single nucleotide polymorphisms (SNPs), in the deoxyribonuclease I (DNase I) gene which remarkably reduce or abolish the activity are assumed to be substantially responsible for the genetic backgrounds determining susceptibility to autoimmune dysfunction. Here, we evaluated many genetic variants, including missense and nonsense SNPs, and indel (inframe) variants in the gene, potentially implicated in autoimmune diseases as functional variants resulting in altered activity levels. Eighteen missense and 7 nonsense SNPs, and 9 indel (inframe) variants were found to result in loss of function and disappearance of DNase I activity. Furthermore, considering the positions in the DNase I protein corresponding to the various nonsense SNPs, all of the other nonsense SNPs and frameshift variants registered in the Ensembl database (https://asia.ensembl.org) appear likely to exert a pathogenetic effect through loss of the activity. Accordingly, a total of 60 genetic variants in the DNase 1 gene (DNASE1) inducing abolishment or marked reduction of the DNase I activity could be identified as genetic risk factors for autoimmunity, irrespective of how sparsely they were distributed in the population. It was noteworthy that SNP p.Gln244Arg, reportedly associated with autoimmunity and reducing the activity to about half of that of the wild type, and SNP p.Arg107Gly, abolishing the activity completely, were distributed worldwide and in African populations at the polymorphic level, respectively. On the other hand, with regard to copy number variations in DNASE1 where loss of copy leads to a reduction of the in vivo enzyme activity, only 2 diploid copy numbers were distributed in Japanese and German populations, demonstrating no loss of copy. These exhaustive data for genetic variants in DNASE1 resulting in loss or marked reduction of the DNase I activity are highly informative when considering genetic predisposition leading to autoimmune dysfunction.

GSVD- and tensor GSVD-uncovered patterns of DNA copy-number alterations predict adenocarcinomas survival in general and in response to platinum

More than a quarter of lung, uterine, and ovarian adenocarcinoma (LUAD, USEC, and OV) tumors are resistant to platinum drugs. Only recently and only in OV, patterns of copy-number alterations that predict survival in response to platinum were discovered, and only by using the tensor GSVD to compare Agilent microarray platform-matched profiles of patient-matched normal and primary tumor DNA. Here, we use the GSVD to compare whole-genome sequencing (WGS) and Affymetrix microarray profiles of patient-matched normal and primary LUAD, USEC, and OV tumor DNA. First, the GSVD uncovers patterns similar to one Agilent OV pattern, where a loss of most of the chromosome arm 6p combined with a gain of 12p encode for transformation. Like the Agilent OV pattern, the WGS LUAD and Affymetrix LUAD, USEC, and OV patterns are correlated with shorter survival, in general and in response to platinum. Like the tensor GSVD, the GSVD separates these tumor-exclusive genotypes from experimental inconsistencies. Second, by identifying the shorter survival phenotypes among the WGS- and Affymetrix-profiled tumors, the Agilent pattern proves to be a technology-independent predictor of survival, independent also of the best other indicator at diagnosis, i.e., stage. Third, like no other indicator, the pattern predicts the overall survival of OV patients experiencing progression-free survival, in general and in response to platinum. We conclude that comparative spectral decompositions, such as the GSVD and tensor GSVD, underlie a mathematically universal description of the relationships between a primary tumor's genotype and a patient's overall survival phenotype, which other methods miss.

Integrative multi-omic analysis identifies new drivers and pathways in molecularly distinct subtypes of ALS

Amyotrophic lateral sclerosis (ALS) is an incurable and fatal neurodegenerative disease. Increasing the chances of success for future clinical strategies requires more in-depth knowledge of the molecular basis underlying disease heterogeneity. We recently laid the foundation for a molecular taxonomy of ALS by whole-genome expression profiling of motor cortex from sporadic ALS (SALS) patients. Here, we analyzed copy number variants (CNVs) occurring in the same patients, by using a customized exon-centered comparative genomic hybridization array (aCGH) covering a large panel of ALS-related genes. A large number of novel and known disease-associated CNVs were detected in SALS samples, including several subgroup-specific loci, suggestive of a great divergence of two subgroups at the molecular level. Integrative analysis of copy number profiles with their associated transcriptomic data revealed subtype-specific genomic perturbations and candidate driver genes positively correlated with transcriptional signatures, suggesting a strong interaction between genomic and transcriptomic events in ALS pathogenesis. The functional analysis confirmed our previous pathway-based characterization of SALS subtypes and identified 24 potential candidates for genomic-based patient stratification. To our knowledge, this is the first comprehensive "omics" analysis of molecular events characterizing SALS pathology, providing a road map to facilitate genome-guided personalized diagnosis and treatments for this devastating disease.

Copy number variants in autism spectrum disorders

In recent years, there has been an explosive increase in genetic studies related to autism spectrum disorder (ASD). This implicated the accumulation of a large amount of molecular data that may be used to verify various hypotheses and models developed to explore the complex genetic component of ASD. Several lines of evidence support the view that structural genomic variation contributes to the pathogenesis of ASD. The introduction of more sophisticated techniques for whole-genome screening, including array comparative genome hybridization and high-resolution single nucleotide polymorphism analysis, has allowed to identify an increasing number of ASD susceptibility loci. Copy number variants (CNVs) are the most common type of structural variation in the human genome and are considered important contributors to the pathogenesis of neurodevelopmental disorders, including ASD. In this review, we describe the accumulated evidence concerning the genetic events associated with ASD, and summarize current knowledge about the clinical relevance of CNVs in these disorders.

Familial segregation of a 5q15-q21.2 deletion associated with facial dysmorphism and speech delay

We report a two-generation family with four females harboring an 8.5Mb heterozygous deletion of 5q15-q21.2 who present with dysmorphic craniofacial features and speech delay. We hypothesize haploinsufficiency of CHD1 to be contributing to the clinical features observed in this family.

Low genetic heterogeneity of copy number variations (CNVs) in the genes encoding the human deoxyribonucleases 1-like 3 and II potentially relevant to autoimmunity

Deoxyribonucleases (DNases) might play a role in prevention of autoimmune conditions such as systemic lupus erythematosus through clearance of cell debris resulting from apoptosis and/or necrosis. Previous studies have suggested that variations in the in vivo activities of DNases I-like 3(1L3) and II have an impact on autoimmune-related conditions. The genes for these DNases are known to show copy number variations (CNVs) whereby copy loss leads to a reduction of the in vivo activities of the enzymes, thereby possibly affecting the pathophysiological background of autoimmune diseases. Using a simple newly developed quantitative real-time PCR method, we investigated the distributions of the CNVs for DNASE1L3 and DNASE2 in Japanese and German populations. It was found that only 2 diploid copy numbers for all of these DNASE CNVs was distributed in both of the study populations; no copy loss or gain was evident for any of the autoimmune-related DNase genes. Therefore, it was demonstrated that these human autoimmune-related DNase genes show low genetic diversity of CNVs resulting in alterations of the in vivo levels of DNase activity.

Fine-Scale Characterization of Genomic Structural Variation in the Human Genome Reveals Adaptive and Biomedically Relevant Hotspots

Genomic structural variants (SVs) are distributed nonrandomly across the human genome. The "hotspots" of SVs have been implicated in evolutionary innovations, as well as medical conditions. However, the evolutionary and biomedical features of these hotspots remain incompletely understood. Here, we analyzed data from 2,504 genomes to construct a refined map of 1,148 SV hotspots in human genomes. We confirmed that segmental duplication-related nonallelic homologous recombination is an important mechanistic driver of SV hotspot formation. However, to our surprise, we also found that a majority of SVs in hotspots do not form through such recombination-based mechanisms, suggesting diverse mechanistic and selective forces shaping hotspots. Indeed, our evolutionary analyses showed that the majority of SV hotspots are within gene-poor regions and evolve under relaxed negative selection or neutrality. However, we still found a small subset of SV hotspots harboring genes that are enriched for anthropologically crucial functions and evolve under geography-specific and balancing adaptive forces. These include two independent hotspots on different chromosomes affecting alpha and beta hemoglobin gene clusters. Biomedically, we found that the SV hotspots coincide with breakpoints of clinically relevant, large de novo SVs, significantly more often than genome-wide expectations. For example, we showed that the breakpoints of multiple large SVs, which lead to idiopathic short stature, coincide with SV hotspots. Therefore, the mutational instability in SV hotpots likely enables chromosomal breaks that lead to pathogenic structural variation formations. Overall, our study contributes to a better understanding of the mutational and adaptive landscape of the genome.

Insights into genetics, human biology and disease gleaned from family based genomic studies

Identifying genes and variants contributing to rare disease phenotypes and Mendelian conditions informs biology and medicine, yet potential phenotypic consequences for variation of >75% of the ~20,000 annotated genes in the human genome are lacking. Technical advances to assess rare variation genome-wide, particularly exome sequencing (ES), enabled establishment in the United States of the National Institutes of Health (NIH)-supported Centers for Mendelian Genomics (CMGs) and have facilitated collaborative studies resulting in novel "disease gene" discoveries. Pedigree-based genomic studies and rare variant analyses in families with suspected Mendelian conditions have led to the elucidation of hundreds of novel disease genes and highlighted the impact of de novo mutational events, somatic variation underlying nononcologic traits, incompletely penetrant alleles, phenotypes with high locus heterogeneity, and multilocus pathogenic variation. Herein, we highlight CMG collaborative discoveries that have contributed to understanding both rare and common diseases and discuss opportunities for future discovery in single-locus Mendelian disorder genomics. Phenotypic annotation of all human genes; development of bioinformatic tools and analytic methods; exploration of non-Mendelian modes of inheritance including reduced penetrance, multilocus variation, and oligogenic inheritance; construction of allelic series at a locus; enhanced data sharing worldwide; and integration with clinical genomics are explored. Realizing the full contribution of rare disease research to functional annotation of the human genome, and further illuminating human biology and health, will lay the foundation for the Precision Medicine Initiative.

TBX6-associated congenital scoliosis (TACS) as a clinically distinguishable subtype of congenital scoliosis: further evidence supporting the compound inheritance and TBX6 gene dosage model

PURPOSE

To characterize clinically measurable endophenotypes, implicating the TBX6 compound inheritance model.

METHODS

Patients with congenital scoliosis (CS) from China(N = 345, cohort 1), Japan (N = 142, cohort 2), and the United States (N = 10, cohort 3) were studied. Clinically measurable endophenotypes were compared according to the TBX6 genotypes. A mouse model for Tbx6 compound inheritance (N = 52) was investigated by micro computed tomography (micro-CT). A clinical diagnostic algorithm (TACScore) was developed to assist in clinical recognition of TBX6-associated CS (TACS).

RESULTS

In cohort 1, TACS patients (N = 33) were significantly younger at onset than the remaining CS patients (P = 0.02), presented with one or more hemivertebrae/butterfly vertebrae (P = 4.9 × 10^‒8), and exhibited vertebral malformations involving the lower part of the spine (T8-S5, P = 4.4 × 10^‒3); observations were confirmed in two replication cohorts. Simple rib anomalies were prevalent in TACS patients (P = 3.1 × 10^‒7), while intraspinal anomalies were uncommon (P = 7.0 × 10^‒7). A clinically usable TACScore was developed with an area under the curve (AUC) of 0.9 (P = 1.6 × 10^‒15). A Tbx6^{-/mh (mild-hypomorphic)} mouse model supported that a gene dosage effect underlies the TACS phenotype.

CONCLUSION

TACS is a clinically distinguishable entity with consistent clinically measurable endophenotypes. The type and distribution of vertebral column abnormalities in TBX6/Tbx6 compound inheritance implicate subtle perturbations in gene dosage as a cause of spine developmental birth defects responsible for about 10% of CS.

Guidelines for DNA recombination and repair studies: Cellular assays of DNA repair pathways

Understanding the plasticity of genomes has been greatly aided by assays for recombination, repair and mutagenesis. These assays have been developed in microbial systems that provide the advantages of genetic and molecular reporters that can readily be manipulated. Cellular assays comprise genetic, molecular, and cytological reporters. The assays are powerful tools but each comes with its particular advantages and limitations. Here the most commonly used assays are reviewed, discussed, and presented as the guidelines for future studies.

Genomic Structural Variations Within Five Continental Populations of Drosophila melanogaster

Chromosomal structural variations (SV) including insertions, deletions, inversions, and translocations occur within the genome and can have a significant effect on organismal phenotype. Some of these effects are caused by structural variations containing genes. Large structural variations represent a significant amount of the genetic diversity within a population. We used a global sampling of Drosophila melanogaster (Ithaca, Zimbabwe, Beijing, Tasmania, and Netherlands) to represent diverse populations within the species. We used long-read sequencing and optical mapping technologies to identify SVs in these genomes. Among the five lines examined, we found an average of 2,928 structural variants within these genomes. These structural variations varied greatly in size and location, included many exonic regions, and could impact adaptation and genomic evolution.

Mathematically universal and biologically consistent astrocytoma genotype encodes for transformation and predicts survival phenotype

DNA alterations have been observed in astrocytoma for decades. A copy-number genotype predictive of a survival phenotype was only discovered by using the generalized singular value decomposition (GSVD) formulated as a comparative spectral decomposition. Here, we use the GSVD to compare whole-genome sequencing (WGS) profiles of patient-matched astrocytoma and normal DNA. First, the GSVD uncovers a genome-wide pattern of copy-number alterations, which is bounded by patterns recently uncovered by the GSVDs of microarray-profiled patient-matched glioblastoma (GBM) and, separately, lower-grade astrocytoma and normal genomes. Like the microarray patterns, the WGS pattern is correlated with an approximately one-year median survival time. By filling in gaps in the microarray patterns, the WGS pattern reveals that this biologically consistent genotype encodes for transformation via the Notch together with the Ras and Shh pathways. Second, like the GSVDs of the microarray profiles, the GSVD of the WGS profiles separates the tumor-exclusive pattern from normal copy-number variations and experimental inconsistencies. These include the WGS technology-specific effects of guanine-cytosine content variations across the genomes that are correlated with experimental batches. Third, by identifying the biologically consistent phenotype among the WGS-profiled tumors, the GBM pattern proves to be a technology-independent predictor of survival and response to chemotherapy and radiation, statistically better than the patient's age and tumor's grade, the best other indicators, and MGMT promoter methylation and IDH1 mutation. We conclude that by using the complex structure of the data, comparative spectral decompositions underlie a mathematically universal description of the genotype-phenotype relations in cancer that other methods miss.

Novel applications of array comparative genomic hybridization in molecular diagnostics

INTRODUCTION

In 2004, the implementation of array comparative genomic hybridization (array comparative genome hybridization [CGH]) into clinical practice marked a new milestone for genetic diagnosis. Array CGH and single-nucleotide polymorphism (SNP) arrays enable genome-wide detection of copy number changes in a high resolution, and therefore microarray has been recognized as the first-tier test for patients with intellectual disability or multiple congenital anomalies, and has also been applied prenatally for detection of clinically relevant copy number variations in the fetus. Area covered: In this review, the authors summarize the evolution of array CGH technology from their diagnostic laboratory, highlighting exonic SNP arrays developed in the past decade which detect small intragenic copy number changes as well as large DNA segments for the region of heterozygosity. The applications of array CGH to human diseases with different modes of inheritance with the emphasis on autosomal recessive disorders are discussed. Expert commentary: An exonic array is a powerful and most efficient clinical tool in detecting genome wide small copy number variants in both dominant and recessive disorders. However, whole-genome sequencing may become the single integrated platform for detection of copy number changes, single-nucleotide changes as well as balanced chromosomal rearrangements in the near future.

TRKA expression and NTRK1 gene copy number across solid tumours

AIMS

Neurotrophic Tropomyosin Kinase Receptor 1 (NTRK1) gene encodes for the protein Tropomyosin-related kinase A (TRKA). Deregulated activity of TRKA has been shown to have oncogenic potential. We present here the results of an immunohistochemical (IHC) observational cohort study of TRKA expression together with gene copy number (GCN) assessment in various solid tumours.

METHODS

Formalin-fixed, paraffin-embedded consecutive samples of different tumour types were tested for TRKA expression. Samples showing TRKA IHC staining in at least 10% of cells were analysed by fluorescence in situ hybridisation to assess NTRK1 gene rearrangements and/or individual GCN gain. All patients underwent this molecular assessment within the phase I ALKA-001 clinical trial.

RESULTS

1043 samples were tested and annotation for histology was available in 1023. Most of the samples were colorectal adenocarcinoma (CRC) (n=550, 52.7%) and lung adenocarcinoma (n=312, 29.9%). 24 samples (2.3%) were biliary tract carcinoma (BTC). Overall, 17 (1.6%) samples were characterised by TRKA IHC expression (four weak, eight moderate, five strong): 9/17 lung adenocarcinoma, 3/17 CRC, 3/17 BTC, 1/17 thyroid cancer and 1/17 cancer of unknown primary. Of these, 1/17 with strong TRKA IHC staining displayed NTRK1 gene rearrangement and 15/17 NTRK1 GCN gain by FISH. No correlation was found between intensity of TRKA IHC staining and number of copies of NTRK1.

CONCLUSIONS

TRKA expression can be found in 1.6% of solid tumours and can be paralleled by NTRK1 gene rearrangements or mostly GCN gain. The prognostic and translational therapeutic impact of the latter remains to be established.

Non Random Distribution of DMD Deletion Breakpoints and Implication of Double Strand Breaks Repair and Replication Error Repair Mechanisms

BACKGROUND

Dystrophinopathies are mostly caused by copy number variations, especially deletions, in the dystrophin gene (DMD). Despite the large size of the gene, deletions do not occur randomly but mainly in two hot spots, the main one involving exons 45 to 55. The underlying mechanisms are complex and implicate two main mechanisms: Non-homologous end joining (NHEJ) and micro-homology mediated replication-dependent recombination (MMRDR).

OBJECTIVE

Our goals were to assess the distribution of intronic breakpoints (BPs) in the genomic sequence of the main hot spot of deletions within DMD gene and to search for specific sequences at or near to BPs that might promote BP occurrence or be associated with DNA break repair.

METHODS

Using comparative genomic hybridization microarray, 57 deletions within the intron 44 to 55 region were mapped. Moreover, 21 junction fragments were sequenced to search for specific sequences.

RESULTS

Non-randomly distributed BPs were found in introns 44, 47, 48, 49 and 53 and 50% of BPs clustered within genomic regions of less than 700bp. Repeated elements (REs), known to promote gene rearrangement via several mechanisms, were present in the vicinity of 90% of clustered BPs and less frequently (72%) close to scattered BPs, illustrating the important role of such elements in the occurrence of DMD deletions. Palindromic and TTTAAA sequences, which also promote DNA instability, were identified at fragment junctions in 20% and 5% of cases, respectively. Micro-homologies (76%) and insertions or deletions of small sequences were frequently found at BP junctions.

CONCLUSIONS

Our results illustrate, in a large series of patients, the important role of RE and other genomic features in DNA breaks, and the involvement of different mechanisms in DMD gene deletions: Mainly replication error repair mechanisms, but also NHEJ and potentially aberrant firing of replication origins. A combination of these mechanisms may also be possible.