Evolution of immune genes is associated with the Black Death.

Infectious diseases are among the strongest selective pressures driving human evolution1,2. This includes the single greatest mortality event in recorded history, the first outbreak of the second pandemic of plague, commonly called the Black Death, which was caused by the bacterium Yersinia pestis3. This pandemic devastated Afro-Eurasia, killing up to 30-50% of the population4. To identify loci that may have been under selection during the Black Death, we characterized genetic variation around immune-related genes from 206 ancient DNA extracts, stemming from two different European populations before, during and after the Black Death. Immune loci are strongly enriched for highly differentiated sites relative to a set of non-immune loci, suggesting positive selection. We identify 245 variants that are highly differentiated within the London dataset, four of which were replicated in an independent cohort from Denmark, and represent the strongest candidates for positive selection. The selected allele for one of these variants, rs2549794, is associated with the production of a full-length (versus truncated) ERAP2 transcript, variation in cytokine response to Y. pestis and increased ability to control intracellular Y. pestis in macrophages. Finally, we show that protective variants overlap with alleles that are today associated with increased susceptibility to autoimmune diseases, providing empirical evidence for the role played by past pandemics in shaping present-day susceptibility to disease.

Helical coiling of metaphase chromatids.

Chromatids of mitotic chromosomes were suggested to coil into a helix in early cytological studies and this assumption was recently supported by chromosome conformation capture (3C) sequencing. Still, direct differential visualization of a condensed chromatin fibre confirming the helical model was lacking. Here, we combined Hi-C analysis of purified metaphase chromosomes, biopolymer modelling and spatial structured illumination microscopy of large fluorescently labeled chromosome segments to reveal the chromonema - a helically-wound, 400 nm thick chromatin thread forming barley mitotic chromatids. Chromatin from adjacent turns of the helix intermingles due to the stochastic positioning of chromatin loops inside the chromonema. Helical turn size varies along chromosome length, correlating with chromatin density. Constraints on the observable dimensions of sister chromatid exchanges further supports the helical chromonema model.

Limitation of current probe design for oligo-cross-FISH, exemplified by chromosome evolution studies in duckweeds.

Duckweeds represent a small, free-floating aquatic family (Lemnaceae) of the monocot order Alismatales with the fastest growth rate among flowering plants. They comprise five genera (Spirodela, Landoltia, Lemna, Wolffiella, and Wolffia) varying in genome size and chromosome number. Spirodela polyrhiza had the first sequenced duckweed genome. Cytogenetic maps are available for both species of the genus Spirodela (S. polyrhiza and S. intermedia). However, elucidation of chromosome homeology and evolutionary chromosome rearrangements by cross-FISH using Spirodela BAC probes to species of other duckweed genera has not been successful so far. We investigated the potential of chromosome-specific oligo-FISH probes to address these topics. We designed oligo-FISH probes specific for one S. intermedia and one S. polyrhiza chromosome (Fig. 1a). Our results show that these oligo-probes cross-hybridize with the homeologous regions of the other congeneric species, but are not suitable to uncover chromosomal homeology across duckweeds genera. This is most likely due to too low sequence similarity between the investigated genera and/or too low probe density on the target genomes. Finally, we suggest genus-specific design of oligo-probes to elucidate chromosome evolution across duckweed genera.

Whole-chromosome paints in maize reveal rearrangements, nuclear domains, and chromosomal relationships.

Whole-chromosome painting probes were developed for each of the 10 chromosomes of maize by producing amplifiable libraries of unique sequences of oligonucleotides that can generate labeled probes through transcription reactions. These paints allow identification of individual homologous chromosomes for many applications as demonstrated in somatic root tip metaphase cells, in the pachytene stage of meiosis, and in interphase nuclei. Several chromosomal aberrations were examined as proof of concept for study of various rearrangements using probes that cover the entire chromosome and that label diverse varieties. The relationship of the supernumerary B chromosome and the normal chromosomes was examined with the finding that there is no detectable homology between any of the normal A chromosomes and the B chromosome. Combined with other chromosome-labeling techniques, a complete set of whole-chromosome oligonucleotide paints lays the foundation for future studies of the structure, organization, and evolution of genomes.

Genetic resiliency and the Black Death: No apparent loss of mitogenomic diversity due to the Black Death in medieval London and Denmark.

OBJECTIVES: In the 14th century AD, medieval Europe was severely affected by the Great European Famine as well as repeated bouts of disease, including the Black Death, causing major demographic shifts. This high volatility led to increased mobility and migration due to new labor and economic opportunities, as evidenced by documentary and stable isotope data. This study uses ancient DNA (aDNA) isolated from skeletal remains to examine whether evidence for large-scale population movement can be gleaned from the complete mitochondrial genomes of 264 medieval individuals from England (London) and Denmark. MATERIALS AND METHODS: Using a novel library-conserving approach to targeted capture, we recovered 264 full mitochondrial genomes from the petrous portion of the temporal bones and teeth and compared genetic diversity across the medieval period within and between English (London) and Danish populations and with contemporary populations through population pairwise ΦST analysis. RESULTS: We find no evidence of significant differences in genetic diversity spatially or temporally in our dataset, yet there is a high degree of haplotype diversity in our medieval samples with little exact sequence sharing. DISCUSSION: The mitochondrial genomes of both medieval Londoners and medieval Danes suggest high mitochondrial diversity before, during and after the Black Death. While our mitochondrial genomic data lack geographically correlated signals, these data could be the result of high, continual female migration before and after the Black Death or may simply indicate a large female effective population size unaffected by the upheaval of the medieval period. Either scenario suggests a genetic resiliency in areas of northwestern medieval Europe.

Eggshell palaeogenomics: Palaeognath evolutionary history revealed through ancient nuclear and mitochondrial DNA from Madagascan elephant bird (Aepyornis sp.) eggshell.

Palaeognaths, the sister group of all other living birds (neognaths), were once considered to be vicariant relics from the breakup of the Gondwanan supercontinent. However, recent molecular studies instead argue for dispersal of volant ancestors across marine barriers. Resolving this debate hinges upon accurately reconstructing their evolutionary relationships and dating their divergences, which often relies on phylogenetic information from extinct relatives and nuclear genomes. Mitogenomes from the extinct elephant birds of Madagascar have helped inform the palaeognath phylogeny; however, nuclear information has remained unavailable. Here, we use ancient DNA (aDNA) extracted from fossil eggshell, together with target enrichment and next-generation sequencing techniques, to reconstruct an additional new mitogenome from Aepyornis sp. with 33.5X coverage. We also recover the first elephant bird nuclear aDNA, represented by 12,500bp of exonic information. While we confirm that elephant birds are sister taxa to the kiwi, our data suggests that, like neognaths, palaeognaths underwent an explosive radiation between 69 and 52Ma-well after the break-up of Gondwana, and more rapidly than previously estimated from mitochondrial data alone. These results further support the idea that ratites primarily diversified immediately following the Cretaceous-Palaeogene mass extinction and convergently evolved flightlessness. Our study reinforces the importance of including information from the nuclear genome of extinct taxa for recovering deep evolutionary relationships. Furthermore, with approximately 3% endogenous aDNA retrieved, avian eggshell can be a valuable substrate for recovering high quality aDNA. We suggest that elephant bird whole genome recovery is ultimately achievable, and will provide future insights into the evolution these birds.

Ancient whole genome enrichment using baits built from modern DNA.

We report metrics from complete genome capture of nuclear DNA from extinct mammoths using biotinylated RNAs transcribed from an Asian elephant DNA extract. Enrichment of the nuclear genome ranged from 1.06- to 18.65-fold, to an apparent maximum threshold of ∼80% on-target. This projects an order of magnitude less costly complete genome sequencing from long-dead organisms, even when a reference genome is unavailable for bait design.

Versatile design and synthesis platform for visualizing genomes with Oligopaint FISH probes.

A host of observations demonstrating the relationship between nuclear architecture and processes such as gene expression have led to a number of new technologies for interrogating chromosome positioning. Whereas some of these technologies reconstruct intermolecular interactions, others have enhanced our ability to visualize chromosomes in situ. Here, we describe an oligonucleotide- and PCR-based strategy for fluorescence in situ hybridization (FISH) and a bioinformatic platform that enables this technology to be extended to any organism whose genome has been sequenced. The oligonucleotide probes are renewable, highly efficient, and able to robustly label chromosomes in cell culture, fixed tissues, and metaphase spreads. Our method gives researchers precise control over the sequences they target and allows for single and multicolor imaging of regions ranging from tens of kilobases to megabases with the same basic protocol. We anticipate this technology will lead to an enhanced ability to visualize interphase and metaphase chromosomes.

Detection and characterization of human pathogenic viruses circulating in community wastewater using multi target microarrays and polymerase chain reaction.

Sewage pollution remains the most significant source of human waterborne pathogens. This study describes the detection and characterization of human enteric viruses in community wastewaters using cell culture coupled with multiple target microarrays (with a total of 780 unique probes targeting 27 different groups of both DNA and RNA viruses) and polymerase chain reaction (PCR) assays. Over a 13-month sampling period, RNA viruses (astroviruses and enteroviruses) were more frequently detected compared to DNA viruses (adenoviruses, particularly type 41 and BK polyomavirus). Overall, many more viruses were shed during the winter months (December-February) compared to the summer months. Exploration of the multiple types of enteric viruses particularly in winter months identified much more significant prevalence of key viral pathogens associated with sewage pollution of the water environment than previously realized and seasonal disinfection used in some parts of the world may lead to a seeding of ambient waters. Molecular characterization of pathogenic viruses in community wastewater will improve the understanding of the potential risk of waterborne disease transmission of viral pathogens.

Reply to Barton et al: signatures of natural selection during the Black Death.

Barton et al.1 raise several statistical concerns regarding our original analyses2 that highlight the challenge of inferring natural selection using ancient genomic data. We show here that these concerns have limited impact on our original conclusions. Specifically, we recover the same signature of enrichment for high FST values at the immune loci relative to putatively neutral sites after switching the allele frequency estimation method to a maximum likelihood approach, filtering to only consider known human variants, and down-sampling our data to the same mean coverage across sites. Furthermore, using permutations, we show that the rs2549794 variant near ERAP2 continues to emerge as the strongest candidate for selection (p = 1.2×10-5), falling below the Bonferroni-corrected significance threshold recommended by Barton et al. Importantly, the evidence for selection on ERAP2 is further supported by functional data demonstrating the impact of the ERAP2 genotype on the immune response to Y. pestis and by epidemiological data from an independent group showing that the putatively selected allele during the Black Death protects against severe respiratory infection in contemporary populations.

Evolution combined with genomic study elucidates genetic bases of isobutanol tolerance in Escherichia coli.

BACKGROUND: Isobutanol is a promising next-generation biofuel with demonstrated high yield microbial production, but the toxicity of this molecule reduces fermentation volumetric productivity and final titer. Organic solvent tolerance is a complex, multigenic phenotype that has been recalcitrant to rational engineering approaches. We apply experimental evolution followed by genome resequencing and a gene expression study to elucidate genetic bases of adaptation to exogenous isobutanol stress. RESULTS: The adaptations acquired in our evolved lineages exhibit antagonistic pleiotropy between minimal and rich medium, and appear to be specific to the effects of longer chain alcohols. By examining genotypic adaptation in multiple independent lineages, we find evidence of parallel evolution in marC, hfq, mdh, acrAB, gatYZABCD, and rph genes. Many isobutanol tolerant lineages show reduced RpoS activity, perhaps related to mutations in hfq or acrAB. Consistent with the complex, multigenic nature of solvent tolerance, we observe adaptations in a diversity of cellular processes. Many adaptations appear to involve epistasis between different mutations, implying a rugged fitness landscape for isobutanol tolerance. We observe a trend of evolution targeting post-transcriptional regulation and high centrality nodes of biochemical networks. Collectively, the genotypic adaptations we observe suggest mechanisms of adaptation to isobutanol stress based on remodeling the cell envelope and surprisingly, stress response attenuation. CONCLUSIONS: We have discovered a set of genotypic adaptations that confer increased tolerance to exogenous isobutanol stress. Our results are immediately useful to further efforts to engineer more isobutanol tolerant host strains of E. coli for isobutanol production. We suggest that rpoS and post-transcriptional regulators, such as hfq, RNA helicases, and sRNAs may be interesting mutagenesis targets for future global phenotype engineering.

OligoArrayDb: pangenomic oligonucleotide microarray probe sets database.

OligoArrayDb is a comprehensive database containing pangenomic oligonucleotide microarray probe sets designed for most of the sequenced genomes that are not covered by commercial catalog arrays. The availability of probe sequences, associated with custom microarray fabrication services offered by many companies and cores presents the unequalled possibility to perform microarray experiments on most of the sequenced organisms. OligoArrayDb contains more than 2.8 probes per gene in average for more than 600 organisms, mostly archaea and bacteria strains available from public database. On average, 98% of the annotated genes have at least one probe which is predicted to be specific to its intended target in >94% of the cases. OligoArrayDb is weekly updated as new sequenced genomes become available. Probe sequences, in addition to a comprehensive set of annotations can be downloaded from this database. OligoArrayDb is publicly accessible online at http://berry.engin.umich.edu/oligoarraydb.

Meiotic crossovers characterized by haplotype-specific chromosome painting in maize.

Meiotic crossovers (COs) play a critical role in generating genetic variation and maintaining faithful segregation of homologous chromosomes during meiosis. We develop a haplotype-specific fluorescence in situ hybridization (FISH) technique that allows visualization of COs directly on metaphase chromosomes. Oligonucleotides (oligos) specific to chromosome 10 of maize inbreds B73 and Mo17, respectively, are synthesized and labeled as FISH probes. The parental and recombinant chromosome 10 in B73 x Mo17 F1 hybrids and F2 progenies can be unambiguously identified by haplotype-specific FISH. Analysis of 58 F2 plants reveals lack of COs in the entire proximal half of chromosome 10. However, we detect COs located in regions very close to the centromere in recombinant inbred lines from an intermated B73 x Mo17 population, suggesting effective accumulation of COs in recombination-suppressed chromosomal regions through intermating and the potential to generate favorable allelic combinations of genes residing in these regions.

In situ-synthesized virulence and marker gene biochip for detection of bacterial pathogens in water.

Pathogen detection tools with high reliability are needed for various applications, including food and water safety and clinical diagnostics. In this study, we designed and validated an in situ-synthesized biochip for detection of 12 microbial pathogens, including a suite of pathogens relevant to water safety. To enhance the reliability of presence/absence calls, probes were designed for multiple virulence and marker genes (VMGs) of each pathogen, and each VMG was targeted by an average of 17 probes. Hybridization of the biochip with amplicon mixtures demonstrated that 95% of the initially designed probes behaved as predicted in terms of positive/negative signals. The probes were further validated using DNA obtained from three different types of water samples and spiked with pathogen genomic DNA at decreasing relative abundance. Excellent specificity for making presence/absence calls was observed by using a cutoff of 0.5 for the positive fraction (i.e., the fraction of probes yielding a positive signal for a given VMG). A split multiplex PCR design for simultaneous amplification of the VMGs resulted in a detection limit of between 0.1 and 0.01% relative abundance, depending on the type of pathogen and the VMG. Thermodynamic analysis of the hybridization patterns obtained with DNA from the different water samples demonstrated that probes with a hybridization Gibbs free energy of approximately -19.3 kcal/mol provided the best trade-off between sensitivity and specificity. The developed biochip may be used to detect the described bacterial pathogens in water samples when parallel and specific detection is required.

On-chip non-equilibrium dissociation curves and dissociation rate constants as methods to assess specificity of oligonucleotide probes.

Nucleic acid hybridization serves as backbone for many high-throughput systems for detection, expression analysis, comparative genomics and re-sequencing. Specificity of hybridization between probes and intended targets is always critical. Approaches to ensure and evaluate specificity include use of mismatch probes, obtaining dissociation curves rather than single temperature hybridizations, and comparative hybridizations. In this study, we quantify effects of mismatch type and position on intensity of hybridization signals and provide a new approach based on dissociation rate constants to evaluate specificity of hybridized signals in complex target mixtures. Using an extensive set of 18mer oligonucleotide probes on an in situ synthesized biochip platform, we demonstrate that mismatches in the center of the probe are more discriminating than mismatches toward the extremities of the probe and mismatches toward the attached end are less discriminating than those toward the loose end. The observed destabilizing effect of a mismatch type agreed in general with predictions using the nearest neighbor model. Use of a new parameter, specific dissociation temperature (T(d-w), temperature of maximum specific dissociation rate constant), obtained from probe-target duplex dissociation profiles considerably improved the evaluation of specificity. These results have broad implications for hybridization data obtained from complex mixtures of nucleic acids.

An open-access long oligonucleotide microarray resource for analysis of the human and mouse transcriptomes.

Two collections of oligonucleotides have been designed for preparing pangenomic human and mouse microarrays. A total of 148,993 and 121,703 oligonucleotides were designed against human and mouse transcripts. Quality scores were created in order to select 25,342 human and 24,109 mouse oligonucleotides. They correspond to: (i) a BLAST-specificity score; (ii) the number of expressed sequence tags matching each probe; (iii) the distance to the 3' end of the target mRNA. Scores were also used to compare in silico the two microarrays with commercial microarrays. The sets described here, called RNG/MRC collections, appear at least as specific and sensitive as those from the commercial platforms. The RNG/MRC collections have now been used by an Anglo-French consortium to distribute more than 3500 microarrays to the academic community. Ad hoc identification of tissue-specific transcripts and a approximately 80% correlation with hybridizations performed on Affymetrix GeneChiptrade mark suggest that the RNG/MRC microarrays perform well. This work provides a comprehensive open resource for investigators working on human and mouse transcriptomes, as well as a generic method to generate new microarray collections in other organisms. All information related to these probes, as well as additional information about commercial microarrays have been stored in a freely-accessible database called MEDIANTE.

Comparative Oligo-FISH Mapping: An Efficient and Powerful Methodology To Reveal Karyotypic and Chromosomal Evolution.

Developing the karyotype of a eukaryotic species relies on identification of individual chromosomes, which has been a major challenge for most nonmodel plant and animal species. We developed a novel chromosome identification system by selecting and labeling oligonucleotides (oligos) located in specific regions on every chromosome. We selected a set of 54,672 oligos (45 nt) based on single copy DNA sequences in the potato genome. These oligos generated 26 distinct FISH signals that can be used as a "bar code" or "banding pattern" to uniquely label each of the 12 chromosomes from both diploid and polyploid (4× and 6×) potato species. Remarkably, the same bar code can be used to identify the 12 homeologous chromosomes among distantly related Solanum species, including tomato and eggplant. Accurate karyotypes based on individually identified chromosomes were established in six Solanum species that have diverged for >15 MY. These six species have maintained a similar karyotype; however, modifications to the FISH signal bar code led to the discovery of two reciprocal chromosomal translocations in Solanum etuberosum and S. caripense We also validated these translocations by oligo-based chromosome painting. We demonstrate that the oligo-based FISH techniques are powerful new tools for chromosome identification and karyotyping research, especially for nonmodel plant species.

The evolutionary history of dogs in the Americas.

Dogs were present in the Americas before the arrival of European colonists, but the origin and fate of these precontact dogs are largely unknown. We sequenced 71 mitochondrial and 7 nuclear genomes from ancient North American and Siberian dogs from time frames spanning ~9000 years. Our analysis indicates that American dogs were not derived from North American wolves. Instead, American dogs form a monophyletic lineage that likely originated in Siberia and dispersed into the Americas alongside people. After the arrival of Europeans, native American dogs almost completely disappeared, leaving a minimal genetic legacy in modern dog populations. The closest detectable extant lineage to precontact American dogs is the canine transmissible venereal tumor, a contagious cancer clone derived from an individual dog that lived up to 8000 years ago.

Cytophobic surface modification of microfluidic arrays for in situ parallel peptide synthesis and cell adhesion assays.

A combination of PEG-based surface passivation techniques and spatially addressable SPPS (solid-phase peptide synthesis) was used to demonstrate a highly specific cell-peptide adhesion assay on a microfluidic platform. The surface of a silicon-glass microchip was modified to form a mixed self-assembled monolayer that presented PEG moieties interspersed with reactive amino terminals. The PEG provided biomolecular inertness and the reactive amino groups were used for consequent peptide synthesis. The cytophobicity of the surface was characterized by on-chip fluorescent binding assays and was found to be resistant to nonspecific attachment of cells and proteins. An integrated system for parallel peptide synthesis on this reactive amino surface was developed using photogenerated acid chemistry and digital microlithography. A constant synthesis efficiency of >98% was observed for up to 7mer peptides. To demonstrate specific cell adhesion on these synthetic peptide arrays, variations of a 7mer cell binding peptide that binds to murine B lymphoma cells were synthesized. Sequence-specific binding was observed on incubation with fluorescently labeled, intact murine B lymphoma cells, and key residues for binding were identified by deletional analysis.

Diagnostic microarray for 14 water and foodborne pathogens using a flatbed scanner.

Parallel detection approaches are of interest to many researchers interested in identifying multiple water and foodborne pathogens simultaneously. Availability and cost-effectiveness are two key factors determining the usefulness of such approaches for laboratories with limited resources. In this study, we developed and validated a high-density microarray for simultaneous screening of 14 bacterial pathogens using an approach that employs gold labeling with silver enhancement (GLS) protocol. In total, 8887 probes (50-mer) were designed using an in-house database of virulence and marker genes (VMGs), and synthesized in quadruplicate on glass slides using an in-situ synthesis technology. Target VMG amplicons were obtained using multiplex polymerase chain reaction (PCR), labeled with biotin, and hybridized to the microarray. The signals generated after gold deposition and silver enhancement, were quantified using a flatbed scanner having 2-µm resolution. Data analysis indicated that reliable presence/absence calls could be made, if: i) over four probes were used per gene, ii) the signal-to-noise ratio (SNR) cutoff was greater than or equal to two, and iii) the positive fraction (PF), i.e., number of probes with SNR≥2 for a given VMG was greater than 0.75. Hybridization of the array with blind samples resulted in 100% correct calls, and no false positive. Because amplicons were obtained by multiplex PCR, sensitivity of this method is similar to PCR. This assay is an inexpensive and reliable technique for high throughput screening of multiple pathogens.