Single-Molecule Studies on PolySUMO Proteins Reveal Their Mechanical Flexibility

Proteins with β-sandwich and β-grasp topologies are resistant to mechanical unfolding as shown by single-molecule force spectroscopy studies. Their high mechanical stability has generally been associated with the mechanical clamp geometry present at the termini. However, there is also evidence for the importance of interactions other than the mechanical clamp in providing mechanical stability, which needs to be tested thoroughly. Here, we report the mechanical unfolding properties of ubiquitin-like proteins (SUMO1 and SUMO2) and their comparison with those of ubiquitin. Although ubiquitin and SUMOs have similar size and structural topology, they differ in their sequences and structural contacts, making them ideal candidates to understand the variations in the mechanical stability of a given protein topology. We observe a two-state unfolding pathway for SUMO1 and SUMO2, similar to that of ubiquitin. Nevertheless, the unfolding forces of SUMO1 (∼130 pN) and SUMO2 (∼120 pN) are lower than that of ubiquitin (∼190 pN) at a pulling speed of 400 nm/s, indicating their lower mechanical stability. The mechanical stabilities of SUMO proteins and ubiquitin are well correlated with the number of interresidue contacts present in their structures. From pulling speed-dependent mechanical unfolding experiments and Monte Carlo simulations, we find that the unfolding potential widths of SUMO1 (∼0.51 nm) and SUMO2 (∼0.33 nm) are much larger than that of ubiquitin (∼0.19 nm), indicating that SUMO1 is six times and SUMO2 is three times mechanically more flexible than ubiquitin. These findings might also be important in understanding the functional differences between ubiquitin and SUMOs.     

Integration of New Genes into Cellular Networks,and Their Structural Maturation

It has been recently discovered that new genes can originate de novo from noncoding DNA, and several biological traits including expression or sequence composition form a continuum from noncoding sequences to conserved genes. In this article, using yeast genes I test whether the integration of new genes into cellular networks and their structural maturation shows such a continuum by analyzing their changes with gene age. I show that 1) The number of regulatory, protein–protein, and genetic interactions increases continuously with gene age, although with very different rates. New regulatory interactions emerge rapidly within a few million years, while the number of protein–protein and genetic interactions increases slowly, with a rate of 2–2.25 × 10⁻⁸/year and 4.8 × 10⁻⁸/year, respectively. 2) Gene essentiality evolves relatively quickly: the youngest essential genes appear in proto-genes ∼14 MY old. 3) In contrast to interactions, the secondary structure of proteins and their robustness to mutations indicate that new genes face a bottleneck in their evolution: proto-genes are characterized by high β-strand content, high aggregation propensity, and low robustness against mutations, while conserved genes are characterized by lower strand content and higher stability, most likely due to the higher probability of gene loss among young genes and accumulation of neutral mutations.     

Quantifying separation and similarity in a Saccharomyces cerevisiae metapopulation

Eukaryotic microbes are key ecosystem drivers; however, we have little theory and few data elucidating the processes influencing their observed population patterns. Here we provide an in-depth quantitative analysis of population separation and similarity in the yeast Saccharomyces cerevisiae with the aim of providing a more detailed account of the population processes occurring in microbes. Over 10 000 individual isolates were collected from native plants, vineyards and spontaneous ferments of fruit from six major regions spanning 1000 km across New Zealand. From these, hundreds of S. cerevisiae genotypes were obtained, and using a suite of analytical methods we provide comprehensive quantitative estimates for both population structure and rates of gene flow or migration. No genetic differentiation was detected within geographic regions, even between populations inhabiting native forests and vineyards. We do, however, reveal a picture of national population structure at scales above ∼100 km with distinctive populations in the more remote Nelson and Central Otago regions primarily contributing to this. In addition, differential degrees of connectivity between regional populations are observed and correlate with the movement of fruit by the New Zealand wine industry. This suggests some anthropogenic influence on these observed population patterns.     

Inferring the Genetic Basis of Sex Determination from the Genome of a Dioecious Nightshade

Dissecting the genetic mechanisms underlying dioecy (i.e., separate female and male individuals) is critical for understanding the evolution of this pervasive reproductive strategy. Nonetheless, the genetic basis of sex determination remains unclear in many cases, especially in systems where dioecy has arisen recently. Within the economically important plant genus Solanum (∼2,000 species), dioecy is thought to have evolved independently at least 4 times across roughly 20 species. Here, we generate the first genome sequence of a dioecious Solanum and use it to ascertain the genetic basis of sex determination in this species. We de novo assembled and annotated the genome of Solanum appendiculatum (assembly size: ∼750 Mb scaffold N50: 0.92 Mb; ∼35,000 genes), identified sex-specific sequences and their locations in the genome, and inferred that males in this species are the heterogametic sex. We also analyzed gene expression patterns in floral tissues of males and females, finding approximately 100 genes that are differentially expressed between the sexes. These analyses, together with observed patterns of gene-family evolution specific to S. appendiculatum, consistently implicate a suite of genes from the regulatory network controlling pectin degradation and modification in the expression of sex. Furthermore, the genome of a species with a relatively young sex-determination system provides the foundational resources for future studies on the independent evolution of dioecy in this clade.     

Association Mapping across Numerous Traits Reveals Patterns of Functional Variation in Maize

Phenotypic variation in natural populations results from a combination of genetic effects, environmental effects, and gene-by-environment interactions. Despite the vast amount of genomic data becoming available, many pressing questions remain about the nature of genetic mutations that underlie functional variation. We present the results of combining genome-wide association analysis of 41 different phenotypes in ∼5,000 inbred maize lines to analyze patterns of high-resolution genetic association among of 28.9 million single-nucleotide polymorphisms (SNPs) and ∼800,000 copy-number variants (CNVs). We show that genic and intergenic regions have opposite patterns of enrichment, minor allele frequencies, and effect sizes, implying tradeoffs among the probability that a given polymorphism will have an effect, the detectable size of that effect, and its frequency in the population. We also find that genes tagged by GWAS are enriched for regulatory functions and are ∼50% more likely to have a paralog than expected by chance, indicating that gene regulation and gene duplication are strong drivers of phenotypic variation. These results will likely apply to many other organisms, especially ones with large and complex genomes like maize.     

Prediction of Structure of Human WNT-CRD (FZD) Complex for Computational Drug Repurposing

The observed genetic alterations of various extracellular and intracellular WNT (Wingless, Int-1 proto-oncogene) signaling components can result in an increase or decrease in gene expression, and hence can be obstructed proficiently. These genetics target sites may include the prevention of WNT-FZD (Frizzled) binding, destruction of β-catenin and formation of Axin, APC and GSK-3β complex. Hence, the localized targeting of these interacting partners can help in devising novel inhibitors against WNT signaling. Our present study is an extension of our previous work, in which we proposed the co-regulated expression pattern of the WNT gene cluster (WNT-1, WNT-6, WNT-10A and WNT-10B) in human breast carcinoma. We present here the computationally modeled three dimensional structure of human WNT-1 in complex with the FZD-1 CRD (Cysteine Rich Domain) receptor. The dimeric cysteine-rich domain was found to fit into the evolutionarily conserved U-shaped groove of WNT protein. The two ends of the U- shaped cleft contain N-terminal and C-terminal hydrophobic residues, thus providing a strong hydrophobic moiety for the frizzled receptor and serving as the largest binding pocket for WNT-FZD interaction. Detailed structural analysis of this cleft revealed a maximum atomic distance of ∼28 Å at the surface, narrowing down to ∼17 Å and again increasing up to ∼27 Å at the bottom. Altogether, structural prediction analysis of WNT proteins was performed to reveal newer details about post-translational modification sites and to map the novel pharmacophore models for potent WNT inhibitors.     

Giga-Pixel Lensfree Holographic Microscopy and Tomography Using Color Image Sensors

We report Giga-pixel lensfree holographic microscopy and tomography using color sensor-arrays such as CMOS imagers that exhibit Bayer color filter patterns. Without physically removing these color filters coated on the sensor chip, we synthesize pixel super-resolved lensfree holograms, which are then reconstructed to achieve ∼350 nm lateral resolution, corresponding to a numerical aperture of ∼0.8, across a field-of-view of ∼20.5 mm². This constitutes a digital image with ∼0.7 Billion effective pixels in both amplitude and phase channels (i.e., ∼1.4 Giga-pixels total). Furthermore, by changing the illumination angle (e.g., ±50°) and scanning a partially-coherent light source across two orthogonal axes, super-resolved images of the same specimen from different viewing angles are created, which are then digitally combined to synthesize tomographic images of the object. Using this dual-axis lensfree tomographic imager running on a color sensor-chip, we achieve a 3D spatial resolution of ∼0.35 µm×0.35 µm×∼2 µm, in x, y and z, respectively, creating an effective voxel size of ∼0.03 µm³ across a sample volume of ∼5 mm³, which is equivalent to >150 Billion voxels. We demonstrate the proof-of-concept of this lensfree optical tomographic microscopy platform on a color CMOS image sensor by creating tomograms of micro-particles as well as a wild-type C. elegans nematode.     

Evolutionary Genomics of Structural Variation in Asian Rice (Oryza sativa) Domestication

Structural variants (SVs) are a largely unstudied feature of plant genome evolution, despite the fact that SVs contribute substantially to phenotypes. In this study, we discovered SVs across a population sample of 347 high-coverage, resequenced genomes of Asian rice (Oryza sativa) and its wild ancestor (O. rufipogon). In addition to this short-read data set, we also inferred SVs from whole-genome assemblies and long-read data. Comparisons among data sets revealed different features of genome variability. For example, genome alignment identified a large (∼4.3 Mb) inversion in indica rice varieties relative to japonica varieties, and long-read analyses suggest that ∼9% of genes from the outgroup (O. longistaminata) are hemizygous. We focused, however, on the resequencing sample to investigate the population genomics of SVs. Clustering analyses with SVs recapitulated the rice cultivar groups that were also inferred from SNPs. However, the site-frequency spectrum of each SV type—which included inversions, duplications, deletions, translocations, and mobile element insertions—was skewed toward lower frequency variants than synonymous SNPs, suggesting that SVs may be predominantly deleterious. Among transposable elements, SINE and mariner insertions were found at especially low frequency. We also used SVs to study domestication by contrasting between rice and O. rufipogon. Cultivated genomes contained ∼25% more derived SVs and mobile element insertions than O. rufipogon, indicating that SVs contribute to the cost of domestication in rice. Peaks of SV divergence were enriched for known domestication genes, but we also detected hundreds of genes gained and lost during domestication, some of which were enriched for traits of agronomic interest.     

Whole-Genome Resequencing of Worldwide Wild and Domestic Sheep Elucidates Genetic Diversity,Introgression, and Agronomically Important Loci

Domestic sheep and their wild relatives harbor substantial genetic variants that can form the backbone of molecular breeding, but their genome landscapes remain understudied. Here, we present a comprehensive genome resource for wild ovine species, landraces and improved breeds of domestic sheep, comprising high-coverage (∼16.10×) whole genomes of 810 samples from 7 wild species and 158 diverse domestic populations. We detected, in total, ∼121.2 million single nucleotide polymorphisms, ∼61 million of which are novel. Some display significant (P < 0.001) differences in frequency between wild and domestic species, or are private to continent-wide or individual sheep populations. Retained or introgressed wild gene variants in domestic populations have contributed to local adaptation, such as the variation in the HBB associated with plateau adaptation. We identified novel and previously reported targets of selection on morphological and agronomic traits such as stature, horn, tail configuration, and wool fineness. We explored the genetic basis of wool fineness and unveiled a novel mutation (chr25: T7,068,586C) in the 3′-UTR of IRF2BP2 as plausible causal variant for fleece fiber diameter. We reconstructed prehistorical migrations from the Near Eastern domestication center to South-and-Southeast Asia and found two main waves of migrations across the Eurasian Steppe and the Iranian Plateau in the Early and Late Bronze Ages. Our findings refine our understanding of genome variation as shaped by continental migrations, introgression, adaptation, and selection of sheep.     

Development of a Site-Directed Integration Plasmid for Heterologous Gene Expression in Mycoplasma gallisepticum

Deciphering the molecular basis of the interactions between the parasite Mycoplasma gallisepticum and its avian hosts suffers from the lack of genetic tools available for the pathogen. In the absence of well established methods for targeted disruption of relevant M. gallisepticum genes, we started to develop suicide vectors and equipped them with a short fragment of M. gallisepticum origin or replication (oriC _MG). We failed to create a disruption vector, although by adding a further short fragment of the M. gallisepticum tufB upstream region we created a “Trojan horse” plasmid. This is fully integrated into the genomic DNA of M. gallisepticum, always at the same site, oriC _MG, and is able to carry and express any gene of interest in the genetic background of M. gallisepticum. Successful expression of a heterologous gene was shown with the lacZ gene of E. coli. When used for gene complementation or expression of hybrid genes in M. gallisepticum, a site-specific combined integration/expression vector constitutes an improvement on randomly integrating transposons, which might have unexpected effects on the expression of chromosomal genes.     

A Common 16p11.2 Inversion Underlies the Joint Susceptibility to Asthma and Obesity

The prevalence of asthma and obesity is increasing worldwide, and obesity is a well-documented risk factor for asthma. The mechanisms underlying this association and parallel time trends remain largely unknown but genetic factors may be involved. Here, we report on a common ∼0.45 Mb genomic inversion at 16p11.2 that can be accurately genotyped via SNP array data. We show that the inversion allele protects against the joint occurrence of asthma and obesity in five large independent studies (combined sample size of 317 cases and 543 controls drawn from a total of 5,809 samples; combined OR = 0.48, p = 5.5 × 10⁻⁶). Allele frequencies show remarkable worldwide population stratification, ranging from 10% in East Africa to 49% in Northern Europe, consistent with discordant and extreme genetic drifts or adaptive selections after human migration out of Africa. Inversion alleles strongly correlate with expression levels of neighboring genes, especially TUFM (p = 3.0 × 10⁻⁴⁰) that encodes a mitochondrial protein regulator of energy balance and inhibitor of type 1 interferon, and other candidates for asthma (IL27) and obesity (APOB48R and SH2B1). Therefore, by affecting gene expression, the ∼0.45 Mb 16p11.2 inversion provides a genetic basis for the joint susceptibility to asthma and obesity, with a population attributable risk of 39.7%. Differential mitochondrial function and basal energy balance of inversion alleles might also underlie the potential selection signature that led to their uneven distribution in world populations.     

Geometrical Membrane Curvature as an Allosteric Regulator of Membrane Protein Structure and Function

Transmembrane proteins are embedded in cellular membranes of varied lipid composition and geometrical curvature. Here, we studied for the first time the allosteric effect of geometrical membrane curvature on transmembrane protein structure and function. We used single-channel optical analysis of the prototypic transmembrane β-barrel α-hemolysin (α-HL) reconstituted on immobilized single small unilamellar liposomes of different diameter and therefore curvature. Our data demonstrate that physiologically abundant geometrical membrane curvatures can enforce a dramatic allosteric regulation (1000-fold inhibition) of α-HL permeability. High membrane curvatures (1/diameter ∼1/40 nm⁻¹) compressed the effective pore diameter of α-HL from 14.2 ± 0.8 Å to 11.4 ± 0.6 Å. This reduction in effective pore area (∼40%) when combined with the area compressibility of α-HL revealed an effective membrane tension of ∼50 mN/m and a curvature-imposed protein deformation energy of ∼7 k_BT. Such substantial energies have been shown to conformationally activate, or unfold, β-barrel and α-helical transmembrane proteins, suggesting that membrane curvature could likely regulate allosterically the structure and function of transmembrane proteins in general.     

Low levels of realized seed and pollen gene flow and strong spatial genetic structure in a small,isolated and fragmented population of the tropical tree Copaifera langsdorffii Desf

Over the past century, the Brazilian Atlantic forest has been reduced to small, isolated fragments of forest. Reproductive isolation theories predict a loss of genetic diversity and increases in inbreeding and spatial genetic structure (SGS) in such populations. We analysed eight microsatellite loci to investigate the pollen and seed dispersal patterns, genetic diversity, inbreeding and SGS of the tropical tree Copaifera langsdorffii in a small (4.8 ha), isolated population. All 112 adult trees and 128 seedlings found in the stand were sampled, mapped and genotyped. Seedlings had significantly lower levels of genetic diversity (A=16.5±0.45, mean±95% s.e.; H_e=0.838±0.006) than did adult trees (A=23.2±0.81; H_e=0.893±0.030). Parentage analysis did not indicate any seed immigration (m_seeds=0) and the pollen immigration rate was very low (m_pollen=0.047). The average distance of realized pollen dispersal within the stand was 94 m, with 81% of the pollen travelling <150 m. A significant negative correlation was found between the frequency and distance of pollen dispersal (r=−0.79, P<0.01), indicating that short-distance pollinations were more frequent. A significant SGS for both adults (∼50 m) and seedlings (∼20 m) was also found, indicating that most of the seeds were dispersed over short distances. The results suggested that the spatial isolation of populations by habitat fragmentation can restrict seed and pollen gene flow, increase SGS and affect the genetic diversity of future generations.     

Investigation of PAS and CNBH domain interactions in hERG channels and effects of long-QT syndrome-causing mutations with surface plasmon resonance

Human ether-á-go-go-related gene (hERG) channels are key regulators of cardiac repolarization, neuronal excitability, and tumorigenesis. hERG channels contain N-terminal Per-Arnt-Sim (PAS) and C-terminal cyclic nucleotide-binding homology (CNBH) domains with many long-QT syndrome (LQTS)-causing mutations located at the interface between these domains. Despite the importance of PAS/CNBH domain interactions, little is known about their affinity. Here, we used the surface plasmon resonance (SPR) technique to investigate interactions between isolated PAS and CNBH domains and the effects of LQTS-causing mutations R20G, N33T, and E58D, located at the PAS/CNBH domain interface, on these interactions. We determined that the affinity of the PAS/CNBH domain interactions was ∼1.4 μM. R20G and E58D mutations had little effect on the domain interaction affinity, while N33T abolished the domain interactions. Interestingly, mutations in the intrinsic ligand, a conserved stretch of amino acids occupying the beta-roll cavity in the CNBH domain, had little effect on the affinity of PAS/CNBH domain interactions. Additionally, we determined that the isolated PAS domains formed oligomers with an interaction affinity of ∼1.6 μM. Coexpression of the isolated PAS domains with the full-length hERG channels or addition of the purified PAS protein inhibited hERG currents. These PAS/PAS interactions can have important implications for hERG function in normal and pathological conditions associated with increased surface density of channels or interaction with other PAS-domain-containing proteins. Taken together, our study provides the first account of the binding affinities for wild-type and mutant hERG PAS and CNBH domains and highlights the potential functional significance of PAS/PAS domain interactions.     

Rates and Genomic Consequences of Spontaneous Mutational Events in Drosophila melanogaster

Because spontaneous mutation is the source of all genetic diversity, measuring mutation rates can reveal how natural selection drives patterns of variation within and between species. We sequenced eight genomes produced by a mutation-accumulation experiment in Drosophila melanogaster. Our analysis reveals that point mutation and small indel rates vary significantly between the two different genetic backgrounds examined. We also find evidence that ∼2% of mutational events affect multiple closely spaced nucleotides. Unlike previous similar experiments, we were able to estimate genome-wide rates of large deletions and tandem duplications. These results suggest that, at least in inbred lines like those examined here, mutational pressures may result in net growth rather than contraction of the Drosophila genome. By comparing our mutation rate estimates to polymorphism data, we are able to estimate the fraction of new mutations that are eliminated by purifying selection. These results suggest that ∼99% of duplications and deletions are deleterious—making them 10 times more likely to be removed by selection than nonsynonymous mutations. Our results illuminate not only the rates of new small- and large-scale mutations, but also the selective forces that they encounter once they arise.     

Detection of Salmonella spp. Using a Generic and Differential FRET-PCR

To facilitate the detection of Salmonella and to be able to rapidly and conveniently determine the species/subspecies present, we developed and tested a generic and differential FRET-PCR targeting their tetrathionate reductase response regulator gene. The differential pan-Salmonella FRET-PCR we developed successfully detected seven plasmids that contained partial sequences of S. bongori and the six S. enterica subspecies. The detection limit varied from ∼5 copies of target gene/per PCR reaction for S. enterica enterica to ∼200 for S. bongori. Melting curve analysis demonstrated a T _m of ∼68°C for S. enterica enterica, ∼62.5°C for S. enterica houtenae and S. enterica diarizonae, ∼57°C for S. enterica indica, and ∼54°C for S. bongori, S. enterica salamae and S. enterica arizonae. The differential pan-Salmonella FRET-PCR also detected and determined the subspecies of 4 reference strains and 47 Salmonella isolated from clinically ill birds or pigs. Finally, we found it could directly detect and differentiate Salmonella in feline (5/50 positive; 10%; one S. enterica salamae and 4 S. enterica enterica) and canine feces (15/114 positive; 13.2%; all S. enterica enterica). The differential pan-Salmonella FRET-PCR failed to react with 96 non-Salmonella bacterial strains. Our experiments show the differential pan-Salmonella FRET-PCR we developed is a rapid, sensitive and specific method to detect and differentiate Salmonella.     

The Evolutionary Potential of Phenotypic Mutations

Errors in protein synthesis, so-called phenotypic mutations, are orders-of-magnitude more frequent than genetic mutations. Here, we provide direct evidence that alternative protein forms and phenotypic variability derived from translational errors paved the path to genetic, evolutionary adaptations via gene duplication. We explored the evolutionary origins of Saccharomyces cerevisiae IDP3 - an NADP-dependent isocitrate dehydrogenase mediating fatty acids ß-oxidation in the peroxisome. Following the yeast whole genome duplication, IDP3 diverged from a cytosolic ancestral gene by acquisition of a C-terminal peroxisomal targeting signal. We discovered that the pre-duplicated cytosolic IDPs are partially localized to the peroxisome owing to +1 translational frameshifts that bypass the stop codon and unveil cryptic peroxisomal targeting signals within the 3’-UTR. Exploring putative cryptic signals in all 3’-UTRs of yeast genomes, we found that other enzymes related to NADPH production such as pyruvate carboxylase 1 (PYC1) might be prone to peroxisomal localization via cryptic signals. Using laboratory evolution we found that these translational frameshifts are rapidly imprinted via genetic single base deletions occurring within the very same gene location. Further, as exemplified here, the sequences that promote translational frameshifts are also more prone to genetic deletions. Thus, genotypes conferring higher phenotypic variability not only meet immediate challenges by unveiling cryptic 3’-UTR sequences, but also boost the potential for future genetic adaptations.