Divergent sensory and immune gene evolution in sea turtles with contrasting demographic and life histories.

Sea turtles represent an ancient lineage of marine vertebrates that evolved from terrestrial ancestors over 100 Mya. The genomic basis of the unique physiological and ecological traits enabling these species to thrive in diverse marine habitats remains largely unknown. Additionally, many populations have drastically declined due to anthropogenic activities over the past two centuries, and their recovery is a high global conservation priority. We generated and analyzed high-quality reference genomes for the leatherback (Dermochelys coriacea) and green (Chelonia mydas) turtles, representing the two extant sea turtle families. These genomes are highly syntenic and homologous, but localized regions of noncollinearity were associated with higher copy numbers of immune, zinc-finger, and olfactory receptor (OR) genes in green turtles, with ORs related to waterborne odorants greatly expanded in green turtles. Our findings suggest that divergent evolution of these key gene families may underlie immunological and sensory adaptations assisting navigation, occupancy of neritic versus pelagic environments, and diet specialization. Reduced collinearity was especially prevalent in microchromosomes, with greater gene content, heterozygosity, and genetic distances between species, supporting their critical role in vertebrate evolutionary adaptation. Finally, diversity and demographic histories starkly contrasted between species, indicating that leatherback turtles have had a low yet stable effective population size, exhibit extremely low diversity compared with other reptiles, and harbor a higher genetic load compared with green turtles, reinforcing concern over their persistence under future climate scenarios. These genomes provide invaluable resources for advancing our understanding of evolution and conservation best practices in an imperiled vertebrate lineage.

Draft de novo Genome Assembly of the Elusive Jaguarundi, Puma yagouaroundi.

The Puma lineage within the family Felidae consists of 3 species that last shared a common ancestor around 4.9 million years ago. Whole-genome sequences of 2 species from the lineage were previously reported: the cheetah (Acinonyx jubatus) and the mountain lion (Puma concolor). The present report describes a whole-genome assembly of the remaining species, the jaguarundi (Puma yagouaroundi). We sequenced the genome of a male jaguarundi with 10X Genomics linked reads and assembled the whole-genome sequence. The assembled genome contains a series of scaffolds that reach the length of chromosome arms and is similar in scaffold contiguity to the genome assemblies of cheetah and puma, with a contig N50 = 100.2 kbp and a scaffold N50 = 49.27 Mbp. We assessed the assembled sequence of the jaguarundi genome using BUSCO, aligned reads of the sequenced individual and another published female jaguarundi to the assembled genome, annotated protein-coding genes, repeats, genomic variants and their effects with respect to the protein-coding genes, and analyzed differences of the 2 jaguarundis from the reference mitochondrial genome. The jaguarundi genome assembly and its annotation were compared in quality, variants, and features to the previously reported genome assemblies of puma and cheetah. Computational analyzes used in the study were implemented in transparent and reproducible way to allow their further reuse and modification.

Diagnostic utility of telomere length testing in a hospital-based setting.

Telomere length (TL) predicts the onset of cellular senescence in vitro but the diagnostic utility of TL measurement in clinical settings is not fully known. We tested the value of TL measurement by flow cytometry and FISH (flowFISH) in patients with mutations in telomerase and telomere maintenance genes. TL had a discrete and reproducible normal range with definable upper and lower boundaries. While TL above the 50th age-adjusted percentile had a 100% negative predictive value for clinically relevant mutations, the lower threshold in mutation carriers was age-dependent, and adult mutation carriers often overlapped with the lowest decile of controls. The extent of telomere shortening correlated with the age at diagnosis as well as the short telomere syndrome phenotype. Extremely short TL caused bone marrow failure and immunodeficiency in children and young adults, while milder defects manifested as pulmonary fibrosis-emphysema in adults. We prospectively examined whether TL altered treatment decisions for newly diagnosed idiopathic bone marrow failure patients and found abnormally short TL enriched for patients with mutations in some inherited bone marrow failure genes, such as RUNX1, in addition to telomerase and telomere maintenance genes. The result was actionable, altering the choice of treatment regimen and/or hematopoietic stem cell donor in one-fourth of the cases (9 of 38, 24%). We conclude that TL measurement by flowFISH, when used for targeted clinical indications and in limited settings, can influence treatment decisions in ways that improve outcome.

Characterization of the Rat Osteosarcoma Cell Line UMR-106 by Long-Read Technologies Identifies a Large Block of Amplified Genes Associated with Human Disease.

BACKGROUND/OBJECTIVES: The rat osteosarcoma cell line UMR-106 is widely used for the study of bone cancer biology but it has not been well characterized with modern genomic methods. METHODS: To better understand the biology of UMR-106 cells we used a combination of optical genome mapping (OGM), long-read sequencing nanopore sequencing and RNA sequencing.The UMR-106 genome was compared to a strain-matched Sprague-Dawley rat for variants associated with human osteosarcoma while expression data were contrasted with a public osteoblast dataset. RESULTS: Using the COSMIC database to identify the most affected genes in human osteosarcomas we found somatic mutations in Tp53 and H3f3a. OGM identified a relatively small number of differences between the cell line and a strain-matched control animal but did detect a ~45 Mb block of amplification that included Myc on chromosome 7 which was confirmed by long-read sequencing. The amplified region showed several blocks of non-contiguous rearranged sequence implying complex rearrangements during their formation and included 14 genes reported as biomarkers in human osteosarcoma, many of which also showed increased transcription. A comparison of 5mC methylation from the nanopore reads of tumor and control samples identified genes with distinct differences including the OS marker Cdkn2a. CONCLUSIONS: This dataset illustrates the value of long DNA methods for the characterization of cell lines and how inter-species analysis can inform us about the genetic nature underlying mutations that underpin specific tumor types. The data should be a valuable resource for investigators studying osteosarcoma, in general, and specifically the UMR-106 model.

Genomic Markers Associated with Cytomegalovirus DNAemia in Kidney Transplant Recipients.

Human cytomegalovirus (CMV) is a major pathogen after solid organ transplantation, leading to high morbidity and mortality. Transplantation from a CMV-seropositive donor to a CMV-seronegative recipient (D+/R-) is associated with high risk of CMV disease. However, that risk is not uniform, suggesting a role for host factors in immune control of CMV. To identify host genetic factors that control CMV DNAemia post transplantation, we performed a whole-exome association study in two cohorts of D+/R- kidney transplant recipients. Quantitative CMV DNA was measured for at least one year following transplantation. Several CMV-protective single-nucleotide polymorphisms (SNPs) were identified in the first cohort (72 patients) but were not reproducible in the second cohort (126 patients). A meta-analysis of both cohorts revealed several SNPs that were significantly associated with protection from CMV DNAemia. The copy number variation of several genes was significantly different between recipients with and without CMV DNAemia. Amongst patients with CMV DNAemia in the second cohort, several variants of interest (p < 5 × 10-5), the most common of which was NLRC5, were associated with peak viral load. We provide new predictive genetic markers for protection of CMV DNAemia. These markers should be validated in larger cohorts.

A Chromosome-Length Assembly of the Hawaiian Monk Seal (Neomonachus schauinslandi): A History of "Genetic Purging" and Genomic Stability.

The Hawaiian monk seal (HMS) is the single extant species of tropical earless seals of the genus Neomonachus. The species survived a severe bottleneck in the late 19th century and experienced subsequent population declines until becoming the subject of a NOAA-led species recovery effort beginning in 1976 when the population was fewer than 1000 animals. Like other recovering species, the Hawaiian monk seal has been reported to have reduced genetic heterogeneity due to the bottleneck and subsequent inbreeding. Here, we report a chromosomal reference assembly for a male animal produced using a variety of methods. The final assembly consisted of 16 autosomes, an X, and portions of the Y chromosomes. We compared variants in this animal to other HMS and to a frequently sequenced human sample, confirming about 12% of the variation seen in man. To confirm that the reference animal was representative of the HMS, we compared his sequence to that of 10 other individuals and noted similarly low variation in all. Variation in the major histocompatibility (MHC) genes was nearly absent compared to the orthologous human loci. Demographic analysis predicts that Hawaiian monk seals have had a long history of small populations preceding the bottleneck, and their current low levels of heterozygosity may indicate specialization to a stable environment. When we compared our reference assembly to that of other species, we observed significant conservation of chromosomal architecture with other pinnipeds, especially other phocids. This reference should be a useful tool for future evolutionary studies as well as the long-term management of this species.

A novel canis lupus familiaris reference genome improves variant resolution for use in breed-specific GWAS.

Reference genome fidelity is critically important for genome wide association studies, yet most vary widely from the study population. A typical whole genome sequencing approach implies short-read technologies resulting in fragmented assemblies with regions of ambiguity. Further information is lost by economic necessity when genotyping populations, as lower resolution technologies such as genotyping arrays are commonly used. Here, we present a phased reference genome for Canis lupus familiaris using high molecular weight DNA-sequencing technologies. We tested wet laboratory and bioinformatic approaches to demonstrate a minimum workflow to generate the 2.4 gigabase genome for a Labrador Retriever. The de novo assembly required eight Oxford Nanopore R9.4 flowcells (∼23X depth) and running a 10X Genomics library on the equivalent of one lane of an Illumina NovaSeq S1 flowcell (∼88X depth), bringing the cost of generating a nearly complete reference genome to less than $10K (USD). Mapping of short-read data from 10 Labrador Retrievers against this reference resulted in 1% more aligned reads versus the current reference (CanFam3.1, P < 0.001), and a 15% reduction of variant calls, increasing the chance of identifying true, low-effect size variants in a genome-wide association studies. We believe that by incorporating the cost to produce a full genome assembly into any large-scale genotyping project, an investigator can improve study power, decrease costs, and optimize the overall scientific value of their study.

Chromosomal-level genome assembly of the scimitar-horned oryx: Insights into diversity and demography of a species extinct in the wild.

Captive populations provide a valuable insurance against extinctions in the wild. However, they are also vulnerable to the negative impacts of inbreeding, selection and drift. Genetic information is therefore considered a critical aspect of conservation management. Recent developments in sequencing technologies have the potential to improve the outcomes of management programmes; however, the transfer of these approaches to applied conservation has been slow. The scimitar-horned oryx (Oryx dammah) is a North African antelope that has been extinct in the wild since the early 1980s and is the focus of a large-scale and long-term reintroduction project. To enable the selection of suitable founder individuals, facilitate post-release monitoring and improve captive breeding management, comprehensive genomic resources are required. Here, we used 10X Chromium sequencing together with Hi-C contact mapping to develop a chromosomal-level genome assembly for the species. The resulting assembly contained 29 chromosomes with a scaffold N50 of 100.4 Mb, and displayed strong chromosomal synteny with the cattle genome. Using resequencing data from six additional individuals, we demonstrated relatively high genetic diversity in the scimitar-horned oryx compared to other mammals, despite it having experienced a strong founding event in captivity. Additionally, the level of diversity across populations varied according to management strategy. Finally, we uncovered a dynamic demographic history that coincided with periods of climate variation during the Pleistocene. Overall, our study provides a clear example of how genomic data can uncover valuable insights into captive populations and contributes important resources to guide future management decisions of an endangered species.

Characterization of the genomic architecture and mutational spectrum of a small cell prostate carcinoma.

We present the use of a series of laboratory, analytical and interpretation methods to investigate personalized cancer care for a case of small cell prostate carcinoma (SCPC), a rare and aggressive tumor with poor prognosis, for which the underlying genomic architecture and mutational spectrum has not been well characterized. We performed both SNP genotyping and exome sequencing of a Virchow node metastasis from a patient with SCPC. A variety of methods were used to analyze and interpret the tumor genome for copy number variation, loss of heterozygosity (LOH), somatic mosaicism and mutations in genes from known cancer pathways. The combination of genotyping and exome sequencing approaches provided more information than either technique alone. The results showed widespread evidence of copy number changes involving most chromosomes including the possible loss of both alleles of CDKN1B (p27/Kip1). LOH was observed for the regions encompassing the tumor suppressors TP53, RB1, and CHD1. Predicted damaging somatic mutations were observed in the retained TP53 and RB1 alleles. Mutations in other genes that may be functionally relevant were noted, especially the recently reported high confidence cancer drivers FOXA1 and CCAR1. The disruption of multiple cancer drivers underscores why SCPC may be such a difficult cancer to manage.

Identification of an HMGB3 frameshift mutation in a family with an X-linked colobomatous microphthalmia syndrome using whole-genome and X-exome sequencing.

IMPORTANCE: Microphthalmias are rare disorders whose genetic bases are not fully understood. HMGB3 is a new candidate gene for X-linked forms of this disease. OBJECTIVE: To identify the causative gene in a pedigree with an X-linked colobomatous microphthalmos phenotype. DESIGN, SETTING, AND PARTICIPANTS: Whole-genome sequencing and chromosome X-exome-targeted sequencing were performed at the High Throughput Sequencing Laboratory of the Genetic Resources Core Facility at the Johns Hopkins University School of Medicine on the DNA of the male proband and informatically filtered to identify rare variants. Polymerase chain reaction and Sanger sequencing were used to confirm the variant in the proband and the carrier status of his mother. Thirteen unrelated male patients with a similar phenotype were also screened. MAIN OUTCOMES AND MEASURES: Whole-genome and X-exome sequencing to identify a frameshift variant in HMGB3. RESULTS: A 2-base pair frameshift insertion (c.477_478insTA, coding for p.Lys161Ilefs*54) in the HGMB3 gene was found in the proband and his carrier mother but not in the unrelated patients. The mutation, confirmed by 3 orthogonal methods, alters an evolutionarily conserved region of the HMGB3 protein from a negatively charged polyglutamic acid tract to a positively charged arginine-rich motif that is likely to interfere with normal protein function. CONCLUSIONS AND RELEVANCE: In this family, microphthalmia, microcephaly, intellectual disability, and short stature are associated with a mutation on the X chromosome in the HMGB3 gene. HMGB3 should be considered when performing genetic studies of patients with similar phenotypes.

Mutations in Alström protein impair terminal differentiation of cardiomyocytes.

Cardiomyocyte cell division and replication in mammals proceed through embryonic development and abruptly decline soon after birth. The process governing cardiomyocyte cell cycle arrest is poorly understood. Here we carry out whole-exome sequencing in an infant with evidence of persistent postnatal cardiomyocyte replication to determine the genetic risk factors. We identify compound heterozygous ALMS1 mutations in the proband, and confirm their presence in her affected sibling, one copy inherited from each heterozygous parent. Next, we recognize homozygous or compound heterozygous truncating mutations in ALMS1 in four other children with high levels of postnatal cardiomyocyte proliferation. Alms1 mRNA knockdown increases multiple markers of proliferation in cardiomyocytes, the percentage of cardiomyocytes in G2/M phases, and the number of cardiomyocytes by 10% in cultured cells. Homozygous Alms1-mutant mice have increased cardiomyocyte proliferation at 2 weeks postnatal compared with wild-type littermates. We conclude that deficiency of Alström protein impairs postnatal cardiomyocyte cell cycle arrest.

Mutations in the TGF-ß repressor SKI cause Shprintzen-Goldberg syndrome with aortic aneurysm.

Elevated transforming growth factor (TGF)-ß signaling has been implicated in the pathogenesis of syndromic presentations of aortic aneurysm, including Marfan syndrome (MFS) and Loeys-Dietz syndrome (LDS). However, the location and character of many of the causal mutations in LDS intuitively imply diminished TGF-ß signaling. Taken together, these data have engendered controversy regarding the specific role of TGF-ß in disease pathogenesis. Shprintzen-Goldberg syndrome (SGS) has considerable phenotypic overlap with MFS and LDS, including aortic aneurysm. We identified causative variation in ten individuals with SGS in the proto-oncogene SKI, a known repressor of TGF-ß activity. Cultured dermal fibroblasts from affected individuals showed enhanced activation of TGF-ß signaling cascades and higher expression of TGF-ß-responsive genes relative to control cells. Morpholino-induced silencing of SKI paralogs in zebrafish recapitulated abnormalities seen in humans with SGS. These data support the conclusions that increased TGF-ß signaling is the mechanism underlying SGS and that high signaling contributes to multiple syndromic presentations of aortic aneurysm.