The incidence of malignant melanoma continues to increase each year with poor prognosis for survival in many relapse cases. To reverse this trend, whole body response measures are needed to discover collaborative paths to primary and secondary malignancy. Several species of fish provide excellent melanoma models because fish and human melanocytes both appear in the epidermis, and fish and human pigment cell tumors share conserved gene expression signatures. For the first time, we have examined the whole body transcriptome response to invasive melanoma as a prelude to using transcriptome profiling to screen for drugs in a medaka (Oryzias latipes) model. We generated RNA-seq data from whole body RNA isolates for controls and melanoma fish. After testing for differential expression, 396 genes had significantly different expression (adjusted p-value <0.02) in the whole body transcriptome between melanoma and control fish; 379 of these genes were matched to human orthologs with 233 having annotated human gene symbols and 14 matched genes that contain putative deleterious variants in human melanoma at varying levels of recurrence. A detailed canonical pathway evaluation for significant enrichment showed the top scoring pathway to be antigen presentation but also included the expected melanocyte development and pigmentation signaling pathway. Results revealed a profound down-regulation of genes involved in the immune response, especially the innate immune system. We hypothesize that the developing melanoma actively suppresses the immune system responses of the body in reacting to the invasive malignancy, and that this mal-adaptive response contributes to disease progression, a result that suggests our whole-body transcriptomic approach merits further use. In these findings, we also observed novel genes not yet identified in human melanoma expression studies and uncovered known and new candidate drug targets for further testing in this malignant melanoma medaka model.

Genome-wide recombination is essential for genome stability, evolution, and speciation. Mouse Tex11, an X-linked meiosis-specific gene, promotes meiotic recombination and chromosomal synapsis. Here, we report that TEX11 is mutated in infertile men with non-obstructive azoospermia and that an analogous mutation in the mouse impairs meiosis. Genetic screening of a large cohort of idiopathic infertile men reveals that TEX11 mutations, including frameshift and splicing acceptor site mutations, cause infertility in 1% of azoospermic men. Functional evaluation of three analogous human TEX11 missense mutations in transgenic mouse models identified one mutation (V748A) as a potential infertility allele and found two mutations non-causative. In the mouse model, an intronless autosomal Tex11 transgene functionally substitutes for the X-linked Tex11 gene, providing genetic evidence for the X-to-autosomal retrotransposition evolution phenomenon. Furthermore, we find that TEX11 protein levels modulate genome-wide recombination rates in both sexes. These studies indicate that TEX11 alleles affecting expression level or substituting single amino acids may contribute to variations in recombination rates between sexes and among individuals in humans.

Gibbons are small arboreal apes that display an accelerated rate of evolutionary chromosomal rearrangement and occupy a key node in the primate phylogeny between Old World monkeys and great apes. Here we present the assembly and analysis of a northern white-cheeked gibbon (Nomascus leucogenys) genome. We describe the propensity for a gibbon-specific retrotransposon (LAVA) to insert into chromosome segregation genes and alter transcription by providing a premature termination site, suggesting a possible molecular mechanism for the genome plasticity of the gibbon lineage. We further show that the gibbon genera (Nomascus, Hylobates, Hoolock and Symphalangus) experienced a near-instantaneous radiation ∼5 million years ago, coincident with major geographical changes in southeast Asia that caused cycles of habitat compression and expansion. Finally, we identify signatures of positive selection in genes important for forelimb development (TBX5) and connective tissues (COL1A1) that may have been involved in the adaptation of gibbons to their arboreal habitat.

Natural populations subjected to strong environmental selection pressures offer a window into the genetic underpinnings of evolutionary change. Cavefish populations, Astyanax mexicanus (Teleostei: Characiphysi), exhibit repeated, independent evolution for a variety of traits including eye degeneration, pigment loss, increased size and number of taste buds and mechanosensory organs, and shifts in many behavioural traits. Surface and cave forms are interfertile making this system amenable to genetic interrogation; however, lack of a reference genome has hampered efforts to identify genes responsible for changes in cave forms of A. mexicanus. Here we present the first de novo genome assembly for Astyanax mexicanus cavefish, contrast repeat elements to other teleost genomes, identify candidate genes underlying quantitative trait loci (QTL), and assay these candidate genes for potential functional and expression differences. We expect the cavefish genome to advance understanding of the evolutionary process, as well as, analogous human disease including retinal dysfunction.

The hookworm Necator americanus is the predominant soil-transmitted human parasite. Adult worms feed on blood in the small intestine, causing iron-deficiency anemia, malnutrition, growth and development stunting in children, and severe morbidity and mortality during pregnancy in women. We report sequencing and assembly of the N. americanus genome (244 Mb, 19,151 genes). Characterization of this first hookworm genome sequence identified genes orchestrating the hookworm's invasion of the human host, genes involved in blood feeding and development, and genes encoding proteins that represent new potential drug targets against hookworms. N. americanus has undergone a considerable and unique expansion of immunomodulator proteins, some of which we highlight as potential treatments against inflammatory diseases. We also used a protein microarray to demonstrate a postgenomic application of the hookworm genome sequence. This genome provides an invaluable resource to boost ongoing efforts toward fundamental and applied postgenomic research, including the development of new methods to control hookworm and human immunological diseases.

Tarsiers are phylogenetically located between the most basal strepsirrhines and the most derived anthropoid primates. While they share morphological features with both groups, they also possess uncommon primate characteristics, rendering their evolutionary history somewhat obscure. To investigate the molecular basis of such attributes, we present here a new genome assembly of the Philippine tarsier (Tarsius syrichta), and provide extended analyses of the genome and detailed history of transposable element insertion events. We describe the silencing of Alu monomers on the lineage leading to anthropoids, and recognize an unexpected abundance of long terminal repeat-derived and LINE1-mobilized transposed elements (Tarsius interspersed elements; TINEs). For the first time in mammals, we identify a complete mitochondrial genome insertion within the nuclear genome, then reveal tarsier-specific, positive gene selection and posit population size changes over time. The genomic resources and analyses presented here will aid efforts to more fully understand the ancient characteristics of primate genomes.

A higher prevalence of inherited disorders among companion animals are often rooted in their historical restricted artificial selection for a variety of observed phenotypes that eventually decreased genetic diversity. Cats have been afflicted with many inherited diseases due to domestication and intense breed selection. Advances in sequencing technology have generated a more comprehensive way to access genetic information from an individual, allowing identification of putative disease-causing variants and in practice a means to avoid their spread and thus better pedigree management. We examine variants in three domestic shorthair cats and then calculated overall genetic diversity to extrapolate the benefits of this data for breeding programs within a feline colony.

Amplicons--large, nearly identical repeats in direct or inverted orientation--are abundant in the male-specific region of the human Y chromosome (MSY) and provide targets for intrachromosomal non-allelic homologous recombination (NAHR). Thus far, NAHR events resulting in deletions, duplications, inversions, or isodicentric chromosomes have been reported only for amplicon pairs located exclusively on the short arm (Yp) or the long arm (Yq). Here we report our finding of four men with Y chromosomes that evidently formed by intrachromosomal NAHR between inverted repeat pairs comprising one amplicon on Yp and one amplicon on Yq. In two men with spermatogenic failure, sister-chromatid crossing-over resulted in pseudoisoYp chromosome formation and loss of distal Yq. In two men with normal spermatogenesis, intrachromatid crossing-over generated pericentric inversions. These findings highlight the recombinogenic nature of the MSY, as intrachromosomal NAHR occurs for nearly all Y-chromosome amplicon pairs, even those located on opposing chromosome arms.

One of the most disputed issues in primate evolution and thus of our own primate roots, is the phylogenetic position of the Southeast Asian tarsier. While much molecular data indicate a basal place in the primate tree shared with strepsirrhines (prosimian monophyly hypothesis), data also exist supporting either an earlier divergence in primates (tarsier-first hypothesis) or a close relationship with anthropoid primates (Haplorrhini hypothesis). The use of retroposon insertions embedded in the Tarsius genome afforded us the unique opportunity to directly test all three hypotheses via three pairwise genome alignments. From millions of retroposons, we found 104 perfect orthologous insertions in both tarsiers and anthropoids to the exclusion of strepsirrhines, providing conflict-free evidence for the Haplorrhini hypothesis, and none supporting either of the other two positions. Thus, tarsiers are clearly the sister group to anthropoids in the clade Haplorrhini.

Gorillas are humans' closest living relatives after chimpanzees, and are of comparable importance for the study of human origins and evolution. Here we present the assembly and analysis of a genome sequence for the western lowland gorilla, and compare the whole genomes of all extant great ape genera. We propose a synthesis of genetic and fossil evidence consistent with placing the human-chimpanzee and human-chimpanzee-gorilla speciation events at approximately 6 and 10 million years ago. In 30% of the genome, gorilla is closer to human or chimpanzee than the latter are to each other; this is rarer around coding genes, indicating pervasive selection throughout great ape evolution, and has functional consequences in gene expression. A comparison of protein coding genes reveals approximately 500 genes showing accelerated evolution on each of the gorilla, human and chimpanzee lineages, and evidence for parallel acceleration, particularly of genes involved in hearing. We also compare the western and eastern gorilla species, estimating an average sequence divergence time 1.75 million years ago, but with evidence for more recent genetic exchange and a population bottleneck in the eastern species. The use of the genome sequence in these and future analyses will promote a deeper understanding of great ape biology and evolution.

The degradome -the complete repertoire of proteases in an organism- is involved in multiple key biological and pathological processes. Previous studies in several organisms have yielded sets of curated protease sequences which may be used to characterize the degradome in a novel genome by similarity. Differences between degradomes can then be related to physiological traits of the species under study. Therefore, the sequencing of the zebra finch genome allows the comparison between the degradomes of mammals and birds and may help to understand the biological peculiarities of the zebra finch.

Although the mammalian X and Y chromosomes evolved from a single pair of autosomes, they are highly differentiated: the Y chromosome is dramatically smaller than the X and has lost most of its genes. The surviving genes are a specialized set with extraordinary evolutionary longevity. Most mammalian lineages have experienced delayed, or relatively recent, loss of at least one conserved Y-linked gene. An extreme example of this phenomenon is in the Japanese spiny rat, where the Y chromosome has disappeared altogether. In this species, many Y-linked genes were rescued by transposition to new genomic locations, but until our work presented here, this has been considered an isolated case.

Little is known about how human Y-Chromosome gene expression directly contributes to differences between XX (female) and XY (male) individuals in nonreproductive tissues. Here, we analyzed quantitative profiles of Y-Chromosome gene expression across 36 human tissues from hundreds of individuals. Although it is often said that Y-Chromosome genes are lowly expressed outside the testis, we report many instances of elevated Y-Chromosome gene expression in a nonreproductive tissue. A notable example is EIF1AY, which encodes eukaryotic translation initiation factor 1A Y-linked, together with its X-linked homolog EIF1AX Evolutionary loss of a Y-linked microRNA target site enabled up-regulation of EIF1AY, but not of EIF1AX, in the heart. Consequently, this essential translation initiation factor is nearly twice as abundant in male as in female heart tissue at the protein level. Divergence between the X and Y Chromosomes in regulatory sequence can therefore lead to tissue-specific Y-Chromosome-driven sex biases in expression of critical, dosage-sensitive regulatory genes.

Non-small-cell lung cancer (NSCLC) accounts for most cancer-related deaths worldwide. Liquid biopsy by a blood draw to detect circulating tumor cells (CTCs) is a tool for molecular profiling of cancer using single-cell and next-generation sequencing (NGS) technologies. The aim of the study was to identify somatic variants in single CTCs isolated from NSCLC patients by targeted NGS. Thirty-one subjects (20 NSCLC patients, 11 smokers without cancer) were enrolled for blood draws (7.5 mL). CTCs were identified by immunofluorescence, individually retrieved, and DNA-extracted. Targeted NGS was performed to detect somatic variants (single-nucleotide variants (SNVs) and insertions/deletions (Indels)) across 65 oncogenes and tumor suppressor genes. Cancer-associated variants were classified using OncoKB database. NSCLC patients had significantly higher CTC counts than control smokers (p = 0.0132; Mann-Whitney test). Analyzing 23 CTCs and 13 white blood cells across seven patients revealed a total of 644 somatic variants that occurred in all CTCs within the same subject, ranging from 1 to 137 per patient. The highest number of variants detected in ≥1 CTC within a patient was 441. A total of 18/65 (27.7%) genes were highly mutated. Mutations with oncogenic impact were identified in functional domains of seven oncogenes/tumor suppressor genes (NF1, PTCH1, TP53, SMARCB1, SMAD4, KRAS, and ERBB2). Single CTC-targeted NGS detects heterogeneous and shared mutational signatures within and between NSCLC patients. CTC single-cell genomics have potential for integration in NSCLC precision oncology.

The Amazon molly is a unique clonal fish species that originated from an interspecies hybrid between Poecilia species P. mexicana and P. latipinna It reproduces by gynogenesis, which eliminates paternal genomic contribution to offspring. An earlier study showed that Amazon molly shows biallelic expression for a large portion of the genome, leading to two main questions: (1) Are the allelic expression patterns from the initial hybridization event stabilized or changed during establishment of the asexual species and its further evolution? (2) Is allelic expression biased toward one parental allele a stochastic or adaptive process? To answer these questions, the allelic expression of P. formosa siblings was assessed to investigate intra- and inter-cohort allelic expression variability. For comparison, interspecies hybrids between P. mexicana and P. latipinna were produced in the laboratory to represent the P. formosa ancestor. We have identified inter-cohort and intra-cohort variation in parental allelic expression. The existence of inter-cohort divergence suggests functional P. formosa allelic expression patterns do not simply reflect the atavistic situation of the first interspecies hybrid but potentially result from long-term selection of transcriptional fitness. In addition, clonal fish show a transcriptional trend representing minimal intra-clonal variability in allelic expression patterns compared to the corresponding hybrids. The intra-clonal similarity in gene expression translates to sophisticated genetic functional regulation at the individuum level. These findings suggest the parental alleles inherited by P. formosa form tightly regulated genetic networks that lead to a stable transcriptomic landscape within clonal individuals.

Genetically resistant or susceptible chickens to Marek's disease (MD) have been widely used models to identify the molecular determinants of these phenotypes. However, these prior studies lacked the basic identification and understanding of immune cell types that could be translated toward improved MD control. To gain insights into specific immune cell types and their responses to Marek's disease virus (MDV) infection, we used single-cell RNA sequencing (scRNAseq) on splenic cells from MD resistant and susceptible birds. In total, 14,378 cells formed clusters that identified various immune cell types. Lymphocytes, specifically T cell subtypes, were the most abundant with significant proportional changes in some subtypes upon infection. The largest number of differentially expressed genes (DEG) response was seen in granulocytes, while macrophage DEGs differed in directionality by subtype and line. Among the most DEG in almost all immune cell types were granzyme and granulysin, both associated with cell-perforating processes. Protein interactive network analyses revealed multiple overlapping canonical pathways within both lymphoid and myeloid cell lineages. This initial estimation of the chicken immune cell type landscape and its accompanying response will greatly aid efforts in identifying specific cell types and improving our knowledge of host response to viral infection.

It is widely believed that at least three nonoverlapping regions of the human Y chromosome-AZFa, AZFb, and AZFc ("azoospermia factors" a, b, and c)-are essential for normal spermatogenesis. These intervals are defined by interstitial Y-chromosome deletions that impair or extinguish spermatogenesis. Deletion breakpoints, mechanisms, and lengths, as well as inventories of affected genes, have been elucidated for deletions of AZFa and of AZFc but not for deletions of AZFb or of AZFb plus AZFc. We studied three deletions of AZFb and eight deletions of AZFb plus AZFc, as assayed by the STSs defining these intervals. Guided by Y-chromosome sequence, we localized breakpoints precisely and were able to sequence nine of the deletion junctions. Homologous recombination can explain seven of these deletions but not the remaining two. This fact and our discovery of breakpoint hotspots suggest that factors in addition to homology underlie these deletions. The deletions previously thought to define AZFb were found to extend from palindrome P5 to the proximal arm of palindrome P1, 1.5 Mb within AZFc. Thus, they do not define a genomic region separate from AZFc. We also found that the deletions of AZFb plus AZFc, as assayed by standard STSs heretofore available, in fact extend from P5 to the distal arm of P1 and spare distal AZFc. Both classes of deletions are massive: P5/proximal-P1 deletions encompass up to 6.2 Mb and remove 32 genes and transcripts; P5/distal-P1 deletions encompass up to 7.7 Mb and remove 42 genes and transcripts. To our knowledge, these are the largest of all human interstitial deletions for which deletion junctions and complete intervening sequence are available. The restriction of the associated phenotype to spermatogenic failure indicates the remarkable functional specialization of the affected regions of the Y chromosome.

Comparative genomics continues illuminating amniote genome evolution, but for many lineages our understanding remains incomplete. Here, we refine the assembly (CPI 3.0.3 NCBI AHGY00000000.2) and develop a cytogenetic map of the painted turtle (Chrysemys picta-CPI) genome, the first in turtles and in vertebrates with temperature-dependent sex determination. A comparison of turtle genomes with those of chicken, selected nonavian reptiles, and human revealed shared and novel genomic features, such as numerous chromosomal rearrangements. The largest conserved syntenic blocks between birds and turtles exist in four macrochromosomes, whereas rearrangements were evident in these and other chromosomes, disproving that turtles and birds retain fully conserved macrochromosomes for greater than 300 Myr. C-banding revealed large heterochromatic blocks in the centromeric region of only few chromosomes. The nucleolar-organizing region (NOR) mapped to a single CPI microchromosome, whereas in some turtles and lizards the NOR maps to nonhomologous sex-chromosomes, thus revealing independent translocations of the NOR in various reptilian lineages. There was no evidence for recent chromosomal fusions as interstitial telomeric-DNA was absent. Some repeat elements (CR1-like, Gypsy) were enriched in the centromeres of five chromosomes, whereas others were widespread in the CPI genome. Bacterial artificial chromosome (BAC) clones were hybridized to 18 of the 25 CPI chromosomes and anchored to a G-banded ideogram. Several CPI sex-determining genes mapped to five chromosomes, and homology was detected between yet other CPI autosomes and the globally nonhomologous sex chromosomes of chicken, other turtles, and squamates, underscoring the independent evolution of vertebrate sex-determining mechanisms.

High-resolution genetic and physical maps are invaluable tools for building accurate genome assemblies, and interpreting results of genome-wide association studies (GWAS). Previous genetic and physical maps anchored good quality draft assemblies of the domestic cat genome, enabling the discovery of numerous genes underlying hereditary disease and phenotypes of interest to the biomedical science and breeding communities. However, these maps lacked sufficient marker density to order thousands of shorter scaffolds in earlier assemblies, which instead relied heavily on comparative mapping with related species. A high-resolution map would aid in validating and ordering chromosome scaffolds from existing and new genome assemblies. Here, we describe a high-resolution genetic linkage map of the domestic cat genome based on genotyping 453 domestic cats from several multi-generational pedigrees on the Illumina 63K SNP array. The final maps include 58,055 SNP markers placed relative to 6637 markers with unique positions, distributed across all autosomes and the X chromosome. Our final sex-averaged maps span a total autosomal length of 4464 cM, the longest described linkage map for any mammal, confirming length estimates from a previous microsatellite-based map. The linkage map was used to order and orient the scaffolds from a substantially more contiguous domestic cat genome assembly (Felis catus v8.0), which incorporated ∼20 × coverage of Illumina fragment reads. The new genome assembly shows substantial improvements in contiguity, with a nearly fourfold increase in N50 scaffold size to 18 Mb. We use this map to report probable structural errors in previous maps and assemblies, and to describe features of the recombination landscape, including a massive (∼50 Mb) recombination desert (of virtually zero recombination) on the X chromosome that parallels a similar desert on the porcine X chromosome in both size and physical location.

The human Y chromosome began to evolve from an autosome hundreds of millions of years ago, acquiring a sex-determining function and undergoing a series of inversions that suppressed crossing over with the X chromosome. Little is known about the recent evolution of the Y chromosome because only the human Y chromosome has been fully sequenced. Prevailing theories hold that Y chromosomes evolve by gene loss, the pace of which slows over time, eventually leading to a paucity of genes, and stasis. These theories have been buttressed by partial sequence data from newly emergent plant and animal Y chromosomes, but they have not been tested in older, highly evolved Y chromosomes such as that of humans. Here we finished sequencing of the male-specific region of the Y chromosome (MSY) in our closest living relative, the chimpanzee, achieving levels of accuracy and completion previously reached for the human MSY. By comparing the MSYs of the two species we show that they differ radically in sequence structure and gene content, indicating rapid evolution during the past 6 million years. The chimpanzee MSY contains twice as many massive palindromes as the human MSY, yet it has lost large fractions of the MSY protein-coding genes and gene families present in the last common ancestor. We suggest that the extraordinary divergence of the chimpanzee and human MSYs was driven by four synergistic factors: the prominent role of the MSY in sperm production, 'genetic hitchhiking' effects in the absence of meiotic crossing over, frequent ectopic recombination within the MSY, and species differences in mating behaviour. Although genetic decay may be the principal dynamic in the evolution of newly emergent Y chromosomes, wholesale renovation is the paramount theme in the continuing evolution of chimpanzee, human and perhaps other older MSYs.