Lack of Dosage Balance and Incomplete Dosage Compensation in the ZZ/ZW Gila Monster (Heloderma suspectum) Revealed by De Novo Genome Assembly

Reptiles exhibit a variety of modes of sex determination, including both temperature-dependent and genetic mechanisms. Among those species with genetic sex determination, sex chromosomes of varying heterogamety (XX/XY and ZZ/ZW) have been observed with different degrees of differentiation. Karyotype studies have demonstrated that Gila monsters (Heloderma suspectum) have ZZ/ZW sex determination and this system is likely homologous to the ZZ/ZW system in the Komodo dragon (Varanus komodoensis), but little else is known about their sex chromosomes. Here, we report the assembly and analysis of the Gila monster genome. We generated a de novo draft genome assembly for a male using 10X Genomics technology. We further generated and analyzed short-read whole genome sequencing and whole transcriptome sequencing data for three males and three females. By comparing female and male genomic data, we identified four putative Z chromosome scaffolds. These putative Z chromosome scaffolds are homologous to Z-linked scaffolds identified in the Komodo dragon. Further, by analyzing RNAseq data, we observed evidence of incomplete dosage compensation between the Gila monster Z chromosome and autosomes and a lack of balance in Z-linked expression between the sexes. In particular, we observe lower expression of the Z in females (ZW) than males (ZZ) on a global basis, though we find evidence suggesting local gene-by-gene compensation. This pattern has been observed in most other ZZ/ZW systems studied to date and may represent a general pattern for female heterogamety in vertebrates.

Genomic and demographic processes differentially influence genetic variation across the human X chromosome

Many forces influence genetic variation across the genome including mutation, recombination, selection, and demography. Increased mutation and recombination both lead to increases in genetic diversity in a region-specific manner, while complex demographic patterns shape patterns of diversity on a more global scale. While these processes act across the entire genome, the X chromosome is particularly interesting because it contains several distinct regions that are subject to different combinations and strengths of these forces: the pseudoautosomal regions (PARs) and the X-transposed region (XTR). The X chromosome thus can serve as a unique model for studying how genetic and demographic forces act in different contexts to shape patterns of observed variation. We therefore sought to explore diversity, divergence, and linkage disequilibrium in each region of the X chromosome using genomic data from 26 human populations. Across populations, we find that both diversity and substitution rate are consistently elevated in PAR1 and the XTR compared to the rest of the X chromosome. In contrast, linkage disequilibrium is lowest in PAR1, consistent with the high recombination rate in this region, and highest in the region of the X chromosome that does not recombine in males. However, linkage disequilibrium in the XTR is intermediate between PAR1 and the autosomes, and much lower than the non-recombining X. Finally, in addition to these global patterns, we also observed variation in ratios of X versus autosomal diversity consistent with population-specific evolutionary history as well. While our results were generally consistent with previous work, two unexpected observations emerged. First, our results suggest that the XTR does not behave like the rest of the recombining X and may need to be evaluated separately in future studies. Second, the different regions of the X chromosome appear to exhibit unique patterns of linked selection across different human populations. Together, our results highlight profound regional differences across the X chromosome, simultaneously making it an ideal system for exploring the action of evolutionary forces as well as necessitating its careful consideration and treatment in genomic analyses.

Sequence-based machine learning reveals 3D genome differences between bonobos and chimpanzees

Phenotypic divergence between closely related species, including bonobos and chimpanzees (genus Pan), is largely driven by variation in gene regulation. The 3D structure of the genome mediates gene expression; however, genome folding differences in Pan are not well understood. Here, we apply machine learning to predict genome-wide 3D genome contact maps from DNA sequence for 56 bonobos and chimpanzees, encompassing all five extant lineages. We use a pairwise approach to estimate 3D divergence between individuals from the resulting contact maps in 4,420 1 Mb genomic windows. While most pairs were similar, ∼17% were predicted to be substantially divergent in genome folding. The most dissimilar maps were largely driven by single individuals with rare variants that produce unique 3D genome folding in a region. We also identified 89 genomic windows where bonobo and chimpanzee contact maps substantially diverged, including several windows harboring genes associated with traits implicated in Pan phenotypic divergence. We used in silico mutagenesis to identify 51 3D-modifying variants in these bonobo-chimpanzee divergent windows, finding that 34 or 66.67% induce genome folding changes via CTCF binding motif disruption. Our results reveal 3D genome variation at the population-level and identify genomic regions where changes in 3D folding may contribute to phenotypic differences in our closest living relatives.

Incomplete dosage balance and dosage compensation in the ZZ/ZW Gila monster (Heloderma suspectum) revealed by de novo genome assembly

Reptiles exhibit a variety of modes of sex determination, including both temperature-dependent and genetic mechanisms. Among those species with genetic sex determination, sex chromosomes of varying heterogamety (XX/XY and ZZ/ZW) have been observed with different degrees of differentiation. Karyotype studies have demonstrated that Gila monsters (Heloderma suspectum) have ZZ/ZW sex determination and this system is likely homologous to the ZZ/ZW system in the Komodo dragon (Varanus komodoensis), but little else is known about their sex chromosomes. Here, we report the assembly and analysis of the Gila monster genome. We generated a de novo draft genome assembly for a male using 10X Genomics technology. We further generated and analyzed short-read whole genome sequencing and whole transcriptome sequencing data for three males and three females. By comparing female and male genomic data, we identified four putative Z-chromosome scaffolds. These putative Z-chromosome scaffolds are homologous to Z-linked scaffolds identified in the Komodo dragon. Further, by analyzing RNAseq data, we observed evidence of incomplete dosage compensation between the Gila monster Z chromosome and autosomes and a lack of balance in Z-linked expression between the sexes. In particular, we observe lower expression of the Z in females (ZW) than males (ZZ) on a global basis, though we find evidence suggesting local gene-by-gene compensation. This pattern has been observed in most other ZZ/ZW systems studied to date and may represent a general pattern for female heterogamety in vertebrates.

High-throughput phenotyping methods for quantifying hair fiber morphology

Quantifying the continuous variation in human scalp hair morphology is of interest to anthropologists, geneticists, dermatologists and forensic scientists, but existing methods for studying hair form are time-consuming and not widely used. Here, we present a high-throughput sample preparation protocol for the imaging of both longitudinal (curvature) and cross-sectional scalp hair morphology. Additionally, we describe and validate a new Python package designed to process longitudinal and cross-sectional hair images, segment them, and provide measurements of interest. Lastly, we apply our methods to an admixed African-European sample (n = 140), demonstrating the benefit of quantifying hair morphology over classification, and providing evidence that the relationship between cross-sectional morphology and curvature may be an artefact of population stratification rather than a causal link.

Comparative genomic analysis of sifakas (Propithecus) reveals selection for folivory and high heterozygosity despite endangered status

Sifakas (genus Propithecus) are critically endangered, large-bodied diurnal lemurs that eat leaf-based diets and show corresponding anatomical and microbial adaptations to folivory. We report on the genome assembly of Coquerel's sifaka (P. coquereli) and the resequenced genomes of Verreaux's (P. verreauxi), the golden-crowned (P. tattersalli), and the diademed (P. diadema) sifakas. We find high heterozygosity in all sifakas compared with other primates and endangered mammals. Demographic reconstructions nevertheless suggest declines in effective population size beginning before human arrival on Madagascar. Comparative genomic analyses indicate pervasive accelerated evolution in the ancestral sifaka lineage affecting genes in several complementary pathways relevant to folivory, including nutrient absorption and xenobiotic and fatty acid metabolism. Sifakas show convergent evolution at the level of the pathway, gene family, gene, and amino acid substitution with other folivores. Although sifakas have relatively generalized diets, the physiological challenges of habitual folivory likely led to strong selection.

Urine as a high-quality source of host genomic DNA from wild populations

The ability to generate genomic data from wild animal populations has the potential to give unprecedented insight into the population history and dynamics of species in their natural habitats. However, for many species, it is impossible legally, ethically or logistically to obtain tissue samples of quality sufficient for genomic analyses. In this study we evaluate the success of multiple sources of genetic material (faeces, urine, dentin and dental calculus) and several capture methods (shotgun, whole-genome, exome) in generating genome-scale data in wild eastern chimpanzees (Pan troglodytes schweinfurthii) from Gombe National Park, Tanzania. We found that urine harbours significantly more host DNA than other sources, leading to broader and deeper coverage across the genome. Urine also exhibited a lower rate of allelic dropout. We found exome sequencing to be far more successful than both shotgun sequencing and whole-genome capture at generating usable data from low-quality samples such as faeces and dental calculus. These results highlight urine as a promising and untapped source of DNA that can be noninvasively collected from wild populations of many species.

Comparing the selective landscape of TLR7 and TLR8 across primates reveals unique sites under positive selection in Alouatta

Among primates, susceptibility to yellow fever (YFV), a single-stranded (ss) RNA virus, ranges from complete resistance to high susceptibility. Howler monkeys (genus Alouatta) are the most susceptible to YFV. In order to identify Alouatta-specific genetic factors that may be responsible for their susceptibility, we collected skin samples from howler monkey museum specimens of the species A. caraya and A. guariba clamitans. We compared the rate of nonsynonymous to synonymous (dN/dS) changes of Toll-like receptor (TLR) 7 and TLR8, the two genes responsible for detecting all ssRNA viruses, across the Primate order. Overall, we found that the TLR7 gene is under stronger purifying selection in howler monkeys compared to other New World and Old World primates, but TLR8 is under the same selective pressure. When we evaluated dN/dS at each codon, we found six codons under positive selection in Alouatta TLR8 and two codons under positive selection in TLR7. The changes in TLR7 are unique to A. guariba clamitans and are found in functionally important regions likely to affect detection of ssRNA viruses by TLR7/TLR8, as well as downstream signaling. These amino acid differences in A. guariba clamitans may play a role in YFV susceptibility. These results have implications for identifying genetic factors affecting YFV susceptibility in primates.

Positively selected variants in functionally important regions of TLR7 in Alouatta guariba clamitans with yellow fever virus exposure in Northern Argentina

BACKGROUND

In 2007-2009, a major yellow fever virus (YFV) outbreak in Northern Argentina decimated the local howler monkey (Alouatta) population.

AIMS

To evaluate whether the surviving howler monkeys possess advantageous genetic variants inherited from monkeys alive prior to the YFV outbreak, we explored the relationship between Toll-like receptor (TLR) 7 and TLR8 gene variation and YFV susceptibility.

METHODS

We used samples from Alouatta individuals in Misiones, Argentina alive before the YFV outbreak, individuals that died during the outbreak, and individuals that survived the outbreak and are alive today. We measured genetic divergence between Alouatta YFV exposure groups and evaluated Alouatta-specific substitutions for functional consequences.

RESULTS

We did not find different allele frequencies in the post-YFV exposure Alouatta group compared to the pre-exposure group. We identified three nonsynonymous variants in TLR7 in Alouatta guariba clamitans. Two of these substitutions are under positive selection in functionally important regions of the gene.

DISCUSSION AND CONCLUSIONS

Our results did not indicate that surviving howler monkey spossess advantageous genetic variants at greater frequency than those alive before the YFV outbreak. However, the positively selected unique coding differences in A. guariba clamitans are in the region important in pathogen detection which may affect YFV resistance. Morework is necessary to fully explore this hypothesis.

Discovery of a New TLR Gene and Gene Expansion Event through Improved Desert Tortoise Genome Assembly with Chromosome-Scale Scaffolds

Toll-like receptors (TLRs) are a complex family of innate immune genes that are well characterized in mammals and birds but less well understood in nonavian sauropsids (reptiles). The advent of highly contiguous draft genomes of nonmodel organisms enables study of such gene families through analysis of synteny and sequence identity. Here, we analyze TLR genes from the genomes of 22 tetrapod species. Findings reveal a TLR8 gene expansion in crocodilians and turtles (TLR8B), and a second duplication (TLR8C) specifically within turtles, followed by pseudogenization of that gene in the nonfreshwater species (desert tortoise and green sea turtle). Additionally, the Mojave desert tortoise (Gopherus agassizii) has a stop codon in TLR8B (TLR8-1) that is polymorphic among conspecifics. Revised orthology further reveals a new TLR homolog, TLR21-like, which is exclusive to lizards, snakes, turtles, and crocodilians. These analyses were made possible by a new draft genome assembly of the desert tortoise (gopAga2.0), which used chromatin-based assembly to yield draft chromosomal scaffolds (L50 = 26 scaffolds, N50 = 28.36 Mb, longest scaffold = 107 Mb) and an enhanced de novo genome annotation with 25,469 genes. Our three-step approach to orthology curation and comparative analysis of TLR genes shows what new insights are possible using genome assemblies with chromosome-scale scaffolds that permit integration of synteny conservation data.

Virus Discovery in Desert Tortoise Fecal Samples: Novel Circular Single-Stranded DNA Viruses

The Sonoran Desert tortoise Gopherus morafkai is adapted to the desert, and plays an important ecological role in this environment. There is limited information on the viral diversity associated with tortoises (family Testudinidae), and to date no DNA virus has been identified associated with these animals. This study aimed to assess the diversity of DNA viruses associated with the Sonoran Desert tortoise by sampling their fecal matter. A viral metagenomics approach was used to identify the DNA viruses in fecal samples from wild Sonoran Desert tortoises in Arizona, USA. In total, 156 novel single-stranded DNA viruses were identified from 40 fecal samples. Those belonged to two known viral families, the Genomoviridae (n = 27) and Microviridae (n = 119). In addition, 10 genomes were recovered that belong to the unclassified group of circular-replication associated protein encoding single-stranded (CRESS) DNA virus and five circular molecules encoding viral-like proteins.

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.

Evolution of Dosage Compensation in Anolis carolinensis, a Reptile with XX/XY Chromosomal Sex Determination

In species with highly heteromorphic sex chromosomes, the degradation of one of the sex chromosomes will result in unequal gene expression between the sexes (e.g. between XX females and XY males) and between the sex chromosomes and the autosomes. Dosage compensation is a process whereby genes on the sex chromosomes achieve equal gene expression. We compared genome-wide levels of transcription between males and females, and between the X chromosome and the autosomes in the green anole, Anolis carolinensis. We present evidence for dosage compensation between the sexes, and between the sex chromosomes and the autosomes. When dividing the X chromosome into regions based on linkage groups, we discovered that genes in the first reported X-linked region, anole linkage group b (LGb), exhibit complete dosage compensation, although the rest of the X-linked genes exhibit incomplete dosage compensation. Our data further suggest that the mechanism of this dosage compensation is upregulation of the X chromosome in males. We report that approximately 10% of coding genes, most of which are on the autosomes, are differentially expressed between males and females. In addition, genes on the X chromosome exhibited higher ratios of nonsynonymous to synonymous substitution than autosomal genes, consistent with the fast-X effect. Our results from the green anole add an additional observation of dosage compensation in a species with XX/XY sex determination.

Genomic signatures of sex-biased demography: progress and prospects

Sex-biased demographic events have played a crucial role in shaping human history. Many of these processes affect genetic variation and can therefore leave detectable signatures in the genome because autosomal, X-linked, Y-linked, and mitochondrial DNA inheritance differ between sexes. Here, we discuss how sex-biased processes shape patterns of genetic diversity across the genome, review recent genomic evidence for sex-biased demography in modern human populations, and suggest directions for future research.

The origin of snakes: revealing the ecology, behavior, and evolutionary history of early snakes using genomics, phenomics, and the fossil record

Background

The highly derived morphology and astounding diversity of snakes has long inspired debate regarding the ecological and evolutionary origin of both the snake total-group (Pan-Serpentes) and crown snakes (Serpentes). Although speculation abounds on the ecology, behavior, and provenance of the earliest snakes, a rigorous, clade-wide analysis of snake origins has yet to be attempted, in part due to a dearth of adequate paleontological data on early stem snakes. Here, we present the first comprehensive analytical reconstruction of the ancestor of crown snakes and the ancestor of the snake total-group, as inferred using multiple methods of ancestral state reconstruction. We use a combined-data approach that includes new information from the fossil record on extinct crown snakes, new data on the anatomy of the stem snakes Najash rionegrina, Dinilysia patagonica, and Coniophis precedens, and a deeper understanding of the distribution of phenotypic apomorphies among the major clades of fossil and Recent snakes. Additionally, we infer time-calibrated phylogenies using both new ‘tip-dating’ and traditional node-based approaches, providing new insights on temporal patterns in the early evolutionary history of snakes.

Results

Comprehensive ancestral state reconstructions reveal that both the ancestor of crown snakes and the ancestor of total-group snakes were nocturnal, widely foraging, non-constricting stealth hunters. They likely consumed soft-bodied vertebrate and invertebrate prey that was subequal to head size, and occupied terrestrial settings in warm, well-watered, and well-vegetated environments. The snake total-group – approximated by the Coniophis node – is inferred to have originated on land during the middle Early Cretaceous (~128.5 Ma), with the crown-group following about 20 million years later, during the Albian stage. Our inferred divergence dates provide strong evidence for a major radiation of henophidian snake diversity in the wake of the Cretaceous-Paleogene (K-Pg) mass extinction, clarifying the pattern and timing of the extant snake radiation. Although the snake crown-group most likely arose on the supercontinent of Gondwana, our results suggest the possibility that the snake total-group originated on Laurasia.

Conclusions

Our study provides new insights into when, where, and how snakes originated, and presents the most complete picture of the early evolution of snakes to date. More broadly, we demonstrate the striking influence of including fossils and phenotypic data in combined analyses aimed at both phylogenetic topology inference and ancestral state reconstruction.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0358-5) contains supplementary material, which is available to authorized users.

Selective insectivory at Toro-Semliki, Uganda: comparative analyses suggest no 'savanna' chimpanzee pattern

Chimpanzee (Pan troglodytes) insectivory across Africa is ubiquitous. Insects provide a significant nutritional payoff and may be important for chimpanzees in dry, open habitats with narrow diets. We tested this hypothesis at Semliki, Uganda, a long-term dry study site. We evaluated prospects for insectivory by measuring insect abundance along de novo transects and trails, monitoring social insect colonies, and surveying available raw materials for elementary technology. We determined the frequency and nature of insectivory through behavioral observation and fecal analysis. We then compared our results with those from 15 other long-term chimpanzee study sites using a cluster analysis. We found that Semliki chimpanzees are one of the most insectivorous populations studied to date in terms of frequency of consumption, but they are very selective in their insectivory, regularly consuming only weaver ants (Oecophylla longinoda) and honey and bees from hives of Apis mellifera. This selectivity obtains despite having a full range of typical prey species available in harvestable quantities. We suggest that Semliki chimpanzees may face ecological time constraints and therefore bias their predation toward prey taxa that can be quickly consumed. Geographical proximity correlated with the results of the cluster analysis, while rainfall, a relatively gross measure of environment, did not. Because broad taxonomic groups of insects were used in analyses, prey availability was unlikely to have a strong effect on this pattern. Instead, we suggest that transmission of cultural knowledge may play a role in determining chimpanzee prey selection across Africa. Further study is needed to test these hypotheses.

The apolipoprotein E (APOE) gene appears functionally monomorphic in chimpanzees (Pan troglodytes)

Background

The human apolipoprotein E (APOE) gene is polymorphic, with three primary alleles (E2, E3, E4) that differ at two key non-synonymous sites. These alleles are functionally different in how they bind to lipoproteins, and this genetic variation is associated with phenotypic variation for several medical traits, including cholesterol levels, cardiovascular health, Alzheimer’s disease risk, and longevity. The relative frequencies of these alleles vary across human populations, and the evolution and maintenance of this diversity is much debated. Previous studies comparing human and chimpanzee APOE sequences found that the chimpanzee sequence is most similar to the human E4 allele, although the resulting chimpanzee protein might function like the protein coded for by the human E3 allele. However, these studies have used sequence data from a single chimpanzee and do not consider whether chimpanzees, like humans, show intra-specific and subspecific variation at this locus.

Methodology and Principal Findings

To examine potential intraspecific variation, we sequenced the APOE gene of 32 chimpanzees. This sample included 20 captive individuals representing the western subspecies (P. troglodytes verus) and 12 wild individuals representing the eastern subspecies (P. t. schweinfurthii). Variation in our resulting sequences was limited to one non-coding, intronic SNP, which showed fixed differences between the two subspecies. We also compared APOE sequences for all available ape genera and fossil hominins. The bonobo APOE protein is identical to that of the chimpanzee, and the Denisovan APOE exhibits all four human-specific, non-synonymous changes and appears functionally similar to the human E4 allele.

Conclusions

We found no coding variation within and between chimpanzee populations, suggesting that the maintenance of functionally diverse APOE polymorphisms is a unique feature of human evolution.

Reversal of midazolam-induced laryngospasm with flumazenil

Midazolam is a commonly used benzodiazepine that is ideal for conscious sedation during a variety of procedures. Laryngospasm is listed by the manufacturer as a rare side effect of midazolam. Flumazenil is a competitive inhibitor of the benzodiazepines that is available for reversal at the end of such a procedure or when too much sedation has been achieved. We present a case of a 61-year-old man who was to undergo conscious sedation before cardioversion. Shortly after receiving midazolam, he developed laryngospasm and resultant respiratory distress; the laryngospasm was reversed by a small dose of flumazenil. The symptoms recurred after approximately 25 minutes and were again reversed with flumazenil. Neither midazolam-induced laryngospasm nor its reversal with flumazenil have previously been reported. A brief discussion regarding laryngospasm and the use of flumazenil are included.