SAMBLASTER: fast duplicate marking and structural variant read extraction

Theoretical characterisation of strand cross-correlation in ChIP-seq

BACKGROUND

Strand cross-correlation profiles are used for both peak calling pre-analysis and quality control (QC) in chromatin immunoprecipitation followed by sequencing (ChIP-seq) analysis. Despite its potential for robust and accurate assessments of signal-to-noise ratio (S/N) because of its peak calling independence, it remains unclear what aspects of quality such strand cross-correlation profiles actually measure.

RESULTS

We introduced a simple model to simulate the mapped read-density of ChIP-seq and then derived the theoretical maximum and minimum of cross-correlation coefficients between strands. The results suggest that the maximum coefficient of typical ChIP-seq samples is directly proportional to the number of total mapped reads and the square of the ratio of signal reads, and inversely proportional to the number of peaks and the length of read-enriched regions. Simulation analysis supported our results and evaluation using 790 ChIP-seq data obtained from the public database demonstrated high consistency between calculated cross-correlation coefficients and estimated coefficients based on the theoretical relations and peak calling results. In addition, we found that the mappability-bias-correction improved sensitivity, enabling differentiation of maximum coefficients from the noise level. Based on these insights, we proposed virtual S/N (VSN), a novel peak call-free metric for S/N assessment. We also developed PyMaSC, a tool to calculate strand cross-correlation and VSN efficiently. VSN achieved most consistent S/N estimation for various ChIP targets and sequencing read depths. Furthermore, we demonstrated that a combination of VSN and pre-existing peak calling results enable the estimation of the numbers of detectable peaks for posterior experiments and assess peak calling results.

CONCLUSIONS

We present the first theoretical insights into the strand cross-correlation, and the results reveal the potential and the limitations of strand cross-correlation analysis. Our quality assessment framework using VSN provides peak call-independent QC and will help in the evaluation of peak call analysis in ChIP-seq experiments.

Molecular basis of a new ovine model for human 3M syndrome-2

Background

Brachygnathia, cardiomegaly and renal hypoplasia syndrome (BCRHS, OMIA 001595–9940) is a previously reported recessively inherited disorder in Australian Poll Merino/Merino sheep. Affected lambs are stillborn with various congenital defects as reflected in the name of the disease, as well as short stature, a short and broad cranium, a small thoracic cavity, thin ribs and brachysternum. The BCRHS phenotype shows similarity to certain human short stature syndromes, in particular the human 3M syndrome-2. Here we report the identification of a likely disease-causing variant and propose an ovine model for human 3M syndrome-2.

Results

Eight positional candidate genes were identified among the 39 genes in the approximately 1 Mb interval to which the disease was mapped previously. Obscurin like cytoskeletal adaptor 1 (OBSL1) was selected as a strong positional candidate gene based on gene function and the resulting phenotypes observed in humans with mutations in this gene. Whole genome sequencing of an affected lamb (BCRHS3) identified a likely causal variant ENSOARG00000020239:g.220472248delC within OBSL1. Sanger sequencing of seven affected, six obligate carrier, two phenotypically unaffected animals from the original flock and one unrelated control animal validated the variant. A genotyping assay was developed to genotype 583 animals from the original flock, giving an estimated allele frequency of 5%.

Conclusions

The identification of a likely disease-causing variant resulting in a frameshift (p.(Val573Trpfs*119)) in the OBSL1 protein has enabled improved breeding management of the implicated flock. The opportunity for an ovine model for human 3M syndrome and ensuing therapeutic research is promising given the availability of carrier ram semen for BCRHS.

Cancer-specific CTCF binding facilitates oncogenic transcriptional dysregulation

BACKGROUND

The three-dimensional genome organization is critical for gene regulation and can malfunction in diseases like cancer. As a key regulator of genome organization, CCCTC-binding factor (CTCF) has been characterized as a DNA-binding protein with important functions in maintaining the topological structure of chromatin and inducing DNA looping. Among the prolific binding sites in the genome, several events with altered CTCF occupancy have been reported as associated with effects in physiology or disease. However, hitherto there is no comprehensive survey of genome-wide CTCF binding patterns across different human cancers.

RESULTS

To dissect functions of CTCF binding, we systematically analyze over 700 CTCF ChIP-seq profiles across human tissues and cancers and identify cancer-specific CTCF binding patterns in six cancer types. We show that cancer-specific lost and gained CTCF binding events are associated with altered chromatin interactions, partially with DNA methylation changes, and rarely with sequence mutations. While lost bindings primarily occur near gene promoters, most gained CTCF binding events exhibit enhancer activities and are induced by oncogenic transcription factors. We validate these findings in T cell acute lymphoblastic leukemia cell lines and patient samples and show that oncogenic NOTCH1 induces specific CTCF binding and they cooperatively activate expression of target genes, indicating transcriptional condensation phenomena.

CONCLUSIONS

Specific CTCF binding events occur in human cancers. Cancer-specific CTCF binding can be induced by other transcription factors to regulate oncogenic gene expression. Our results substantiate CTCF binding alteration as a functional epigenomic signature of cancer.

Two Novel Pathogenic Variants Confirm RMND1 Causative Role in Perrault Syndrome with Renal Involvement

RMND1 (required for meiotic nuclear division 1 homolog) pathogenic variants are known to cause combined oxidative phosphorylation deficiency (COXPD11), a severe multisystem disorder. In one patient, a homozygous RMND1 pathogenic variant, with an established role in COXPD11, was associated with a Perrault-like syndrome. We performed a thorough clinical investigation and applied a targeted multigene hearing loss panel to reveal the cause of hearing loss, ovarian dysfunction (two cardinal features of Perrault syndrome) and chronic kidney disease in two adult female siblings. Two compound heterozygous missense variants, c.583G>A (p.Gly195Arg) and c.818A>C (p.Tyr273Ser), not previously associated with disease, were identified in RMND1 in both patients, and their segregation with disease was confirmed in family members. The patients have no neurological or intellectual impairment, and nephrological evaluation predicts a benign course of kidney disease. Our study presents the mildest, so far reported, RMND1-related phenotype and delivers the first independent confirmation that RMND1 is causally involved in the development of Perrault syndrome with renal involvement. This highlights the importance of including RMND1 to the list of Perrault syndrome causative factors and provides new insight into the clinical manifestation of RMND1 deficiency.

Parallel Genetic Origin of Foot Feathering in Birds

Understanding the genetic basis of similar phenotypes shared between lineages is a long-lasting research interest. Even though animal evolution offers many examples of parallelism, for many phenotypes little is known about the underlying genes and mutations. We here use a combination of whole-genome sequencing, expression analyses, and comparative genomics to study the parallel genetic origin of ptilopody (Pti) in chicken. Ptilopody (or foot feathering) is a polygenic trait that can be observed in domesticated and wild avian species and is characterized by the partial or complete development of feathers on the ankle and feet. In domesticated birds, ptilopody is easily selected to fixation, though extensive variation in the type and level of feather development is often observed. By means of a genome-wide association analysis, we identified two genomic regions associated with ptilopody. At one of the loci, we identified a 17-kb deletion affecting PITX1 expression, a gene known to encode a transcription regulator of hindlimb identity and development. Similarly to pigeon, at the second loci, we observed ectopic expression of TBX5, a gene involved in forelimb identity and a key determinant of foot feather development. We also observed that the trait evolved only once as foot-feathered birds share the same haplotype upstream TBX5. Our findings indicate that in chicken and pigeon ptilopody is determined by the same set of genes that affect similar molecular pathways. Our study confirms that ptilopody has evolved through parallel evolution in chicken and pigeon.

A Spontaneous Aggressive ERα+ Mammary Tumor Model Is Driven by Kras Activation

The NRL-PRL murine model, defined by mammary-selective transgenic rat prolactin ligand rPrl expression, establishes spontaneous ER+ mammary tumors in nulliparous females, mimicking the association between elevated prolactin (PRL) and risk for development of ER+ breast cancer in postmenopausal women. Whole-genome and exome sequencing in a discovery cohort (n = 5) of end-stage tumors revealed canonical activating mutations and copy number amplifications of Kras. The frequent mutations in this pathway were validated in an extension cohort, identifying activating Ras alterations in 79% of tumors (23 of 29). Transcriptome analyses over the course of oncogenesis revealed marked alterations associated with Ras activity in established tumors compared with preneoplastic tissues; in cell-intrinsic processes associated with mitosis, cell adhesion, and invasion; as well as in the surrounding tumor environment. These genomic analyses suggest that PRL induces a selective bottleneck for spontaneous Ras-driven tumors that may model a subset of aggressive clinical ER+ breast cancers.

Extensive Homoplasy but No Evidence of Convergent Evolution of Repeat Numbers at MIRU Loci in Modern Mycobacterium tuberculosis Lineages

More human deaths have been attributable to Mycobacterium tuberculosis than any other pathogen, and the epidemic is sustained by ongoing transmission. Various typing schemes have been developed to identify strain-specific differences and track transmission dynamics in affected communities, with recent introduction of whole genome sequencing providing the most accurate assessment. Mycobacterial interspersed repetitive unit (MIRU) typing is a family of variable number tandem repeat schemes that have been widely used to study the molecular epidemiology of M. tuberculosis. MIRU typing was used in most well-resourced settings to perform routine molecular epidemiology. Instances of MIRU homoplasy have been observed in comparison with sequence-based phylogenies, limiting its discriminatory value. A fundamental question is whether the observed homoplasy arises purely through stochastic processes, or whether there is evidence of natural selection. We compared repeat numbers at 24 MIRU loci with a whole genome sequence-based phylogeny of 245 isolates representing three modern M. tuberculosis lineages. This analysis demonstrated extensive homoplasy of repeat numbers, but did not detect any evidence of natural selection of repeat numbers, at least since the ancestral branching of the three modern lineages of M. tuberculosis. In addition, we observed good sensitivity but poor specificity and positive predictive values of MIRU-24 to detect clusters of recent transmission, as defined by whole-genome single nucleotide polymorphism analysis. These findings provide mechanistic insight, and support a transition away from VNTR-based typing toward sequence-based typing schemes for both research and public health purposes.

Multiplexed Quantification of Four Neuroblastoma DNA Targets in a Single Droplet Digital PCR Reaction

The detection and characterization of cell-free DNA (cfDNA) in peripheral blood from neuroblastoma patients may serve as a minimally invasive approach to liquid biopsy. Major challenges in the analysis of cfDNA purified from blood samples are small sample volumes and low cfDNA concentrations. Droplet digital PCR (ddPCR) is a technology suitable for analyzing low levels of cfDNA. Reported here are two quadruplexed ddPCR assay protocols that reliably quantify MYCN and ALK copy numbers in a single reaction together with the two reference genes, NAGK and AFF3, and accurately estimate ALK^F1174L (exon 23 position 3522, C>A) and ALK^R1275Q (exon 25 position 3824, G>A) mutant allele fractions using cfDNA as input. The separation of positive and negative droplets was optimized for detecting two targets in each ddPCR fluorescence channel by the adjustment of the probe and primer concentrations of each target molecule. The quadruplexed assays were validated using a panel of 10 neuroblastoma cell lines and paired blood plasma and primary neuroblastoma samples from nine patients. Accuracy and sensitivity thresholds in quadruplexed assays corresponded well with those from the respective duplexed assays. Presented are two robust quadruplexed ddPCR protocols applicable in the routine clinical setting and that require only minimal plasma volumes for the assessment of MYCN and ALK oncogene status.

Sequencing identifies multiple early introductions of SARS-CoV-2 to the New York City Region

Effective public response to a pandemic relies upon accurate measurement of the extent and dynamics of an outbreak. Viral genome sequencing has emerged as a powerful approach to link seemingly unrelated cases, and large-scale sequencing surveillance can inform on critical epidemiological parameters. Here, we report the analysis of 864 SARS-CoV-2 sequences from cases in the New York City metropolitan area during the COVID-19 outbreak in Spring 2020. The majority of cases had no recent travel history or known exposure, and genetically linked cases were spread throughout the region. Comparison to global viral sequences showed that early transmission was most linked to cases from Europe. Our data are consistent with numerous seeds from multiple sources and a prolonged period of unrecognized community spreading. This work highlights the complementary role of genomic surveillance in addition to traditional epidemiological indicators.

A FACS-based approach to obtain viable eosinophils from human adipose tissue

Eosinophils have been widely investigated in asthma and allergic diseases. More recently, new insights into the biology of these cells has illustrated eosinophils contribute to homeostatic functions in health such as regulation of adipose tissue glucose metabolism. Human translational studies are limited by the difficulty of obtaining cells taken directly from their tissue environment, relying instead on eosinophils isolated from peripheral blood. Isolation techniques for tissue-derived eosinophils can result in unwanted cell or ribonuclease activation, leading to poor cell viability or RNA quality, which may impair analysis of effector activities of these cells. Here we demonstrate a technique to obtain eosinophils from human adipose tissue samples for the purpose of downstream molecular analysis. From as little as 2 g of intact human adipose tissue, greater than 10⁴ eosinophils were purified by fluorescence-activated cell sorting (FACS) protocol resulting in ≥ 99% purity and ≥ 95% viable eosinophils. We demonstrated that the isolated eosinophils could undergo epigenetic analysis to determine differences in DNA methylation in various settings. Here we focused on comparing eosinophils isolated from human peripheral blood vs human adipose tissue. Our results open the door to future mechanistic investigations to better understand the role of tissue resident eosinophils in different context.

A genome-wide scan for candidate lethal variants in Thoroughbred horses

Domestic animal populations are often characterised by high rates of inbreeding and low effective population sizes due to selective breeding practices. These practices can result in otherwise rare recessive deleterious alleles drifting to high frequencies, resulting in reduced fertility rates. This study aimed to identify potential recessive lethal haplotypes in the Thoroughbred horse breed, a closed population that has been selectively bred for racing performance. In this study, we identified a haplotype in the LY49B gene that shows strong evidence of being homozygous lethal, despite having high frequencies of heterozygotes in Thoroughbreds and other domestic horse breeds. Variant analysis of whole-genome sequence data identified two SNPs in the 3'UTR of the LY49B gene that may result in loss of function. Analysis of transcriptomic data from equine embryonic tissue revealed that LY49B is expressed in the trophoblast during placentation stage of development. These findings suggest that LY49B may have an essential, but as yet unknown function in the implantation stage of equine development. Further investigation of this region may allow for the development of a genetic test to improve fertility rates in horse populations. Identification of other lethal variants could assist in improving natural levels of fertility in horse populations.

A comparison of humans and baboons suggests germline mutation rates do not track cell divisions

In humans, most germline mutations are inherited from the father. This observation has been widely interpreted as reflecting the replication errors that accrue during spermatogenesis. If so, the male bias in mutation should be substantially lower in a closely related species with similar rates of spermatogonial stem cell divisions but a shorter mean age of reproduction. To test this hypothesis, we resequenced two 3-4 generation nuclear families (totaling 29 individuals) of olive baboons (Papio anubis), who reproduce at approximately 10 years of age on average, and analyzed the data in parallel with three 3-generation human pedigrees (26 individuals). We estimated a mutation rate per generation in baboons of 0.57×10-8 per base pair, approximately half that of humans. Strikingly, however, the degree of male bias in germline mutations is approximately 4:1, similar to that of humans-indeed, a similar male bias is seen across mammals that reproduce months, years, or decades after birth. These results mirror the finding in humans that the male mutation bias is stable with parental ages and cast further doubt on the assumption that germline mutations track cell divisions. Our mutation rate estimates for baboons raise a further puzzle, suggesting a divergence time between apes and Old World monkeys of 65 million years, too old to be consistent with the fossil record; reconciling them now requires not only a slowdown of the mutation rate per generation in humans but also in baboons.

Assessment of Pre-Analytical Sample Handling Conditions for Comprehensive Liquid Biopsy Analysis

Liquid biopsies as a minimally invasive approach have the potential to revolutionize molecular diagnostics. Yet, although protocols for sample handling and the isolation of circulating tumor DNA (ctDNA) are numerous, comprehensive guidelines for diagnostics and research considering all aspects of real-life multicenter clinical studies are currently not available. These include limitations in sample volume, transport, and blood collection tubes. We tested the impact of commonly used (EDTA and heparin) and specialized blood collection tubes and storage conditions on the yield and purity of cell-free DNA for the application in down-stream analysis. Moreover, we evaluated the feasibility of a combined workflow for ctDNA and tumor cell genomic testing and parallel flow cytometric analysis of leukocytes. For genomic analyses, EDTA tubes showed good results if stored for a maximum of 4 hours at room temperature or for up to 24 hours when stored at 4°C. Spike-in experiments revealed that EDTA tubes in combination with density gradient centrifugation allowed the parallel isolation of ctDNA, leukocytes, and low amounts of tumor cells (0.1%) and their immunophenotyping by flow cytometry and down-stream genomic analysis by whole genome sequencing. In conclusion, adhering to time and temperature limits allows the use of routine EDTA blood samples for liquid biopsy analyses. We further provide a workflow enabling the parallel analysis of cell-free and cellular features for disease monitoring and for clonal evolution studies.

Characterizing the dynamic and functional DNA methylation landscape in the developing human cortex

DNA methylation (DNAm) is a key epigenetic regulator of gene expression across development. The developing prenatal brain is a highly dynamic tissue, but our understanding of key drivers of epigenetic variability across development is limited. We, therefore, assessed genomic methylation at over 39 million sites in the prenatal cortex using whole-genome bisulfite sequencing and found loci and regions in which methylation levels are dynamic across development. We saw that DNAm at these loci was associated with nearby gene expression and enriched for enhancer chromatin states in prenatal brain tissue. Additionally, these loci were enriched for genes associated with neuropsychiatric disorders and genes involved with neurogenesis. We also found autosomal differences in DNAm between the sexes during prenatal development, though these have less clear functional consequences. We lastly confirmed that the dynamic methylation at this critical period is specifically CpG methylation, with generally low levels of CpH methylation. Our findings provide detailed insight into prenatal brain development as well as clues to the pathogenesis of psychiatric traits seen later in life.

Rapid Whole-Exome Sequencing as a Diagnostic Tool in a Neonatal/Pediatric Intensive Care Unit

Genetic disorders are the leading cause of infant morbidity and mortality. Due to the large number of genetic diseases, molecular and phenotype heterogeneity and often severe course, these diseases remain undiagnosed. In infants with a suspected acute monogenic disease, rapid whole-exome sequencing (R-WES) can be successfully performed. R-WES (singletons) was performed in 18 unrelated infants with a severe and/or progressing disease with the suspicion of genetic origin hospitalized in an Intensive Care Unit (ICU). Blood samples were also collected from the parents. The results from the R-WES were available after 5-14 days. A conclusive genetic diagnosis was obtained in 13 children, corresponding to an overall diagnostic yield of 72.2%. For nine patients, R-WES was used as a first-tier test. Eight patients were diagnosed with inborn errors of metabolism, mainly mitochondrial diseases. In two patients, the disease was possibly caused by variants in genes which so far have not been associated with human disease (NARS1 and DCAF5). R-WES proved to be an effective diagnostic tool for critically ill infants in ICUs suspected of having a genetic disorder. It also should be considered as a first-tier test after precise clinical description. The quickly obtained diagnosis impacts patient's medical management, and families can receive genetic counseling.

Molecular Characterization of Astrocytoma Progression Towards Secondary Glioblastomas Utilizing Patient-Matched Tumor Pairs

Astrocytomas are primary human brain tumors including diffuse or anaplastic astrocytomas that develop towards secondary glioblastomas over time. However, only little is known about molecular alterations that drive this progression. We measured multi-omics profiles of patient-matched astrocytoma pairs of initial and recurrent tumors from 22 patients to identify molecular alterations associated with tumor progression. Gene copy number profiles formed three major subcluters, but more than half of the patient-matched astrocytoma pairs differed in their gene copy number profiles like astrocytomas from different patients. Chromosome 10 deletions were not observed for diffuse astrocytomas, but occurred in corresponding recurrent tumors. Gene expression profiles formed three other major subclusters and patient-matched expression profiles were much more heterogeneous than their copy number profiles. Still, recurrent tumors showed a strong tendency to switch to the mesenchymal subtype. The direct progression of diffuse astrocytomas to secondary glioblastomas showed the largest number of transcriptional changes. Astrocytoma progression groups were further distinguished by signaling pathway expression signatures affecting cell division, interaction and differentiation. As expected, IDH1 was most frequently mutated closely followed by TP53, but also MUC4 involved in the regulation of apoptosis and proliferation was frequently mutated. Astrocytoma progression groups differed in their mutation frequencies of these three genes. Overall, patient-matched astrocytomas can differ substantially within and between patients, but still molecular signatures associated with the progression to secondary glioblastomas exist and should be analyzed for their potential clinical relevance in future studies.

The prognostic value of a Methylome-based Malignancy Density Scoring System to predict recurrence risk in early-stage Lung Adenocarcinoma

Current NCCN guidelines do not recommend the use of adjuvant chemotherapy for stage IA lung adenocarcinoma patients with R0 surgery. However, 25% to 40% of patients with stage IA disease experience recurrence. Stratifying patients according to the recurrence risk may tailor adjuvant therapy and surveillance imaging for those with a higher risk. However, prognostic markers are often identified by comparing high-risk and low-risk cases which might introduce bias due to the widespread interpatient heterogeneity. Here, we developed a scoring system quantifying the degree of field cancerization in adjacent normal tissues and revealed its association with disease-free survival (DFS).

Methods: We recruited a cohort of 44 patients with resected stage IA lung adenocarcinoma who did not receive adjuvant therapy. Both tumor and adjacent normal tissues were obtained from each patient and subjected to capture-based targeted genomic and epigenomic profiling. A novel methylome-based scoring system namely malignancy density ratio (MD ratio) was developed based on 39 patients by comparing tumor and corresponding adjacent normal tissues of each patient. A MD score was then obtained by Wald statistics. The correlations of MD ratio, MD score, and genomic features with clinical outcome were investigated.

Results: Patients with a high-risk MD ratio showed a significantly shorter postsurgical DFS compared with those with a low-risk MD ratio (HR=4.47, P=0.01). The MD ratio was not associated with T stage (P=1), tumor cell fraction (P=0.748) nor inflammatory status (p=0.548). Patients with a high-risk MD score also demonstrated an inferior DFS (HR=4.69, P=0.039). In addition, multivariate analysis revealed EGFR 19 del (HR=5.39, P=0.012) and MD score (HR= 7.90, P=0.01) were independent prognostic markers.

Conclusion: The novel methylome-based scoring system, developed by comparing the signatures between tumor and corresponding adjacent normal tissues of individual patients, largely minimizes the bias of interpatient heterogeneity and reveals a robust prognostic value in patients with resected lung adenocarcinoma.

Genetic Architecture of Parkinson's Disease in the Indian Population: Harnessing Genetic Diversity to Address Critical Gaps in Parkinson's Disease Research

Over the past two decades, our understanding of Parkinson's disease (PD) has been gleaned from the discoveries made in familial and/or sporadic forms of PD in the Caucasian population. The transferability and the clinical utility of genetic discoveries to other ethnically diverse populations are unknown. The Indian population has been under-represented in PD research. The Genetic Architecture of PD in India (GAP-India) project aims to develop one of the largest clinical/genomic bio-bank for PD in India. Specifically, GAP-India project aims to: (1) develop a pan-Indian deeply phenotyped clinical repository of Indian PD patients; (2) perform whole-genome sequencing in 500 PD samples to catalog Indian genetic variability and to develop an Indian PD map for the scientific community; (3) perform a genome-wide association study to identify novel loci for PD and (4) develop a user-friendly web-portal to disseminate results for the scientific community. Our "hub-spoke" model follows an integrative approach to develop a pan-Indian outreach to develop a comprehensive cohort for PD research in India. The alignment of standard operating procedures for recruiting patients and collecting biospecimens with international standards ensures harmonization of data/bio-specimen collection at the beginning and also ensures stringent quality control parameters for sample processing. Data sharing and protection policies follow the guidelines established by local and national authorities.We are currently in the recruitment phase targeting recruitment of 10,200 PD patients and 10,200 healthy volunteers by the end of 2020. GAP-India project after its completion will fill a critical gap that exists in PD research and will contribute a comprehensive genetic catalog of the Indian PD population to identify novel targets for PD.

Hi-C Identifies Complex Genomic Rearrangements and TAD-Shuffling in Developmental Diseases

Genome-wide analysis methods, such as array comparative genomic hybridization (CGH) and whole-genome sequencing (WGS), have greatly advanced the identification of structural variants (SVs) in the human genome. However, even with standard high-throughput sequencing techniques, complex rearrangements with multiple breakpoints are often difficult to resolve, and predicting their effects on gene expression and phenotype remains a challenge. Here, we address these problems by using high-throughput chromosome conformation capture (Hi-C) generated from cultured cells of nine individuals with developmental disorders (DDs). Three individuals had previously been identified as harboring duplications at the SOX9 locus and six had been identified with translocations. Hi-C resolved the positions of the duplications and was instructive in interpreting their distinct pathogenic effects, including the formation of new topologically associating domains (neo-TADs). Hi-C was very sensitive in detecting translocations, and it revealed previously unrecognized complex rearrangements at the breakpoints. In several cases, we observed the formation of fused-TADs promoting ectopic enhancer-promoter interactions that were likely to be involved in the disease pathology. In summary, we show that Hi-C is a sensible method for the detection of complex SVs in a clinical setting. The results help interpret the possible pathogenic effects of the SVs in individuals with DDs.

A chromosome-scale draft sequence of the Canada fleabane genome

BACKGROUND

Due to the accessibility of underlying technologies the 'Omics', in particular genomics, are becoming commonplace in several fields of research, including the study of agricultural pests. The weed community is starting to embrace these approaches; genome sequences have been made available in the past years, with several other sequencing projects underway, as promoted by the International Weed Genome Consortium. Chromosome-scale sequences are essential to fully exploit the power of genetics and genomics.

RESULTS

We report such an assembly for Conyza canadensis, an important agricultural weed. Third-generation sequencing technology was used to create a genome assembly of 426 megabases, of which nine chromosome-scale scaffolds cover more than 98% of the entire assembled sequence. As this weed was the first to be identified with glyphosate resistance, and since we do not have a firm handle on the genetic mechanisms responsible for several herbicide resistances in the species, the genome sequence was annotated with genes known to be associated with herbicide resistance. A high number of ABC-type transporters, cytochrome P450 and glycosyltransferases (159, 352 and 181, respectively) were identified among the list of ab initio predicted genes.

CONCLUSION

As C. canadensis has a small genome that is syntenic with other Asteraceaes, has a short life cycle and is relatively easy to cross, it has the potential to become a model weed species and, with the chromosome-scale genome sequence, contribute to a paradigm shift in the way non-target site resistance is studied. © 2020 Her Majesty the Queen in Right of CanadaPest Management Science © 2020 Society of Chemical Industry.

A de novo chromosome-level genome assembly of Coregonus sp. "Balchen": One representative of the Swiss Alpine whitefish radiation

Salmonids are of particular interest to evolutionary biologists due to their incredible diversity of life-history strategies and the speed at which many salmonid species have diversified. In Switzerland alone, over 30 species of Alpine whitefish from the subfamily Coregoninae have evolved since the last glacial maximum, with species exhibiting a diverse range of morphological and behavioural phenotypes. This, combined with the whole genome duplication which occurred in the ancestor of all salmonids, makes the Alpine whitefish radiation a particularly interesting system in which to study the genetic basis of adaptation and speciation and the impacts of ploidy changes and subsequent rediploidization on genome evolution. Although well-curated genome assemblies exist for many species within Salmonidae, genomic resources for the subfamily Coregoninae are lacking. To assemble a whitefish reference genome, we carried out PacBio sequencing from one wild-caught Coregonus sp. "Balchen" from Lake Thun to ~90× coverage. PacBio reads were assembled independently using three different assemblers, falcon, canu and wtdbg2 and subsequently scaffolded with additional Hi-C data. All three assemblies were highly contiguous, had strong synteny to a previously published Coregonus linkage map, and when mapping additional short-read data to each of the assemblies, coverage was fairly even across most chromosome-scale scaffolds. Here, we present the first de novo genome assembly for the Salmonid subfamily Coregoninae. The final 2.2-Gb wtdbg2 assembly included 40 scaffolds, an N50 of 51.9 Mb and was 93.3% complete for BUSCOs. The assembly consisted of ~52% transposable elements and contained 44,525 genes.

Pinpointing the Genomic Localizations of Chromatin-Associated Proteins: The Yesterday, Today, and Tomorrow of ChIP-seq

Chromatin-associated proteins are instrumental for controlling spatiotemporal gene expression. Determining where these proteins bind across the genome is critical for understanding gene regulation. A widely used technique at present is ChIP-seq, which leverages chromatin fragmentation, antibody-mediated enrichment, next-generation sequencing, and data analysis to uncover the genomic sequences and patterns of protein-DNA interactions. In this article, we will provide an overview of how ChIP-seq was developed, the key elements of the experimentation and data analysis pipeline, and the recent variations that push the boundaries of precision and cell number requirements. We will also briefly discuss how future development of ChIP-seq may further advance our understanding of chromatin biology. © 2019 by John Wiley & Sons, Inc.

Acquired resistance to DZNep-mediated apoptosis is associated with copy number gains of AHCY in a B-cell lymphoma model

BACKGROUND

Enhancer of zeste homolog 2 (EZH2) is considered an important driver of tumor development and progression by its histone modifying capabilities. Inhibition of EZH2 activity is thought to be a potent treatment option for eligible cancer patients with an aberrant EZH2 expression profile, thus the indirect EZH2 inhibitor 3-Deazaneplanocin A (DZNep) is currently under evaluation for its clinical utility. Although DZNep blocks proliferation and induces apoptosis in different tumor types including lymphomas, acquired resistance to DZNep may limit its clinical application.

METHODS

To investigate possible mechanisms of acquired DZNep resistance in B-cell lymphomas, we generated a DZNep-resistant clone from a previously DZNep-sensitive B-cell lymphoma cell line by long-term treatment with increasing concentrations of DZNep (ranging from 200 to 2000 nM) and compared the molecular profiles of resistant and wild-type clones. This comparison was done using molecular techniques such as flow cytometry, copy number variation assay (OncoScan and TaqMan assays), fluorescence in situ hybridization, Western blot, immunohistochemistry and metabolomics analysis.

RESULTS

Whole exome sequencing did not indicate the acquisition of biologically meaningful single nucleotide variants. Analysis of copy number alterations, however, demonstrated among other acquired imbalances an amplification (about 30 times) of the S-adenosyl-L-homocysteine hydrolase (AHCY) gene in the resistant clone. AHCY is a direct target of DZNep and is critically involved in the biological methylation process, where it catalyzes the reversible hydrolysis of S-adenosyl-L-homocysteine to L-homocysteine and adenosine. The amplification of the AHCY gene is paralleled by strong overexpression of AHCY at both the transcriptional and protein level, and persists upon culturing the resistant clone in a DZNep-free medium.

CONCLUSIONS

This study reveals one possible molecular mechanism how B-cell lymphomas can acquire resistance to DZNep, and proposes AHCY as a potential biomarker for investigation during the administration of EZH2-targeted therapy with DZNep.

Phased Diploid Genome Sequence for the Fast-Growing Microalga Picochlorum celeri

Picochlorum celeri is a fast-growing marine microalga with high biomass productivity. Here, we report the use of PacBio sequencing to assemble the phased diploid genome of P. celeri.

ATAC-seq identifies thousands of extrachromosomal circular DNA in cancer and cell lines

Extrachromosomal circular DNAs (eccDNAs) are somatically mosaic and contribute to intercellular heterogeneity in normal and tumor cells. Because short eccDNAs are poorly chromatinized, we hypothesized that they are sequenced by tagmentation in ATAC-seq experiments without any enrichment of circular DNA. Indeed, ATAC-seq identified thousands of eccDNAs in cell lines that were validated by inverse PCR and by metaphase FISH. ATAC-seq in gliomas and glioblastomas identify hundreds of eccDNAs, including one containing the well-known EGFR gene amplicon from chr7. More than 18,000 eccDNAs, many carrying known cancer driver genes, are identified in a pan-cancer analysis of ATAC-seq libraries from 23 tumor types. Somatically mosaic eccDNAs are identified by ATAC-seq even before amplification is recognized by genome-wide copy number variation measurements. Thus, ATAC-seq is a sensitive method to detect eccDNA present in a tumor at the pre-amplification stage and can be used to predict resistance to therapy.

Investigating APOE, APP-Aβ metabolism genes and Alzheimer's disease GWAS hits in brain small vessel ischemic disease

Alzheimer’s disease and small vessel ischemic disease frequently co-exist in the aging brain. However, pathogenic links between these 2 disorders are yet to be identified. Therefore we used Taqman genotyping, exome and RNA sequencing to investigate Alzheimer’s disease known pathogenic variants and pathways: APOE ε4 allele, APP-Aβ metabolism and late-onset Alzheimer’s disease main genome-wide association loci (APOE, BIN1, CD33, MS4A6A, CD2AP, PICALM, CLU, CR1, EPHA1, ABCA7) in 96 early-onset small vessel ischemic disease Caucasian patients and 368 elderly neuropathologically proven controls (HEX database) and in a mouse model of cerebral hypoperfusion. Only a minority of patients (29%) carried APOE ε4 allele. We did not detect any pathogenic mutation in APP, PSEN1 and PSEN2 and report a burden of truncating mutations in APP-Aß degradation genes. The single-variant association test identified 3 common variants with a likely protective effect on small vessel ischemic disease (0.54>OR > 0.32, adj. p-value <0.05) (EPHA1 p.M900V and p.V160A and CD33 p.A14V). Moreover, 5/17 APP-Aß catabolism genes were significantly upregulated (LogFC > 1, adj. p-val<0.05) together with Apoe, Ms4a cluster and Cd33 during brain hypoperfusion and their overexpression correlated with the ischemic lesion size. Finally, the detection of Aβ oligomers in the hypoperfused hippocampus supported the link between brain ischemia and Alzheimer’s disease pathology.

An Improved Phenotype-Driven Tool for Rare Mendelian Variant Prioritization: Benchmarking Exomiser on Real Patient Whole-Exome Data

Next-generation sequencing has revolutionized rare disease diagnostics, but many patients remain without a molecular diagnosis, particularly because many candidate variants usually survive despite strict filtering. Exomiser was launched in 2014 as a Java tool that performs an integrative analysis of patients' sequencing data and their phenotypes encoded with Human Phenotype Ontology (HPO) terms. It prioritizes variants by leveraging information on variant frequency, predicted pathogenicity, and gene-phenotype associations derived from human diseases, model organisms, and protein-protein interactions. Early published releases of Exomiser were able to prioritize disease-causative variants as top candidates in up to 97% of simulated whole-exomes. The size of the tested real patient datasets published so far are very limited. Here, we present the latest Exomiser version 12.0.1 with many new features. We assessed the performance using a set of 134 whole-exomes from patients with a range of rare retinal diseases and known molecular diagnosis. Using default settings, Exomiser ranked the correct diagnosed variants as the top candidate in 74% of the dataset and top 5 in 94%; not using the patients' HPO profiles (i.e., variant-only analysis) decreased the performance to 3% and 27%, respectively. In conclusion, Exomiser is an effective support tool for rare Mendelian phenotype-driven variant prioritization.

Identification of a novel archaea virus, detected in hydrocarbon polluted Hungarian and Canadian samples

Metagenomics is a helpful tool for the analysis of unculturable organisms and viruses. Viruses that target bacteria and archaea play important roles in the microbial diversity of various ecosystems. Here we show that Methanosarcina virus MV (MetMV), the second Methanosarcina sp. virus with a completely determined genome, is characteristic of hydrocarbon pollution in environmental (soil and water) samples. It was highly abundant in Hungarian hydrocarbon polluted samples and its genome was also present in the NCBI SRA database containing reads from hydrocarbon polluted samples collected in Canada, indicating the stability of its niche and the marker feature of this virus. MetMV, as the only currently identified marker virus for pollution in environmental samples, could contribute to the understanding of the complicated network of prokaryotes and their viruses driving the decomposition of environmental pollutants.

Natural Selection Shapes Variation in Genome-wide Recombination Rate in Drosophila pseudoobscura

While recombination is widely recognized to be a key modulator of numerous evolutionary phenomena, we have a poor understanding of how recombination rate itself varies and evolves within a species. Here, we performed a comprehensive study of recombination rate (rate of meiotic crossing over) in two natural populations of Drosophila pseudoobscura from Utah and Arizona, USA. We used an amplicon sequencing approach to obtain high-quality genotypes in approximately 8,000 individual backcrossed offspring (17 mapping populations with roughly 530 individuals each), for which we then quantified crossovers. Interestingly, variation in recombination rate within and between populations largely manifested as differences in genome-wide recombination rate rather than remodeling of the local recombination landscape. Comparing populations, we discovered individuals from the Utah population displayed on average 8% higher crossover rates than the Arizona population, a statistically significant difference. Using a Q_ST-F_ST analysis, we found that this difference in crossover rate was dramatically higher than expected under neutrality, indicating that this difference may have been driven by natural selection. Finally, using a combination of short- and long-read whole-genome sequencing, we found no significant association between crossover rate and structural variation at the 200-400 kb scale. Our results demonstrate that (1) there is abundant variation in genome-wide crossover rate in natural populations, (2) at the 200-400 kb scale, recombination rate appears to vary largely genome-wide, rather than in specific intervals, and (3) interpopulation differences in recombination rate may be the result of local adaptation.

Nuclear-mitochondrial DNA segments resemble paternally inherited mitochondrial DNA in humans

Several strands of evidence question the dogma that human mitochondrial DNA (mtDNA) is inherited exclusively down the maternal line, most recently in three families where several individuals harbored a 'heteroplasmic haplotype' consistent with biparental transmission. Here we report a similar genetic signature in 7 of 11,035 trios, with allelic fractions of 5-25%, implying biparental inheritance of mtDNA in 0.06% of offspring. However, analysing the nuclear whole genome sequence, we observe likely large rare or unique nuclear-mitochondrial DNA segments (mega-NUMTs) transmitted from the father in all 7 families. Independently detecting mega-NUMTs in 0.13% of fathers, we see autosomal transmission of the haplotype. Finally, we show the haplotype allele fraction can be explained by complex concatenated mtDNA-derived sequences rearranged within the nuclear genome. We conclude that rare cryptic mega-NUMTs can resemble paternally mtDNA heteroplasmy, but find no evidence of paternal transmission of mtDNA in humans.

Biallelic variants in KYNU cause a multisystemic syndrome with hand hyperphalangism

Catel-Manzke syndrome is characterized by the combination of Pierre Robin sequence and radial deviation, shortening as well as clinodactyly of the index fingers, due to an accessory ossification center. Mutations in TGDS have been identified as one cause of Catel-Manzke syndrome, but cannot be found as causative in every patient with the clinical diagnosis. We performed a chromosome microarray and/or exome sequencing in three patients with hand hyperphalangism, heart defect, short stature, and mild to severe developmental delay, all of whom were initially given a clinical diagnosis of Catel-Manzke syndrome. In one patient, we detected a large deletion of exons 1-8 and the missense variant c.1282C > T (p.Arg428Trp) in KYNU (NM_003937.2), whereas homozygous missense variants in KYNU were found in the other two patients (c.989G > A (p.Arg330Gln) and c.326G > C (p.Trp109Ser)). Plasma and urine metabolomic analysis of two patients indicated a block along the tryptophan catabolic pathway and urine organic acid analysis showed excretion of xanthurenic acid. Biallelic loss-of-function mutations in KYNU were recently described as a cause of NAD deficiency with vertebral, cardiac, renal and limb defects; however, no hand hyperphalangism was described in those patients, and Catel-Manzke syndrome was not discussed as a differential diagnosis. In conclusion, we present unrelated patients identified with biallelic variants in KYNU leading to kynureninase deficiency and xanthurenic aciduria as a very likely cause of their hyperphalangism, heart defect, short stature, and developmental delay. We suggest performance of urine organic acid analysis in patients with suspected Catel-Manzke syndrome, particularly in those with cardiac or vertebral defects or without mutations in TGDS.

An Annotated Chromosome-Level Reference Genome of the Red-Eared Slider Turtle (Trachemys scripta elegans)

Among vertebrates, turtles have many unique characteristics providing biologists with opportunities to study novel evolutionary innovations and processes. We present here a high-quality, partially phased, and chromosome-level Red-Eared Slider (Trachemys scripta elegans, TSE) genome as a reference for future research on turtle and tetrapod evolution. This TSE assembly is 2.269 Gb in length, has one of the highest scaffold N50 and N90 values of any published turtle genome to date (N50 = 129.68 Mb and N90 = 19 Mb), and has a total of 28,415 annotated genes. We introduce synteny analyses using BUSCO single-copy orthologs, which reveal two chromosome fusion events accounting for differences in chromosome counts between emydids and other cryptodire turtles and reveal many fission/fusion events for birds, crocodiles, and snakes relative to TSE. This annotated chromosome-level genome will provide an important reference genome for future studies on turtle, vertebrate, and chromosome evolution.

KCNN2 mutation in autosomal-dominant tremulous myoclonus-dystonia

BACKGROUND AND PURPOSE

Despite recent advances in neurogenetics that have facilitated the identification of a number of dystonia genes, many familial dystonia syndromes remain without known cause. The aim of the study was to identify the cause of autosomal dominant tremulous myoclonus-dystonia in a UK kindred with affected individuals in three generations.

METHODS

Known genetic causes of myoclonus-dystonia were excluded. We combined clinical and electrophysiological phenotyping with whole-exome sequencing and Sanger sequencing to identify candidate causal variants in a family with tremulous myoclonus-dystonia.

RESULTS

The core phenotype consisted of childhood-onset dystonia predominantly affecting hands and neck, with a fast tremor with superimposed myoclonus and, in some individuals, subtle cerebellar signs. We identified a novel missense variant in potassium calcium-activated channel subfamily N member 2 (KCNN2) [NM_021614:c.1112G>A:p.(Gly371Glu)], which was the only variant that we were able to identify as segregating with the phenotype over three generations. This variant, which is absent from the most recent version of gnomAD, was predicted to be deleterious by SIFT and PolyPhen-2 and had an overall CADD score of 29.7.

CONCLUSIONS

KCNN2, a member of the KCNN family of potassium channel genes, is highly conserved across species and in humans is highly expressed in the brain, particularly the cerebellum. KCNN2 mutations have never been described as pathological in human disease, but are recognized abnormalities in two rodent models of fast, jerky tremor. Segregation, absence of the variant in the normal population and in-silico prediction of a deleterious effect together with animal models compatible with the clinical phenotype are all in line with KCNN2 mutations being a plausible cause underlying myoclonus-dystonia.

Comparative and population genomics approaches reveal the basis of adaptation to deserts in a small rodent

Organisms that live in deserts offer the opportunity to investigate how species adapt to environmental conditions that are lethal to most plants and animals. In the hot deserts of North America, high temperatures and lack of water are conspicuous challenges for organisms living there. The cactus mouse (Peromyscus eremicus) displays several adaptations to these conditions, including low metabolic rate, heat tolerance, and the ability to maintain homeostasis under extreme dehydration. To investigate the genomic basis of desert adaptation in cactus mice, we built a chromosome‐level genome assembly and resequenced 26 additional cactus mouse genomes from two locations in southern California (USA). Using these data, we integrated comparative, population, and functional genomic approaches. We identified 16 gene families exhibiting significant contractions or expansions in the cactus mouse compared to 17 other Myodontine rodent genomes, and found 232 sites across the genome associated with selective sweeps. Functional annotations of candidate gene families and selective sweeps revealed a pervasive signature of selection at genes involved in the synthesis and degradation of proteins, consistent with the evolution of cellular mechanisms to cope with protein denaturation caused by thermal and hyperosmotic stress. Other strong candidate genes included receptors for bitter taste, suggesting a dietary shift towards chemically defended desert plants and insects, and a growth factor involved in lipid metabolism, potentially involved in prevention of dehydration. Understanding how species adapted to deserts will provide an important foundation for predicting future evolutionary responses to increasing temperatures, droughts and desertification in the cactus mouse and other species.

Fast analysis of scATAC-seq data using a predefined set of genomic regions

Background: Analysis of scATAC-seq data has been recently scaled to thousands of cells. While processing of other types of single cell data was boosted by the implementation of alignment-free techniques, pipelines available to process scATAC-seq data still require large computational resources. We propose here an approach based on pseudoalignment, which reduces the execution times and hardware needs at little cost for precision. Methods: Public data for 10k PBMC were downloaded from 10x Genomics web site. Reads were aligned to various references derived from DNase I Hypersensitive Sites (DHS) using kallisto and quantified with bustools. We compared our results with the ones publicly available derived by cellranger-atac. We subsequently tested our approach on scATAC-seq data for K562 cell line. Results: We found that kallisto does not introduce biases in quantification of known peaks; cells groups identified are consistent with the ones identified from standard method. We also found that cell identification is robust when analysis is performed using DHS-derived reference in place of de novo identification of ATAC peaks. Lastly, we found that our approach is suitable for reliable quantification of gene activity based on scATAC-seq signal, thus allows for efficient labelling of cell groups based on marker genes. Conclusions: Analysis of scATAC-seq data by means of kallisto produces results in line with standard pipelines while being considerably faster; using a set of known DHS sites as reference does not affect the ability to characterize the cell populations.

Gene clustering and copy number variation in alkaloid metabolic pathways of opium poppy

Genes in plant secondary metabolic pathways enable biosynthesis of a range of medically and industrially important compounds, and are often clustered on chromosomes. Here, we study genomic clustering in the benzylisoquinoline alkaloid (BIA) pathway in opium poppy (Papaver somniferum), exploring relationships between gene expression, copy number variation, and metabolite production. We use Hi-C to improve the existing draft genome assembly, yielding chromosome-scale scaffolds that include 35 previously unanchored BIA genes. We find that co-expression of BIA genes increases within clusters and identify candidates with unknown function based on clustering and covariation in expression and alkaloid production. Copy number variation in critical BIA genes correlates with stark differences in alkaloid production, linking noscapine production with an 11-gene deletion, and increased thebaine/decreased morphine production with deletion of a T6ODM cluster. Our results show that the opium poppy genome is still dynamically evolving in ways that contribute to medically and industrially important phenotypes.

Hybridization and introgression drive genome evolution of Dutch elm disease pathogens

Hybridization and the resulting introgression can drive the success of invasive species via the rapid acquisition of adaptive traits. The Dutch elm disease pandemics in the past 100 years were caused by three fungal lineages with permeable reproductive barriers: Ophiostoma ulmi, Ophiostoma novo-ulmi subspecies novo-ulmi and Ophiostoma novo-ulmi subspecies americana. Using whole-genome sequences and growth phenotyping of a worldwide collection of isolates, we show that introgression has been the main driver of genomic diversity and that it impacted fitness-related traits. Introgressions contain genes involved in host-pathogen interactions and reproduction. Introgressed isolates have enhanced growth rate at high temperature and produce different necrosis sizes on an in vivo model for pathogenicity. In addition, lineages diverge in many pathogenicity-associated genes and exhibit differential mycelial growth in the presence of a proxy of a host defence compound, implying an important role of host trees in the molecular and functional differentiation of these pathogens.

Genetic basis of neurodevelopmental disorders in 103 Jordanian families

We recruited 103 families from Jordan with neurodevelopmental disorders (NDD) and patterns of inheritance mostly suggestive of autosomal recessive inheritance. In each family, we investigated at least one affected individual using exome sequencing and an in-house diagnostic variant interpretation pipeline including a search for copy number variation. This approach led us to identify the likely molecular defect in established disease genes in 37 families. We could identify 25 pathogenic nonsense and 11 missense variants as well as 3 pathogenic copy number variants and 1 repeat expansion. Notably, 11 of the disease-causal variants occurred de novo. In addition, we prioritized a homozygous frameshift variant in PUS3 in two sisters with intellectual disability. To our knowledge, PUS3 has been postulated only recently as a candidate disease gene for intellectual disability in a single family with three affected siblings. Our findings provide additional evidence to establish loss of PUS3 function as a cause of intellectual disability.

Assessment of Potential Clinical Role for Exome Sequencing in Schizophrenia

BACKGROUND

There is increasing evidence that certain genetic variants increase the risk of schizophrenia and other neurodevelopmental disorders. Exome sequencing has been shown to have a high diagnostic yield for developmental disability and testing for copy number variants has been advocated for schizophrenia. The diagnostic yield for exome sequencing in schizophrenia is unknown.

METHOD

A sample of 591 exome-sequenced schizophrenia cases and their parents were screened for disruptive and damaging variants in autosomal genes listed in the Genomics England panels for intellectual disability and other neurological disorders.

RESULTS

Previously reported disruptive de novo variants were noted in SETD1A, POGZ, SCN2A, and ZMYND11. Although the loss of function of ZMYND11 is a recognized cause of intellectual disability, it has not previously been noted as a risk factor for schizophrenia. A damaging de novo variant of uncertain significance was noted in NRXN1. A previously reported homozygous damaging variant in BLM is predicted to cause Bloom syndrome in 1 case and 1 case was homozygous for a damaging variant in MCPH1, a result of uncertain significance. There were more than 400 disruptive and damaging variants in the target genes in cases but similar numbers were seen among untransmitted parental alleles and none appeared to be clinically significant.

CONCLUSIONS

The diagnostic yield from exome sequencing in schizophrenia is low. Disruptive and damaging variants seen in known neuropsychiatric genes should not be automatically assumed to have an etiological role if observed in a patient with schizophrenia.

Genomics of sorghum local adaptation to a parasitic plant

Host-parasite coevolution can maintain high levels of genetic diversity in traits involved in species interactions. In many systems, host traits exploited by parasites are constrained by use in other functions, leading to complex selective pressures across space and time. Here, we study genome-wide variation in the staple crop Sorghum bicolor (L.) Moench and its association with the parasitic weed Striga hermonthica (Delile) Benth., a major constraint to food security in Africa. We hypothesize that geographic selection mosaics across gradients of parasite occurrence maintain genetic diversity in sorghum landrace resistance. Suggesting a role in local adaptation to parasite pressure, multiple independent loss-of-function alleles at sorghum LOW GERMINATION STIMULANT 1 (LGS1) are broadly distributed among African landraces and geographically associated with S. hermonthica occurrence. However, low frequency of these alleles within S. hermonthica-prone regions and their absence elsewhere implicate potential trade-offs restricting their fixation. LGS1 is thought to cause resistance by changing stereochemistry of strigolactones, hormones that control plant architecture and below-ground signaling to mycorrhizae and are required to stimulate parasite germination. Consistent with trade-offs, we find signatures of balancing selection surrounding LGS1 and other candidates from analysis of genome-wide associations with parasite distribution. Experiments with CRISPR-Cas9-edited sorghum further indicate that the benefit of LGS1-mediated resistance strongly depends on parasite genotype and abiotic environment and comes at the cost of reduced photosystem gene expression. Our study demonstrates long-term maintenance of diversity in host resistance genes across smallholder agroecosystems, providing a valuable comparison to both industrial farming systems and natural communities.

Unmasking Retinitis Pigmentosa complex cases by a whole genome sequencing algorithm based on open-access tools: hidden recessive inheritance and potential oligogenic variants

Background

Retinitis Pigmentosa (RP) is a clinically and genetically heterogeneous disorder that results in inherited blindness. Despite the large number of genes identified, only ~ 60% of cases receive a genetic diagnosis using targeted-sequencing. The aim of this study was to design a whole genome sequencing (WGS) based approach to increase the diagnostic yield of complex Retinitis Pigmentosa cases.

Methods

WGS was conducted in three family members, belonging to one large apparent autosomal dominant RP family that remained unsolved by previous studies, using Illumina TruSeq library preparation kit and Illumina HiSeq X platform. Variant annotation, filtering and prioritization were performed using a number of open-access tools and public databases. Sanger sequencing of candidate variants was conducted in the extended family members.

Results

We have developed and optimized an algorithm, based on the combination of different open-access tools, for variant prioritization of WGS data which allowed us to reduce significantly the number of likely causative variants pending to be manually assessed and segregated. Following this algorithm, four heterozygous variants in one autosomal recessive gene (USH2A) were identified, segregating in pairs in the affected members. Additionally, two pathogenic alleles in ADGRV1 and PDZD7 could be contributing to the phenotype in one patient.

Conclusions

The optimization of a diagnostic algorithm for WGS data analysis, accompanied by a hypothesis-free approach, have allowed us to unmask the genetic cause of the disease in one large RP family, as well as to reassign its inheritance pattern which implies differences in the clinical management of these cases. These results contribute to increasing the number of cases with apparently dominant inheritance that carry causal mutations in recessive genes, as well as the possible involvement of various genes in the pathogenesis of RP in one patient. Moreover, our WGS-analysis approach, based on open-access tools, can easily be implemented by other researchers and clinicians to improve the diagnostic yield of additional patients with inherited retinal dystrophies.

Synaptonemal Complex-Deficient Drosophila melanogaster Females Exhibit Rare DSB Repair Events, Recurrent Copy-Number Variation, and an Increased Rate of de Novo Transposable Element Movement

Genetic stability depends on the maintenance of a variety of chromosome structures and the precise repair of DNA breaks. During meiosis, programmed double-strand breaks (DSBs) made in prophase I are normally repaired as gene conversions or crossovers. DSBs can also be made by other mechanisms, such as the movement of transposable elements (TEs), which must also be resolved. Incorrect repair of these DNA lesions can lead to mutations, copy-number changes, translocations, and/or aneuploid gametes. In Drosophila melanogaster, as in most organisms, meiotic DSB repair occurs in the presence of a rapidly evolving multiprotein structure called the synaptonemal complex (SC). Here, whole-genome sequencing is used to investigate the fate of meiotic DSBs in D. melanogaster mutant females lacking functional SC, to assay for de novo CNV formation, and to examine the role of the SC in transposable element movement in flies. The data indicate that, in the absence of SC, copy-number variation still occurs and meiotic DSB repair by gene conversion occurs infrequently. Remarkably, an 856-kilobase de novo CNV was observed in two unrelated individuals of different genetic backgrounds and was identical to a CNV recovered in a previous wild-type study, suggesting that recurrent formation of large CNVs occurs in Drosophila. In addition, the rate of novel TE insertion was markedly higher than wild type in one of two SC mutants tested, suggesting that SC proteins may contribute to the regulation of TE movement and insertion in the genome. Overall, this study provides novel insight into the role that the SC plays in genome stability and provides clues as to why the sequence, but not structure, of SC proteins is rapidly evolving.

Genome-wide analysis of HOXC4 and HOXC6 regulated genes and binding sites in prostate cancer cells

Aberrant expression of HOXC6 and HOXC4 is commonly detected in prostate cancer. The high expression of these transcription factors is associated with aggressive prostate cancer and can predict cancer recurrence after treatment. Thus, HOXC4 and HOXC6 are clinically relevant biomarkers of aggressive prostate cancer. However, the molecular mechanisms by which these HOXC genes contribute to prostate cancer is not yet understood. To begin to address the role of HOXC4 and HOXC6 in prostate cancer, we performed RNA-seq analyses before and after siRNA-mediated knockdown of HOXC4 and/or HOXC6 and also performed ChIP-seq to identify genomic binding sites for both of these transcription factors. Our studies demonstrate that HOXC4 and HOXC6 co-localize with HOXB13, FOXA1 and AR, three transcription factors previously shown to contribute to the development of prostate cancer. We suggest that the aberrantly upregulated HOXC4 and HOXC6 proteins may compete with HOXB13 for binding sites, thus altering the prostate transcriptome. This competition model may be applicable to many different human cancers that display increased expression of a HOX transcription factor.

Generation and characterization of a mitotane-resistant adrenocortical cell line

Mitotane is the only drug approved for the therapy of adrenocortical carcinoma (ACC). Its clinical use is limited by the occurrence of relapse during therapy. To investigate the underlying mechanisms in vitro, we here generated mitotane-resistant cell lines. After long-term pulsed treatment of HAC-15 human adrenocortical carcinoma cells with 70 µM mitotane, we isolated monoclonal cell populations of treated cells and controls and assessed their respective mitotane sensitivities by MTT assay. We performed exome sequencing and electron microscopy, conducted gene expression microarray analysis and determined intracellular lipid concentrations in the presence and absence of mitotane. Clonal cell lines established after pulsed treatment were resistant to mitotane (IC50 of 102.2 ± 7.3 µM (n = 12) vs 39.4 ± 6.2 µM (n = 6) in controls (biological replicates, mean ± s.d., P = 0.0001)). Unlike nonresistant clones, resistant clones maintained normal mitochondrial and nucleolar morphology during mitotane treatment. Resistant clones largely shared structural and single nucleotide variants, suggesting a common cell of origin. Resistance depended, in part, on extracellular lipoproteins and was associated with alterations in intracellular lipid homeostasis, including levels of free cholesterol, as well as decreased steroid production. By gene expression analysis, resistant cells showed profound alterations in pathways including steroid metabolism and transport, apoptosis, cell growth and Wnt signaling. These studies establish an in vitro model of mitotane resistance in ACC and point to underlying molecular mechanisms. They may enable future studies to overcome resistance in vitro and improve ACC treatment in vivo.

Genome-wide discovery of DNA variants in cucumber somaclonal lines

The emergence of somaclonal variability in in vitro cultures is undesirable during micropropagation, but this phenomenon may be a source of genetic variability sought by breeders. The main factors that affect the appearance of variability are known, but the exact mechanism has not yet been determined. In this paper, we used next-generation sequencing and comparative genomics to study changes in the genomes of cucumber lines resulting from in vitro regeneration and somaclonal mutation in comparison to a reference, the highly inbred B10 line. The total number of obtained polymorphisms differed between the three somaclonal lines S1, S2 and S3, with 8369, 7591 and 44510, respectively. Polymorphisms occurred most frequently in non-coding regions and were mainly SNPs. High-impact changes accounted for 1%-3% of all polymorphisms and most often caused an open reading frame shift. Functional analysis of genes affected by high impact variants showed that they were related to transport, biosynthetic processes, nucleotide-containing compounds and cellular protein modification processes. The obtained results indicated significant factors affecting somaclonal variability and the appearance of changes in the genome, and demonstrated a lack of dependence between phenotype and the number of genomic polymorphisms.

Oral Cavity Squamous Cell Carcinoma Xenografts Retain Complex Genotypes and Intertumor Molecular Heterogeneity

Herein, we report an oral cavity squamous cell carcinoma (OCSCC) patient-derived xenograft (PDX) platform, with genomic annotation useful for co-clinical trial and mechanistic studies. Genomic analysis included whole-exome sequencing (WES) and transcriptome sequencing (RNA-seq) on 16 tumors and matched PDXs and additional whole-genome sequencing (WGS) on 9 of these pairs as a representative subset of a larger OCSCC PDX repository (n = 63). In 12 models with high purity, more than 90% of variants detected in the tumor were retained in the matched PDX. The genomic landscape across these PDXs reflected OCSCC molecular heterogeneity, including previously described basal, mesenchymal, and classical molecular subtypes. To demonstrate the integration of PDXs into a clinical trial framework, we show that pharmacological intervention in PDXs parallels clinical response and extends patient data. Together, these data describe a repository of OCSCC-specific PDXs and illustrate conservation of primary tumor genotypes, intratumoral heterogeneity, and co-clinical trial application.

Campbell et al. report the genomic fidelity of patient-derived xenograft models from oral cavity squamous cell carcinomas. These models conserve the mutation and expression profile of their matched tumors, validating their use for co-clinical trial and mechanistic studies.

Identifying, understanding, and correcting technical artifacts on the sex chromosomes in next-generation sequencing data

Background

Mammalian X and Y chromosomes share a common evolutionary origin and retain regions of high sequence similarity. Similar sequence content can confound the mapping of short next-generation sequencing reads to a reference genome. It is therefore possible that the presence of both sex chromosomes in a reference genome can cause technical artifacts in genomic data and affect downstream analyses and applications. Understanding this problem is critical for medical genomics and population genomic inference.

Results

Here, we characterize how sequence homology can affect analyses on the sex chromosomes and present XYalign, a new tool that (1) facilitates the inference of sex chromosome complement from next-generation sequencing data; (2) corrects erroneous read mapping on the sex chromosomes; and (3) tabulates and visualizes important metrics for quality control such as mapping quality, sequencing depth, and allele balance. We find that sequence homology affects read mapping on the sex chromosomes and this has downstream effects on variant calling. However, we show that XYalign can correct mismapping, resulting in more accurate variant calling. We also show how metrics output by XYalign can be used to identify XX and XY individuals across diverse sequencing experiments, including low- and high-coverage whole-genome sequencing, and exome sequencing. Finally, we discuss how the flexibility of the XYalign framework can be leveraged for other uses including the identification of aneuploidy on the autosomes. XYalign is available open source under the GNU General Public License (version 3).

Conclusions

Sex chromsome sequence homology causes the mismapping of short reads, which in turn affects downstream analyses. XYalign provides a reproducible framework to correct mismapping and improve variant calling on the sex chromsomes.

Establishment of a novel human cell line retaining the characteristics of the original pancreatic adenocarcinoma, and evaluation of MEK as a therapeutic target

Pancreatic cancer is a lethal solid malignancy with limited therapeutic options. The development of novel therapeutic drugs requires adequate new cell line models. A new pancreatic cancer cell line, designated PDXPC1, was established from one pancreatic ductal adenocarcinoma (PDAC) patient-derived xenograft. The PDXPC1 cells were stably cultured for >2 years and had a stable short tandem repeat profile. The PDXPC1 cell line retained the key mutations of the primary tumor, along with the epithelial origin and other important protein expression. The PDXPC1 cells induced rapid in vivo tumor growth, both subcutaneously and orthotopically, in a mouse model with an elevated CA199 level. The PDXPC1 cells showed weak growth, invasion and migration potency compared to another pancreatic cancer cell line, but were relatively resistant to multiple anti-cancer drugs. Interestingly, the MEK inhibitor trametinib significantly inhibited the proliferation of PDXPC1 cells, and not that of Panc-1 cells, by inactivating MEK/ERK/MYC signaling and activating the apoptotic pathway via Bcl-2 degradation. In conclusion, the PDXPC1 cell line, capturing the major characteristics of the primary tumor, may be a suitable tool for studying the underlying mechanisms of chemo-resistance in PDAC and developing new targeted therapeutic options.

Assessment of Tumor Mutational Burden in Pediatric Tumors by Real-Life Whole-Exome Sequencing and In Silico Simulation of Targeted Gene Panels: How the Choice of Method Could Affect the Clinical Decision?

Background: Tumor mutational burden (TMB) is an emerging genomic biomarker in cancer that has been associated with improved response to immune checkpoint inhibitors (ICIs) in adult cancers. It was described that variability in TMB assessment is introduced by different laboratory techniques and various settings of bioinformatic pipelines. In pediatric oncology, no study has been published describing this variability so far. Methods: In our study, we performed whole exome sequencing (WES, both germline and somatic) and calculated TMB in 106 patients with high-risk/recurrent pediatric solid tumors of 28 distinct cancer types. Subsequently, we used WES data for TMB calculation using an in silico approach simulating two The Food and Drug Administration (FDA)-approved/authorized comprehensive genomic panels for cancer. Results: We describe a strong correlation between WES-based and panel-based TMBs; however, we show that this high correlation is significantly affected by inclusion of only a few hypermutated cases. In the series of nine cases, we determined TMB in two sequentially collected tumor tissue specimens and observed an increase in TMB along with tumor progression. Furthermore, we evaluated the extent to which potential ICI indication could be affected by variability in techniques and bioinformatic pipelines used for TMB assessment. We confirmed that this technological variability could significantly affect ICI indication in pediatric cancer patients; however, this significance decreases with the increasing cut-off values. Conclusions: For the first time in pediatric oncology, we assessed the reliability of TMB estimation across multiple pediatric cancer types using real-life WES and in silico analysis of two major targeted gene panels and confirmed a significant technological variability to be introduced by different laboratory techniques and various settings of bioinformatic pipelines.

Plasma Circulating Tumor DNA and Clonal Hematopoiesis in Metastatic Renal Cell Carcinoma

BACKGROUND

There is a lack of molecularly-informed biomarkers for patients with metastatic renal cell carcinoma (RCC). Plasma cell-free DNA (cfDNA) sequencing is a minimally-invasive alternative to tissue for profiling the genome in other cancers but relevance in metastatic RCC remains unclear.

MATERIALS AND METHODS

Whole blood was collected from 55 patients with metastatic RCC. Plasma cfDNA and leukocyte DNA were subjected to targeted sequencing across 981 cancer genes. Matched tumor tissue from 14 patients was analyzed.

RESULTS

Thirty-three percent of patients had evidence for RCC-derived circulating tumor DNA (ctDNA), significantly lower than patients with metastatic prostate or bladder cancer analyzed using the same approach. Among ctDNA-positive patients, ctDNA fraction averaged only 3.9% and showed no strong association with clinical variables. In these patients, the most commonly mutated genes were VHL, BAP1, and PBRM1, and matched tissue concordance was 77%. Evidence of somatic expansions unrelated to RCC, such as clonal hematopoiesis of indeterminate potential, were detected in 43% of patients. Pathogenic germline mutations in DNA repair genes were detected in 11% of patients. CtDNA-positive patients had shorter overall survival and progression-free survival on first-line therapy. Patients with evidence of clonal hematopoiesis of indeterminate potential had an intermediate prognosis compared with ctDNA-positive and -negative patients.

CONCLUSIONS

CfDNA sequencing enables straightforward characterization of the somatic RCC genome in a minority of patients with metastatic RCC. Owing to low ctDNA abundance, and the presence of non-RCC derived somatic clones in circulation, cfDNA sequencing may not be a simple pan-patient alternative to tissue biopsy in metastatic RCC.