Sources of gene expression variation in a globally diverse human cohort.

Genetic variation influencing gene expression and splicing is a key source of phenotypic diversity. Though invaluable, studies investigating these links in humans have been strongly biased toward participants of European ancestries, diminishing generalizability and hindering evolutionary research. To address these limitations, we developed MAGE, an open-access RNA-seq data set of lymphoblastoid cell lines from 731 individuals from the 1000 Genomes Project spread across 5 continental groups and 26 populations. Most variation in gene expression (92%) and splicing (95%) was distributed within versus between populations, mirroring variation in DNA sequence. We mapped associations between genetic variants and expression and splicing of nearby genes (cis-eQTLs and cis-sQTLs, respective), identifying >15,000 putatively causal eQTLs and >16,000 putatively causal sQTLs that are enriched for relevant epigenomic signatures. These include 1310 eQTLs and 1657 sQTLs that are largely private to previously underrepresented populations. Our data further indicate that the magnitude and direction of causal eQTL effects are highly consistent across populations and that apparent "population-specific" effects observed in previous studies were largely driven by low resolution or additional independent eQTLs of the same genes that were not detected. Together, our study expands understanding of gene expression diversity across human populations and provides an inclusive resource for studying the evolution and function of human genomes.

Familial Clonal Hematopoiesis in a Long Telomere Syndrome.

BACKGROUND: Telomere shortening is a well-characterized cellular aging mechanism, and short telomere syndromes cause age-related disease. However, whether long telomere length is advantageous is poorly understood. METHODS: We examined the clinical and molecular features of aging and cancer in persons carrying heterozygous loss-of-function mutations in the telomere-related gene POT1 and noncarrier relatives. RESULTS: A total of 17 POT1 mutation carriers and 21 noncarrier relatives were initially included in the study, and a validation cohort of 6 additional mutation carriers was subsequently recruited. A majority of the POT1 mutation carriers with telomere length evaluated (9 of 13) had long telomeres (>99th percentile). POT1 mutation carriers had a range of benign and malignant neoplasms involving epithelial, mesenchymal, and neuronal tissues in addition to B- and T-cell lymphoma and myeloid cancers. Five of 18 POT1 mutation carriers (28%) had T-cell clonality, and 8 of 12 (67%) had clonal hematopoiesis of indeterminate potential. A predisposition to clonal hematopoiesis had an autosomal dominant pattern of inheritance, as well as penetrance that increased with age; somatic DNMT3A and JAK2 hotspot mutations were common. These and other somatic driver mutations probably arose in the first decades of life, and their lineages secondarily accumulated a higher mutation burden characterized by a clocklike signature. Successive generations showed genetic anticipation (i.e., an increasingly early onset of disease). In contrast to noncarrier relatives, who had the typical telomere shortening with age, POT1 mutation carriers maintained telomere length over the course of 2 years. CONCLUSIONS: POT1 mutations associated with long telomere length conferred a predisposition to a familial clonal hematopoiesis syndrome that was associated with a range of benign and malignant solid neoplasms. The risk of these phenotypes was mediated by extended cellular longevity and by the capacity to maintain telomeres over time. (Funded by the National Institutes of Health and others.).

Unexpected Distribution of Chitin and Chitin Synthase across Soft-Bodied Cnidarians.

Cnidarians are commonly recognized as sea jellies, corals, or complex colonies such as the Portuguese man-of-war. While some cnidarians possess rigid internal calcareous skeletons (e.g., corals), many are soft-bodied. Intriguingly, genes coding for the chitin-biosynthetic enzyme, chitin synthase (CHS), were recently identified in the model anemone Nematostella vectensis, a species lacking hard structures. Here we report the prevalence and diversity of CHS across Cnidaria and show that cnidarian chitin synthase genes display diverse protein domain organizations. We found that CHS is expressed in cnidarian species and/or developmental stages with no reported chitinous or rigid morphological structures. Chitin affinity histochemistry indicates that chitin is present in soft tissues of some scyphozoan and hydrozoan medusae. To further elucidate the biology of chitin in cnidarian soft tissues, we focused on CHS expression in N. vectensis. Spatial expression data show that three CHS orthologs are differentially expressed in Nematostella embryos and larvae during development, suggesting that chitin has an integral role in the biology of this species. Understanding how a non-bilaterian lineage such as Cnidaria employs chitin may provide new insight into hitherto unknown functions of polysaccharides in animals, as well as their role in the evolution of biological novelty.

Induced Immune Reaction in the Acorn Worm, Saccoglossus kowalevskii, Informs the Evolution of Antiviral Immunity.

Evolutionary perspectives on the deployment of immune factors following infection have been shaped by studies on a limited number of biomedical model systems with a heavy emphasis on vertebrate species. Although their contributions to contemporary immunology cannot be understated, a broader phylogenetic perspective is needed to understand the evolution of immune systems across Metazoa. In our study, we leverage differential gene expression analyses to identify genes implicated in the antiviral immune response of the acorn worm hemichordate, Saccoglossus kowalevskii, and place them in the context of immunity evolution within deuterostomes-the animal clade composed of chordates, hemichordates, and echinoderms. Following acute exposure to the synthetic viral double-stranded RNA analog, poly(I:C), we show that S. kowalevskii responds by regulating the transcription of genes associated with canonical innate immunity signaling pathways (e.g., nuclear factor κB and interferon regulatory factor signaling) and metabolic processes (e.g., lipid metabolism), as well as many genes without clear evidence of orthology with those of model species. Aggregated across all experimental time point contrasts, we identify 423 genes that are differentially expressed in response to poly(I:C). We also identify 147 genes with altered temporal patterns of expression in response to immune challenge. By characterizing the molecular toolkit involved in hemichordate antiviral immunity, our findings provide vital evolutionary context for understanding the origins of immune systems within Deuterostomia.

TIAMMAt: Leveraging Biodiversity to Revise Protein Domain Models, Evidence from Innate Immunity.

Sequence annotation is fundamental for studying the evolution of protein families, particularly when working with nonmodel species. Given the rapid, ever-increasing number of species receiving high-quality genome sequencing, accurate domain modeling that is representative of species diversity is crucial for understanding protein family sequence evolution and their inferred function(s). Here, we describe a bioinformatic tool called Taxon-Informed Adjustment of Markov Model Attributes (TIAMMAt) which revises domain profile hidden Markov models (HMMs) by incorporating homologous domain sequences from underrepresented and nonmodel species. Using innate immunity pathways as a case study, we show that revising profile HMM parameters to directly account for variation in homologs among underrepresented species provides valuable insight into the evolution of protein families. Following adjustment by TIAMMAt, domain profile HMMs exhibit changes in their per-site amino acid state emission probabilities and insertion/deletion probabilities while maintaining the overall structure of the consensus sequence. Our results show that domain revision can heavily impact evolutionary interpretations for some families (i.e., NLR's NACHT domain), whereas impact on other domains (e.g., rel homology domain and interferon regulatory factor domains) is minimal due to high levels of sequence conservation across the sampled phylogenetic depth (i.e., Metazoa). Importantly, TIAMMAt revises target domain models to reflect homologous sequence variation using the taxonomic distribution under consideration by the user. TIAMMAt's flexibility to revise any subset of the Pfam database using a user-defined taxonomic pool will make it a valuable tool for future protein evolution studies, particularly when incorporating (or focusing) on nonmodel species.

Genomic adaptations to chemosymbiosis in the deep-sea seep-dwelling tubeworm Lamellibrachia luymesi.

BACKGROUND: Symbiotic relationships between microbes and their hosts are widespread and diverse, often providing protection or nutrients, and may be either obligate or facultative. However, the genetic mechanisms allowing organisms to maintain host-symbiont associations at the molecular level are still mostly unknown, and in the case of bacterial-animal associations, most genetic studies have focused on adaptations and mechanisms of the bacterial partner. The gutless tubeworms (Siboglinidae, Annelida) are obligate hosts of chemoautotrophic endosymbionts (except for Osedax which houses heterotrophic Oceanospirillales), which rely on the sulfide-oxidizing symbionts for nutrition and growth. Whereas several siboglinid endosymbiont genomes have been characterized, genomes of hosts and their adaptations to this symbiosis remain unexplored. RESULTS: Here, we present and characterize adaptations of the cold seep-dwelling tubeworm Lamellibrachia luymesi, one of the longest-lived solitary invertebrates. We sequenced the worm's ~ 688-Mb haploid genome with an overall completeness of ~ 95% and discovered that L. luymesi lacks many genes essential in amino acid biosynthesis, obligating them to products provided by symbionts. Interestingly, the host is known to carry hydrogen sulfide to thiotrophic endosymbionts using hemoglobin. We also found an expansion of hemoglobin B1 genes, many of which possess a free cysteine residue which is hypothesized to function in sulfide binding. Contrary to previous analyses, the sulfide binding mediated by zinc ions is not conserved across tubeworms. Thus, the sulfide-binding mechanisms in sibgolinids need to be further explored, and B1 globins might play a more important role than previously thought. Our comparative analyses also suggest the Toll-like receptor pathway may be essential for tolerance/sensitivity to symbionts and pathogens. Several genes related to the worm's unique life history which are known to play important roles in apoptosis, cell proliferation, and aging were also identified. Last, molecular clock analyses based on phylogenomic data suggest modern siboglinid diversity originated in 267 mya (± 70 my) support previous hypotheses indicating a Late Mesozoic or Cenozoic origins of approximately 50-126 mya for vestimentiferans. CONCLUSIONS: Here, we elucidate several specific adaptations along various molecular pathways that link phenome to genome to improve understanding of holobiont evolution. Our findings of adaptation in genomic mechanisms to reducing environments likely extend to other chemosynthetic symbiotic systems.

Mitogenomics Reveals a Novel Genetic Code in Hemichordata.

The diverse array of codon reassignments demonstrate that the genetic code is not universal in nature. Exploring mechanisms underlying codon reassignment is critical for understanding the evolution of the genetic code during translation. Hemichordata, comprising worm-like Enteropneusta and colonial filter-feeding Pterobranchia, is the sister taxon of echinoderms and is more distantly related to chordates. However, only a few hemichordate mitochondrial genomes have been sequenced, hindering our understanding of mitochondrial genome evolution within Deuterostomia. In this study, we sequenced four mitochondrial genomes and two transcriptomes, including representatives of both major hemichordate lineages and analyzed together with public available data. Contrary to the current understanding of the mitochondrial genetic code in hemichordates, our comparative analyses suggest that UAA encodes Tyr instead of a "Stop" codon in the pterobranch lineage Cephalodiscidae. We also predict that AAA encodes Lys in pterobranch and enteropneust mitochondrial genomes, contradicting the previous assumption that hemichordates share the same genetic code with echinoderms for which AAA encodes Asn. Thus, we propose a new mitochondrial genetic code for Cephalodiscus and a revised code for enteropneusts. Moreover, our phylogenetic analyses are largely consistent with previous phylogenomic studies. The only exception is the phylogenetic position of the enteropneust Stereobalanus, whose placement as sister to all other described enteropneusts. With broader taxonomic sampling, we provide evidence that evolution of mitochondrial gene order and genetic codes in Hemichordata are more dynamic than previously thought and these findings provide insights into mitochondrial genome evolution within this clade.

Phylogenomics offers resolution of major tunicate relationships.

Tunicata, a diverse clade of approximately 3000 described species of marine, filter-feeding chordates, is of great interest to researchers because tunicates are the closest living relatives of vertebrates and they facilitate comparative studies of our own biology. The group also includes numerous invasive species that cause considerable economic damage and some species of tunicates are edible. Despite their diversity and importance, relationships among major lineages of Tunicata are not completely resolved. Here, we supplemented public data with transcriptomes from seven species spanning the diversity of Tunicata and conducted phylogenomic analyses on data sets of up to 798 genes. Sensitivity analyses were employed to examine the influences of reducing compositional heterogeneity and branch-length heterogeneity. All analyses maximally supported a monophyletic Tunicata within Olfactores (Vertebrata + Tunicata). Within Tunicata, all analyses recovered Appendicularia sister to the rest of Tunicata and confirmed (with maximal support) that Thaliacea is nested within Ascidiacea. Stolidobranchia is the sister taxon to all other tunicates except Appendicularia. In most analyses, phlebobranch tunicates were recovered paraphyletic with respect to Aplousobranchia. Support for this topology varied but was strong in some cases. However, when only the 50 best genes based on compositional heterogeneity were analysed, we recovered Phlebobranchia and Aplousobranchia reciprocally monophyletic with strong support, consistent with most traditional morphology-based hypotheses. Examination of internode certainty also cast doubt on results of phlebobranch paraphyly, which may be due to limited taxon sampling. Taken together, these results provide a higher-level phylogenetic framework for our closest living invertebrate relatives.

Toll-like receptor pathway evolution in deuterostomes.

Animals have evolved an array of pattern-recognition receptor families essential for recognizing conserved molecular motifs characteristic of pathogenic microbes. One such family is the Toll-like receptors (TLRs). On pathogen binding, TLRs initiate specialized cytokine signaling catered to the class of invading pathogen. This signaling is pivotal for activating adaptive immunity in vertebrates, suggesting a close evolutionary relationship between innate and adaptive immune systems. Despite significant advances toward understanding TLR-facilitated immunity in vertebrates, knowledge of TLR pathway evolution in other deuterostomes is limited. By analyzing genomes and transcriptomes across 37 deuterostome taxa, we shed light on the evolution and diversity of TLR pathway signaling elements. Here, we show that the deuterostome ancestor possessed a molecular toolkit homologous to that which drives canonical MYD88-dependent TLR signaling in contemporary mammalian lineages. We also provide evidence that TLR3-facilitated antiviral signaling predates the origin of its TCAM1 dependence recognized in the vertebrates. SARM1, a negative regulator of TCAM1-dependent pathways in vertebrates, was also found to be present across all major deuterostome lineages despite the apparent absence of TCAM1 in invertebrate deuterostomes. Whether the presence of SARM1 is the result of its role in immunity regulation, neuron physiology, or a function of both is unclear. Additionally, Bayesian phylogenetic analyses corroborate several lineage-specific TLR gene expansions in urchins and cephalochordates. Importantly, our results underscore the need to sample across taxonomic groups to understand evolutionary patterns of the innate immunity foundation on which complex immunological novelties arose.

The Global Diversity of Hemichordata.

Phylum Hemichordata, composed of worm-like Enteropneusta and colonial Pterobranchia, has been reported to only contain about 100 species. However, recent studies of hemichordate phylogeny and taxonomy suggest the species number has been largely underestimated. One issue is that species must be described by experts, and historically few taxonomists have studied this group of marine invertebrates. Despite this previous lack of coverage, interest in hemichordates has piqued in the past couple of decades, as they are critical to understanding the evolution of chordates-as acorn worms likely resemble the deuterostome ancestor more closely than any other extant animal. This review provides an overview of our current knowledge of hemichordates, focusing specifically on their global biodiversity, geographic distribution, and taxonomy. Using information available in the World Register of Marine Species and published literature, we assembled a list of 130 described, extant species. The majority (83%) of these species are enteropneusts, and more taxonomic descriptions are forthcoming. Ptychoderidae contained the greatest number of species (41 species), closely followed by Harrimaniidae (40 species), of the recognized hemichordate families. Hemichordates are found throughout the world's oceans, with the highest reported numbers by regions with marine labs and diligent taxonomic efforts (e.g. North Pacific and North Atlantic). Pterobranchs are abundant in Antarctica, but have also been found at lower latitudes. We consider this a baseline report and expect new species of Hemichordata will continue to be discovered and described as new marine habitats are characterized and explored.