Discovery and genotyping of structural variation from long-read haploid genome sequence data

In an effort to more fully understand the full spectrum of human genetic variation, we generated deep single-molecule, real-time (SMRT) sequencing data from two haploid human genomes. By using an assembly-based approach (SMRT-SV), we systematically assessed each genome independently for structural variants (SVs) and indels resolving the sequence structure of 461,553 genetic variants from 2 bp to 28 kbp in length. We find that &gt;89% of these variants have been missed as part of analysis of the 1000 Genomes Project even after adjusting for more common variants (MAF &gt; 1%). We estimate that this theoretical human diploid differs by as much as ∼16 Mbp with respect to the human reference, with long-read sequencing data providing a fivefold increase in sensitivity for genetic variants ranging in size from 7 bp to 1 kbp compared with short-read sequence data. Although a large fraction of genetic variants were not detected by short-read approaches, once the alternate allele is sequence-resolved, we show that 61% of SVs can be genotyped in short-read sequence data sets with high accuracy. Uncoupling discovery from genotyping thus allows for the majority of this missed common variation to be genotyped in the human population. Interestingly, when we repeat SV detection on a pseudodiploid genome constructed in silico by merging the two haploids, we find that ∼59% of the heterozygous SVs are no longer detected by SMRT-SV. These results indicate that haploid resolution of long-read sequencing data will significantly increase sensitivity of SV detection.</jats:p

Impact of animal strain on gene expression in a rat model of acute cardiac rejection

<p>Abstract</p> <p>Background</p> <p>The expression levels of many genes show wide natural variation among strains or populations. This study investigated the potential for animal strain-related genotypic differences to confound gene expression profiles in acute cellular rejection (ACR). Using a rat heart transplant model and 2 different rat strains (Dark Agouti, and Brown Norway), microarrays were performed on native hearts, transplanted hearts, and peripheral blood mononuclear cells (PBMC).</p> <p>Results</p> <p>In heart tissue, strain alone affected the expression of only 33 probesets while rejection affected the expression of 1368 probesets (FDR 10% and FC ≥ 3). Only 13 genes were affected by both strain and rejection, which was < 1% (13/1368) of all probesets differentially expressed in ACR. However, for PBMC, strain alone affected 265 probesets (FDR 10% and FC ≥ 3) and the addition of ACR had little further effect. Pathway analysis of these differentially expressed strain effect genes connected them with immune response, cell motility and cell death, functional themes that overlap with those related to ACR. After accounting for animal strain, additional analysis identified 30 PBMC candidate genes potentially associated with ACR.</p> <p>Conclusion</p> <p>In ACR, genetic background has a large impact on the transcriptome of immune cells, but not heart tissue. Gene expression studies of ACR should avoid study designs that require cross strain comparisons between leukocytes.</p

Improved assembly and variant detection of a haploid human genome using single-molecule, high-fidelity long reads

The sequence and assembly of human genomes using long-read sequencing technologies has revolutionized our understanding of structural variation and genome organization. We compared the accuracy, continuity, and gene annotation of genome assemblies generated from either high-fidelity (HiFi) or continuous long-read (CLR) datasets from the same complete hydatidiform mole human genome. We find that the HiFi sequence data assemble an additional 10% of duplicated regions and more accurately represent the structure of tandem repeats, as validated with orthogonal analyses. As a result, an additional 5 Mbp of pericentromeric sequences are recovered in the HiFi assembly, resulting in a 2.5-fold increase in the NG50 within 1 Mbp of the centromere (HiFi 480.6 kbp, CLR 191.5 kbp). Additionally, the HiFi genome assembly was generated in significantly less time with fewer computational resources than the CLR assembly. Although the HiFi assembly has significantly improved continuity and accuracy in many complex regions of the genome, it still falls short of the assembly of centromeric DNA and the largest regions of segmental duplication using existing assemblers. Despite these shortcomings, our results suggest that HiFi may be the most effective standalone technology for de novo assembly of human genomes

COMPUTER SIMULATIONS OF POSSIBLE FUTURES FOR TWO FLOCKS OF WHOOPING CRANES

We conducted computer simulations using the program VORTEX (version 7) to project population sizes, growth rates, genetic diversity, and probabilities of extinction over the next 100 years for 2 flocks of whooping cranes (Grus americana), the Aransas/Wood Buffalo population and the experimental Florida population. Standard runs based on best estimates of demographic. genetic, and environmental parameter values were used as a baseline to which several alternative scenarios were compared. Results generally supported the conclusion of the earlier Population Viability Assessment (Mirande et al. 1991) that the AransaslWood Buffalo population will continue to grow steadily with less than a 1 % probability of extinction. It was noted, however, that a combination of negative factors such as shrinking habitat and increased probabilities of catastrophes accompanied by increased mortality rates could put this population at risk. Results for the Florida population were less optimistic. The standard run produced a population growth rate (r) of only 0.0026 for the next 100 years, and this shifted down to -0.0001 over a 200-year time frame. Adult mortality in this flock would have to be about 20% lower than the predicted value (10%) in order to raise growth rates to above r = 0.02. Amount and duration of supplementation of the Florida flock had minimal impacts on the long-tenn growth rate of the flock. It is the enduring rates of mortality, breeding, and disease risk that will have major effects on this population. For example, if disease risks tum out to be greater than the best-estimate scenario, this population could face a relatively high risk of extinction (17%). The formula for success in Florida is lower adult mortality, lower age of first breeding, lower disease risk, and higher productivity than the best-guess estimates. Fortunately, there are some potential management interventions (e.g., predator control, vaccines and health monitoring, selective introductions to balance the sex ratio of the flock) that may be able to push the odds in favor of success

Anticancer effects of phytocannabinoids used with chemotherapy in leukaemia cells can be improved by altering the sequence of their administration

Phytocannabinoids possess anticancer activity when used alone, and a number have also been shown to combine favourably with each other in vitro in leukaemia cells to generate improved activity. We have investigated the effect of pairing cannabinoids and assessed their anticancer activity in cell line models. Those most effective were then used with the common anti-leukaemia drugs cytarabine and vincristine, and the effects of this combination therapy on cell death studied in vitro. Results show a number of cannabinoids could be paired together to generate an effect superior to that achieved if the components were used individually. For example, in HL60 cells, the IC50 values at 48 h for cannabidiol (CBD) and tetrahydrocannabinol (THC) when used alone were 8 and 13 µM, respectively; however, if used together, it was 4 µM. Median-effect analysis confirmed the benefit of using cannabinoids in pairs, with calculated combination indices being <1 in a number of cases. The most efficacious cannabinoid-pairs subsequently synergised further when combined with the chemotherapy agents, and were also able to sensitise leukaemia cells to their cytotoxic effects. The sequence of administration of these drugs was important though; using cannabinoids after chemotherapy resulted in greater induction of apoptosis, whilst this was the opposite when the schedule of administration was reversed. Our results suggest that when certain cannabinoids are paired together, the resulting product can be combined synergistically with common anti-leukaemia drugs allowing the dose of the cytotoxic agents to be dramatically reduced yet still remain efficacious. Nevertheless, the sequence of drug administration is crucial to the success of these triple combinations and should be considered when planning such treatments

An evolutionary driver of interspersed segmental duplications in primates

Background The complex interspersed pattern of segmental duplications in humans is responsible for rearrangements associated with neurodevelopmental disease, including the emergence of novel genes important in human brain evolution. We investigate the evolution of LCR16a, a putative driver of this phenomenon that encodes one of the most rapidly evolving human–ape gene families, nuclear pore interacting protein (NPIP). Results Comparative analysis shows that LCR16a has independently expanded in five primate lineages over the last 35 million years of primate evolution. The expansions are associated with independent lineage-specific segmental duplications flanking LCR16a leading to the emergence of large interspersed duplication blocks at non-orthologous chromosomal locations in each primate lineage. The intron-exon structure of the NPIP gene family has changed dramatically throughout primate evolution with different branches showing characteristic gene models yet maintaining an open reading frame. In the African ape lineage, we detect signatures of positive selection that occurred after a transition to more ubiquitous expression among great ape tissues when compared to Old World and New World monkeys. Mouse transgenic experiments from baboon and human genomic loci confirm these expression differences and suggest that the broader ape expression pattern arose due to mutational changes that emerged in cis. Conclusions LCR16a promotes serial interspersed duplications and creates hotspots of genomic instability that appear to be an ancient property of primate genomes. Dramatic changes to NPIP gene structure and altered tissue expression preceded major bouts of positive selection in the African ape lineage, suggestive of a gene undergoing strong adaptive evolution

Recurrent structural variation, clustered sites of selection, and disease risk for the complement factor H (CFH) gene family

Data deposition: The data reported in this paper have been deposited as a National Center for Biotechnology Information BioProject (accession no. PRJNA401648). Author contributions: S.C. and E.E.E. designed research; S.C., C.B., L.H., K.P., K.M.M., M.S., A.E.W., V.D., T.A.G.-L., and R.K.W. performed research; S.C., J.H., C.B., L.H., K.P., K.M.M., M.S., A.E.W., V.D., F.G., A.J.R., R.H.G., T.A.G.-L., R.K.W., B.H.F.W., P.N.B., R.A., and E.E.E. contributed new reagents/analytic tools; S.C., B.J.N., J.H., and E.E.E. analyzed data; and S.C., B.J.N., and E.E.E. wrote the paper.Peer reviewedPublisher PD

A high-quality bonobo genome refines the analysis of hominid evolution

The divergence of chimpanzee and bonobo provides one of the few examples of recent hominid speciation1,2. Here we describe a fully annotated, high-quality bonobo genome assembly, which was constructed without guidance from reference genomes by applying a multiplatform genomics approach. We generate a bonobo genome assembly in which more than 98% of genes are completely annotated and 99% of the gaps are closed, including the resolution of about half of the segmental duplications and almost all of the full-length mobile elements. We compare the bonobo genome to those of other great apes1,3,4,5 and identify more than 5,569 fixed structural variants that specifically distinguish the bonobo and chimpanzee lineages. We focus on genes that have been lost, changed in structure or expanded in the last few million years of bonobo evolution. We produce a high-resolution map of incomplete lineage sorting and estimate that around 5.1% of the human genome is genetically closer to chimpanzee or bonobo and that more than 36.5% of the genome shows incomplete lineage sorting if we consider a deeper phylogeny including gorilla and orangutan. We also show that 26% of the segments of incomplete lineage sorting between human and chimpanzee or human and bonobo are non-randomly distributed and that genes within these clustered segments show significant excess of amino acid replacement compared to the rest of the genome

Allorecognition in the Tasmanian Devil (Sarcophilus harrisii), an Endangered Marsupial Species with Limited Genetic Diversity

Tasmanian devils (Sarcophilus harrisii) are on the verge of extinction due to a transmissible cancer, devil facial tumour disease (DFTD). This tumour is an allograft that is transmitted between individuals without immune recognition of the tumour cells. The mechanism to explain this lack of immune recognition and acceptance is not well understood. It has been hypothesized that lack of genetic diversity at the Major Histocompatibility Complex (MHC) allowed the tumour cells to grow in genetically similar hosts without evoking an immune response to alloantigens. We conducted mixed lymphocyte reactions and skin grafts to measure functional MHC diversity in the Tasmanian devil population. The limited MHC diversity was sufficient to produce measurable mixed lymphocyte reactions. There was a wide range of responses, from low or no reaction to relatively strong responses. The highest responses occurred when lymphocytes from devils from the east of Tasmania were mixed with lymphocytes from devils from the west of Tasmania. All of the five successful skin allografts were rejected within 14 days after surgery, even though little or no MHC I and II mismatches were found. Extensive T-cell infiltration characterised the immune rejection. We conclude that Tasmanian devils are capable of allogeneic rejection. Consequently, a lack of functional allorecognition mechanisms in the devil population does not explain the transmission of a contagious cancer