Mobile Element Insertions Research Articles

Detection of trait-associated structural variations using short-read sequencing

Genomic structural variation (SV) affects genetic and phenotypic characteristics in diverse organisms, but the lack of reliable methods to detect SV has hindered genetic analysis. We developed a computational algorithm (MOPline) that includes missing call recovery combined with high-confidence SV call selection and genotyping using short-read whole-genome sequencing (WGS) data. Using 3,672 high-coverage WGS datasets, MOPline stably detected ∼16,000 SVs per individual, which is over ∼1.7-3.3-fold higher than previous large-scale projects while exhibiting a comparable level of statistical quality metrics. We imputed SVs from 181,622 Japanese individuals for 42 diseases and 60 quantitative traits. A genome-wide association study with the imputed SVs revealed 41 top-ranked or nearly top-ranked genome-wide significant SVs, including 8 exonic SVs with 5 novel associations and enriched mobile element insertions. This study demonstrates that short-read WGS data can be used to identify rare and common SVs associated with a variety of traits.

Diagnostic rate of germline pathogenic variants in pancreatic ductal adenocarcinoma patients using whole genome sequencing.

Identification of germline pathogenic variants in cancer patients is critical for treatment planning, genetic counseling, and health policymaking. However, previous estimates of the prevalence of germline etiology of pancreatic ductal adenocarcinoma (PDAC) were biased because they were based only on sequencing data of protein-coding regions of known PDAC candidate genes. To determine the percentage of patients with PDAC carrying germline pathogenic variants, we enrolled the inpatients from the digestive health clinics, hematology and oncology clinics, and surgical clinics of a single tertiary medical center in Taiwan for whole genome sequencing (WGS) analysis of genomic DNA. The virtual gene panel of 750 genes comprised PDAC candidate genes and those listed in the COSMIC Cancer Gene Census. The genetic variant types under investigation included single nucleotide substitutions, small indels, structural variants, and mobile element insertions (MEIs). In 8 of 24 (33.3%) patients with PDAC, we identified pathogenic/likely pathogenic variants, including single nucleotide substitutions and small indels in ATM, BRCA1, BRCA2, POLQ, SPINK1 and CASP8, as well as structural variants in CDC25C and USP44. We identified additional patients carrying variants that could potentially affect splicing. This cohort study demonstrates that an extensive analysis of the abundant information yielded by the WGS approach can uncover many pathogenic variants that could be missed by traditional panel-based or whole exome sequencing-based approaches. The percentage of patients with PDAC carrying germline variants might be much higher than previously expected.

Genome graphs detect human polymorphisms in active epigenomic state during influenza infection

Genetic variants, including mobile element insertions (MEIs), are known to impact the epigenome. We hypothesized that genome graphs, which encapsulate genetic diversity, could reveal missing epigenomic signals. To test this, we sequenced the epigenome of monocyte-derived macrophages from 35 ancestrally diverse individuals before and after influenza infection, allowing us to investigate the role of MEIs in immunity. We characterized genetic variants and MEIs using linked reads and built a genome graph. Mapping epigenetic data revealed 2.3%-3% novel peaks for H3K4me1, H3K27ac chromatin immunoprecipitation sequencing (ChIP-seq), and ATAC-seq. Additionally, the use of a genome graph modified some quantitative trait loci estimates and revealed 375 polymorphic MEIs in an active epigenomic state. Among these is an AluYh3 polymorphism whose chromatin state changed after infection and was associated with the expression of TRIM25, a gene that restricts influenza RNA synthesis. Our results demonstrate that graph genomes can reveal regulatory regions that would have been overlooked by other approaches.

A Novel Alu Element Insertion in ATM Induces Exon Skipping in Suspected HBOC Patients

The vast majority of patients at risk of hereditary breast and/or ovarian cancer (HBOC) syndrome remain without a molecular diagnosis after routine genetic testing. One type of genomic alteration that is commonly missed by diagnostic pipelines is mobile element insertions (MEIs). Here, we reanalyzed multigene panel data from suspected HBOC patients using the MEI detection tool Mobster. A novel Alu element insertion in ATM intron 54 (ATM:c.8010+30_8010+31insAluYa5) was identified as a potential contributing factor in seven patients. Transcript analysis of patient-derived RNA from three heterozygous carriers revealed exon 54 skipping in 38% of total ATM transcripts. To manifest the direct association between the Alu element insertion and the aberrant splice pattern, HEK293T and MCF7 cells were transfected with wild-type or Alu element-carrying minigene constructs. On average, 77% of plasmid-derived transcripts lacked exon 54 in the presence of the Alu element insertion compared to only 4.7% of transcripts expressed by the wild-type minigene. These results strongly suggest ATM:c.8010+30_8010+31insAluYa5 as the main driver of ATM exon 54 skipping. Since this exon loss is predicted to cause a frameshift and a premature stop codon, mutant transcripts are unlikely to translate into functional proteins. Based on its estimated frequency of up to 0.05% in control populations, we propose to consider ATM:c.8010+30_8010+31insAluYa5 in suspected HBOC patients and to clarify its role in carcinogenesis through future epidemiological and functional analyses. Generally, the implementation of MEI detection tools in diagnostic sequencing pipelines could increase the diagnostic yield, as MEIs are likely underestimated contributors to genetic diseases.

Variant spectrum of von Hippel-Lindau disease and its genomic heterogeneity in Japan.

Von Hippel-Lindau (VHL) disease is an autosomal dominant, inherited syndrome with variants in the VHL gene, causing predisposition to multi-organ neoplasms with vessel abnormality. Germline variants in VHL can be detected in 80-90% of patients clinically diagnosed with VHL disease. Here, we summarize the results of genetic tests for 206 Japanese VHL families, and elucidate the molecular mechanisms of VHL disease, especially in variant-negative unsolved cases. Of the 206 families, genetic diagnosis was positive in 175 families (85%), including 134 families (65%) diagnosed by exon sequencing (15 novel variants) and 41 (20%) diagnosed by multiplex ligation-dependent probe amplification (MLPA) (one novel variant). The deleterious variants were significantly enriched in VHL disease Type 1. Interestingly, five synonymous or non-synonymous variants within exon 2 caused exon 2 skipping, which is the first report of exon 2 skipping caused by several missense variants. Whole genome and target deep sequencing analysis were performed for 22 unsolved cases with no variant identified and found three cases with VHL mosaicism (variant allele frequency: 2.5-22%), one with mobile element insertion in the VHL promoter region, and two with a pathogenic variant of BAP1 or SDHB. The variants associated with VHL disease are heterogeneous, and for more accuracy of the genetic diagnosis of VHL disease, comprehensive genome and DNA/RNA analyses are required to detect VHL mosaicism, complicated structure variants and other related gene variants.

Diversity and prevalence of type VI secretion system effectors in clinical Pseudomonas aeruginosa isolates.

Pseudomonas aeruginosa is an opportunistic pathogen and a major driver of morbidity and mortality in people with Cystic Fibrosis (CF). The Type VI secretion system (T6SS) is a molecular nanomachine that translocates effectors across the bacterial membrane into target cells or the extracellular environment enabling intermicrobial interaction. P. aeruginosa encodes three T6SS clusters, the H1-, H2- and H3-T6SS, and numerous orphan islands. Genetic diversity of T6SS-associated effectors in P. aeruginosa has been noted in reference strains but has yet to be explored in clinical isolates. Here, we perform a comprehensive bioinformatic analysis of the pangenome and T6SS effector genes in 52 high-quality clinical P. aeruginosa genomes isolated from CF patients and housed in the Personalised Approach to P. aeruginosa strain repository. We confirm that the clinical CF isolate pangenome is open and principally made up of accessory and unique genes that may provide strain-specific advantages. We observed genetic variability in some effector/immunity encoding genes and show that several well-characterised vgrG and PAAR islands are absent from numerous isolates. Our analysis shows clear evidence of disruption to T6SS genomic loci through transposon, prophage, and mobile genetic element insertions. We identified an orphan vgrG island in P. aeruginosa strain PAK and five clinical isolates using in silico analysis which we denote vgrG7, predicting a gene within this cluster to encode a Tle2 lipase family effector. Close comparison of T6SS loci in clinical isolates compared to reference P. aeruginosa strain PAO1 revealed the presence of genes encoding eight new T6SS effectors with the following putative functions: cytidine deaminase, lipase, metallopeptidase, NADase, and pyocin. Finally, the prevalence of characterised and putative T6SS effectors were assessed in 532 publicly available P. aeruginosa genomes, which suggests the existence of accessory effectors. Our in silico study of the P. aeruginosa T6SS exposes a level of genetic diversity at T6SS genomic loci not seen to date within P. aeruginosa, particularly in CF isolates. As understanding the effector repertoire is key to identifying the targets of T6SSs and its efficacy, this comprehensive analysis provides a path for future experimental characterisation of these mediators of intermicrobial competition and host manipulation.

The Diagnostic Yield and Difficulties of Utilizing Soft-clipped Read Clusters Encountered in Clinical Exome Sequencing.

Despite the wide use of next generation sequencing, there are still many difficulties in detecting structural variants. A split read is one of the clues of structural variants and is represented as a soft-clipped read in the raw sequencing data. Considering that most of the breakpoints of structural variants reside in non-coding regions, split read information has not been routinely used in exome sequencing or targeted panel sequencing. Recently, SCRAMble, a software capable of detecting mobile element insertion (MEI) and deletion based on soft-clipped read clusters (SCRCs), was shown to provide an additional diagnostic yield of 0.03 - 0.25%. SCRAMble is the only software that can be used for exome sequencing or targeted panel sequencing to detect structural variants based on SCRC information. The aim of present study was to establish a working procedure of utilizing SCRC information using SCRAMble in clinical exome sequencing and to assess its diagnostic yield. Raw sequencing data of clinical exome sequencing were retrospectively analyzed using SCRAMble to search MEIs and deletions. SCRAMble software was installed according to the manufacturer's instructions and default parameters except for one, mei-score, which was adjusted for sensitivity, were used. RefSeq gene annotation was performed for both MEI and deletion calls using ANNOVAR. Blacklist-based filtering was used to reduce candidate MEI/deletion calls. Clinical relevance was manually evaluated for the remaining variant calls. One diagnostic MEI, which is a founder variant in East Asia, was detected in two cases (2/266, 0.75%). In addition, two diagnostic deletions, which had been previously detected in depth-of-coverage (DOC)-based copy number variant (CNV) callers, were detected (2/266, 0.75%). Base-level breakpoints that could not be derived by the DOC-based callers were identified for these two deletions using SCRAMble. Most SCRCs were repetitive among cases and blacklist-based filtering reduced candidate MEI/deletion calls by 49.5 - 94.5%, leaving a considerable number of variant calls to be manually validated. SCRC screening in exome or targeted panel sequencing may provide additional diagnostic yield either by pathogenic MEI detection or reassurance of deletions identified by DOC-based CNV callers. Development of an efficient filtering algorithm is warranted.

Whole F9 gene sequencing identified deep intronic variations in genetically unresolved hemophilia B patients

BackgroundThe disease-causative variant remains unidentified in approximately 0.5% to 2% of hemophilia B patients using conventional genetic investigations, and F9 deep intronic variations could be responsible for these phenotypes. ObjectivesThis study aimed to characterize deep intronic variants in hemophilia B patients for whom genetic investigations failed. MethodsWe performed whole F9 sequencing in 17 genetically unsolved hemophilia B patients. The pathogenic impact of the candidate variants identified was studied using both in silico analysis (MaxEntScan and spliceAI) and minigene assay. ResultsIn total, 9 candidate variants were identified in 15 patients; 7 were deep intronic substitutions and 2 corresponded to insertions of mobile elements. The most frequent variants found were c.278-1806A>C and the association of c.278-1244A>G and c.392-864T>G, identified in 4 and 6 unrelated individuals, respectively. In silico analysis predicted splicing impact for 4 substitutions (c.278-1806A>C, c.392-864T>G, c.724-2385G>T, c.723+4297T>A). Minigene assay showed a deleterious splicing impact for these 4 substitutions and also for the c.278-1786_278-1785insLINE. In the end, 5 variants were classified as likely pathogenic using the American College of Medical Genetics and Genomics guidelines, and 4 as of unknown significance. Thus, the hemophilia B-causing variant was identified in 13/17 (76%) families. ConclusionWe elucidated the causing defect in three-quarters of the families included in this study, and we reported new F9 deep intronic variants that can cause hemophilia B.

A second look at exome sequencing data: detecting mobile elements insertion in a rare disease cohort.

About 0.3% of all variants are due to de novo mobile element insertions (MEIs). The massive development of next-generation sequencing has made it possible to identify MEIs on a large scale. We analyzed exome sequencing (ES) data from 3232 individuals (2410 probands) with developmental and/or neurological abnormalities, with MELT, a tool designed to identify MEIs. The results were filtered by frequency, impacted region and gene function. Following phenotype comparison, two candidates were identified in two unrelated probands. The first mobile element (ME) was found in a patient referred for poikilodermia. A homozygous insertion was identified in the FERMT1 gene involved in Kindler syndrome. RNA study confirmed its pathological impact on splicing. The second ME was a de novo Alu insertion in the GRIN2B gene involved in intellectual disability, and detected in a patient with a developmental disorder. The frequency of de novo exonic MEIs in our study is concordant with previous studies on ES data. This project, which aimed to identify pathological MEIs in the coding sequence of genes, confirms that including detection of MEs in the ES pipeline can increase the diagnostic rate. This work provides additional evidence that ES could be used alone as a diagnostic exam.

Detection and Characterization of a De Novo Alu Retrotransposition Event Causing NKX2-1-Related Disorder.

Heterozygous NKX2-1 loss-of-function variants cause combinations of hyperkinetic movement disorders (MDs, particularly childhood-onset chorea), pulmonary dysfunction, and hypothyroidism. Mobile element insertions (MEIs) are potential disease-causing structural variants whose detection in routine diagnostics remains challenging. To establish the molecular diagnosis of two first-degree relatives with clinically suspected NKX2-1-related disorder who had negative NKX2-1 Sanger (SS), whole-exome (WES), and whole-genome (WGS) sequencing. The proband's WES was analyzed for MEIs. A candidate MEI in NKX2-1 underwent optimized SS after plasmid cloning. Functional studies exploring NKX2-1 haploinsufficiency at RNA and protein levels were performed. A 347-bp AluYa5 insertion with a 65-bp poly-A tail followed by a 16-bp duplication of the pre-insertion wild-type sequence in exon 3 of NKX2-1 (ENST00000354822.7:c.556_557insAlu541_556dup) segregated with the disease phenotype. We identified a de novo exonic AluYa5 insertion causing NKX2-1-related disorder in SS/WES/WGS-negative cases, suggesting that MEI analysis of short-read sequencing data or targeted long-read sequencing could unmask the molecular diagnosis of unsolved MD cases. © 2022 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Comparative genomics reveals new insights into the evolution of the IncA and IncC family of plasmids.

Incompatibility groups IncA and IncC plasmids are of great concern due to their ability to disseminate antibiotic resistance in bacteria via conjugative transfer. A deep understanding of their genomic structures and evolutionary characteristics is of great significance for improving our knowledge about its multidrug-resistance evolution and dissemination. However, current knowledge of their backbone structure, features of core functional modules and the characteristics of variable regions is based on a few plasmids, which highlights the need for a comprehensive systematic study. The present study thoroughly compared and analysed 678 IncA and IncC plasmid genomes. We found that their core functional genes were occasionally deficient and sometimes existed as multiple functional copies/multiple families, which resulted in much diversity. The phylogeny of 13 core functional genes corresponded well to the plasmid subtypes. The conjugative transfer system gained diverse complexity and exhibited many previously unnoticed types with multiple combinations. The insertion of mobile genetic elements (MGEs) in plasmids varied between types and was present in 4 insertion spots in different types of plasmids with certain types of transposons, integrons and insertion sequences. The impact of gene duplication, deletion, the insertion of MGEs, genome rearrangement and recombination resulted in the complex dynamic variable backbone of IncA and IncC plasmids. And IncA and IncC plasmids were more complex than their closest relative SXT/R391 integrative conjugative elements (ICEs), which included nearly all of the diversity of SXT/R391 in key systems. Our work demonstrated a global and systematic view of the IncA and IncC plasmids and provides many new insights into their genome evolution.

Owl Monkey Alu Insertion Polymorphisms and Aotus Phylogenetics.

Owl monkeys (genus Aotus), or "night monkeys" are platyrrhine primates in the Aotidae family. Early taxonomy only recognized one species, Aotus trivirgatus, until 1983, when Hershkovitz proposed nine unique species designations, classified into red-necked and gray-necked species groups based predominately on pelage coloration. Recent studies questioned this conventional separation of the genus and proposed designations based on the geographical location of wild populations. Alu retrotransposons are a class of mobile element insertion (MEI) widely used to study primate phylogenetics. A scaffold-level genome assembly for one Aotus species, Aotus nancymaae [Anan_2.0], facilitated large-scale ascertainment of nearly 2000 young lineage-specific Alu insertions. This study provides candidate oligonucleotides for locus-specific PCR assays for over 1350 of these elements. For 314 Alu elements across four taxa with multiple specimens, PCR analyses identified 159 insertion polymorphisms, including 21 grouping A. nancymaae and Aotus azarae (red-necked species) as sister taxa, with Aotus vociferans and A. trivirgatus (gray-necked) being more basal. DNA sequencing identified five novel Alu elements from three different taxa. The Alu datasets reported in this study will assist in species identification and provide a valuable resource for Aotus phylogenetics, population genetics and conservation strategies when applied to wild populations.

Genomics and epigenetics guided identification of tissue-specific genomic safe harbors

BackgroundGenomic safe harbors are regions of the genome that can maintain transgene expression without disrupting the function of host cells. Genomic safe harbors play an increasingly important role in improving the efficiency and safety of genome engineering. However, limited safe harbors have been identified.ResultsHere, we develop a framework to facilitate searches for genomic safe harbors by integrating information from polymorphic mobile element insertions that naturally occur in human populations, epigenomic signatures, and 3D chromatin organization. By applying our framework to polymorphic mobile element insertions identified in the 1000 Genomes project and the Genotype-Tissue Expression (GTEx) project, we identify 19 candidate safe harbors in blood cells and 5 in brain cells. For three candidate sites in blood, we demonstrate the stable expression of transgene without disrupting nearby genes in host erythroid cells. We also develop a computer program, Genomics and Epigenetic Guided Safe Harbor mapper (GEG-SH mapper), for knowledge-based tissue-specific genomic safe harbor selection.ConclusionsOur study provides a new knowledge-based framework to identify tissue-specific genomic safe harbors. In combination with the fast-growing genome engineering technologies, our approach has the potential to improve the overall safety and efficiency of gene and cell-based therapy in the near future.

Allele-biased expression of the bovine APOB gene associated with the cholesterol deficiency defect suggests cis-regulatory enhancer effects of the LTR retrotransposon insertion

The insertion of an endogenous retroviral long terminal repeat (LTR) sequence into the bovine apolipoprotein B (APOB) gene is causal to the inherited genetic defect cholesterol deficiency (CD) observed in neonatal and young calves. Affected calves suffer from developmental abnormalities, symptoms of incurable diarrhoea and often die within weeks to a few months after birth. Neither the detailed effects of the LTR insertion on APOB expression profile nor the specific mode of inheritance nor detailed phenotypic consequences of the mutation are undisputed. In our study, we analysed German Holstein dairy heifers at the peak of hepatic metabolic load and exposed to an additional pathogen challenge for clinical, metabolic and hepatic transcriptome differences between wild type (CDF) and heterozygote carriers of the mutation (CDC). Our data revealed that a divergent allele-biased expression pattern of the APOB gene in heterozygous CDC animals leads to a tenfold higher expression of exons upstream and a decreased expression of exons downstream of the LTR insertion compared to expression levels of CDF animals. This expression pattern could be a result of enhancer activity induced by the LTR insertion, in addition to a previously reported artificial polyadenylation signal. Thus, our data support a regulatory potential of mobile element insertions. With regard to the phenotype generated by the LTR insertion, heterozygote CDC carriers display significantly differential hepatic expression of genes involved in cholesterol biosynthesis and lipid metabolism. Phenotypically, CDC carriers show a significantly affected lipomobilization compared to wild type animals. These results reject a completely recessive mode of inheritance for the CD defect, which should be considered for selection decisions in the affected population. Exemplarily, our results illustrate the regulatory impact of mobile element insertions not only on specific host target gene expression but also on global transcriptome profiles with subsequent biological, functional and phenotypic consequences in a natural in-vivo model of a non-model mammalian organism.

Accumulation of endosymbiont genomes in an insect autosome followed by endosymbiont replacement

Eukaryotic genomes can acquire bacterial DNA via lateral gene transfer (LGT).1 A prominent source of LGT is Wolbachia,2 a widespread endosymbiont of arthropods and nematodes that is transmitted maternally through female germline cells.3,4 The DNA transfer from the Wolbachia endosymbiont wAna to Drosophila ananassae is extensive5-7 and has been localized to chromosome 4, contributing to chromosome expansion in this lineage.6 As has happened frequently with claims of bacteria-to-eukaryote LGT, the contribution of wAna transfers to the expanded size of D.ananassae chromosome 4 has been specifically contested8 owing to an assembly where Wolbachia sequences were classified as contaminants and removed.9 Here, long-read sequencing with DNA from a Wolbachia-cured line enabled assembly of 4.9 Mbp of nuclear Wolbachia transfers (nuwts) in D.ananassae and a 24-kbp nuclear mitochondrial transfer. The nuwts are <8,000 years old in at least two locations in chromosome 4 with at least one whole-genome integration followed by rapid extensive duplication of most of the genome with regions that have up to 10 copies. The genes in nuwts are accumulating small indels and mobile element insertions. Among the highly duplicated genes are cifA and cifB, two genes associated with Wolbachia-mediated Drosophila cytoplasmic incompatibility. The wAna strain that was the source of nuwts was subsequently replaced by a different wAna endosymbiont. Direct RNA Nanopore sequencing of Wolbachia-cured lines identified nuwt transcripts, including spliced transcripts, but functionality, if any, remains elusive.

Refinements and considerations for trio whole-genome sequence analysis when investigating Mendelian diseases presenting in early childhood

SummaryTo facilitate early deployment of whole-genome sequencing (WGS) for severely ill children, a standardized pipeline for WGS analysis with timely turnaround and primary care pediatric uptake is needed. We developed a bioinformatics pipeline for comprehensive gene-agnostic trio WGS analysis of children suspected of having an undiagnosed monogenic disease that included detection and interpretation of primary genetic mechanisms of disease, including SNVs/indels, CNVs/SVs, uniparental disomy (UPD), imprinted genes, short tandem repeat expansions, mobile element insertions, SMN1/2 copy number calling, and mitochondrial genome variants. We assessed primary care practitioner experience and competence in a large cohort of 521 families (comprising 90% WGS trios). Children were identified by primary practitioners for recruitment, and we used the UK index of multiple deprivation to confirm lack of patient socio-economic status ascertainment bias. Of the 521 children sequenced, 176 (34%) received molecular diagnoses, with rates as high as 45% for neurology clinics. Twenty-three of the diagnosed cases (13%) required bespoke methods beyond routine SNV/CNV analysis. In our multidisciplinary clinician user experience assessment, both pediatricians and clinical geneticists expressed strong support for rapid WGS early in the care pathway, but requested further training in determining patient selection, consenting, and variant interpretation. Rapid trio WGS provides an efficacious single-pass screening test for children when deployed by primary practitioners in clinical settings that carry high a priori risk for rare pediatric disease presentations.

Rearrangement processes and structural variations show evidence of selection in oesophageal adenocarcinomas

Oesophageal adenocarcinoma (OAC) provides an ideal case study to characterize large-scale rearrangements. Using whole genome short-read sequencing of 383 cases, for which 214 had matched whole transcriptomes, we observed structural variations (SV) with a predominance of deletions, tandem duplications and inter-chromosome junctions that could be identified as LINE-1 mobile element (ME) insertions. Complex clusters of rearrangements resembling breakage-fusion-bridge cycles or extrachromosomal circular DNA accounted for 22% of complex SVs affecting known oncogenes. Counting SV events affecting known driver genes substantially increased the recurrence rates of these drivers. After excluding fragile sites, we identified 51 candidate new drivers in genomic regions disrupted by SVs, including ETV5, KAT6B and CLTC. RUNX1 was the most recurrently altered gene (24%), with many deletions inactivating the RUNT domain but preserved the reading frame, suggesting an altered protein product. These findings underscore the importance of identification of SV events in OAC with implications for targeted therapies.

Recently Integrated Alu Elements in Capuchin Monkeys: A Resource for Cebus/Sapajus Genomics.

Capuchins are platyrrhines (monkeys found in the Americas) within the Cebidae family. For most of their taxonomic history, the two main morphological types of capuchins, gracile (untufted) and robust (tufted), were assigned to a single genus, Cebus. Further, all tufted capuchins were assigned to a single species, Cebus apella, despite broad geographic ranges spanning Central and northern South America. In 2012, tufted capuchins were assigned to their genus, Sapajus, with eight currently recognized species and five Cebus species, although these numbers are still under debate. Alu retrotransposons are a class of mobile element insertion (MEI) widely used to study primate phylogenetics. However, Alu elements have rarely been used to study capuchins. Recent genome-level assemblies for capuchins (Cebus imitator; [Cebus_imitator_1.0] and Sapajus apella [GSC_monkey_1.0]) facilitated large scale ascertainment of young lineage-specific Alu insertions. Reported here are 1607 capuchin specific and 678 Sapajus specific Alu insertions along with candidate oligonucleotides for locus-specific PCR assays for many elements. PCR analyses identified 104 genus level and 51 species level Alu insertion polymorphisms. The Alu datasets reported in this study provide a valuable resource that will assist in the classification of archival samples lacking phenotypic data and for the study of capuchin phylogenetic relationships.

Detection of a pathogenic Alu element insertion in PALB2 gene from targeted NGS diagnostic data.

Despite routine analysis of a large panel of genes, pathogenic variants are only detected in approximately 20% of families with hereditary breast and/or ovarian cancer. Mobile element insertions (MEI) are known to cause genetic diseases in humans, but remain challenging to detect. Retrospective analysis of targeted next-generation sequencing (NGS) data from 359 patients was performed using a dedicated MEI detection pipeline. We detected one MEI in exon 9 of the PALB2 gene in a woman with a family history of breast cancer. The pathogenic variant, c.2872_2888delins114AluL2, disrupts the PALB2 coding sequence and leads to the production of a truncated protein, p.(Gln958Valfs*38). This is the first report of a pathogenic MEI in PALB2. This study illustrates that MEI analysis may help to improve molecular diagnostic yield and can be performed from targeted NGS data used for routine diagnosis.

Whole-genome sequencing of 1,171 elderly admixed individuals from Brazil

As whole-genome sequencing (WGS) becomes the gold standard tool for studying population genomics and medical applications, data on diverse non-European and admixed individuals are still scarce. Here, we present a high-coverage WGS dataset of 1,171 highly admixed elderly Brazilians from a census-based cohort, providing over 76 million variants, of which ~2 million are absent from large public databases. WGS enables identification of ~2,000 previously undescribed mobile element insertions without previous description, nearly 5 Mb of genomic segments absent from the human genome reference, and over 140 alleles from HLA genes absent from public resources. We reclassify and curate pathogenicity assertions for nearly four hundred variants in genes associated with dominantly-inherited Mendelian disorders and calculate the incidence for selected recessive disorders, demonstrating the clinical usefulness of the present study. Finally, we observe that whole-genome and HLA imputation could be significantly improved compared to available datasets since rare variation represents the largest proportion of input from WGS. These results demonstrate that even smaller sample sizes of underrepresented populations bring relevant data for genomic studies, especially when exploring analyses allowed only by WGS.