SciCombinator

The recent genealogical history of human populations is a complex mosaic formed by individual migration, large-scale population movements, and other demographic events. Population genomics datasets can provide a window into this recent history, as rare traces of recent shared genetic ancestry are detectable due to long segments of shared genomic material. We make use of genomic data for 2,257 Europeans (in the Population Reference Sample [POPRES] dataset) to conduct one of the first surveys of recent genealogical ancestry over the past 3,000 years at a continental scale. We detected 1.9 million shared long genomic segments, and used the lengths of these to infer the distribution of shared ancestors across time and geography. We find that a pair of modern Europeans living in neighboring populations share around 2-12 genetic common ancestors from the last 1,500 years, and upwards of 100 genetic ancestors from the previous 1,000 years. These numbers drop off exponentially with geographic distance, but since these genetic ancestors are a tiny fraction of common genealogical ancestors, individuals from opposite ends of Europe are still expected to share millions of common genealogical ancestors over the last 1,000 years. There is also substantial regional variation in the number of shared genetic ancestors. For example, there are especially high numbers of common ancestors shared between many eastern populations that date roughly to the migration period (which includes the Slavic and Hunnic expansions into that region). Some of the lowest levels of common ancestry are seen in the Italian and Iberian peninsulas, which may indicate different effects of historical population expansions in these areas and/or more stably structured populations. Population genomic datasets have considerable power to uncover recent demographic history, and will allow a much fuller picture of the close genealogical kinship of individuals across the world.

Although most eukaryotes reproduce sexually at some moment of their life cycle, as much as a fifth of fungal species were thought to reproduce exclusively asexually. Nevertheless, recent studies have revealed the occurrence of sex in some of these supposedly asexual species. For industrially relevant fungi, for which inoculums are produced by clonal-subcultures since decades, the potentiality for sex is of great interest for strain improvement strategies. Here, we investigated the sexual capability of the fungus Penicillium roqueforti, used as starter for blue cheese production. We present indirect evidence suggesting that recombination could be occurring in this species. The screening of a large sample of strains isolated from diverse substrates throughout the world revealed the existence of individuals of both mating types, even in the very same cheese. The MAT genes, involved in fungal sexual compatibility, appeared to evolve under purifying selection, suggesting that they are still functional. The examination of the recently sequenced genome of the FM 164 cheese strain enabled the identification of the most important genes known to be involved in meiosis, which were found to be highly conserved. Linkage disequilibria were not significant among three of the six marker pairs and 11 out of the 16 possible allelic combinations were found in the dataset. Finally, the detection of signatures of repeat induced point mutations (RIP) in repeated sequences and transposable elements reinforces the conclusion that P. roqueforti underwent more or less recent sex events. In this species of high industrial importance, the induction of a sexual cycle would open the possibility of generating new genotypes that would be extremely useful to diversify cheese products.

Recent advances in whole-genome sequencing have brought the vision of personal genomics and genomic medicine closer to reality. However, current methods lack clinical accuracy and the ability to describe the context (haplotypes) in which genome variants co-occur in a cost-effective manner. Here we describe a low-cost DNA sequencing and haplotyping process, long fragment read (LFR) technology, which is similar to sequencing long single DNA molecules without cloning or separation of metaphase chromosomes. In this study, ten LFR libraries were made using only ∼100 picograms of human DNA per sample. Up to 97% of the heterozygous single nucleotide variants were assembled into long haplotype contigs. Removal of false positive single nucleotide variants not phased by multiple LFR haplotypes resulted in a final genome error rate of 1 in 10 megabases. Cost-effective and accurate genome sequencing and haplotyping from 10-20 human cells, as demonstrated here, will enable comprehensive genetic studies and diverse clinical applications.

The unicellular eukaryote Tetrahymena thermophila has seven mating types. Cells can mate only when they recognize cells of a different mating type as non-self. As a ciliate, Tetrahymena separates its germline and soma into two nuclei. During growth the somatic nucleus is responsible for all gene transcription while the germline nucleus remains silent. During mating, a new somatic nucleus is differentiated from a germline nucleus and mating type is decided by a stochastic process. We report here that the somatic mating type locus contains a pair of genes arranged head-to-head. Each gene encodes a mating type-specific segment and a transmembrane domain that is shared by all mating types. Somatic gene knockouts showed both genes are required for efficient non-self recognition and successful mating, as assessed by pair formation and progeny production. The germline mating type locus consists of a tandem array of incomplete gene pairs representing each potential mating type. During mating, a complete new gene pair is assembled at the somatic mating type locus; the incomplete genes of one gene pair are completed by joining to gene segments at each end of germline array. All other germline gene pairs are deleted in the process. These programmed DNA rearrangements make this a fascinating system of mating type determination.

The evolutionary importance of hybridization and introgression has long been debated. Hybrids are usually rare and unfit, but even infrequent hybridization can aid adaptation by transferring beneficial traits between species. Here we use genomic tools to investigate introgression in Heliconius, a rapidly radiating genus of neotropical butterflies widely used in studies of ecology, behaviour, mimicry and speciation. We sequenced the genome of Heliconius melpomene and compared it with other taxa to investigate chromosomal evolution in Lepidoptera and gene flow among multiple Heliconius species and races. Among 12,669 predicted genes, biologically important expansions of families of chemosensory and Hox genes are particularly noteworthy. Chromosomal organization has remained broadly conserved since the Cretaceous period, when butterflies split from the Bombyx (silkmoth) lineage. Using genomic resequencing, we show hybrid exchange of genes between three co-mimics, Heliconius melpomene, Heliconius timareta and Heliconius elevatus, especially at two genomic regions that control mimicry pattern. We infer that closely related Heliconius species exchange protective colour-pattern genes promiscuously, implying that hybridization has an important role in adaptive radiation.

Inter-individual variation in facial shape is one of the most noticeable phenotypes in humans, and it is clearly under genetic regulation; however, almost nothing is known about the genetic basis of normal human facial morphology. We therefore conducted a genome-wide association study for facial shape phenotypes in multiple discovery and replication cohorts, considering almost ten thousand individuals of European descent from several countries. Phenotyping of facial shape features was based on landmark data obtained from three-dimensional head magnetic resonance images (MRIs) and two-dimensional portrait images. We identified five independent genetic loci associated with different facial phenotypes, suggesting the involvement of five candidate genes-PRDM16, PAX3, TP63, C5orf50, and COL17A1-in the determination of the human face. Three of them have been implicated previously in vertebrate craniofacial development and disease, and the remaining two genes potentially represent novel players in the molecular networks governing facial development. Our finding at PAX3 influencing the position of the nasion replicates a recent GWAS of facial features. In addition to the reported GWA findings, we established links between common DNA variants previously associated with NSCL/P at 2p21, 8q24, 13q31, and 17q22 and normal facial-shape variations based on a candidate gene approach. Overall our study implies that DNA variants in genes essential for craniofacial development contribute with relatively small effect size to the spectrum of normal variation in human facial morphology. This observation has important consequences for future studies aiming to identify more genes involved in the human facial morphology, as well as for potential applications of DNA prediction of facial shape such as in future forensic applications.

Obesity is associated with vitamin D deficiency, and both are areas of active public health concern. We explored the causality and direction of the relationship between body mass index (BMI) and 25-hydroxyvitamin D [25(OH)D] using genetic markers as instrumental variables (IVs) in bi-directional Mendelian randomization (MR) analysis.

BACKGROUND: The only known albino gorilla, named Snowflake, was a male wild born individual from Equatorial Guinea who lived at the Barcelona Zoo for almost 40 years. He was diagnosed with non-syndromic oculocutaneous albinism, i.e. white hair, light eyes, pink skin, photophobia and reduced visual acuity. Despite previous efforts to explain the genetic cause, this is still unknown. Here, we study the genetic cause of his albinism and making use of whole genome sequencing data we find a higher inbreeding coefficient compared to other gorillas. RESULTS: We successfully identified the causal genetic variant for Snowflake’s albinism, a non-synonymous single nucleotide variant located in a transmembrane region of SLC45A2. This transporter is known to be involved in oculocutaneous albinism type 4 (OCA4) in humans. We provide experimental evidence that shows that this amino acid replacement alters the membrane spanning capability of this transmembrane region. Finally, we provide a comprehensive study of genome-wide patterns of autozygogosity revealing that Snowflake’s parents were related, being this the first report of inbreeding in a wild born Western lowland gorilla. CONCLUSIONS: In this study we demonstrate how the use of whole genome sequencing can be extended to link genotype and phenotype in non-model organisms and it can be a powerful tool in conservation genetics (e.g., inbreeding and genetic diversity) with the expected decrease in sequencing cost.

Understanding the evolutionary history of microbial pathogens is critical for mitigating the impacts of emerging infectious diseases on economically and ecologically important host species. We used a genome resequencing approach to resolve the evolutionary history of an important microbial pathogen, the chytrid Batrachochytrium dendrobatidis (Bd), which has been implicated in amphibian declines worldwide. We sequenced the genomes of 29 isolates of Bd from around the world, with an emphasis on North, Central, and South America because of the devastating effect that Bd has had on amphibian populations in the New World. We found a substantial amount of evolutionary complexity in Bd with deep phylogenetic diversity that predates observed global amphibian declines. By investigating the entire genome, we found that even the most recently evolved Bd clade (termed the global panzootic lineage) contained more genetic variation than previously reported. We also found dramatic differences among isolates and among genomic regions in chromosomal copy number and patterns of heterozygosity, suggesting complex and heterogeneous genome dynamics. Finally, we report evidence for selection acting on the Bd genome, supporting the hypothesis that protease genes are important in evolutionary transitions in this group. Bd is considered an emerging pathogen because of its recent effects on amphibians, but our data indicate that it has a complex evolutionary history that predates recent disease outbreaks. Therefore, it is important to consider the contemporary effects of Bd in a broader evolutionary context and identify specific mechanisms that may have led to shifts in virulence in this system.

Evolutionary novelties have been important in the history of life, but their origins are usually difficult to examine in detail. We previously described the evolution of a novel trait, aerobic citrate utilization (Cit(+)), in an experimental population of Escherichia coli. Here we analyse genome sequences to investigate the history and genetic basis of this trait. At least three distinct clades coexisted for more than 10,000 generations before its emergence. The Cit(+) trait originated in one clade by a tandem duplication that captured an aerobically expressed promoter for the expression of a previously silent citrate transporter. The clades varied in their propensity to evolve this novel trait, although genotypes able to do so existed in all three clades, implying that multiple potentiating mutations arose during the population’s history. Our findings illustrate the importance of promoter capture and altered gene regulation in mediating the exaptation events that often underlie evolutionary innovations.

Data from the 1000 genomes project (1KGP) and Complete Genomics (CG) have dramatically increased the numbers of known genetic variants and challenge several assumptions about the reference genome and its uses in both clinical and research settings. Specifically, 34% of published array-based GWAS studies for a variety of diseases utilize probes that overlap unanticipated single nucleotide polymorphisms (SNPs), indels, or structural variants. Linkage disequilibrium (LD) block length depends on the numbers of markers used, and the mean LD block size decreases from 16 kb to 7 kb,when HapMap-based calculations are compared to blocks computed from1KGP data. Additionally, when 1KGP and CG variants are compared, 19% of the single nucleotide variants (SNVs) reported from common genomes are unique to one dataset; likely a result of differences in data collection methodology, alignment of reads to the reference genome, and variant-calling algorithms. Together these observations indicate that current research resources and informatics methods do not adequately account for the high level of variation that already exists in the human population and significant efforts are needed to create resources that can accurately assess personal genomics for health, disease, and predict treatment outcomes.

To correlate the variable clinical features of oestrogen-receptor-positive breast cancer with somatic alterations, we studied pretreatment tumour biopsies accrued from patients in two studies of neoadjuvant aromatase inhibitor therapy by massively parallel sequencing and analysis. Eighteen significantly mutated genes were identified, including five genes (RUNX1, CBFB, MYH9, MLL3 and SF3B1) previously linked to haematopoietic disorders. Mutant MAP3K1 was associated with luminal A status, low-grade histology and low proliferation rates, whereas mutant TP53 was associated with the opposite pattern. Moreover, mutant GATA3 correlated with suppression of proliferation upon aromatase inhibitor treatment. Pathway analysis demonstrated that mutations in MAP2K4, a MAP3K1 substrate, produced similar perturbations as MAP3K1 loss. Distinct phenotypes in oestrogen-receptor-positive breast cancer are associated with specific patterns of somatic mutations that map into cellular pathways linked to tumour biology, but most recurrent mutations are relatively infrequent. Prospective clinical trials based on these findings will require comprehensive genome sequencing.

ABSTRACT A novel human coronavirus (HCoV-EMC/2012) was isolated from a man with acute pneumonia and renal failure in June 2012. This report describes the complete genome sequence, genome organization, and expression strategy of HCoV-EMC/2012 and its relation with known coronaviruses. The genome contains 30,119 nucleotides and contains at least 10 predicted open reading frames, 9 of which are predicted to be expressed from a nested set of seven subgenomic mRNAs. Phylogenetic analysis of the replicase gene of coronaviruses with completely sequenced genomes showed that HCoV-EMC/2012 is most closely related to Tylonycteris bat coronavirus HKU4 (BtCoV-HKU4) and Pipistrellus bat coronavirus HKU5 (BtCoV-HKU5), which prototype two species in lineage C of the genus Betacoronavirus. In accordance with the guidelines of the International Committee on Taxonomy of Viruses, and in view of the 75% and 77% amino acid sequence identity in 7 conserved replicase domains with BtCoV-HKU4 and BtCoV-HKU5, respectively, we propose that HCoV-EMC/2012 prototypes a novel species in the genus Betacoronavirus. HCoV-EMC/2012 may be most closely related to a coronavirus detected in Pipistrellus pipistrellus in The Netherlands, but because only a short sequence from the most conserved part of the RNA-dependent RNA polymerase-encoding region of the genome was reported for this bat virus, its genetic distance from HCoV-EMC remains uncertain. HCoV-EMC/2012 is the sixth coronavirus known to infect humans and the first human virus within betacoronavirus lineage C. IMPORTANCE Coronaviruses are capable of infecting humans and many animal species. Most infections caused by human coronaviruses are relatively mild. However, the outbreak of severe acute respiratory syndrome (SARS) caused by SARS-CoV in 2002 to 2003 and the fatal infection of a human by HCoV-EMC/2012 in 2012 show that coronaviruses are able to cause severe, sometimes fatal disease in humans. We have determined the complete genome of HCoV-EMC/2012 using an unbiased virus discovery approach involving next-generation sequencing techniques, which enabled subsequent state-of-the-art bioinformatics, phylogenetics, and taxonomic analyses. By establishing its complete genome sequence, HCoV-EMC/2012 was characterized as a new genotype which is closely related to bat coronaviruses that are distant from SARS-CoV. We expect that this information will be vital to rapid advancement of both clinical and vital research on this emerging pathogen.

Many groups, including our own, have proposed the use of DNA methylation profiles as biomarkers for various disease states. While much research has been done identifying DNA methylation signatures in cancer vs. normal etc., we still lack sufficient knowledge of the role that differential methylation plays during normal cellular differentiation and tissue specification. We also need thorough, genome level studies to determine the meaning of methylation of individual CpG dinucleotides in terms of gene expression.

ABSTRACT Horizontal gene transfer (HGT) is largely responsible for increasing the incidence of antibiotic-resistant infections worldwide. While studies have focused on HGT in vivo, this work investigates whether the ability of pathogens to persist in the environment, particularly on touch surfaces, may also play an important role. Escherichia coli, virulent clone ST131, and Klebsiella pneumoniae harboring extended-spectrum-β-lactamase (ESBL) bla(CTX-M-15) and metallo-β-lactamase bla(NDM-1), respectively, exhibited prolonged survival on stainless steel, with approximately 10(4) viable cells remaining from an inoculum of 10(7) CFU per cm(2) after 1 month at 21°C. HGT of bla to an antibiotic-sensitive but azide-resistant recipient E. coli strain occurred on stainless steel dry touch surfaces and in suspension but not on dry copper. The conjugation frequency was approximately 10 to 50 times greater and occurred immediately, and resulting transconjugants were more stable with ESBL E. coli as the donor cell than with K. pneumoniae, but bla(NDM-1) transfer increased with time. Transconjugants also exhibited the same resistance profile as the donor, suggesting multiple gene transfer. Rapid death, inhibition of respiration, and destruction of genomic and plasmid DNA of both pathogens occurred on copper alloys accompanied by a reduction in bla copy number. Naked E. coli DNA degraded on copper at 21°C and 37°C but slowly at 4°C, suggesting a direct role for the metal. Persistence of viable pathogenic bacteria on touch surfaces may not only increase the risk of infection transmission but may also contribute to the spread of antibiotic resistance by HGT. The use of copper alloys as antimicrobial touch surfaces may help reduce infection and HGT. IMPORTANCE Horizontal gene transfer (HGT) conferring resistance to many classes of antimicrobials has resulted in a worldwide epidemic of nosocomial and community infections caused by multidrug-resistant microorganisms, leading to suggestions that we are in effect returning to the preantibiotic era. While studies have focused on HGT in vivo, this work investigates whether the ability of pathogens to persist in the environment, particularly on touch surfaces, may also play an important role. Here we show prolonged (several-week) survival of multidrug-resistant Escherichia coli and Klebsiella pneumoniae on stainless steel surfaces. Plasmid-mediated HGT of β-lactamase genes to an azide-resistant recipient E. coli strain occurred when the donor and recipient cells were mixed together on stainless steel and in suspension but not on copper surfaces. In addition, rapid death of both antibiotic-resistant strains and destruction of plasmid and genomic DNA were observed on copper and copper alloy surfaces, which could be useful in the prevention of infection spread and gene transfer.

Understanding of soil processes is essential for addressing the global issues of food security, disease transmission and climate change. However, techniques for observing soil biology are lacking. We present a heterogeneous, porous, transparent substrate for in situ 3D imaging of living plants and root-associated microorganisms using particles of the transparent polymer, Nafion, and a solution with matching optical properties. Minerals and fluorescent dyes were adsorbed onto the Nafion particles for nutrient supply and imaging of pore size and geometry. Plant growth in transparent soil was similar to that in soil. We imaged colonization of lettuce roots by the human bacterial pathogen Escherichia coli O157:H7 showing micro-colony development. Micro-colonies may contribute to bacterial survival in soil. Transparent soil has applications in root biology, crop genetics and soil microbiology.

Pest and pathogen losses jeopardise global food security and ever since the 19(th) century Irish famine, potato late blight has exemplified this threat. The causal oomycete pathogen, Phytophthora infestans, undergoes major population shifts in agricultural systems via the successive emergence and migration of asexual lineages. The phenotypic and genotypic bases of these selective sweeps are largely unknown but management strategies need to adapt to reflect the changing pathogen population. Here, we used molecular markers to document the emergence of a lineage, termed 13_A2, in the European P. infestans population, and its rapid displacement of other lineages to exceed 75% of the pathogen population across Great Britain in less than three years. We show that isolates of the 13_A2 lineage are among the most aggressive on cultivated potatoes, outcompete other aggressive lineages in the field, and overcome previously effective forms of plant host resistance. Genome analyses of a 13_A2 isolate revealed extensive genetic and expression polymorphisms particularly in effector genes. Copy number variations, gene gains and losses, amino-acid replacements and changes in expression patterns of disease effector genes within the 13_A2 isolate likely contribute to enhanced virulence and aggressiveness to drive this population displacement. Importantly, 13_A2 isolates carry intact and in planta induced Avrblb1, Avrblb2 and Avrvnt1 effector genes that trigger resistance in potato lines carrying the corresponding R immune receptor genes Rpi-blb1, Rpi-blb2, and Rpi-vnt1.1. These findings point towards a strategy for deploying genetic resistance to mitigate the impact of the 13_A2 lineage and illustrate how pathogen population monitoring, combined with genome analysis, informs the management of devastating disease epidemics.

Spermatogenesis is a complex process reliant upon interactions between germ cells (GC) and supporting somatic cells. Testicular Sertoli cells (SC) support GCs during maturation through physical attachment, the provision of nutrients, and protection from immunological attack. This role is facilitated by an active cytoskeleton of parallel microtubule arrays that permit transport of nutrients to GCs, as well as translocation of spermatids through the seminiferous epithelium during maturation. It is well established that chemical perturbation of SC microtubule remodelling leads to premature GC exfoliation demonstrating that microtubule remodelling is an essential component of male fertility, yet the genes responsible for this process remain unknown. Using a random ENU mutagenesis approach, we have identified a novel mouse line displaying male-specific infertility, due to a point mutation in the highly conserved ATPase domain of the novel KATANIN p60-related microtubule severing protein Katanin p60 subunit A-like1 (KATNAL1). We demonstrate that Katnal1 is expressed in testicular Sertoli cells (SC) from 15.5 days post-coitum (dpc) and that, consistent with chemical disruption models, loss of function of KATNAL1 leads to male-specific infertility through disruption of SC microtubule dynamics and premature exfoliation of spermatids from the seminiferous epithelium. The identification of KATNAL1 as an essential regulator of male fertility provides a significant novel entry point into advancing our understanding of how SC microtubule dynamics promotes male fertility. Such information will have resonance both for future treatment of male fertility and the development of non-hormonal male contraceptives.

A recent slew of ENCODE Consortium publications, specifically the article signed by all Consortium members, put forward the idea that more than 80% of the human genome is functional. This claim flies in the face of current estimates according to which the fraction of the genome that is evolutionarily conserved through purifying selection is under 10%. Thus, according to the ENCODE Consortium, a biological function can be maintained indefinitely without selection, which implies that at least 80 - 10 = 70% of the genome is perfectly invulnerable to deleterious mutations, either because no mutation can ever occur in these “functional” regions, or because no mutation in these regions can ever be deleterious. This absurd conclusion was reached through various means, chiefly (1) by employing the seldom used “causal role” definition of biological function and then applying it inconsistently to different biochemical properties, (2) by committing a logical fallacy known as “affirming the consequent,” (3) by failing to appreciate the crucial difference between “junk DNA” and “garbage DNA,” (4) by using analytical methods that yield biased errors and inflate estimates of functionality, (5) by favoring statistical sensitivity over specificity, and (6) by emphasizing statistical significance rather than the magnitude of the effect. Here, we detail the many logical and methodological transgressions involved in assigning functionality to almost every nucleotide in the human genome. The ENCODE results were predicted by one of its authors to necessitate the rewriting of textbooks. We agree, many textbooks dealing with marketing, mass-media hype, and public relations may well have to be rewritten.

Why evolvability appears to have increased over evolutionary time is an important unresolved biological question. Unlike most candidate explanations, this paper proposes that increasing evolvability can result without any pressure to adapt. The insight is that if evolvability is heritable, then an unbiased drifting process across genotypes can still create a distribution of phenotypes biased towards evolvability, because evolvable organisms diffuse more quickly through the space of possible phenotypes. Furthermore, because phenotypic divergence often correlates with founding niches, niche founders may on average be more evolvable, which through population growth provides a genotypic bias towards evolvability. Interestingly, the combination of these two mechanisms can lead to increasing evolvability without any pressure to out-compete other organisms, as demonstrated through experiments with a series of simulated models. Thus rather than from pressure to adapt, evolvability may inevitably result from any drift through genotypic space combined with evolution’s passive tendency to accumulate niches.