How structural variants shape avian phenotypes: Lessons from model systems

Despite receiving significant recent attention, the relevance of structural variation (SV) in driving phenotypic diversity remains understudied, although recent advances in long-read sequencing, bioinformatics and pangenomic approaches have enhanced SV detection. We review the role of SVs in shaping phenotypes in avian model systems, and identify some general patterns in SV type, length and their associated traits. We found that most of the avian SVs so far identified are short indels in chickens, which are frequently associated with changes in body weight and plumage colouration. Overall, we found that relatively short SVs are more frequently detected, likely due to a combination of their prevalence compared to large SVs, and a detection bias, stemming primarily from the widespread use of short-read sequencing and associated analytical methods. SVs most commonly involve non-coding regions, especially introns, and when patterns of inheritance were reported, SVs associated primarily with dominant discrete traits. We summarise several examples of phenotypic convergence across different species, mediated by different SVs in the same or different genes and different types of changes in the same gene that can lead to various phenotypes. Complex rearrangements and supergenes, which can simultaneously affect and link several genes, tend to have pleiotropic phenotypic effects. Additionally, SVs commonly co-occur with single-nucleotide polymorphisms, highlighting the need to consider all types of genetic changes to understand the basis of phenotypic traits. We end by summarising expectations for when long-read technologies become commonly implemented in non-model birds, likely leading to an increase in SV discovery and characterisation. The growing interest in this subject suggests an increase in our understanding of the phenotypic effects of SVs in upcoming years.

Estrogen receptor alpha mutations, truncations, heterodimers, and therapies

Annual breast cancer (BCa) deaths have declined since its apex in 1989 concomitant with widespread adoption of hormone therapies that target estrogen receptor alpha (ERα), the prominent nuclear receptor expressed in ∼80% of BCa. However, up to ∼50% of ER + patients with high-risk disease experience post endocrine therapy relapse and metastasis to distant organs. The vast majority of BCa mortality occurs in this setting, highlighting the inadequacy of current therapies. Genomic abnormalities to ESR1, the gene encoding ERα, emerge under prolonged selective pressure to enable endocrine therapy resistance. These genetic lesions include focal gene amplifications, hotspot missense mutations in the ligand binding domain, truncations, fusions, and complex interactions with other nuclear receptors. Tumor cells utilize aberrant ERα activity to proliferate, spread and evade therapy in BCa as well as other cancers. Cutting edge studies on ERα structural and transcriptional relationships are being harnessed to produce new therapies that have shown benefits in patients with ESR1 hotspot mutations. In this review we discuss the history of ERα, current research unlocking unknown aspects of ERα signaling including the structural basis for receptor antagonism, and future directions of ESR1 investigation. In addition, we discuss the development of endocrine therapies from their inception to present day and survey new avenues of drug development to improve pharmaceutical profiles, targeting, and efficacy.

Olfactory receptor genes and chromosome 11 structural aberrations: Players or spectators?

The largest multi-gene family in metazoans is the family of olfactory receptor (OR) genes. Human ORs are organized in clusters over most chromosomes and seem to include >0.1% the human genome. Because 369 out of 856 OR genes are mapped on chromosome 11 (HSA11), we sought to determine whether they mediate structural rearrangements involving this chromosome. To this aim, we analyzed 220 specimens collected during diagnostic procedures involving structural rearrangements of chromosome 11. A total of 222 chromosomal abnormalities were included, consisting of inversions, deletions, translocations, duplications, and one insertion, detected by conventional chromosome analysis and/or fluorescence in situ hybridization (FISH) and array comparative genomic hybridization (array-CGH). We verified by bioinformatics and statistical approaches the occurrence of breakpoints in cytobands with or without OR genes. We found that OR genes are not involved in chromosome 11 reciprocal translocations, suggesting that different DNA motifs and mechanisms based on homology or non-homology recombination can cause chromosome 11 structural alterations. We also considered the proximity between the chromosomal territories of chromosome 11 and its partner chromosomes involved in the translocations by using the deposited Hi-C data concerning the possible occurrence of chromosome interactions. Interestingly, most of the breakpoints are located in regions highly involved in chromosome interactions. Further studies should be carried out to confirm the potential role of chromosome territories' proximity in promoting genome structural variation, so fundamental in our understanding of the molecular basis of medical genetics and evolutionary genetics.

Whole-Genome Sequencing Analysis of Male Breast Cancer Unveils Novel Structural Events and Potential Therapeutic Targets

The molecular characterization of male breast cancer (MaBC) has received limited attention in research, mostly because of its low incidence rate, accounting for only 0.5% to 1% of all reported cases of breast cancer each year. Managing MaBC presents significant challenges, with most treatment protocols being adapted from those developed for female breast cancer. Utilizing whole-genome sequencing (WGS) and state-of-the-art analyses, the genomic features of 10 MaBC cases (n = 10) were delineated and correlated with clinical and histopathologic characteristics. Using fluorescence in situ hybridization, an additional cohort of 18 patients was interrogated to supplement WGS findings. The genomic landscape of MaBC uncovered significant genetic alterations that could influence diagnosis and treatment. We found common somatic mutations in key driver genes, such as FAT1, GATA3, SMARCA4, and ARID2. Our study also mapped out structural variants that impact cancer-associated genes, such as ARID1A, ESR1, GATA3, NTRK1, and NF1. Using a WGS-based classifier, homologous recombination deficiency (HRD) was identified in 2 cases, both presenting with deleterious variants in BRCA2. Noteworthy was the observation of FGFR1 amplification in 21% of cases. Altogether, we identified at least 1 potential therapeutic target in 8 of the 10 cases, including high tumor mutational burden, FGFR1 amplification, and HRD. Our study is the first WGS characterization of MaBC, which uncovered potentially relevant variants, including structural events in cancer genes, HRD signatures, and germline pathogenic mutations. Our results demonstrate unique genetic markers and potential treatment targets in MaBC, thereby underlining the necessity of tailoring treatment strategies for this understudied patient population. These WGS-based findings add to the growing knowledge of MaBC genomics and highlight the need to expand research on this type of cancer.

Leveraging multiple approaches for the detection of pathogenic deep intronic variants in developmental and epileptic encephalopathies: A case report

About 50% of individuals with developmental and epileptic encephalopathies (DEEs) are unsolved following genetic testing. Deep intronic variants, defined as >100 bp from exon-intron junctions, contribute to disease by affecting the splicing of mRNAs in clinically relevant genes. Identifying deep intronic pathogenic variants is challenging and resource intensive, and interpretation is difficult due to limited functional annotations. We aimed to identify deep intronic variants in individuals suspected to have unsolved single gene DEEs. In a research cohort of unsolved cases of DEEs, we searched for children with a DEE syndrome predominantly caused by variants in specific genes in >80% of described cases. We identified two children with Dravet syndrome and one individual with classic lissencephaly. Multiple sequencing and bioinformatics strategies were employed to interrogate intronic regions in SCN1A and PAFAH1B1. A novel de novo deep intronic 12 kb deletion in PAFAH1B1 was identified in the individual with lissencephaly. We showed experimentally that the deletion disrupts mRNA splicing, which results in partial intron retention after exon 2 and disruption of the highly conserved LisH motif. We demonstrate that targeted interrogation of deep intronic regions using multiple genomics technologies, coupled with functional analysis, can reveal hidden causes of unsolved monogenic DEE syndromes. PLAIN LANGUAGE SUMMARY: Deep intronic variants can cause disease by affecting the splicing of mRNAs in clinically relevant genes. A deep intronic deletion that caused abnormal splicing of the PAFAH1B1 gene was identified in a patient with classic lissencephaly. Our findings reinforce that targeted interrogation of deep intronic regions and functional analysis can reveal hidden causes of unsolved epilepsy syndromes.

Germline structural variation globally impacts the cancer transcriptome including disease-relevant genes

Germline variation and somatic alterations contribute to the molecular profile of cancers. We combine RNA with whole genome sequencing across 1,218 cancer patients to determine the extent germline structural variants (SVs) impact expression of nearby genes. For hundreds of genes, recurrent and common germline SV breakpoints within 100 kb associate with increased or decreased expression in tumors spanning various tissues of origin. A significant fraction of germline SV expression associations involves duplication of intergenic enhancers or 3' UTR disruption. Genes altered by both somatic and germline SVs include ATRX and CEBPA. Genes essential in cancer cell lines include BARD1 and IRS2. Genes with both expression and germline SV breakpoint patterns associated with patient survival include GCLM. Our results capture a class of phenotypic variation at work in the disease setting, including genes with cancer roles. Specific germline SVs represent potential cancer risk variants for genetic testing, including those involving genes with targeting implications.

On- and off-target effects of paired CRISPR-Cas nickase in primary human cells

Undesired on- and off-target effects of CRISPR-Cas nucleases remain a challenge in genome editing. While the use of Cas9 nickases has been shown to minimize off-target mutagenesis, their use in therapeutic genome editing has been hampered by a lack of efficacy. To overcome this limitation, we and others have developed double-nickase-based strategies to generate staggered DNA double-strand breaks to mediate gene disruption or gene correction with high efficiency. However, the impact of paired single-strand nicks on genome integrity has remained largely unexplored. Here, we developed a novel CAST-seq pipeline, dual CAST, to characterize chromosomal aberrations induced by paired CRISPR-Cas9 nickases at three different loci in primary keratinocytes derived from patients with epidermolysis bullosa. While targeting COL7A1, COL17A1, or LAMA3 with Cas9 nucleases caused previously undescribed chromosomal rearrangements, no chromosomal translocations were detected following paired-nickase editing. While the double-nicking strategy induced large deletions/inversions within a 10 kb region surrounding the target sites at all three loci, similar to the nucleases, the chromosomal on-target aberrations were qualitatively different and included a high proportion of insertions. Taken together, our data indicate that double-nickase approaches combine efficient editing with greatly reduced off-target effects but still leave substantial chromosomal aberrations at on-target sites.

Long read sequencing on its way to the routine diagnostics of genetic diseases

The clinical application of technological progress in the identification of DNA alterations has always led to improvements of diagnostic yields in genetic medicine. At chromosome side, from cytogenetic techniques evaluating number and gross structural defects to genomic microarrays detecting cryptic copy number variants, and at molecular level, from Sanger method studying the nucleotide sequence of single genes to the high-throughput next-generation sequencing (NGS) technologies, resolution and sensitivity progressively increased expanding considerably the range of detectable DNA anomalies and alongside of Mendelian disorders with known genetic causes. However, particular genomic regions (i.e., repetitive and GC-rich sequences) are inefficiently analyzed by standard genetic tests, still relying on laborious, time-consuming and low-sensitive approaches (i.e., southern-blot for repeat expansion or long-PCR for genes with highly homologous pseudogenes), accounting for at least part of the patients with undiagnosed genetic disorders. Third generation sequencing, generating long reads with improved mappability, is more suitable for the detection of structural alterations and defects in hardly accessible genomic regions. Although recently implemented and not yet clinically available, long read sequencing (LRS) technologies have already shown their potential in genetic medicine research that might greatly impact on diagnostic yield and reporting times, through their translation to clinical settings. The main investigated LRS application concerns the identification of structural variants and repeat expansions, probably because techniques for their detection have not evolved as rapidly as those dedicated to single nucleotide variants (SNV) identification: gold standard analyses are karyotyping and microarrays for balanced and unbalanced chromosome rearrangements, respectively, and southern blot and repeat-primed PCR for the amplification and sizing of expanded alleles, impaired by limited resolution and sensitivity that have not been significantly improved by the advent of NGS. Nevertheless, more recently, with the increased accuracy provided by the latest product releases, LRS has been tested also for SNV detection, especially in genes with highly homologous pseudogenes and for haplotype reconstruction to assess the parental origin of alleles with de novo pathogenic variants. We provide a review of relevant recent scientific papers exploring LRS potential in the diagnosis of genetic diseases and its potential future applications in routine genetic testing.

T/myeloid mixed phenotype acute leukaemia harbouring TLX3::BCL11B with TLX3 activation

T/myeloid mixed phenotype acute leukaemia (MPAL) is a rare aggressive acute leukaemia with poorly understood pathogenesis. Herein, we report two cases of T/myeloid MPAL harbouring BCL11B-associated structural variants that activate TLX3 (TLX3::BCL11B-TLX3-activation) by genome sequencing and transcriptomic analyses. Both patients were young males with extramedullary involvement. Cooperative gene alterations characteristic of T/myeloid MPAL and T-lymphoblastic leukaemia (T-ALL) were detected. Both patients achieved initial remission following lineage-matched ALL-based therapy with one patient requiring a lineage-switched myeloid-based therapy. Our study is the first to demonstrate the clinicopathological and genomic features of TLX3::BCL11B-TLX3-activated T/myeloid MPAL and provide insights into leukaemogenesis.

Optical genome mapping identifies structural variants in potentially new cancer predisposition candidate genes in pediatric cancer patients

Genetic predisposition is one of the major risk factors for pediatric cancer, with ~10% of children being carriers of a predisposing germline alteration. It is likely that this is the tip of the iceberg and many children are underdiagnosed, as most of the analysis focuses on single or short nucleotide variants, not considering the full spectrum of DNA alterations. Hence, we applied optical genome mapping (OGM) to our cohort of 34 pediatric cancer patients to perform an unbiased germline screening and analyze the frequency of structural variants (SVs) and their impact on cancer predisposition. All children were clinically highly suspicious for germline alterations (concomitant conditions or congenital anomalies, positive family cancer history, particular cancer type, synchronous or metachronous tumors), but whole exome sequencing (WES) had failed to detect pathogenic variants in cancer predisposing genes. OGM detected a median of 49 rare SVs (range 27-149) per patient. By analysis of 18 patient-parent trios, we identified three de novo SVs. Moreover, we discovered a likely pathogenic deletion of exon 3 in the known cancer predisposition gene BRCA2, and identified a duplication in RPA1, which might represent a new cancer predisposition gene. We conclude that optical genome mapping is a suitable tool for detecting potentially predisposing SVs in addition to WES in pediatric cancer patients.

A graph-based pan-genome of Brassica oleracea provides new insights into its domestication and morphotype diversification

The domestication of Brassica oleracea has resulted in diverse morphological types with distinct patterns of organ development. Here we report a graph-based pan-genome of B. oleracea constructed from high-quality genome assemblies of different morphotypes. The pan-genome harbors over 200 structural variant hotspot regions enriched in auxin- and flowering-related genes. Population genomic analyses revealed that early domestication of B. oleracea focused on leaf or stem development. Gene flows resulting from agricultural practices and variety improvement were detected among different morphotypes. Selective-sweep and pan-genome analyses identified an auxin-responsive small auxin up-regulated RNA gene and a CLAVATA3/ESR-RELATED family gene as crucial players in leaf-stem differentiation during the early stage of B. oleracea domestication and the BoKAN1 gene as instrumental in shaping the leafy heads of cabbage and Brussels sprouts. Our pan-genome and functional analyses further revealed that variations in the BoFLC2 gene play key roles in the divergence of vernalization and flowering characteristics among different morphotypes, and variations in the first intron of BoFLC3 are involved in fine-tuning the flowering process in cauliflower. This study provides a comprehensive understanding of the pan-genome of B. oleracea and sheds light on the domestication and differential organ development of this globally important crop species.

Structural genomic variation in the inbred Scandinavian wolf population contributes to the realized genetic load but is positively affected by immigration

When populations decrease in size and may become isolated, genomic erosion by loss of diversity from genetic drift and accumulation of deleterious mutations is likely an inevitable consequence. In such cases, immigration (genetic rescue) is necessary to restore levels of genetic diversity and counteract inbreeding depression. Recent work in conservation genomics has studied these processes focusing on the genetic diversity of single nucleotide polymorphisms. In contrast, our knowledge about structural genomic variation (insertions, deletions, duplications and inversions) in endangered species is limited. We analysed whole-genome, short-read sequences from 212 wolves from the inbred Scandinavian population and from neighbouring populations in Finland and Russia, and detected >35,000 structural variants (SVs) after stringent quality and genotype frequency filtering; >26,000 high-confidence variants remained after manual curation. The majority of variants were shorter than 1 kb, with a distinct peak in the length distribution of deletions at 190 bp, corresponding to insertion events of SINE/tRNA-Lys elements. The site frequency spectrum of SVs in protein-coding regions was significantly shifted towards rare alleles compared to putatively neutral variants, consistent with purifying selection. The realized genetic load of SVs in protein-coding regions increased with inbreeding levels in the Scandinavian population, but immigration provided a genetic rescue effect by lowering the load and reintroducing ancestral alleles at loci fixed for derived SVs. Our study shows that structural variation comprises a common type of in part deleterious mutations in endangered species and that establishing gene flow is necessary to mitigate the negative consequences of loss of diversity.

Reference-free Structural Variant Detection in Microbiomes via Long-read Coassembly Graphs

Bacterial genome dynamics are vital for understanding the mechanisms underlying microbial adaptation, growth, and their broader impact on host phenotype. Structural variants (SVs), genomic alterations of 10 base pairs or more, play a pivotal role in driving evolutionary processes and maintaining genomic heterogeneity within bacterial populations. While SV detection in isolate genomes is relatively straightforward, metagenomes present broader challenges due to absence of clear reference genomes and presence of mixed strains. In response, our proposed method rhea, forgoes reference genomes and metagenome-assembled genomes (MAGs) by encompassing a single metagenome coassembly graph constructed from all samples in a series. The log fold change in graph coverage between subsequent samples is then calculated to call SVs that are thriving or declining throughout the series. We show rhea to outperform existing methods for SV and horizontal gene transfer (HGT) detection in two simulated mock metagenomes, which is particularly noticeable as the simulated reads diverge from reference genomes and an increase in strain diversity is incorporated. We additionally demonstrate use cases for rhea on series metagenomic data of environmental and fermented food microbiomes to detect specific sequence alterations between subsequent time and temperature samples, suggesting host advantage. Our innovative approach leverages raw read patterns rather than references or MAGs to include all sequencing reads in analysis, and thus provide versatility in studying SVs across diverse and poorly characterized microbial communities for more comprehensive insights into microbial genome dynamics.

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.

Monitoring Genomic Structural Rearrangements Resulting from Gene Editing

The cytogenomics-based methodology of directional genomic hybridization (dGH) enables the detection and quantification of a more comprehensive spectrum of genomic structural variants than any other approach currently available, and importantly, does so on a single-cell basis. Thus, dGH is well-suited for testing and/or validating new advancements in CRISPR-Cas9 gene editing systems. In addition to aberrations detected by traditional cytogenetic approaches, the strand specificity of dGH facilitates detection of otherwise cryptic intra-chromosomal rearrangements, specifically small inversions. As such, dGH represents a powerful, high-resolution approach for the quantitative monitoring of potentially detrimental genomic structural rearrangements resulting from exposure to agents that induce DNA double-strand breaks (DSBs), including restriction endonucleases and ionizing radiations. For intentional genome editing strategies, it is critical that any undesired effects of DSBs induced either by the editing system itself or by mis-repair with other endogenous DSBs are recognized and minimized. In this paper, we discuss the application of dGH for assessing gene editing-associated structural variants and the potential heterogeneity of such rearrangements among cells within an edited population, highlighting its relevance to personalized medicine strategies.

A graph clustering algorithm for detection and genotyping of structural variants from long reads

BACKGROUND

Structural variants (SVs) are genomic polymorphisms defined by their length (>50 bp). The usual types of SVs are deletions, insertions, translocations, inversions, and copy number variants. SV detection and genotyping is fundamental given the role of SVs in phenomena such as phenotypic variation and evolutionary events. Thus, methods to identify SVs using long-read sequencing data have been recently developed.

FINDINGS

We present an accurate and efficient algorithm to predict germline SVs from long-read sequencing data. The algorithm starts collecting evidence (signatures) of SVs from read alignments. Then, signatures are clustered based on a Euclidean graph with coordinates calculated from lengths and genomic positions. Clustering is performed by the DBSCAN algorithm, which provides the advantage of delimiting clusters with high resolution. Clusters are transformed into SVs and a Bayesian model allows to precisely genotype SVs based on their supporting evidence. This algorithm is integrated into the single sample variants detector of the Next Generation Sequencing Experience Platform, which facilitates the integration with other functionalities for genomics analysis. We performed multiple benchmark experiments, including simulation and real data, representing different genome profiles, sequencing technologies (PacBio HiFi, ONT), and read depths.

CONCLUSION

The results show that our approach outperformed state-of-the-art tools on germline SV calling and genotyping, especially at low depths, and in error-prone repetitive regions. We believe this work significantly contributes to the development of bioinformatic strategies to maximize the use of long-read sequencing technologies.

Genomic consequences associated with Agrobacterium-mediated transformation of plants

Attenuated strains of the naturally occurring plant pathogen Agrobacterium tumefaciens can transfer virtually any DNA sequence of interest to model plants and crops. This has made Agrobacterium-mediated transformation (AMT) one of the most commonly used tools in agricultural biotechnology. Understanding AMT, and its functional consequences, is of fundamental importance given that it sits at the intersection of many fundamental fields of study, including plant-microbe interactions, DNA repair/genome stability, and epigenetic regulation of gene expression. Despite extensive research and use of AMT over the last 40 years, the extent of genomic disruption associated with integrating exogenous DNA into plant genomes using this method remains underappreciated. However, new technologies like long-read sequencing make this disruption more apparent, complementing previous findings from multiple research groups that have tackled this question in the past. In this review, we cover progress on the molecular mechanisms involved in Agrobacterium-mediated DNA integration into plant genomes. We also discuss localized mutations at the site of insertion and describe the structure of these DNA insertions, which can range from single copy insertions to large concatemers, consisting of complex DNA originating from different sources. Finally, we discuss the prevalence of large-scale genomic rearrangements associated with the integration of DNA during AMT with examples. Understanding the intended and unintended effects of AMT on genome stability is critical to all plant researchers who use this methodology to generate new genetic variants.

Lipid-Related Domestication Accounts for the Extreme Cold Sensitivity of Semiwild and Tropic Xishuangbanna Cucumber (Cucumis sativus L. var. xishuangbannanesis)

Xishuangbanna (XIS) cucumber (Cucumis sativus L. var. xishuangbannanesis) is a semiwild variety originating from low latitude tropic areas, and therefore shows extreme cold sensitivity and heat tolerance. Here, we mapped the quantitative trait loci (QTLs) that control the cold sensitivity and heat tolerance of XIS cucumber seedlings. Using bulked segregant analysis (BSA), we identified three QTLs (HTT1.1, HTT3.1, and HTT3.2, with a total length of 11.98 Mb) for heat tolerance and two QTLs (LTT6.1 and LTT6.2, with a total length of 8.74 Mb) for cold sensitivity. The QTL LTT6.1 was then narrowed down to a length of 641 kb by using kompetitive allele-specific PCR (KASP) markers. Based on structural variants (SVs) and single-nucleotide polymorphisms (SNPs), we found the LTT6.1 is covered by a high divergent region including a 50 kb deletion in the XIS49 genome, which affects the gene structure of lipase abhydrolase domain containing 6 (ABHD6, Csa_6G032560). Accordingly, there is a very big difference in lipid composition, but not in other osmoprotectants like free amino acids and fatty acids, between XIS49 and cultivated cucumber CL. Moreover, we calculated the composite likelihood ratio (CLR) and identified selective sweeps from 115 resequencing data, and found that lipid- and fatty-acid-related processes are major aspects in the domestication of the XIS group cucumber. LTT6.1 is a particularly special region positioned nearby lipid-related selective sweeps. These studies above suggested that the lipid-related domestication of XIS cucumbers should account for their extreme cold sensitivity.

Structural genomic variation and migratory behavior in a wild songbird

Structural variants (SVs) are a major source of genetic variation; and descriptions in natural populations and connections with phenotypic traits are beginning to accumulate in the literature. We integrated advances in genomic sequencing and animal tracking to begin filling this knowledge gap in the Eurasian blackcap. Specifically, we (a) characterized the genome-wide distribution, frequency, and overall fitness effects of SVs using haplotype-resolved assemblies for 79 birds, and (b) used these SVs to study the genetics of seasonal migration. We detected >15 K SVs. Many SVs overlapped repetitive regions and exhibited evidence of purifying selection suggesting they have overall deleterious effects on fitness. We used estimates of genomic differentiation to identify SVs exhibiting evidence of selection in blackcaps with different migratory strategies. Insertions and deletions dominated the SVs we identified and were associated with genes that are either directly (e.g., regulatory motifs that maintain circadian rhythms) or indirectly (e.g., through immune response) related to migration. We also broke migration down into individual traits (direction, distance, and timing) using existing tracking data and tested if genetic variation at the SVs we identified could account for phenotypic variation at these traits. This was only the case for 1 trait-direction-and 1 specific SV (a deletion on chromosome 27) accounted for much of this variation. Our results highlight the evolutionary importance of SVs in natural populations and provide insight into the genetic basis of seasonal migration.

Genomic Variations Between Fibrolamellar and Conventional Hepatocellular Carcinomas

Aim The aim of this study is to define genomic variations between fibrolamellar hepatocellular carcinoma (FL-HCC) and conventional hepatocellular carcinoma (HCC) Methods This study used the American Association for Cancer Research (AACR) Project GENIE data as a foundational element. Specifically, information about both fibrolamellar and conventional hepatocellular carcinoma was retrieved from this database. Results A total of 719 patients diagnosed with HCC and 52 individuals presenting with FL-HCC underwent thorough analysis. Notably, distinct variances in gene alterations were observed between the two cohorts. Predominantly, the HCC group exhibited frequent occurrences of mutations within the TP53 and CTNNB1 genes. Conversely, DNAJB1 fusion was uniquely identified in FL-HCC cases. Conclusion This study significantly broadens our understanding of the genetic makeup associated with FL-HCC and HCC. It is particularly notable because it reveals clear disparities in gene modifications between FL-HCC and HCC. Further investigation is essential to unravel the functional consequences of these genetic variances. This exploration will aid in the development of targeted therapeutic approaches to enhance the prognosis of patients diagnosed with diverse subtypes of HCC.

Enhancers dysfunction in the 3D genome of cancer cells

Eukaryotic genomes are spatially organized inside the cell nucleus, forming a threedimensional (3D) architecture that allows for spatial separation of nuclear processes and for controlled expression of genes required for cell identity specification and tissue homeostasis. Hence, it is of no surprise that mis-regulation of genome architecture through rearrangements of the linear genome sequence or epigenetic perturbations are often linked to aberrant gene expression programs in tumor cells. Increasing research efforts have shed light into the causes and consequences of alterations of 3D genome organization. In this review, we summarize the current knowledge on how 3D genome architecture is dysregulated in cancer, with a focus on enhancer highjacking events and their contribution to tumorigenesis. Studying the functional effects of genome architecture perturbations on gene expression in cancer offers a unique opportunity for a deeper understanding of tumor biology and sets the basis for the discovery of novel therapeutic targets.

Structural variants and short tandem repeats impact gene expression and splicing in bovine testis tissue

Structural variants (SVs) and short tandem repeats (STRs) are significant sources of genetic variation. However, the impacts of these variants on gene regulation have not been investigated in cattle. Here, we genotyped and characterized 19,408 SVs and 374,821 STRs in 183 bovine genomes and investigated their impact on molecular phenotypes derived from testis transcriptomes. We found that 71% STRs were multiallelic. The vast majority (95%) of STRs and SVs were in intergenic and intronic regions. Only 37% SVs and 40% STRs were in high linkage disequilibrium (LD) (R2 > 0.8) with surrounding SNPs/insertions and deletions (Indels), indicating that SNP-based association testing and genomic prediction are blind to a nonnegligible portion of genetic variation. We showed that both SVs and STRs were more than 2-fold enriched among expression and splicing QTL (e/sQTL) relative to SNPs/Indels and were often associated with differential expression and splicing of multiple genes. Deletions and duplications had larger impacts on splicing and expression than any other type of SV. Exonic duplications predominantly increased gene expression either through alternative splicing or other mechanisms, whereas expression- and splicing-associated STRs primarily resided in intronic regions and exhibited bimodal effects on the molecular phenotypes investigated. Most e/sQTL resided within 100 kb of the affected genes or splicing junctions. We pinpoint candidate causal STRs and SVs associated with the expression of SLC13A4 and TTC7B and alternative splicing of a lncRNA and CAPP1. We provide a catalog of STRs and SVs for taurine cattle and show that these variants contribute substantially to gene expression and splicing variation.

A genotype-phenotype correlation in split-hand/foot malformation type 1: further refinement of the phenotypic subregions within the 7q21.3 locus

Background: Split-hand/foot malformation type 1 (SHFM1) refers to the group of rare congenital limb disorders defined by the absence or hypoplasia of the central rays of the autopods with or without accompanying anomalies, such as hearing loss, craniofacial malformation, and ectodermal dysplasia. Consequently, the condition is characterized by clinical variability that hinders diagnostic and counseling procedures. SHFM1 is caused by pathogenic variants affecting the DLX5/6 genes and/or their tissue-specific enhancers at the 7q21.3 locus. Herein, we report on seven patients from five unrelated Polish families affected by variable symptoms of the SHFM1 spectrum, all harboring 7q21.3 or 7q21.2-q21.3 rearrangements, and provide a genotype-phenotype correlation in the studied cohort. Methods: We applied GTG banding, array-based comparative genomic hybridization (aCGH), and whole-genome sequencing (WGS) in order to identify the causative aberrations in all affected patients. Results: The identified pathogenic structural variants included deletions and/or translocations involving the 7q21.3 locus, i.e., t(7;10)(q21.3;q22.2) and t(7;12)(q21.3;q21.2) in all affected individuals. Interestingly, a sporadic carrier of the latter aberration presented the SHFM1 phenotype with additional features overlapping with Baker-Gordon syndrome (BAGOS), which resulted from the translocation breakpoint at chromosome 12 within the SYT1 gene. Conclusion: Clinical variability of the studied cohort reflects the composition of the DLX5/6 regulatory elements that were dislocated from their target genes by chromosomal rearrangements. The correlation of our data with the previously published observations enabled us to update the phenotypic subregions and regulatory units within the SHFM1 locus. In addition, we present the first case of SHFM1 and BAGOS-like phenotype that resulted from translocation breakpoints at chromosomes 7 and 12, both of which were pathogenic, and consequently, we show the first evidence that BAGOS can also result from the regulatory loss-of-function SYT1 mutations. In this paper, we emphasize the utility of sequence-based approaches in molecular diagnostics of disorders caused by regulatory structural variants.

VariantSurvival: a tool to identify genotype-treatment response

Motivation: For a number of neurological diseases, such as Alzheimer's disease, amyotrophic lateral sclerosis, and many others, certain genes are known to be involved in the disease mechanism. A common question is whether a structural variant in any such gene may be related to drug response in clinical trials and how this relationship can contribute to the lifecycle of drug development. Results: To this end, we introduce VariantSurvival, a tool that identifies changes in survival relative to structural variants within target genes. VariantSurvival matches annotated structural variants with genes that are clinically relevant to neurological diseases. A Cox regression model determines the change in survival between the placebo and clinical trial groups with respect to the number of structural variants in the drug target genes. We demonstrate the functionality of our approach with the exemplary case of the SETX gene. VariantSurvival has a user-friendly and lightweight graphical user interface built on the shiny web application package.

Comprehensive single-cell genome analysis at nucleotide resolution using the PTA Analysis Toolbox

Detection of somatic mutations in single cells has been severely hampered by technical limitations of whole-genome amplification. Novel technologies including primary template-directed amplification (PTA) significantly improved the accuracy of single-cell whole-genome sequencing (WGS) but still generate hundreds of artifacts per amplification reaction. We developed a comprehensive bioinformatic workflow, called the PTA Analysis Toolbox (PTATO), to accurately detect single base substitutions, insertions-deletions (indels), and structural variants in PTA-based WGS data. PTATO includes a machine learning approach and filtering based on recurrence to distinguish PTA artifacts from true mutations with high sensitivity (up to 90%), outperforming existing bioinformatic approaches. Using PTATO, we demonstrate that hematopoietic stem cells of patients with Fanconi anemia, which cannot be analyzed using regular WGS, have normal somatic single base substitution burdens but increased numbers of deletions. Our results show that PTATO enables studying somatic mutagenesis in the genomes of single cells with unprecedented sensitivity and accuracy.

Whole-Genome Sequencing Identified New Structural Variations in the DMD Gene That Cause Duchenne Muscular Dystrophy in Two Girls

Dystrophinopathies are the most common muscle diseases, especially in men. In women, on the other hand, a manifestation of Duchenne muscular dystrophy is rare due to X-chromosomal inheritance. We present two young girls with severe muscle weakness, muscular dystrophies, and creatine kinase (CK) levels exceeding 10,000 U/L. In the skeletal muscle tissues, dystrophin staining reaction showed mosaicism. The almost entirely skewed X-inactivation in both cases supported the possibility of a dystrophinopathy. Despite standard molecular diagnostics (including multiplex ligation-dependent probe amplification (MLPA) and next generation sequencing (NGS) gene panel sequencing), the genetic cause of the girls' conditions remained unknown. However, whole-genome sequencing revealed two reciprocal translocations between their X chromosomes and chromosome 5 and chromosome 19, respectively. In both cases, the breakpoints on the X chromosomes were located directly within the DMD gene (in introns 54 and 7, respectively) and were responsible for the patients' phenotypes. Additional techniques such as Sanger sequencing, conventional karyotyping and fluorescence in situ hybridization (FISH) confirmed the disruption of DMD gene in both patients through translocations. These findings underscore the importance of accurate clinical data combined with histopathological analysis in pinpointing the suspected underlying genetic disorder. Moreover, our study illustrates the viability of whole-genome sequencing as a time-saving and highly effective method for identifying genetic factors responsible for complex genetic constellations in Duchenne muscular dystrophy (DMD).

Long-Read DNA Sequencing: Recent Advances and Remaining Challenges

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

Copy Number Variation on ABCC2-DNMBP Loci Affects the Diversity and Composition of the Fecal Microbiota in Pigs

Genetic variation in the pig genome partially modulates the composition of porcine gut microbial communities. Previous studies have been focused on the association between single nucleotide polymorphisms (SNPs) and the gut microbiota, but little is known about the relationship between structural variants and fecal microbial traits. The main goal of this study was to explore the association between porcine genome copy number variants (CNVs) and the diversity and composition of pig fecal microbiota. For this purpose, we used whole-genome sequencing data to undertake a comprehensive identification of CNVs followed by a genome-wide association analysis between the estimated CNV status and the fecal bacterial diversity in a commercial Duroc pig population. A CNV predicted as gain (DUP) partially harboring ABCC2-DNMBP loci was associated with richness (P = 5.41 × 10-5, false discovery rate [FDR] = 0.022) and Shannon α-diversity (P = 1.42 × 10-4, FDR = 0.057). The in silico predicted gain of copies was validated by real-time quantitative PCR (qPCR), and its segregation, and positive association with the richness and Shannon α-diversity of the porcine fecal bacterial ecosystem was confirmed in an unrelated F1 (Duroc × Iberian) cross. Our results advise the relevance of considering the role of host-genome structural variants as potential modulators of microbial ecosystems and suggest the ABCC2-DNMBP CNV as a host-genetic factor for the modulation of the diversity and composition of the fecal microbiota in pigs. IMPORTANCE A better understanding of the environmental and host factors modulating gut microbiomes is a topic of greatest interest. Recent evidence suggests that genetic variation in the pig genome partially controls the composition of porcine gut microbiota. However, since previous studies have been focused on the association between single nucleotide polymorphisms and the fecal microbiota, little is known about the relationship between other sources of genetic variation, like the structural variants and microbial traits. Here, we identified, experimentally validated, and replicated in an independent population a positive link between the gain of copies of ABCC2-DNMBP loci and the diversity and composition of pig fecal microbiota. Our results advise the relevance of considering the role of host-genome structural variants as putative modulators of microbial ecosystems and open the possibility of implementing novel holobiont-based management strategies in breeding programs for the simultaneous improvement of microbial traits and host performance.

Schizophrenia-associated somatic copy-number variants from 12,834 cases reveal recurrent NRXN1 and ABCB11 disruptions

While germline copy-number variants (CNVs) contribute to schizophrenia (SCZ) risk, the contribution of somatic CNVs (sCNVs)-present in some but not all cells-remains unknown. We identified sCNVs using blood-derived genotype arrays from 12,834 SCZ cases and 11,648 controls, filtering sCNVs at loci recurrently mutated in clonal blood disorders. Likely early-developmental sCNVs were more common in cases (0.91%) than controls (0.51%, p = 2.68e-4), with recurrent somatic deletions of exons 1-5 of the NRXN1 gene in five SCZ cases. Hi-C maps revealed ectopic, allele-specific loops forming between a potential cryptic promoter and non-coding cis-regulatory elements upon 5' deletions in NRXN1. We also observed recurrent intragenic deletions of ABCB11, encoding a transporter implicated in anti-psychotic response, in five treatment-resistant SCZ cases and showed that ABCB11 is specifically enriched in neurons forming mesocortical and mesolimbic dopaminergic projections. Our results indicate potential roles of sCNVs in SCZ risk.

Genome sequencing and comprehensive rare-variant analysis of 465 families with neurodevelopmental disorders

Despite significant progress in unraveling the genetic causes of neurodevelopmental disorders (NDDs), a substantial proportion of individuals with NDDs remain without a genetic diagnosis after microarray and/or exome sequencing. Here, we aimed to assess the power of short-read genome sequencing (GS), complemented with long-read GS, to identify causal variants in participants with NDD from the National Institute for Health and Care Research (NIHR) BioResource project. Short-read GS was conducted on 692 individuals (489 affected and 203 unaffected relatives) from 465 families. Additionally, long-read GS was performed on five affected individuals who had structural variants (SVs) in technically challenging regions, had complex SVs, or required distal variant phasing. Causal variants were identified in 36% of affected individuals (177/489), and a further 23% (112/489) had a variant of uncertain significance after multiple rounds of re-analysis. Among all reported variants, 88% (333/380) were coding nuclear SNVs or insertions and deletions (indels), and the remainder were SVs, non-coding variants, and mitochondrial variants. Furthermore, long-read GS facilitated the resolution of challenging SVs and invalidated variants of difficult interpretation from short-read GS. This study demonstrates the value of short-read GS, complemented with long-read GS, in investigating the genetic causes of NDDs. GS provides a comprehensive and unbiased method of identifying all types of variants throughout the nuclear and mitochondrial genomes in individuals with NDD.

Recent Advances in Pathology: the 2023 Annual Review Issue of The Journal of Pathology

The 2023 Annual Review Issue of The Journal of Pathology, Recent Advances in Pathology, contains 12 invited reviews on topics of current interest in pathology. This year, our subjects include immuno-oncology and computational pathology approaches for diagnostic and research applications in human disease. Reviews on the tissue microenvironment include the effects of apoptotic cell-derived exosomes, how understanding the tumour microenvironment predicts prognosis, and the growing appreciation of the diverse functions of fibroblast subtypes in health and disease. We also include up-to-date reviews of modern aspects of the molecular basis of malignancies, and our final review covers new knowledge of vascular and lymphatic regeneration in cardiac disease. All of the reviews contained in this issue are written by expert groups of authors selected to discuss the recent progress in their particular fields and all articles are freely available online (https://pathsocjournals.onlinelibrary.wiley.com/journal/10969896). © 2023 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

The spectrum and significance of secondary (co-occurring) genetic alterations in sarcomas: the hallmarks of sarcomagenesis

Bone and soft tissue tumors are generally classified into complex karyotype sarcomas versus those with recurrent genetic alterations, often in the form of gene fusions. In this review, we provide an overview of important co-occurring genomic alterations, organized by biological mechanisms and covering a spectrum of genomic alteration types: mutations (single-nucleotide variations or indels) in oncogenes or tumor suppressor genes, copy number alterations, transcriptomic signatures, genomic complexity indices (e.g. CINSARC), and complex genomic structural variants. We discuss the biological and prognostic roles of these so-called secondary or co-occurring alterations, arguing that recognition and detection of these alterations may be significant for our understanding and management of mesenchymal tumors. On a related note, we also discuss major recurrent alterations in so-called complex karyotype sarcomas. These secondary alterations are essential to sarcomagenesis via a variety of mechanisms, such as inactivation of tumor suppressors, activation of proliferative signal transduction, telomere maintenance, and aberrant regulation of epigenomic/chromatin remodeling players. The use of comprehensive genomic profiling, including targeted next-generation sequencing panels or whole-exome sequencing, may be incorporated into clinical workflows to offer more comprehensive, potentially clinically actionable information. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Diallel panel reveals a significant impact of low-frequency genetic variants on gene expression variation in yeast

Unraveling the genetic sources of gene expression variation is essential to better understand the origins of phenotypic diversity in natural populations. Genome-wide association studies identified thousands of variants involved in gene expression variation, however, variants detected only explain part of the heritability. In fact, variants such as low-frequency and structural variants (SVs) are poorly captured in association studies. To assess the impact of these variants on gene expression variation, we explored a half-diallel panel composed of 323 hybrids originated from pairwise crosses of 26 natural Saccharomyces cerevisiae isolates. Using short- and long-read sequencing strategies, we established an exhaustive catalog of single nucleotide polymorphisms (SNPs) and SVs for this panel. Combining this dataset with the transcriptomes of all hybrids, we comprehensively mapped SNPs and SVs associated with gene expression variation. While SVs impact gene expression variation, SNPs exhibit a higher effect size with an overrepresentation of low-frequency variants compared to common ones. These results reinforce the importance of dissecting the heritability of complex traits with a comprehensive catalog of genetic variants at the population level.

Nanopore long-read next-generation sequencing for detection of mitochondrial DNA large-scale deletions

Primary mitochondrial diseases are progressive genetic disorders affecting multiple organs and characterized by mitochondrial dysfunction. These disorders can be caused by mutations in nuclear genes coding proteins with mitochondrial localization or by genetic defects in the mitochondrial genome (mtDNA). The latter include point pathogenic variants and large-scale deletions/rearrangements. MtDNA molecules with the wild type or a variant sequence can exist together in a single cell, a condition known as mtDNA heteroplasmy. MtDNA single point mutations are typically detected by means of Next-Generation Sequencing (NGS) based on short reads which, however, are limited for the identification of structural mtDNA alterations. Recently, new NGS technologies based on long reads have been released, allowing to obtain sequences of several kilobases in length; this approach is suitable for detection of structural alterations affecting the mitochondrial genome. In the present work we illustrate the optimization of two sequencing protocols based on long-read Oxford Nanopore Technology to detect mtDNA structural alterations. This approach presents strong advantages in the analysis of mtDNA compared to both short-read NGS and traditional techniques, potentially becoming the method of choice for genetic studies on mtDNA.

Genome-Wide Discovery of Structural Variants Reveals Distinct Variant Dynamics for Two Closely Related Monilinia Species

Structural variants (SVs) are variants with sizes bigger than 50 base pairs and capable of changing the size, copy number, location, orientation, and sequence content of genomic DNA. Although these variants have been proven to be extensive and involved in many evolutionary processes along the tree of life, there is still insufficient information on many fungal plant pathogens. In this study, the extent of SVs, as well as single nucleotide polymorphisms (SNPs), have been determined for two prominent species of the Monilinia genus (the causal agents of brown rot disease in pome and stone fruits): M. fructicola and M. laxa for the first time. The genomes of M. fructicola were found to be more variant-rich in contrast to M. laxa based on the reference-based variant calling (with a total number of 266.618 and 190.599 SNPs; 1540 and 918 SVs, respectively). The extent, as well as distribution of SVs, presented high conservation within the species and high diversity between the species. Investigation of potential functional effects of characterized variants revealed high potential relevance of SVs. Moreover, the detailed characterization of copy number variations (CNVs) for each isolate revealed that around 0.67% of M. fructicola genomes and 2.06% of M. laxa genomes are CN variables. The variant catalog as well as distinct variant dynamics within and between the species presented in this study opens doors for many further research questions.

Structural genomic variants in thoracic aortic disease

PURPOSE OF REVIEW

Structural genomic variants have emerged as a relevant cause for several disorders, including intellectual disability, neuropsychiatric disorders, cancer and congenital heart disease. In this review, we will discuss the current knowledge about the involvement of structural genomic variants and, in particular, copy number variants in the development of thoracic aortic and aortic valve disease.

RECENT FINDINGS

There is a growing interest in the identification of structural variants in aortopathy. Copy number variants identified in thoracic aortic aneurysms and dissections, bicuspid aortic valve related aortopathy, Williams-Beuren syndrome and Turner syndrome are discussed in detail. Most recently, the first inversion disrupting FBN1 has been reported as a cause for Marfan syndrome.

SUMMARY

During the past 15 years, the knowledge on the role of copy number variants as a cause for aortopathy has grown significantly, which is partially due to the development of novel technologies including next-generation sequencing. Although copy number variants are now often investigated on a routine basis in diagnostic laboratories, more complex structural variants such as inversions, which require the use of whole genome sequencing, are still relatively new to the field of thoracic aortic and aortic valve disease.

The EN-TEx resource of multi-tissue personal epigenomes & variant-impact models

Understanding how genetic variants impact molecular phenotypes is a key goal of functional genomics, currently hindered by reliance on a single haploid reference genome. Here, we present the EN-TEx resource of 1,635 open-access datasets from four donors (∼30 tissues × ∼15 assays). The datasets are mapped to matched, diploid genomes with long-read phasing and structural variants, instantiating a catalog of >1 million allele-specific loci. These loci exhibit coordinated activity along haplotypes and are less conserved than corresponding, non-allele-specific ones. Surprisingly, a deep-learning transformer model can predict the allele-specific activity based only on local nucleotide-sequence context, highlighting the importance of transcription-factor-binding motifs particularly sensitive to variants. Furthermore, combining EN-TEx with existing genome annotations reveals strong associations between allele-specific and GWAS loci. It also enables models for transferring known eQTLs to difficult-to-profile tissues (e.g., from skin to heart). Overall, EN-TEx provides rich data and generalizable models for more accurate personal functional genomics.

Usefulness and Limitations of Multiple Ligation-Dependent Probe Amplification in Antithrombin Deficiency

Multiplex ligation-dependent probe amplification (MLPA) identifies genetic structural variants in SERPINC1 in 5% of cases with antithrombin deficiency (ATD), the most severe congenital thrombophilia. Our aim was to unravel the utility and limitations of MLPA in a large cohort of unrelated patients with ATD (N = 341). MLPA identified 22 structural variants (SVs) causing ATD (6.5%). MLPA did not detect SVs affecting introns (four cases), and the diagnosis was inaccurate in two cases according to long-range PCR or nanopore sequencing. MLPA was used to detect possible hidden SVs in 61 cases with type I deficiency with single nucleotide variations (SNVs) or small insertion/deletion (INDEL). One case had a false deletion of exon 7, as the 29-bp deletion affected an MLPA probe. We evaluated 32 variants affecting MLPA probes: 27 SNVs and 5 small INDELs. In three cases, MLPA gave false-positive results, all diagnosed as deletions of the affected exon: a small INDEL complex, and two SNVs affecting MLPA probes. Our study confirms the utility of MLPA to detect SVs in ATD, but also shows some limitations in detecting intronic SVs. MLPA renders imprecise and false-positive results for genetic defects which affect MLPA probes. Our results encourage the validation of MLPA results.

Genome Size Changes by Duplication, Divergence, and Insertion in Caenorhabditis Worms

Genome size has been measurable since the 1940s but we still do not understand genome size variation. Caenorhabditis nematodes show strong conservation of chromosome number but vary in genome size between closely related species. Androdioecy, where populations are composed of males and self-fertile hermaphrodites, evolved from outcrossing, female-male dioecy, three times in this group. In Caenorhabditis, androdioecious genomes are 10-30% smaller than dioecious species, but in the nematode Pristionchus, androdioecy evolved six times and does not correlate with genome size. Previous hypotheses include genome size evolution through: 1) Deletions and "genome shrinkage" in androdioecious species; 2) Transposable element (TE) expansion and DNA loss through large deletions (the "accordion model"); and 3) Differing TE dynamics in androdioecious and dioecious species. We analyzed nematode genomes and found no evidence for these hypotheses. Instead, nematode genome sizes had strong phylogenetic inertia with increases in a few dioecious species, contradicting the "genome shrinkage" hypothesis. TEs did not explain genome size variation with the exception of the DNA transposon Mutator which was twice as abundant in dioecious genomes. Across short and long evolutionary distances Caenorhabditis genomes evolved through small structural mutations including gene-associated duplications and insertions. Seventy-one protein families had significant, parallel decreases across androdioecious Caenorhabditis including genes involved in the sensory system, regulatory proteins and membrane-associated immune responses. Our results suggest that within a dynamic landscape of frequent small rearrangements in Caenorhabditis, reproductive mode mediates genome evolution by altering the precise fates of individual genes, proteins, and the phenotypes they underlie.

Genomic evidence reveals intraspecific divergence of the hot-spring snake (Thermophis baileyi), an endangered reptile endemic to the Qinghai-Tibet plateau

Understanding how and why species evolve requires knowledge on intraspecific divergence. In this study, we examined intraspecific divergence in the endangered hot-spring snake (Thermophis baileyi), an endemic species on the Qinghai-Tibet Plateau (QTP). Whole-genome resequencing of 58 sampled individuals from 15 populations was performed to identify the drivers of intraspecific divergence and explore the potential roles of genes under selection. Our analyses resolved three groups, with major intergroup admixture occurring in regions of group contact. Divergence probably occurred during the Pleistocene as a result of glacial climatic oscillations, Yadong-Gulu rift, and geothermal fields differentiation, while complex gene flow between group pairs reflected a unique intraspecific divergence pattern on the QTP. Intergroup fixed loci involved selected genes functionally related to divergence and local adaptation, especially adaptation to hot spring microenvironments in different geothermal fields. Analysis of structural variants, genetic diversity, inbreeding, and genetic load indicated that the hot-spring snake population has declined to a low level with decreased diversity, which is important for the conservation management of this endangered species. Our study demonstrated that the integration of demographic history, gene flow, genomic divergence genes, and other information is necessary to distinguish the evolutionary processes involved in speciation.

Genome assembly, structural variants, and genetic differentiation between lake whitefish young species pairs (Coregonus sp.) with long and short reads

Nascent pairs of ecologically differentiated species offer an opportunity to get a better glimpse at the genetic architecture of speciation. Of particular interest is our recent ability to consider a wider range of genomic variants, not only single-nucleotide polymorphisms (SNPs), thanks to long-read sequencing technology. We can now identify structural variants (SVs) such as insertions, deletions and other rearrangements, allowing further insights into the genetic architecture of speciation and how different types of variants are involved in species differentiation. Here, we investigated genomic patterns of differentiation between sympatric species pairs (Dwarf and Normal) belonging to the lake whitefish (Coregonus clupeaformis) species complex. We assembled the first reference genomes for both C. clupeaformis sp. Normal and C. clupeaformis sp. Dwarf, annotated the transposable elements and analysed the genomes in the light of related coregonid species. Next, we used a combination of long- and short-read sequencing to characterize SVs and genotype them at the population scale using genome-graph approaches, showing that SVs cover five times more of the genome than SNPs. We then integrated both SNPs and SVs to investigate the genetic architecture of species differentiation in two different lakes and highlighted an excess of shared outliers of differentiation. In particular, a large fraction of SVs differentiating the two species correspond to insertions or deletions of transposable elements (TEs), suggesting that TE accumulation may represent a key component of genetic divergence between the Dwarf and Normal species. Together, our results suggest that SVs may play an important role in speciation and that, by combining second- and third-generation sequencing, we now have the ability to integrate SVs into speciation genomics.

The role of structural variants in pest adaptation and genome evolution of the Colorado potato beetle, Leptinotarsa decemlineata (Say)

Structural variation has been associated with genetic diversity and adaptation. Despite these observations, it is not clear what their relative importance is for evolution, especially in rapidly adapting species. Here, we examine the significance of structural polymorphisms in pesticide resistance evolution of the agricultural super-pest, the Colorado potato beetle, Leptinotarsa decemlineata. By employing a parent offspring trio sequencing procedure, we develop highly contiguous reference genomes to characterize structural variation. These updated assemblies represent >100-fold improvement of contiguity and include derived pest and ancestral nonpest individuals. We identify >200,000 structural variations, which appear to be nonrandomly distributed across the genome as they co-occur with transposable elements and genes. Structural variations intersect with exons in a large proportion of gene annotations (~20%) that are associated with insecticide resistance (including cytochrome P450s), development, and transcription. To understand the role structural variations play in adaptation, we measure their allele frequencies among an additional 57 individuals using whole genome resequencing data, which represents pest and nonpest populations of North America. Incorporating multiple independent tests to detect the signature of natural selection using SNP data, we identify 14 genes that are probably under positive selection, include structural variations, and SNPs of elevated frequency within the pest lineages. Among these, three are associated with insecticide resistance based on previous research. One of these genes, CYP4g15, is coinduced during insecticide exposure with glycosyltransferase-13, which is a duplicated gene enclosed within a structural variant adjacent to the CYP4g15 genic region. These results demonstrate the significance of structural variations as a genomic feature to describe species history, genetic diversity, and adaptation.

Genomic variation across Chinook salmon populations reveals effects of a duplication on migration alleles and supports fine scale structure

The distribution of ecotypic variation in natural populations is influenced by neutral and adaptive evolutionary forces that are challenging to disentangle. This study provides a high-resolution portrait of genomic variation in Chinook salmon (Oncorhynchus tshawytscha) with emphasis on a region of major effect for ecotypic variation in migration timing. With a filtered data set of ~13 million single nucleotide polymorphisms (SNPs) from low-coverage whole genome resequencing of 53 populations (3566 barcoded individuals), we contrasted patterns of genomic structure within and among major lineages and examined the extent of a selective sweep at a major effect region underlying migration timing (GREB1L/ROCK1). Neutral variation provided support for fine-scale structure of populations, while allele frequency variation in GREB1L/ROCK1 was highly correlated with mean return timing for early and late migrating populations within each of the lineages (r²  = .58-.95; p < .001). However, the extent of selection within the genomic region controlling migration timing was much narrower in one lineage (interior stream-type) compared to the other two major lineages, which corresponded to the breadth of phenotypic variation in migration timing observed among lineages. Evidence of a duplicated block within GREB1L/ROCK1 may be responsible for reduced recombination in this portion of the genome and contributes to phenotypic variation within and across lineages. Lastly, SNP positions across GREB1L/ROCK1 were assessed for their utility in discriminating migration timing among lineages, and we recommend multiple markers nearest the duplication to provide highest accuracy in conservation applications such as those that aim to protect early migrating Chinook salmon. These results highlight the need to investigate variation throughout the genome and the effects of structural variants on ecologically relevant phenotypic variation in natural species.

Assessment of copy number variants in three Brazilian locally adapted cattle breeds using whole-genome re-sequencing data

Further characterization of genetic structural variations should strongly focus on small and endangered local breeds given their role in unraveling genes and structural variants underlying selective pressures and phenotype variation. A comprehensive genome-wide assessment of copy number variations (CNVs) based on whole-genome re-sequencing data was performed on three Brazilian locally adapted cattle breeds (Caracu Caldeano, Crioulo Lageano, and Pantaneiro) using the ARS-UCD1.2 genome assembly. Data from 36 individuals with an average coverage depth of 14.07× per individual was used. A total of 24 945 CNVs were identified distributed among the breeds (Caracu Caldeano = 7285, Crioulo Lageano = 7297, and Pantaneiro = 10 363). Deletion events were 1.75-2.07-fold higher than duplications, and the total length of CNVs is composed mostly of a high number of segments between 10 and 30 kb. CNV regions (CNVRs) are not uniformly scattered throughout the genomes (n = 463), and 105 CNVRs were found overlapping among the studied breeds. Functional annotation of the CNVRs revealed variants with high consequence on protein sequence harboring relevant genes, in which we highlighted the BOLA-DQB, BOLA-DQA5, CD1A, β-defensins, PRG3, and ULBP21 genes. Enrichment analysis based on the gene list retrieved from the CNVRs disclosed over-represented terms (p < 0.01) strongly associated with immunity and cattle resilience to harsh environments. Additionally, QTL associated with body conformation and dairy-related traits were also unveiled within the CNVRs. These results provide better understanding of the selective forces shaping the genome of such cattle breeds and identify traces of natural selection pressures by which these populations have been exposed to challenging environmental conditions.

The Power of Clinical Diagnosis for Deciphering Complex Genetic Mechanisms in Rare Diseases

Complex genetic disease mechanisms, such as structural or non-coding variants, currently pose a substantial difficulty in frontline diagnostic tests. They thus may account for most unsolved rare disease patients regardless of the clinical phenotype. However, the clinical diagnosis can narrow the genetic focus to just a couple of genes for patients with well-established syndromes defined by prominent physical and/or unique biochemical phenotypes, allowing deeper analyses to consider complex genetic origin. Then, clinical-diagnosis-driven genome sequencing strategies may expedite the development of testing and analytical methods to account for complex disease mechanisms as well as to advance functional assays for the confirmation of complex variants, clinical management, and the development of new therapies.

Third-generation genome sequencing implicates medium-sized structural variants in chronic schizophrenia

Background

Schizophrenia (SCZ) is a heterogeneous psychiatric disorder, with significant contribution from genetic factors particularly for chronic cases with negative symptoms and cognitive deficits. To date, Genome Wide Association Studies (GWAS) and exome sequencing have associated SCZ with a number of single nucleotide polymorphisms (SNPs) and copy number variants (CNVs), but there is still missing heritability. Medium-sized structural variants (SVs) are difficult to detect using SNP arrays or second generation sequencing, and may account for part of the missing heritability of SCZ.

Aims and objectives

To identify SVs associated with severe chronic SCZ across the whole genome.

Study design

10 multiplex families with probands suffering from chronic SCZ with negative symptoms and cognitive deficits were recruited, with all their affected members demonstrating uni-lineal inheritance. Control subjects comprised one affected member from the affected lineage, and unaffected members from each paternal and maternal lineage.

Methods

Third generation sequencing was applied to peripheral blood samples from 10 probands and 5 unaffected controls. Bioinformatic tools were used to identify SVs from the long sequencing reads, with confirmation of findings in probands by short-read Illumina sequencing, Sanger sequencing and visual manual validation with Integrated Genome Browser.

Results

In the 10 probands, we identified and validated 88 SVs (mostly in introns and medium-sized), within 79 genes, which were absent in the 5 unaffected control subjects. These 79 genes were enriched in 20 biological pathways which were related to brain development, neuronal migration, neurogenesis, neuronal/synaptic function, learning/memory, and hearing. These identified SVs also showed evidence for enrichment of genes that are highly expressed in the adolescent striatum.

Conclusion

A substantial part of the missing heritability in SCZ may be explained by medium-sized SVs detectable only by third generation sequencing. We have identified a number of such SVs potentially conferring risk for SCZ, which implicate multiple brain-related genes and pathways. In addition to previously-identified pathways involved in SCZ such as neurodevelopment and neuronal/synaptic functioning, we also found novel evidence for enrichment in hearing-related pathways and genes expressed in the adolescent striatum.

Evolutionary genomics: Insights from the invasive European starlings

Two fundamental questions for evolutionary studies are the speed at which evolution occurs, and the way that this evolution may present itself within an organism's genome. Evolutionary studies on invasive populations are poised to tackle some of these pressing questions, including understanding the mechanisms behind rapid adaptation, and how it facilitates population persistence within a novel environment. Investigation of these questions are assisted through recent developments in experimental, sequencing, and analytical protocols; in particular, the growing accessibility of next generation sequencing has enabled a broader range of taxa to be characterised. In this perspective, we discuss recent genetic findings within the invasive European starlings in Australia, and outline some critical next steps within this research system. Further, we use discoveries within this study system to guide discussion of pressing future research directions more generally within the fields of population and evolutionary genetics, including the use of historic specimens, phenotypic data, non-SNP genetic variants (e.g., structural variants), and pan-genomes. In particular, we emphasise the need for exploratory genomics studies across a range of invasive taxa so we can begin understanding broad mechanisms that underpin rapid adaptation in these systems. Understanding how genetic diversity arises and is maintained in a population, and how this contributes to adaptability, requires a deep understanding of how evolution functions at the molecular level, and is of fundamental importance for the future studies and preservation of biodiversity across the globe.

Cas9-Mediated Nanopore Sequencing Enables Precise Characterization of Structural Variants in CCM Genes

Deletions in the CCM1, CCM2, and CCM3 genes are a common cause of familial cerebral cavernous malformations (CCMs). In current molecular genetic laboratories, targeted next-generation sequencing or multiplex ligation-dependent probe amplification are mostly used to identify copy number variants (CNVs). However, both techniques are limited in their ability to specify the breakpoints of CNVs and identify complex structural variants (SVs). To overcome these constraints, we established a targeted Cas9-mediated nanopore sequencing approach for CNV detection with single nucleotide resolution. Using a MinION device, we achieved complete coverage for the CCM genes and determined the exact size of CNVs in positive controls. Long-read sequencing for a CCM1 and CCM2 CNV revealed that the adjacent ANKIB1 and NACAD genes were also partially or completely deleted. In addition, an interchromosomal insertion and an inversion in CCM2 were reliably re-identified by long-read sequencing. The refinement of CNV breakpoints by long-read sequencing enabled fast and inexpensive PCR-based variant confirmation, which is highly desirable to reduce costs in subsequent family analyses. In conclusion, Cas9-mediated nanopore sequencing is a cost-effective and flexible tool for molecular genetic diagnostics which can be easily adapted to various target regions.

The genomic basis of reproductive and migratory behaviour in a polymorphic salmonid

Recent ecotypic differentiation provides unique opportunities to investigate the genomic basis and architecture of local adaptation, while offering insights into how species form and persist. Sockeye salmon (Oncorhynchus nerka) exhibit migratory and resident ("kokanee") ecotypes, which are further distinguished into shore-spawning and stream-spawning reproductive ecotypes. Here, we analysed 36 sockeye (stream-spawning) and kokanee (stream- and shore-spawning) genomes from a system where they co-occur and have recent common ancestry (Okanagan Lake/River in British Columbia, Canada) to investigate the genomic basis of reproductive and migratory behaviour. Examination of the genomic landscape of differentiation, differences in allele frequencies and genotype-phenotype associations revealed three main blocks of sequence differentiation on chromosomes 7, 12 and 20, associated with migratory behaviour, spawning location and spawning timing. Structural variants identified in these same areas suggest they could contribute to ecotypic differentiation directly as causal variants or via maintenance of their genomic architecture through recombination suppression mechanisms. Genes in these regions were related to spatial memory and swimming endurance (SYNGAP, TPM3), as well as eye and brain development (including SIX6), potentially associated with differences in migratory behaviour and visual habitats across spawning locations, respectively. Additional genes (GREB1L, ROCK1) identified here have been associated with timing of migration in other salmonids and could explain variation in timing of O. nerka spawning. Together, these results based on the joint analysis of sequence and structural variation represent a significant advance in our understanding of the genomic landscape of ecotypic differentiation at different stages in the speciation continuum.