Different rates of substitution may produce different phylogenies of the eutherian mammals

Summary In an attempt to resolve some points of branching order in the phylogeny of the eutherian mammals, a phylogenetic analysis of 26 nuclear and 6 mitochondrial genes was undertaken using a maximum likelihood method on a constant rate stochastic model of molecular evolution. Seventeen of the nuclear genes gave a primates/artiodactyls grouping highest support whereas three of the mitochondrial genes found a rodents/artiodactyls grouping to be best supported. The primates/rodents grouping was never the best supported. On the assumption that rodents are indeed an outgroup to primates and artiodactyls and that the latter taxa diverged 70 million years ago, an estimation was made, for each gene, of the time of divergence of the rodent lineage. In most cases such estimates were beyond the limits set by present interpretations of the paleontological record as were many estimates of the divergence time of mouse and rat. These results suggest that, although there is locus variation, the divergent position of the rodent lineage may be an artifact of an elevated rate of nucleotide substitution in this order.     

Protein design of IgG/TCR chimeras for the co-expression of Fab-like moieties within bispecific antibodies

Immunoglobulins and T cell receptors (TCRs) share common sequences and structures. With the goal of creating novel bispecific antibodies (BsAbs), we generated chimeric molecules, denoted IgG_TCRs, where the Fv regions of several antibodies were fused to the constant domains of the α/β TCR. Replacing C_H1 with Cα and C_L with Cβ, respectively, was essential for achieving at least partial heavy chain/light chain assembly. Further optimization of the linker regions between the variable and constant domains, as well as replacement of the large FG loop of Cβ with a canonical β-turn, was necessary to consistently obtain full heavy chain/light chain assembly. The optimized IgG_TCR molecules were evaluated biophysically and shown to maintain the binding properties of their parental antibodies. A few BsAbs were generated by co-expressing native Fabs and IgG_TCR Fabs within the same molecular construct. We demonstrate that the IgG_TCR designs steered each of the light chains within the constructs to specifically pair with their cognate heavy chain counterparts. We did find that even with complete constant domain specificity between the C_H1/C_L and Cα/Cβ domains of the Fabs, strong variable domain interactions can dominate the pairing specificity and induce some mispairing. Overall, the IgG_TCR designs described here are a first step toward the generation of novel BsAbs that may be directed toward the treatment of multi-faceted and complex diseases.     

The SKINT1-Like Gene Is Inactivated in Hominoids But Not in All Primate Species: Implications for the Origin of Dendritic Epidermal T Cells

Dendritic epidermal T cells, which express an invariant Vγ5Vδ1 T-cell receptor and account for 95% of all resident T cells in the mouse epidermis, play a critical role in skin immune surveillance. These γδ T cells are generated by positive selection in the fetal thymus, after which they migrate to the skin. The development of dendritic epidermal T cells is critically dependent on the Skint1 gene expressed specifically in keratinocytes and thymic epithelial cells, suggesting an indispensable role for Skint1 in the selection machinery for specific intraepithelial lymphocytes. Phylogenetically, rodents have functional SKINT1 molecules, but humans and chimpanzees have a SKINT1-like (SKINT1L) gene with multiple inactivating mutations. In the present study, we analyzed SKINT1L sequences in representative primate species and found that all hominoid species have a common inactivating mutation, but that Old World monkeys such as olive baboons, green monkeys, cynomolgus macaques and rhesus macaques have apparently functional SKINT1L sequences, indicating that SKINT1L was inactivated in a common ancestor of hominoids. Interestingly, the epidermis of cynomolgus macaques contained a population of dendritic-shaped γδ T cells expressing a semi-invariant Vγ10/Vδ1 T-cell receptor. However, this population of macaque T cells differed from rodent dendritic epidermal T cells in that their Vγ10/Vδ1 T-cell receptors displayed junctional diversity and expression of Vγ10 was not epidermis-specific. Therefore, macaques do not appear to have rodent-type dendritic epidermal T cells despite having apparently functional SKINT1L. Comprehensive bioinformatics analysis indicates that SKINT1L emerged in an ancestor of placental mammals but was inactivated or lost multiple times in mammalian evolution and that Skint1 arose by gene duplication in a rodent lineage, suggesting that authentic dendritic epidermal T cells are presumably unique to rodents.     

Two Functional Epithelial Sodium Channel Isoforms Are Present in Rodents despite Pronounced Evolutionary Pseudogenization and Exon Fusion

The epithelial sodium channel (ENaC) plays a key role in salt and water homeostasis in tetrapod vertebrates. There are four ENaC subunits (α, β, γ, δ), forming heterotrimeric αβγ- or δβγ-ENaCs. Although the physiology of αβγ-ENaC is well understood, for decades the field has stalled with respect to δβγ-ENaC due to the lack of mammalian model organisms. The SCNN1D gene coding for δ-ENaC was previously believed to be absent in rodents, hindering studies using standard laboratory animals. We analyzed all currently available rodent genomes and discovered that SCNN1D is present in rodents but was independently lost in five rodent lineages, including the Muridae (mice and rats). The independent loss of SCNN1D in rodent lineages may be constrained by phylogeny and taxon-specific adaptation to dry habitats, however habitat aridity does not provide a selection pressure for maintenance of SCNN1D across Rodentia. A fusion of two exons coding for a structurally flexible region in the extracellular domain of δ-ENaC appeared in the Hystricognathi (a group that includes guinea pigs). This conserved pattern evolved at least 41 Ma and represents a new autapomorphic feature for this clade. Exon fusion does not impair functionality of guinea pig (Cavia porcellus) δβγ-ENaC expressed in Xenopus oocytes. Electrophysiological characterization at the whole-cell and single-channel level revealed conserved biophysical features and mechanisms controlling guinea pig αβγ- and δβγ-ENaC function as compared with human orthologs. Guinea pigs therefore represent commercially available mammalian model animals that will help shed light on the physiological function of δ-ENaC.     

Glycogen Synthase Kinase-3β (GSK-3β) Inhibition Enhances Dendritic Cell-based Cancer Vaccine Potency via Suppression of Interferon-γ-induced Indoleamine 2,3-Dioxygenase Expression

Indoleamine 2,3-dioxygenase (IDO) functions as a crucial mediator of tumor-mediated immune tolerance by causing T-cell suppression via tryptophan starvation in a tumor environment. Glycogen synthase kinase-3β (GSK-3β) is also involved in immune and anti-tumor responses. However, the relativity of these proteins has not been as well defined. Here, we found that GSK-3β-dependent IDO expression in the dendritic cell (DC) plays a role in anti-tumor activity via the regulation of CD8⁺ T-cell polarization and cytotoxic T lymphocyte activity. By the inhibition of GSK-3β, attenuated IDO expression and impaired JAK1/2-Stat signaling crucial for IDO expression were observed. Protein kinase Cδ (PKCδ) activity and the interaction between JAK1/2 and Stat3, which are important for IDO expression, were also reduced by GSK-3β inhibition. CD8⁺ T-cell proliferation mediated by OVA-pulsed DC was blocked by interferon (IFN)-γ-induced IDO expression via GSK-3β activity. Specific cytotoxic T lymphocyte activity mediated by OVA-pulsed DC against OVA-expressing EG7 thymoma cells but not OVA-nonexpressing EL4 thymoma cells was also attenuated by the expressed IDO via IFN-γ-induced activation of GSK-3β. Furthermore, tumor growth that was suppressed with OVA-pulsed DC vaccination was restored by IDO-expressing DC via IFN-γ-induced activation of GSK-3β in an OVA-expressing murine EG7 thymoma model. Taken together, DC-based immune response mediated by interferon-γ-induced IDO expression via GSK-3β activity not only regulates CD8⁺ T-cell proliferation and cytotoxic T lymphocyte activity but also modulates OVA-pulsed DC vaccination against EG7 thymoma.     

A New Class of Wheat Gliadin Genes and Proteins

The utility of mining DNA sequence data to understand the structure and expression of cereal prolamin genes is demonstrated by the identification of a new class of wheat prolamins. This previously unrecognized wheat prolamin class, given the name δ-gliadins, is the most direct ortholog of barley γ3-hordeins. Phylogenetic analysis shows that the orthologous δ-gliadins and γ3-hordeins form a distinct prolamin branch that existed separate from the γ-gliadins and γ-hordeins in an ancestral Triticeae prior to the branching of wheat and barley. The expressed δ-gliadins are encoded by a single gene in each of the hexaploid wheat genomes. This single δ-gliadin/γ3-hordein ortholog may be a general feature of the Triticeae tribe since examination of ESTs from three barley cultivars also confirms a single γ3-hordein gene. Analysis of ESTs and cDNAs shows that the genes are expressed in at least five hexaploid wheat cultivars in addition to diploids Triticum monococcum and Aegilops tauschii. The latter two sequences also allow assignment of the δ-gliadin genes to the A and D genomes, respectively, with the third sequence type assumed to be from the B genome. Two wheat cultivars for which there are sufficient ESTs show different patterns of expression, i.e., with cv Chinese Spring expressing the genes from the A and B genomes, while cv Recital has ESTs from the A and D genomes. Genomic sequences of Chinese Spring show that the D genome gene is inactivated by tandem premature stop codons. A fourth δ-gliadin sequence occurs in the D genome of both Chinese Spring and Ae. tauschii, but no ESTs match this sequence and limited genomic sequences indicates a pseudogene containing frame shifts and premature stop codons. Sequencing of BACs covering a 3 Mb region from Ae. tauschii locates the δ-gliadin gene to the complex Gli-1 plus Glu-3 region on chromosome 1.     

Phosphoantigen/IL2 Expansion and Differentiation of Vγ2Vδ2 T Cells Increase Resistance to Tuberculosis in Nonhuman Primates

Dominant Vγ2Vδ2 T-cell subset exist only in primates, and recognize phosphoantigen from selected pathogens including M. tuberculosis(Mtb). In vivo function of Vγ2Vδ2 T cells in tuberculosis remains unknown. We conducted mechanistic studies to determine whether earlier expansion/differentiation of Vγ2Vδ2 T cells during Mtb infection could increase immune resistance to tuberculosis in macaques. Phosphoantigen/IL-2 administration specifically induced major expansion and pulmonary trafficking/accumulation of phosphoantigen-specific Vγ2Vδ2 T cells, significantly reduced Mtb burdens and attenuated tuberculosis lesions in lung tissues compared to saline/BSA or IL-2 controls. Expanded Vγ2Vδ2 T cells differentiated into multifunctional effector subpopulations capable of producing anti-TB cytokines IFNγ, perforin and granulysin, and co-producing perforin/granulysin in lung tissue. Mechanistically, perforin/granulysin-producing Vγ2Vδ2 T cells limited intracellular Mtb growth, and macaque granulysin had Mtb-bactericidal effect, and inhibited intracellular Mtb in presence of perforin. Furthermore, phosphoantigen/IL2-expanded Vγ2Vδ2 T effector cells produced IL-12, and their expansion/differentiation led to enhanced pulmonary responses of peptide-specific CD4+/CD8+ Th1-like cells. These results provide first in vivo evidence implicating that early expansion/differentiation of Vγ2Vδ2 T effector cells during Mtb infection increases resistance to tuberculosis. Thus, data support a rationale for conducting further studies of the γδ T-cell-targeted treatment of established TB, which might ultimately help explore single or adjunctive phosphoantigen expansion of Vγ2Vδ2 T-cell subset as intervention of MDR-tuberculosis or HIV-related tuberculosis.     

Distinctive Patterns of Evolution of the δ-Globin Gene (HBD) in Primates

In most vertebrates, hemoglobin (Hb) is a heterotetramer composed of two dissimilar globin chains, which change during development according to the patterns of expression of α- and β-globin family members. In placental mammals, the β-globin cluster includes three early-expressed genes, ε(HBE)-γ(HBG)-ψβ(HBBP1), and the late expressed genes, δ (HBD) and β (HBB). While HBB encodes the major adult β-globin chain, HBD is weakly expressed or totally silent. Paradoxically, in human populations HBD shows high levels of conservation typical of genes under strong evolutionary constraints, possibly due to a regulatory role in the fetal-to-adult switch unique of Anthropoid primates. In this study, we have performed a comprehensive phylogenetic and comparative analysis of the two adult β-like globin genes in a set of diverse mammalian taxa, focusing on the evolution and functional divergence of HBD in primates. Our analysis revealed that anthropoids are an exception to a general pattern of concerted evolution in placental mammals, showing a high level of sequence conservation at HBD, less frequent and shorter gene conversion events. Moreover, this lineage is unique in the retention of a functional GATA-1 motif, known to be involved in the control of the developmental expression of the β-like globin genes. We further show that not only the mode but also the rate of evolution of the δ-globin gene in higher primates are strictly associated with the fetal/adult β-cluster developmental switch. To gain further insight into the possible functional constraints that have been shaping the evolutionary history of HBD in primates, we calculated dN/dS (ω) ratios under alternative models of gene evolution. Although our results indicate that HBD might have experienced different selective pressures throughout primate evolution, as shown by different ω values between apes and Old World Monkeys + New World Monkeys (0.06 versus 0.43, respectively), these estimates corroborated a constrained evolution for HBD in Anthropoid lineages, which is unlikely to be related to protein function. Collectively, these findings suggest that sequence change at the δ-globin gene has been under strong selective constraints over 65 Myr of primate evolution, likely due to a regulatory role in ontogenic switches of gene expression.     

Evolutionary Paths of the cAMP-Dependent Protein Kinase (PKA) Catalytic Subunits

3′,5′-cyclic adenosine monophosphate (cAMP) dependent protein kinase or protein kinase A (PKA) has served as a prototype for the large family of protein kinases that are crucially important for signal transduction in eukaryotic cells. The PKA catalytic subunits Cα and Cβ, encoded by the two genes PRKACA and PRKACB, respectively, are among the best understood and characterized human kinases. Here we have studied the evolution of this gene family in chordates, arthropods, mollusks and other animals employing probabilistic methods and show that Cα and Cβ arose by duplication of an ancestral PKA catalytic subunit in a common ancestor of vertebrates. The two genes have subsequently been duplicated in teleost fishes. The evolution of the PRKACG retroposon in simians was also investigated. Although the degree of sequence conservation in the PKA Cα/Cβ kinase family is exceptionally high, a small set of signature residues defining Cα and Cβ subfamilies were identified. These conserved residues might be important for functions that are unique to the Cα or Cβ clades. This study also provides a good example of a seemingly simple phylogenetic problem which, due to a very high degree of sequence conservation and corresponding weak phylogenetic signals, combined with problematic nonphylogenetic signals, is nontrivial for state-of-the-art probabilistic phylogenetic methods.     

The Laminin α Chains: Expression, Developmental Transitions, and Chromosomal Locations of α1-5, Identification of Heterotrimeric Laminins 8–11, and Cloning of a Novel α3 Isoform

Laminin trimers composed of α, β, and γ chains are major components of basal laminae (BLs) throughout the body. To date, three α chains (α1–3) have been shown to assemble into at least seven heterotrimers (called laminins 1–7). Genes encoding two additional α chains (α4 and α5) have been cloned, but little is known about their expression, and their protein products have not been identified. Here we generated antisera to recombinant α4 and α5 and used them to identify authentic proteins in tissue extracts. Immunoprecipitation and immunoblotting showed that α4 and α5 assemble into four novel laminin heterotrimers (laminins 8–11: α4β1γ1, α4β2γ1, α5β1γ1, and α5β2γ1, respectively). Using a panel of nucleotide and antibody probes, we surveyed the expression of α1-5 in murine tissues. All five chains were expressed in both embryos and adults, but each was distributed in a distinct pattern at both RNA and protein levels. Overall, α4 and α5 exhibited the broadest patterns of expression, while expression of α1 was the most restricted. Immunohistochemical analysis of kidney, lung, and heart showed that the α chains were confined to extracellular matrix and, with few exceptions, to BLs. All developing and adult BLs examined contained at least one α chain, all α chains were present in multiple BLs, and some BLs contained two or three α chains. Detailed analysis of developing kidney revealed that some individual BLs, including those of the tubule and glomerulus, changed in laminin chain composition as they matured, expressing up to three different α chains and two different β chains in an elaborate and dynamic progression. Interspecific backcross mapping of the five α chain genes revealed that they are distributed on four mouse chromosomes. Finally, we identified a novel full-length α3 isoform encoded by the Lama3 gene, which was previously believed to encode only truncated chains. Together, these results reveal remarkable diversity in BL composition and complexity in BL development.     

Significance of Population Size on the Fixation of Nonsynonymous Mutations in Genes Under Varying Levels of Selection Pressure

Previous studies observed a higher ratio of divergences at nonsynonymous and synonymous sites (ω = d_N/d_S) in species with a small population size compared to that estimated for those with a large population size. Here we examined the theoretical relationship between ω, effective population size (N_e), and selection coefficient (s). Our analysis revealed that when purifying selection is high, ω of species with small N_e is much higher than that of species with large N_e. However the difference between the two ω reduces with the decline in selection pressure (s → 0). We examined this relationship using primate and rodent genes and found that the ω estimated for highly constrained genes of primates was up to 2.9 times higher than that obtained for their orthologous rodent genes. Conversely, for genes under weak purifying selection the ω of primates was only 17% higher than that of rodents. When tissue specificity was used as a proxy for selection pressure we found that the ω of broadly expressed genes of primates was up to 2.1-fold higher than that of their rodent counterparts and this difference was only 27% for tissue specific genes. Since most of the nonsynonymous mutations in constrained or broadly expressed genes are deleterious, fixation of these mutations is influenced by N_e. This results in a higher ω of these genes in primates compared to those from rodents. Conversely, the majority of nonsynonymous mutations in less-constrained or tissue-specific genes are neutral or nearly neutral and therefore fixation of them is largely independent of N_e, which leads to the similarity of ω in primates and rodents.     

Gamma/Delta T-cell functional responses differ after pathogenic human immunodeficiency virus and nonpathogenic simian immunodeficiency virus infections

The objective of this study was to functionally assess gamma/delta (γδ) T cells following pathogenic human immunodeficiency virus (HIV) infection of humans and nonpathogenic simian immunodeficiency virus (SIV) infection of sooty mangabeys. γδ T cells were obtained from peripheral blood samples from patients and sooty mangabeys that exhibited either a CD4-healthy (>200 CD4⁺ T cells/μl blood) or CD4-low (<200 CD4 cells/μl blood) phenotype. Cytokine flow cytometry was utilized to assess production of Th1 cytokines tumor necrosis factor alpha and gamma interferon following ex vivo stimulation with either phorbol myristate acetate/ionomycin or the Vδ2 γδ T-cell receptor agonist isopentenyl pyrophosphate. Sooty mangabeys were observed to have higher percentages of γδ T cells in their peripheral blood than humans did. Following stimulation, γδ T cells from SIV-positive (SIV⁺) mangabeys maintained or increased their ability to express the Th1 cytokines regardless of CD4⁺ T-cell levels. In contrast, HIV-positive (HIV⁺) patients exhibited a decreased percentage of γδ T cells expressing Th1 cytokines following stimulation. This dysfunction is primarily within the Vδ2⁺ γδ T-cell subset which incurred both a decreased overall level in the blood and a reduced Th1 cytokine production. Patients treated with highly active antiretroviral therapy exhibited a partial restoration in their γδ T-cell Th1 cytokine response that was intermediate between the responses of the uninfected and HIV⁺ patients. The SIV⁺ sooty mangabey natural hosts, which do not proceed to clinical AIDS, provide evidence that γδ T-cell dysfunction occurs in HIV⁺ patients and may contribute to HIV disease progression.     

Renal Cell Carcinoma-Infiltrating CD3low Vγ9Vδ1 T Cells Represent Potentially Novel Anti-Tumor Immune Players

Due to the highly immunogenic nature of renal cell carcinoma (RCC), the tumor microenvironment (TME) is enriched with various innate and adaptive immune subsets. In particular, gamma-delta (γδ) T cells can act as potent attractive mediators of adoptive cell transfer immunotherapy because of their unique properties such as non-reliance on major histocompatibility complex expression, their ability to infiltrate human tumors and recognize tumor antigens, relative insensitivity to immune checkpoint molecules, and broad tumor cytotoxicity. Therefore, it is now critical to better characterize human γδ T-cell subsets and their mechanisms in RCCs, especially the stage of differentiation. In this study, we aimed to identify γδ T cells that might have adaptive responses against RCC progression. We characterized γδ T cells in peripheral blood and tumor-infiltrating lymphocytes (TILs) in freshly resected tumor specimens from 20 RCC patients. Furthermore, we performed a gene set enrichment analysis on RNA-sequencing data from The Cancer Genome Atlas (TCGA) derived from normal kidneys and RCC tumors to ascertain the association between γδ T-cell infiltration and anti-cancer immune activity. Notably, RCC-infiltrating CD3^low Vγ9Vδ1 T cells with a terminally differentiated effector memory phenotype with up-regulated activation/exhaustion molecules were newly detected as predominant TILs, and the cytotoxic activity of these cells against RCC was confirmed in vitro. In an additional analysis of the TCGA RCC dataset, γδ T-cell enrichment scores correlated strongly with those for CTLs, Th1 cells, “exhausted” T cells, and M1 macrophages, suggesting active involvement of γδ T cells in anti-tumor rather than pro-tumor activity, and Vδ1 cells were more abundant than Vδ2 or Vδ3 cells in RCC tumor samples. Thus, we posit that Vγ9Vδ1 T cells may represent an excellent candidate for adoptive immunotherapy in RCC patients with a high risk of relapse after surgery.     

The thymic cortical epithelium determines the TCR repertoire of IL-17-producing γδT cells

The thymus provides a specialized microenvironment in which distinct subsets of thymic epithelial cells (TECs) support T-cell development. Here, we describe the significance of cortical TECs (cTECs) in T-cell development, using a newly established mouse model of cTEC deficiency. The deficiency of mature cTECs caused a massive loss of thymic cellularity and impaired the development of αβT cells and invariant natural killer T cells. Unexpectedly, the differentiation of certain γδT-cell subpopulations—interleukin-17-producing Vγ4 and Vγ6 cells—was strongly dysregulated, resulting in the perturbation of γδT-mediated inflammatory responses in peripheral tissues. These findings show that cTECs contribute to the shaping of the TCR repertoire, not only of “conventional” αβT cells but also of inflammatory “innate” γδT cells.