A highly diverse set of novel immunoglobulin-like transcript (NILT) genes in zebrafish indicates a wide range of functions with complex relationships to mammalian receptors.

Multiple novel immunoglobulin-like transcripts (NILTs) have been identified from salmon, trout, and carp. NILTs typically encode activating or inhibitory transmembrane receptors with extracellular immunoglobulin (Ig) domains. Although predicted to provide immune recognition in ray-finned fish, we currently lack a definitive framework of NILT diversity, thereby limiting our predictions for their evolutionary origin and function. In order to better understand the diversity of NILTs and their possible roles in immune function, we identified five NILT loci in the Atlantic salmon (Salmo salar) genome, defined 86 NILT Ig domains within a 3-Mbp region of zebrafish (Danio rerio) chromosome 1, and described 41 NILT Ig domains as part of an alternative haplotype for this same genomic region. We then identified transcripts encoded by 43 different NILT genes which reflect an unprecedented diversity of Ig domain sequences and combinations for a family of non-recombining receptors within a single species. Zebrafish NILTs include a sole putative activating receptor but extensive inhibitory and secreted forms as well as membrane-bound forms with no known signaling motifs. These results reveal a higher level of genetic complexity, interindividual variation, and sequence diversity for NILTs than previously described, suggesting that this gene family likely plays multiple roles in host immunity.

Both Conifer II and Gnetales are characterized by a high frequency of ancient mitochondrial gene transfer to the nuclear genome.

BACKGROUND: Mitochondrial gene transfer/loss is common in land plants, and therefore the fate of missing mitochondrial genes has attracted more and more attention. The gene content of gymnosperm mitochondria varies greatly, supplying a system for studying the evolutionary fate of missing mitochondrial genes. RESULTS: Here, we studied the tempo and pattern of mitochondrial gene transfer/loss in gymnosperms represented by all 13 families, using high-throughput sequencing of both DNA and cDNA. All 41 mitochondrial protein-coding genes were found in cycads, Ginkgo and Pinaceae, whereas multiple mitochondrial genes were absent in Conifer II and Gnetales. In Conifer II, gene transfer from mitochondria to the nucleus followed by loss of the mitochondrial copy was common, but complete loss of a gene in both mitochondrial and nuclear genomes was rare. In contrast, both gene transfer and loss were commonly found in Gnetales. Notably, in Conifer II and Gnetales, the same five mitochondrial genes were transferred to the nuclear genome, and these gene transfer events occurred, respectively, in ancestors of the two lineages. A two-step transfer mechanism (retroprocessing and subsequent DNA-mediated gene transfer) may be responsible for mitochondrial gene transfer in Conifer II and Gnetales. Moreover, the mitochondrial gene content variation is correlated with gene length, GC content, hydrophobicity, and nucleotide substitution rates in land plants. CONCLUSIONS: This study reveals a complete evolutionary scenario for variations of mitochondrial gene transferring in gymnosperms, and the factors responsible for mitochondrial gene content variation in land plants.

Population Genomics Reveals Molecular Determinants of Specialization to Tomato in the Polyphagous Fungal Pathogen Botrytis cinerea in France.

Many fungal plant pathogens encompass multiple populations specialized on different plant species. Understanding the factors underlying pathogen adaptation to their hosts is a major challenge of evolutionary microbiology, and it should help to prevent the emergence of new specialized pathogens on novel hosts. Previous studies have shown that French populations of the gray mold pathogen Botrytis cinerea parasitizing tomato and grapevine are differentiated from each other, and have higher aggressiveness on their host of origin than on other hosts, indicating some degree of host specialization in this polyphagous pathogen. Here, we aimed at identifying the genomic features underlying the specialization of B. cinerea populations to tomato and grapevine. Based on whole genome sequences of 32 isolates, we confirmed the subdivision of B. cinerea pathogens into two genetic clusters on grapevine and another, single cluster on tomato. Levels of genetic variation in the different clusters were similar, suggesting that the tomato-specific cluster has not recently emerged following a bottleneck. Using genome scans for selective sweeps and divergent selection, tests of positive selection based on polymorphism and divergence at synonymous and nonsynonymous sites, and analyses of presence and absence variation, we identified several candidate genes that represent possible determinants of host specialization in the tomato-associated population. This work deepens our understanding of the genomic changes underlying the specialization of fungal pathogen populations.

Transposable elements contribute to dynamic genome content in maize.

Transposable elements (TEs) are ubiquitous components of eukaryotic genomes and can create variation in genome organization and content. Most maize genomes are composed of TEs. We developed an approach to define shared and variable TE insertions across genome assemblies and applied this method to four maize genomes (B73, W22, Mo17 and PH207) with uniform structural annotations of TEs. Among these genomes we identified approximately 400 000 TEs that are polymorphic, encompassing 1.6 Gb of variable TE sequence. These polymorphic TEs include a combination of recent transposition events as well as deletions of older TEs. There are examples of polymorphic TEs within each of the superfamilies of TEs and they are found distributed across the genome, including in regions of recent shared ancestry among individuals. There are many examples of polymorphic TEs within or near maize genes. In addition, there are 2380 gene annotations in the B73 genome that are located within variable TEs, providing evidence for the role of TEs in contributing to the substantial differences in annotated gene content among these genotypes. TEs are highly variable in our survey of four temperate maize genomes, highlighting the major contribution of TEs in driving variation in genome organization and gene content. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://github.com/SNAnderson/maizeTE_variation; https://mcstitzer.github.io/maize_TEs.

Insights into hominin phenotypic and dietary evolution from ancient DNA sequence data.

Nuclear genome sequence data from Neandertals, Denisovans, and archaic anatomically modern humans can be used to complement our understanding of hominin evolutionary biology and ecology through i) direct inference of archaic hominin phenotypes, ii) indirect inference of those phenotypes by identifying the effects of previously-introgressed alleles still present among modern humans, or iii) determining the evolutionary timing of relevant hominin-specific genetic changes. Here we review and reanalyze published Neandertal and Denisovan genome sequence data to illustrate an example of the third approach. Specifically, we infer the timing of five human gene presence/absence changes that may be related to particular hominin-specific dietary changes and discuss these results in the context of our broader reconstructions of hominin evolutionary ecology. We show that pseudogenizing (gene loss) mutations in the TAS2R62 and TAS2R64 bitter taste receptor genes and the MYH16 masticatory myosin gene occurred after the hominin-chimpanzee divergence but before the divergence of the human and Neandertal/Denisovan lineages. The absence of a functional MYH16 protein may explain our relatively reduced jaw muscles; this gene loss may have followed the adoption of cooking behavior. In contrast, salivary amylase gene (AMY1) duplications were not observed in the Neandertal and Denisovan genomes, suggesting a relatively recent origin for the AMY1 copy number gains that are observed in modern humans. Thus, if earlier hominins were consuming large quantities of starch-rich underground storage organs, as previously hypothesized, then they were likely doing so without the digestive benefits of increased salivary amylase production. Our most surprising result was the observation of a heterozygous mutation in the first codon of the TAS2R38 bitter taste receptor gene in the Neandertal individual, which likely would have resulted in a non-functional protein and inter-individual PTC (phenylthiocarbamide) taste sensitivity variation, as also observed in both humans and chimpanzees.

KIR Gene Content Diversity in a Zimbabwean Population: Does KIR2DL2 Have a Role in Protection Against Human Immunodeficiency Virus Infection?

Killer cell immunoglobulin-like receptors (KIRs) mediate natural killer cell function through interaction with their cognate human leukocyte antigen ligands. Thus, KIR gene variants have been implicated in resistance or susceptibility to viral infections. However, research on the role of these variants remains contradictory and inconclusive. In the present study, we investigated KIR gene content diversity and its association with human immunodeficiency virus (HIV) infection in an adult Black Zimbabwean population. Presence or absence of 15 KIR genes was determined in 189 HIV-infected adults and 97 HIV-uninfected blood donors using sequence specific primer polymerase chain reaction. Frequencies of KIR genes, genotypes, and haplotypes were compared between the cases and controls to identify putative associations between KIR gene variants and HIV status. We report in this study the frequencies of 15 KIR genes and 43 KIR genotypes (40 known and 3 novel) among Zimbabweans. Importantly, the frequency of the inhibitory KIR2DL2 gene was significantly higher in the uninfected group (62%) compared to the HIV-infected group (47%) (OR = 0.55, 95% CI: 0.33-0.90, p = 0.019). KIR2DL2/2DL2 homozygosity was also significantly higher in the uninfected group (35%) compared to HIV-infected group (53%) (OR = 0.33, 95% CI: 0.16-0.72, p = 0.005) under a recessive model. We conclude that the KIR2DL2 gene may be involved in protection against HIV infection. It may be possible that inhibitory KIR genes may have an important role to play in HIV acquisition among populations of African origin in whom the activating KIR genes are less frequent compared to among Caucasians.