Whole-exome sequencing (WES) is rapidly evolving into a tool of choice for rapid, and inexpensive identification of molecular genetic lesions within targeted regions of the human genome. While biases in WES coverage of nucleotides in targeted regions are recognized, it is not well understood how repetition of WES improves the interpretation of sequencing results in a clinical diagnostic setting.

Carbapenem-resistant Klebsiella pneumoniae strains are formidable hospital pathogens that pose a serious threat to patients around the globe due to a rising incidence in health care facilities, high mortality rates associated with infection, and potential to spread antibiotic resistance to other bacterial species, such as Escherichia coli Over 6 months in 2011, 17 patients at the National Institutes of Health (NIH) Clinical Center became colonized with a highly virulent, transmissible carbapenem-resistant strain of K. pneumoniae Our real-time genomic sequencing tracked patient-to-patient routes of transmission and informed epidemiologists' actions to monitor and control this outbreak. Two of these patients remained colonized with carbapenemase-producing organisms for at least 2 to 4 years, providing the opportunity to undertake a focused genomic study of long-term colonization with antibiotic-resistant bacteria. Whole-genome sequencing studies shed light on the underlying complex microbial colonization, including mixed or evolving bacterial populations and gain or loss of plasmids. Isolates from NIH patient 15 showed complex plasmid rearrangements, leaving the chromosome and the blaKPC-carrying plasmid intact but rearranging the two other plasmids of this outbreak strain. NIH patient 16 has shown continuous colonization with blaKPC-positive organisms across multiple time points spanning 2011 to 2015. Genomic studies defined a complex pattern of succession and plasmid transmission across two different K. pneumoniae sequence types and an E. coli isolate. These findings demonstrate the utility of genomic methods for understanding strain succession, genome plasticity, and long-term carriage of antibiotic-resistant organisms.

ATAD5, the human ortholog of yeast Elg1, plays a role in PCNA deubiquitination. Since PCNA modification is important to regulate DNA damage bypass, ATAD5 may be important for suppression of genomic instability in mammals in vivo. To test this hypothesis, we generated heterozygous (Atad5(+/m)) mice that were haploinsuffficient for Atad5. Atad5(+/m) mice displayed high levels of genomic instability in vivo, and Atad5(+/m) mouse embryonic fibroblasts (MEFs) exhibited molecular defects in PCNA deubiquitination in response to DNA damage, as well as DNA damage hypersensitivity and high levels of genomic instability, apoptosis, and aneuploidy. Importantly, 90% of haploinsufficient Atad5(+/m) mice developed tumors, including sarcomas, carcinomas, and adenocarcinomas, between 11 and 20 months of age. High levels of genomic alterations were evident in tumors that arose in the Atad5(+/m) mice. Consistent with a role for Atad5 in suppressing tumorigenesis, we also identified somatic mutations of ATAD5 in 4.6% of sporadic human endometrial tumors, including two nonsense mutations that resulted in loss of proper ATAD5 function. Taken together, our findings indicate that loss-of-function mutations in mammalian Atad5 are sufficient to cause genomic instability and tumorigenesis.

Nectins are immunoglobulin superfamily glycoproteins that mediate intercellular adhesion in many vertebrate tissues. Homophilic and heterophilic interactions between nectin family members help mediate tissue patterning. We determined the homophilic binding affinities and heterophilic specificities of all four nectins and the related protein nectin-like 5 (Necl-5) from human and mouse, revealing a range of homophilic interaction strengths and a defined heterophilic specificity pattern. To understand the molecular basis of their adhesion and specificity, we determined the crystal structures of natively glycosylated full ectodomains or adhesive fragments of all four nectins and Necl-5. All of the crystal structures revealed dimeric nectins bound through a stereotyped interface that was previously proposed to represent a cis dimer. However, conservation of this interface and the results of targeted cross-linking experiments showed that this dimer probably represents the adhesive trans interaction. The structure of the dimer provides a simple molecular explanation for the adhesive binding specificity of nectins.

The function of the actin-binding domain of α-catenin, αABD, including its possible role in the direct anchorage of the cadherin-catenin complex to the actin cytoskeleton, has remained uncertain. We identified two point mutations on the αABD surface that interfere with αABD binding to actin and used them to probe the role of α-catenin-actin interactions in adherens junctions. We found that the junctions directly bound to actin via αABD were more dynamic than the junctions bound to actin indirectly through vinculin and that recombinant αABD interacted with cortical actin but not with actin bundles. This interaction resulted in the formation of numerous short-lived cortex-bound αABD clusters. Our data suggest that αABD clustering drives the continuous assembly of transient, actin-associated cadherin-catenin clusters whose disassembly is maintained by actin depolymerization. It appears then that such actin-dependent αABD clustering is a unique molecular mechanism mediating both integrity and reassembly of the cell-cell adhesive interface formed through weak cis- and trans-intercadherin interactions.

The Hawaiian honeycreepers are an avian adaptive radiation containing many endangered and extinct species. They display a dramatic range of phenotypic variation and are a model system for studies of evolution, conservation, disease dynamics and population genetics. Development of a genome-scale resources for this group would augment the quality of research focusing on Hawaiian honeycreepers and facilitate comparative avian genomic research.

DNA methylation is an essential epigenetic process in mammals, intimately involved in gene regulation. Here we address the extent to which genetics, sex, and pregnancy influence genomic DNA methylation by intercrossing 2 inbred mouse strains, C57BL/6N and C3H/HeN, and analyzing DNA methylation in parents and offspring using whole-genome bisulfite sequencing. Differential methylation across genotype is detected at thousands of loci and is preserved on parental alleles in offspring. In comparison of autosomal DNA methylation patterns across sex, hundreds of differentially methylated regions are detected. Comparison of animals with different histories of pregnancy within our study reveals a CpG methylation pattern that is restricted to female animals that had borne offspring. Collectively, our results demonstrate the stability of CpG methylation across generations, clarify the interplay of epigenetics with genetics and sex, and suggest that CpG methylation may serve as an epigenetic record of life events in somatic tissues at loci whose expression is linked to the relevant biology.

Human microbiome studies have revealed the intricate interplay of host immunity and bacterial communities to achieve homeostatic balance. Healthy skin microbial communities are dominated by bacteria with low viral representation1-3, mainly bacteriophage. Specific eukaryotic viruses have been implicated in both common and rare skin diseases, but cataloging skin viral communities has been limited. Alterations in host immunity provide an opportunity to expand our understanding of microbial-host interactions. Primary immunodeficient patients manifest with various viral, bacterial, fungal, and parasitic infections, including skin infections4. Dedicator of cytokinesis 8 (DOCK8) deficiency is a rare primary human immunodeficiency characterized by recurrent cutaneous and systemic infections, as well as atopy and cancer susceptibility5. DOCK8, encoding a guanine nucleotide exchange factor highly expressed in lymphocytes, regulates actin cytoskeleton, which is critical for migration through collagen-dense tissues such as skin6. Analyzing deep metagenomic sequencing data from DOCK8-deficient skin samples demonstrated a notable increase in eukaryotic viral representation and diversity compared with healthy volunteers. De novo assembly approaches identified hundreds of novel human papillomavirus genomes, illuminating microbial dark matter. Expansion of the skin virome in DOCK8-deficient patients underscores the importance of immune surveillance in controlling eukaryotic viral colonization and infection.

Drosophila Dpr (21 paralogs) and DIP proteins (11 paralogs) are cell recognition molecules of the immunoglobulin superfamily (IgSF) that form a complex protein interaction network. DIP and Dpr proteins are expressed in a synaptic layer-specific fashion in the visual system. How interactions between these proteins regulate layer-specific synaptic circuitry is not known. Here we establish that DIP-α and its interacting partners Dpr6 and Dpr10 regulate multiple processes, including arborization within layers, synapse number, layer specificity, and cell survival. We demonstrate that heterophilic binding between Dpr6/10 and DIP-α and homophilic binding between DIP-α proteins promote interactions between processes in vivo. Knockin mutants disrupting the DIP/Dpr binding interface reveal a role for these proteins during normal development, while ectopic expression studies support an instructive role for interactions between DIPs and Dprs in circuit development. These studies support an important role for the DIP/Dpr protein interaction network in regulating cell-type-specific connectivity patterns.

Binding between DIP and Dpr neuronal recognition proteins has been proposed to regulate synaptic connections between lamina and medulla neurons in the Drosophila visual system. Each lamina neuron was previously shown to express many Dprs. Here, we demonstrate, by contrast, that their synaptic partners typically express one or two DIPs, with binding specificities matched to the lamina neuron-expressed Dprs. A deeper understanding of the molecular logic of DIP/Dpr interaction requires quantitative studies on the properties of these proteins. We thus generated a quantitative affinity-based DIP/Dpr interactome for all DIP/Dpr protein family members. This revealed a broad range of affinities and identified homophilic binding for some DIPs and some Dprs. These data, along with full-length ectodomain DIP/Dpr and DIP/DIP crystal structures, led to the identification of molecular determinants of DIP/Dpr specificity. This structural knowledge, along with a comprehensive set of quantitative binding affinities, provides new tools for functional studies in vivo.

Single particle cryo-electron microscopy (cryoEM) is often performed under the assumption that particles are not adsorbed to the air-water interfaces and in thin, vitreous ice. In this study, we performed fiducial-less tomography on over 50 different cryoEM grid/sample preparations to determine the particle distribution within the ice and the overall geometry of the ice in grid holes. Surprisingly, by studying particles in holes in 3D from over 1000 tomograms, we have determined that the vast majority of particles (approximately 90%) are adsorbed to an air-water interface. The implications of this observation are wide-ranging, with potential ramifications regarding protein denaturation, conformational change, and preferred orientation. We also show that fiducial-less cryo-electron tomography on single particle grids may be used to determine ice thickness, optimal single particle collection areas and strategies, particle heterogeneity, and de novo models for template picking and single particle alignment.

The tip link, a filament formed by protocadherin 15 (PCDH15) and cadherin 23, conveys mechanical force from sound waves and head movement to open hair-cell mechanotransduction channels. Tip-link cadherins are thought to have acquired structural features critical for their role in mechanotransduction. Here, we biophysically and structurally characterize the unusual cis-homodimeric architecture of PCDH15. We show that PCDH15 molecules form double-helical assemblies through cis-dimerization interfaces in the extracellular cadherin EC2-EC3 domain region and in a unique membrane-proximal domain. Electron microscopy studies visualize the cis-dimeric PCDH15 assembly and reveal the PCDH15 extracellular domain as a parallel double helix with cis cross-bridges at the two locations we defined. The helical configuration suggests the potential for elasticity through helix winding and unwinding. Functional studies in hair cells show that mutations that perturb PCDH15 dimerization contacts affect mechanotransduction. Together, these data reveal the cis-dimeric architecture of PCDH15 and show that dimerization is critical for sensing mechanical stimuli.

Metastasis is characterized by spreading of neoplastic cells to an organ other than where they originated and is the predominant cause of death among cancer patients. This holds true for melanoma, whose incidence is increasing more rapidly than any other cancer and once disseminated has few therapeutic options. Here we performed whole exome sequencing of two sets of matched normal and metastatic tumor DNAs.

The kallikrein (KLK) gene family comprises the largest uninterrupted locus of serine proteases in the human genome and represents a notable case for studying the evolutionary fate of duplicated genes. In primates, a recent duplication event gave rise to KLK2 and KLK3, both encoding essential proteins for the cascade of seminal plasma liquefaction. We reconstructed the evolutionary history of KLK2 and KLK3 by comparative analysis of the orthologous sequences from 22 primate species, calculated d(N)/d(S) ratios, and addressed the hypothesis of coevolution with their substrates, the semenogelins (SEMG1 and SEMG2). Our findings support the placement of the KLK2-KLK3 duplication in the Catarrhini ancestor and unveil the frequent loss of KLK2 throughout primate evolution by different genomic mechanisms, including unequal crossing-over, deletions, and pseudogenization. We provide evidences for an adaptive evolution of KLK3 toward an expanded enzymatic spectrum, with an effect on the hydrolysis of semen coagulum. Furthermore, we found associations between mating system, the number of SEMG repeat units, and the number of functional KLK2 and KLK3, suggesting complex evolutionary dynamics shaped by reproductive biology.

Although long noncoding RNAs (lncRNAs) are proposed to play essential roles in mammalian neurodevelopment, we know little of their functions from their disruption in vivo. Combining evidence for evolutionary constraint and conserved expression data, we previously identified candidate lncRNAs that might play important and conserved roles in brain function. Here, we demonstrate that the sequence and neuronal transcription of lncRNAs transcribed from the previously uncharacterized Visc locus are conserved across diverse mammals. Consequently, one of these lncRNAs, Visc-2, was selected for targeted deletion in the mouse, and knockout animals were subjected to an extremely detailed anatomical and behavioral characterization. Despite a neurodevelopmental expression pattern of Visc-2 that is highly localized to the cortex and sites of neurogenesis, anomalies in neither cytoarchitecture nor neuroproliferation were identified in knockout mice. In addition, no abnormal motor, sensory, anxiety, or cognitive behavioral phenotypes were observed. These results are important because they contribute to a growing body of evidence that lncRNA loci contribute on average far less to brain and biological functions than protein-coding loci. A high-throughput knockout program focussing on lncRNAs, similar to that currently underway for protein-coding genes, will be required to establish the distribution of their organismal functions.

Bacterial whole-genome sequencing (WGS) of human pathogens has provided unprecedented insights into the evolution of antibiotic resistance. Most studies have focused on identification of resistance mutations, leaving one to speculate on the fate of these mutants once the antibiotic selective pressure is removed. We performed WGS on longitudinal isolates of Acinetobacter baumannii from patients undergoing colistin treatment, and upon subsequent drug withdrawal. In each of the four patients, colistin resistance evolved via mutations at the pmr locus. Upon colistin withdrawal, an ancestral susceptible strain outcompeted resistant isolates in three of the four cases. In the final case, resistance was also lost, but by a compensatory inactivating mutation in the transcriptional regulator of the pmr locus. Notably, this inactivating mutation reduced the probability of reacquiring colistin resistance when subsequently challenged in vitro. On face value, these results supported an in vivo fitness cost preventing the evolution of stable colistin resistance. However, more careful analysis of WGS data identified genomic evidence for stable colistin resistance undetected by clinical microbiological assays. Transcriptional studies validated this genomic hypothesis, showing increased pmr expression of the initial isolate. Moreover, altering the environmental growth conditions of the clinical assay recapitulated the classification as colistin resistant. Additional targeted sequencing revealed that this isolate evolved undetected in a patient undergoing colistin treatment, and was then transmitted to other hospitalized patients, further demonstrating its stability in the absence of colistin. This study provides a unique window into mutational pathways taken in response to antibiotic pressure in vivo, and demonstrates the potential for genome sequence data to predict resistance phenotypes.

Gene clusters containing multiple similar genomic regions in close proximity are of great interest for biomedical studies because of their associations with inherited diseases. However, such regions are difficult to analyze due to their structural complexity and their complicated evolutionary histories, reflecting a variety of large-scale mutational events. In particular, conversion events can mislead inferences about the relationships among these regions, as traced by traditional methods such as construction of phylogenetic trees or multi-species alignments.

Antibodies of the VRC01 class neutralize HIV-1, arise in diverse HIV-1-infected donors, and are potential templates for an effective HIV-1 vaccine. However, the stochastic processes that generate repertoires in each individual of >10(12) antibodies make elicitation of specific antibodies uncertain. Here we determine the ontogeny of the VRC01 class by crystallography and next-generation sequencing. Despite antibody-sequence differences exceeding 50%, antibody-gp120 cocrystal structures reveal VRC01-class recognition to be remarkably similar. B cell transcripts indicate that VRC01-class antibodies require few specific genetic elements, suggesting that naive-B cells with VRC01-class features are generated regularly by recombination. Virtually all of these fail to mature, however, with only a few-likely one-ancestor B cell expanding to form a VRC01-class lineage in each donor. Developmental similarities in multiple donors thus reveal the generation of VRC01-class antibodies to be reproducible in principle, thereby providing a framework for attempts to elicit similar antibodies in the general population.

Whole genome sequence data for small pedigrees has been shown to provide sufficient information to resolve detailed haplotypes in small pedigrees. Using such information, recombinations can be mapped onto chromosomes, compared with the segregation of a disease of interest and used to filter genome sequence variants. We now show that relatively inexpensive SNP array data from small pedigrees can be used in a similar manner to provide a means of identifying regions of interest in exome sequencing projects. We demonstrate that in those situations where one can assume complete penetrance and parental DNA is available, SNP recombination mapping using Boolean logic identifies chromosomal regions identical to those detected by multipoint linkage using microsatellites but with much less computation. We further show that this approach is successful because the probability of a double crossover between informative SNP loci is negligible. Our observations provide a rationale for using SNP arrays and recombination mapping as a rapid and cost-effective means of incorporating chromosome segregation information into exome sequencing projects intended for disease-gene identification.

Many software tools for comparative analysis of genomic sequence data have been released in recent decades. Despite this, it remains challenging to determine evolutionary relationships in gene clusters due to their complex histories involving duplications, deletions, inversions, and conversions. One concept describing these relationships is orthology. Orthologs derive from a common ancestor by speciation, in contrast to paralogs, which derive from duplication. Discriminating orthologs from paralogs is a necessary step in most multispecies sequence analyses, but doing so accurately is impeded by the occurrence of gene conversion events. We propose a refined method of orthology assignment based on two paradigms for interpreting its definition: by genomic context or by sequence content. X-orthology (based on context) traces orthology resulting from speciation and duplication only, while N-orthology (based on content) includes the influence of conversion events. We developed a computational method for automatically mapping both types of orthology on a per-nucleotide basis in gene cluster regions studied by comparative sequencing, and we make this mapping accessible by visualizing the output. All of these steps are incorporated into our newly extended CHAP 2 package. We evaluate our method using both simulated data and real gene clusters (including the well-characterized α-globin and β-globin clusters). We also illustrate use of CHAP 2 by analyzing four more loci: CCL (chemokine ligand), IFN (interferon), CYP2abf (part of cytochrome P450 family 2), and KIR (killer cell immunoglobulin-like receptors). These new methods facilitate and extend our understanding of evolution at these and other loci by adding automated accurate evolutionary inference to the biologist's toolkit. The CHAP 2 package is freely available from http://www.bx.psu.edu/miller_lab.