For many pathogens with environmental stages, or those carried by vectors or intermediate hosts, disease transmission is strongly influenced by pathogen, host, and vector movements across complex landscapes, and thus quantitative measures of movement rate and direction can reveal new opportunities for disease management and intervention. Genetic assignment methods are a set of powerful statistical approaches useful for establishing population membership of individuals. Recent theoretical improvements allow these techniques to be used to cost-effectively estimate the magnitude and direction of key movements in infectious disease systems, revealing important ecological and environmental features that facilitate or limit transmission. Here, we review the theory, statistical framework, and molecular markers that underlie assignment methods, and we critically examine recent applications of assignment tests in infectious disease epidemiology. Research directions that capitalize on use of the techniques are discussed, focusing on key parameters needing study for improved understanding of patterns of disease.

Abstract Re-emergence of schistosomiasis in regions of China where control programs have ceased requires development of molecular-genetic tools to track gene flow and assess genetic diversity of Schistosoma populations. We identified many microsatellite loci in the draft genome of Schistosoma japonicum using defined search criteria and selected a subset for further analysis. From an initial panel of 50 loci, 20 new microsatellites were selected for eventual optimization and application to a panel of worms from endemic areas. All but one of the selected microsatellites contain simple tri-nucleotide repeats. Moderate to high levels of polymorphism were detected. Numbers of alleles ranged from 6 to 14 and observed heterozygosity was always >0.6. The loci reported here will facilitate high resolution population-genetic studies on schistosomes in re-emergent foci.

The interface between the membrane (MS) and cytoplasmic (C) rings of the bacterial flagellar motor couples torque generation to rotation within the membrane. The structure of the C-terminal helices of the integral membrane protein FliF (FliF_C) bound to the N terminal domain of the switch complex protein FliG (FliG_N) reveals that FliG_N folds around FliF_C to produce a topology that closely resembles both the middle and C-terminal domains of FliG. The interface is consistent with solution-state nuclear magnetic resonance, small-angle X-ray scattering, in vivo interaction studies, and cellular motility assays. Co-folding with FliF_C induces substantial conformational changes in FliG_N and suggests that FliF and FliG have the same stoichiometry within the rotor. Modeling the FliF_C:FliG_N complex into cryo-electron microscopy rotor density updates the architecture of the middle and upper switch complex and shows how domain shuffling of a conserved interaction module anchors the cytoplasmic rotor to the membrane.

Phosphorylation of Escherichia coli CheY protein transduces chemoreceptor stimulation to a highly cooperative flagellar motor response. CheY binds to the N-terminal peptide of the FliM motor protein (FliM_N). Constitutively active D13K-Y106W CheY has been an important tool for motor physiology. The crystal structures of CheY and CheY ⋅ FliM_N with and without D13K-Y106W have shown FliM_N-bound CheY contains features of both active and inactive states. We used molecular dynamics (MD) simulations to characterize the CheY conformational landscape accessed by FliM_N and D13K-Y106W. Mutual information measures identified the central features of the long-range CheY allosteric network between D57 phosphorylation site and the FliM_N interface, namely the closure of the α4-β4 hinge and inward rotation of Y- or W106 with W58. We used hydroxy-radical foot printing with mass spectroscopy (XFMS) to track the solvent accessibility of these and other side chains. The solution XFMS oxidation rate correlated with the solvent-accessible area of the crystal structures. The protection of allosteric relay side chains reported by XFMS confirmed the intermediate conformation of the native CheY ⋅ FliM_N complex, the inactive state of free D13K-Y106W CheY, and the MD-based network architecture. We extended the MD analysis to determine temporal coupling and energetics during activation. Coupled aromatic residue rotation was a graded rather than a binary switch, with Y- or W106 side-chain burial correlated with increased FliM_N affinity. Activation entrained CheY fold stabilization to FliM_N affinity. The CheY network could be partitioned into four dynamically coordinated sectors. Residue substitutions mapped to sectors around D57 or the FliM_N interface according to phenotype. FliM_N increased sector size and interactions. These sectors fused between the substituted K13-W106 residues to organize a tightly packed core and novel surfaces that may bind additional sites to explain the cooperative motor response. The community maps provide a more complete description of CheY priming than proposed thus far.

After infection, many factors coordinate the population expansion and differentiation of CD8+ effector and memory T cells. Using data of unparalleled breadth from the Immunological Genome Project, we analyzed the CD8+ T cell transcriptome throughout infection to establish gene-expression signatures and identify putative transcriptional regulators. Notably, we found that the expression of key gene signatures can be used to predict the memory-precursor potential of CD8+ effector cells. Long-lived memory CD8+ cells ultimately expressed a small subset of genes shared by natural killer T and γδ T cells. Although distinct inflammatory milieu and T cell precursor frequencies influenced the differentiation of CD8+ effector and memory populations, core transcriptional signatures were regulated similarly, whether polyclonal or transgenic, and whether responding to bacterial or viral model pathogens. Our results provide insights into the transcriptional regulation that influence memory formation and CD8+ T cell immunity.

Purpose

To describe a case of primary coenzyme Q₁₀ deficiency in a child manifesting as early-onset renal failure, retinal dystrophy, and optic atrophy leading to progressive vision loss.

Methods

Clinical presentation and workup including visual fields, electroretinogram, and optical coherence tomography are presented. Genetic testing was performed.

Results

An eight-year-old female with nephropathy requiring renal transplantation subsequently developed progressive cone-rod dystrophy and optic atrophy. The patient had negative results on a targeted next-generation sequencing retinal dystrophy panel but whole-exome sequencing revealed two variants in COQ2 (likely biallelic), consistent with a diagnosis of primary coenzyme Q₁₀ deficiency.

Conclusions

Primary coenzyme Q₁₀ deficiency is a rare disorder with variable systemic and ocular findings; there is also genetic heterogeneity. Genetic testing aids in the diagnosis of this condition, and variants in the COQ2 and PDSS1 genes appear to have the strongest association with ocular manifestations. Oral supplementation of coenzyme Q₁₀ may slow progression of disease. This case highlights the utility of whole-exome sequencing in the diagnosis of a rare syndromic form of ocular disease and reports a novel phenotypic association for this condition.

Molecular genetic tools are needed to address questions as to the source and dynamics of transmission of the human blood fluke Schistosoma japonicum in regions where human infections have reemerged, and to characterize infrapopulations in individual hosts. The life stage that interests us as a target for collecting genotypic data is the miracidium, a very small larval stage that consequently yields very little DNA for analysis. Here, we report the successful development of a multiplex format permitting genotyping of 17 microsatellite loci in four sequential multiplex reactions using a single miracidium held on a Whatman Classic FTA indicating card. This approach was successful after short storage periods, but after long storage (>4 years), considerable difficulty was encountered in multiplex genotyping, necessitating the use of whole genome amplification (WGA) methods. WGA applied to cards stored for long periods of time resulted in sufficient DNA for accurate and repeatable genotyping. Trials and tests of these methods, as well as application to some field-collected samples, are reported, along with the discussion of the potential insights to be gained from such techniques. These include recognition of sibships among miracidia from a single host, and inference of the minimum number of worm pairs that might be present in a host.

Opisthorchis viverrini is a major public health concern in Southeast Asia. Various reports have suggested that this parasite may represent a species complex, with genetic structure in the region perhaps being dictated by geographical factors and different species of intermediate hosts. We used four microsatellite loci to analyze O. viverrini adult worms originating from six species of cyprinid fish in Thailand and Lao PDR. Two distinct O. viverrini populations were observed. In Ban Phai, Thailand, only one subgroup occurred, hosted by two different fish species. Both subgroups occurred in fish from That Luang, Lao PDR, but were represented to very different degrees among the fish hosts there. Our data suggest that, although geographical separation is more important than fish host specificity in influencing genetic structure, it is possible that two species of Opisthorchis, with little interbreeding, are present near Vientiane in Lao PDR.

Alternative splicing (AS) allows increased diversity and orthogonal regulation of the transcriptional products of mammalian genomes. To assess the distribution and variation of alternative splicing across cell lineages of the immune system, we comprehensively analyzed RNA sequencing and microarray data generated by the Immunological Genome Project Consortium. AS is pervasive: 60% of genes showed frequent AS isoforms in T or B lymphocytes, with 7,599 previously unreported isoforms. Distinct cell specificity was observed, with differential exon skipping in 5% of genes otherwise coexpressed in both B and T cells. The distribution of isoforms was mostly all or none, suggesting on/off switching as a frequent mode of AS regulation in lymphocytes. From the identification of differential exon use in the microarray data, clustering of exon inclusion/exclusion patterns across all Immunological Genome Project cell types showed that ∼70% of AS exons are distributed along a common pattern linked to lineage differentiation and cell cycling. Other AS events distinguished myeloid from lymphoid cells or affected only a small set of exons without clear lineage specificity (e.g., Ptprc). Computational analysis predicted specific associations between AS exons and splicing regulators, which were verified by detection of the hnRPLL/Ptprc connection.

To determine the breadth and underpinning of changes in immunocyte gene expression due to genetic variation in mice, we performed, as part of the Immunological Genome Project, gene expression profiling for CD4(+) T cells and neutrophils purified from 39 inbred strains of the Mouse Phenome Database. Considering both cell types, a large number of transcripts showed significant variation across the inbred strains, with 22% of the transcriptome varying by 2-fold or more. These included 119 loci with apparent complete loss of function, where the corresponding transcript was not expressed in some of the strains, representing a useful resource of "natural knockouts." We identified 1222 cis-expression quantitative trait loci (cis-eQTL) that control some of this variation. Most (60%) cis-eQTLs were shared between T cells and neutrophils, but a significant portion uniquely impacted one of the cell types, suggesting cell type-specific regulatory mechanisms. Using a conditional regression algorithm, we predicted regulatory interactions between transcription factors and potential targets, and we demonstrated that these predictions overlap with regulatory interactions inferred from transcriptional changes during immunocyte differentiation. Finally, comparison of these and parallel data from CD4(+) T cells of healthy humans demonstrated intriguing similarities in variability of a gene's expression: the most variable genes tended to be the same in both species, and there was an overlap in genes subject to strong cis-acting genetic variants. We speculate that this "conservation of variation" reflects a differential constraint on intraspecies variation in expression levels of different genes, either through lower pressure for some genes, or by favoring variability for others.