Search

The analysis of introgression of genomic regions between divergent populations provides an excellent opportunity to determine the genetic basis of reproductive isolation during the early stages of speciation. However, hybridization and subsequent gene flow must be relatively common in order to localize individual loci that resist introgression. In this study, we used next-generation sequencing to study genome-wide patterns of genetic differentiation between two hybridizing subspecies of rabbits (Oryctolagus cuniculus algirus and O. c. cuniculus) that are known to undergo high rates of gene exchange. Our primary objective was to identify specific genes or genomic regions that have resisted introgression and are likely to confer reproductive barriers in natural conditions. On the basis of 326,000 polymorphisms, we found low to moderate overall levels of differentiation between subspecies, and fewer than 200 genomic regions dispersed throughout the genome showing high differentiation consistent with a signature of reduced gene flow. Most differentiated regions were smaller than 200 Kb and contained very few genes. Remarkably, 30 regions were each found to contain a single gene, facilitating the identification of candidate genes underlying reproductive isolation. This gene-level resolution yielded several insights into the genetic basis and architecture of reproductive isolation in rabbits. Regions of high differentiation were enriched on the X-chromosome and near centromeres. Genes lying within differentiated regions were often associated with transcription and epigenetic activities, including chromatin organization, regulation of transcription, and DNA binding. Overall, our results from a naturally hybridizing system share important commonalities with hybrid incompatibility genes identified using laboratory crosses in mice and flies, highlighting general mechanisms underlying the maintenance of reproductive barriers. [ABSTRACT FROM AUTHOR]

Bipolar disorder is a common, heritable mental illness characterized by recurrent episodes of mania and depression. Despite considerable effort to elucidate the genetic underpinnings of bipolar disorder, causative genetic risk factors remain elusive. We conducted a comprehensive genomic analysis of bipolar disorder in a large Old Order Amish pedigree. Microsatellite genotypes and high-density SNP-array genotypes of 388 family members were combined with whole genome sequence data for 50 of these subjects, comprising 18 parent-child trios. This study design permitted evaluation of candidate variants within the context of haplotype structure by resolving the phase in sequenced parent-child trios and by imputation of variants into multiple unsequenced siblings. Non-parametric and parametric linkage analysis of the entire pedigree as well as on smaller clusters of families identified several nominally significant linkage peaks, each of which included dozens of predicted deleterious variants. Close inspection of exonic and regulatory variants in genes under the linkage peaks using family-based association tests revealed additional credible candidate genes for functional studies and further replication in population-based cohorts. However, despite the in-depth genomic characterization of this unique, large and multigenerational pedigree from a genetic isolate, there was no convergence of evidence implicating a particular set of risk loci or common pathways. The striking haplotype and locus heterogeneity we observed has profound implications for the design of studies of bipolar and other related disorders. [ABSTRACT FROM AUTHOR]

Metabolic traits are molecular phenotypes that can drive clinical phenotypes and may predict disease progression. Here, we report results from a metabolome- and genome-wide association study on 1H-NMR urine metabolic profiles. The study was conducted within an untargeted approach, employing a novel method for compound identification. From our discovery cohort of 835 Caucasian individuals who participated in the CoLaus study, we identified 139 suggestively significant (P<5×10−8) and independent associations between single nucleotide polymorphisms (SNP) and metabolome features. Fifty-six of these associations replicated in the TasteSensomics cohort, comprising 601 individuals from São Paulo of vastly diverse ethnic background. They correspond to eleven gene-metabolite associations, six of which had been previously identified in the urine metabolome and three in the serum metabolome. Our key novel findings are the associations of two SNPs with NMR spectral signatures pointing to fucose (rs492602, P = 6.9×10−44) and lysine (rs8101881, P = 1.2×10−33), respectively. Fine-mapping of the first locus pinpointed the FUT2 gene, which encodes a fucosyltransferase enzyme and has previously been associated with Crohn's disease. This implicates fucose as a potential prognostic disease marker, for which there is already published evidence from a mouse model. The second SNP lies within the SLC7A9 gene, rare mutations of which have been linked to severe kidney damage. The replication of previous associations and our new discoveries demonstrate the potential of untargeted metabolomics GWAS to robustly identify molecular disease markers. [ABSTRACT FROM AUTHOR]

Bacteria are believed to play an important role in the fitness and biochemistry of sponges (Porifera). Pseudomonas species (Gammaproteobacteria, Pseudomonadales) are capable of colonizing a broad range of eukaryotic hosts, but knowledge of their diversity and function in freshwater invertebrates is rudimentary. We assessed the diversity, structure and antimicrobial activities of Pseudomonas spp. in the freshwater sponge Ephydatia fluviatilis. Polymerase Chain Reaction – Denaturing Gradient Gel Electrophoresis (PCR-DGGE) fingerprints of the global regulator gene gacA revealed distinct structures between sponge-associated and free-living Pseudomonas communities, unveiling previously unsuspected diversity of these assemblages in freshwater. Community structures varied across E. fluviatilis specimens, yet specific gacA phylotypes could be detected by PCR-DGGE in almost all sponge individuals sampled over two consecutive years. By means of whole-genome fingerprinting, 39 distinct genotypes were found within 90 fluorescent Pseudomonas isolates retrieved from E. fluviatilis. High frequency of in vitro antibacterial (49%), antiprotozoan (35%) and anti-oomycetal (32%) activities was found among these isolates, contrasting less-pronounced basidiomycetal (17%) and ascomycetal (8%) antagonism. Culture extracts of highly predation-resistant isolates rapidly caused complete immobility or lysis of cells of the protozoan Colpoda steinii. Isolates tentatively identified as P. jessenii, P. protegens and P. oryzihabitans showed conspicuous inhibitory traits and correspondence with dominant sponge-associated phylotypes registered by cultivation-independent analysis. Our findings suggest that E. fluviatilis hosts both transient and persistent Pseudomonas symbionts displaying antimicrobial activities of potential ecological and biotechnological value. [ABSTRACT FROM AUTHOR]

AIRmax (maximal inflammation) and AIRmin (minimal inflammation) mice show distinct susceptibilities to pristane-induced arthritis (PIA). The Slc11a1 gene, which regulates macrophage and neutrophil activity, is involved in this infirmity. AIRmaxSS mice homozygous for the non-functional Slc11a1 S (gly169asp) allele obtained by genotype-assisted crosses from AIRmax and AIRmin mice are more susceptible than mice homozygous for the Slc11a1 resistant (R) allele. The present work sought to identify the quantitative trait loci (QTL) regulating PIA and to examine the interactions of these QTL with Slc11a1 alleles in modulating PIA. Mice were given two ip injections of 0.5 mL pristane at 60 day intervals, and the incidence and severity of PIA was scored up to 160 days. Genome-wide linkage studies were performed to search for arthritis QTL in an F2 (AIRmax × AIRmin, n = 290) population. Significant arthritis QTL (LODscore>4) were detected on chromosomes 5 and 8, and suggestive QTL on chromosomes 7, 17 and 19. Global gene expression analyses performed on Affymetrix mouse 1.0 ST bioarrays (27k genes) using RNA from arthritic or control mice paws showed 419 differentially expressed genes between AIRmax and AIRmin mice and demonstrated significantly (P<0.001) over-represented genes related to inflammatory responses and chemotaxis. Up-regulation of the chemokine genes Cxcl1, Cxcl9, Cxcl5, Cxcl13 on chromosome 5 was higher in AIRmaxSS than in the other lines. Macrophage scavenger receptor 1 and hemeoxigenase (decycling) 1 genes on chromosome 8 were also expressed at higher levels in AIRmaxSS mice. Our results show that the gene expression profiles of the two arthritis QTL (on chromosomes 5 and 8) correlate with Slc11a1 alleles, resulting in enhanced AIRmaxSS mice susceptibility to PIA. [ABSTRACT FROM AUTHOR]

Integrating genetic data from families with highly penetrant forms of disease together with genetic data from outbred populations represents a promising strategy to uncover the complete frequency spectrum of risk alleles for complex traits such as rheumatoid arthritis (RA). Here, we demonstrate that rare, low-frequency and common alleles at one gene locus, phospholipase B1 (PLB1), might contribute to risk of RA in a 4-generation consanguineous pedigree (Middle Eastern ancestry) and also in unrelated individuals from the general population (European ancestry). Through identity-by-descent (IBD) mapping and whole-exome sequencing, we identified a non-synonymous c.2263G>C (p.G755R) mutation at the PLB1 gene on 2q23, which significantly co-segregated with RA in family members with a dominant mode of inheritance (P = 0.009). We further evaluated PLB1 variants and risk of RA using a GWAS meta-analysis of 8,875 RA cases and 29,367 controls of European ancestry. We identified significant contributions of two independent non-coding variants near PLB1 with risk of RA (rs116018341 [MAF = 0.042] and rs116541814 [MAF = 0.021], combined P = 3.2×10−6). Finally, we performed deep exon sequencing of PLB1 in 1,088 RA cases and 1,088 controls (European ancestry), and identified suggestive dispersion of rare protein-coding variant frequencies between cases and controls (P = 0.049 for C-alpha test and P = 0.055 for SKAT). Together, these data suggest that PLB1 is a candidate risk gene for RA. Future studies to characterize the full spectrum of genetic risk in the PLB1 genetic locus are warranted. [ABSTRACT FROM AUTHOR]

Genome-wide mapping of transcription factor binding and histone modification reveals complex patterns of interactions. Identifying overlaps in binding patterns by different factors is a major objective of genomic studies, but existing methods to archive large numbers of datasets in a personalised database lack sophistication and utility. Therefore we have developed transcription factor DNA binding site analyser software (BiSA), for archiving of binding regions and easy identification of overlap with or proximity to other regions of interest. Analysis results can be restricted by chromosome or base pair overlap between regions or maximum distance between binding peaks. BiSA is capable of reporting overlapping regions that share common base pairs; regions that are nearby; regions that are not overlapping; and average region sizes. BiSA can identify genes located near binding regions of interest, genomic features near a gene or locus of interest and statistical significance of overlapping regions can also be reported. Overlapping results can be visualized as Venn diagrams. A major strength of BiSA is that it is supported by a comprehensive database of publicly available transcription factor binding sites and histone modifications, which can be directly compared to user data. The documentation and source code are available on http://bisa.sourceforge.net [ABSTRACT FROM AUTHOR]

It is increasingly realised that the molecular clock does not tick at a constant rate. Rather, mitochondrial mutation rates are influenced by factors such as generation length and body mass. This has implications for the use of genetic data in species delimitation. It could be that speciation, as recognised by avian taxonomists, is associated with a certain minimum genetic distance between sister taxa, in which case we would predict no difference in the cytochrome b divergence of sister taxa according to the species' body size or generation time. Alternatively, if what taxonomists recognise as speciation has tended to be associated with the passage of a minimum amount of time since divergence, then there might be less genetic divergence between sister taxa with slower mutation rates, namely those that are heavier and/or with longer generation times. After excluding non-flying species, we analysed a database of over 600 avian sister species pairs, and found that species pairs with longer generation lengths (which tend to be the larger species) showed less cytochrome b divergence. This finding cautions against using any simple unitary criterion of genetic divergence to delimit species. [ABSTRACT FROM AUTHOR]

Campylobacter jejuni is a major cause of bacterial gastroenteritis worldwide, and the capsular polysaccharide (CPS) of this organism is required for persistence and disease. C. jejuni produces over 47 different capsular structures, including a unique O-methyl phosphoramidate (MeOPN) modification present on most C. jejuni isolates. Although the MeOPN structure is rare in nature it has structural similarity to some synthetic pesticides. In this study, we have demonstrated, by whole genome comparisons and high resolution magic angle spinning NMR, that MeOPN modifications are common to several Campylobacter species. Using MeOPN biosynthesis and transferase mutants generated in C. jejuni strain 81–176, we observed that loss of MeOPN from the cell surface correlated with increased invasion of Caco-2 epithelial cells and reduced resistance to killing by human serum. In C. jejuni, the observed serum mediated killing was determined to result primarily from activation of the classical complement pathway. The C. jejuni MeOPN transferase mutant showed similar levels of colonization relative to the wild-type in chickens, but showed a five-fold drop in colonization when co-infected with the wild-type in piglets. In Galleria mellonella waxmoth larvae, the MeOPN transferase mutant was able to kill the insects at wild-type levels. Furthermore, injection of the larvae with MeOPN-linked monosaccharides or CPS purified from the wild-type strain did not result in larval killing, indicating that MeOPN does not have inherent insecticidal activity. [ABSTRACT FROM AUTHOR]

Background: High frequency of physical aggression is the central feature of severe conduct disorder and is associated with a wide range of social, mental and physical health problems. We have previously tested the hypothesis that differential DNA methylation signatures in peripheral T cells are associated with a chronic aggression trajectory in males. Despite the fact that sex differences appear to play a pivotal role in determining the development, magnitude and frequency of aggression, most of previous studies focused on males, so little is known about female chronic physical aggression. We therefore tested here whether or not there is a signature of physical aggression in female DNA methylation and, if there is, how it relates to the signature observed in males. Methodology/Principal Findings: Methylation profiles were created using the method of methylated DNA immunoprecipitation (MeDIP) followed by microarray hybridization and statistical and bioinformatic analyses on T cell DNA obtained from adult women who were found to be on a chronic physical aggression trajectory (CPA) between 6 and 12 years of age compared to women who followed a normal physical aggression trajectory. We confirmed the existence of a well-defined, genome-wide signature of DNA methylation associated with chronic physical aggression in the peripheral T cells of adult females that includes many of the genes similarly associated with physical aggression in the same cell types of adult males. Conclusions: This study in a small number of women presents preliminary evidence for a genome-wide variation in promoter DNA methylation that associates with CPA in women that warrant larger studies for further verification. A significant proportion of these associations were previously observed in men with CPA supporting the hypothesis that the epigenetic signature of early life aggression in females is composed of a component specific to females and another common to both males and females. [ABSTRACT FROM AUTHOR]

Background: Disseminated cancer cells (DCCs) and circulating tumor cells (CTCs) are extremely rare, but comprise the precursors cells of distant metastases or therapy resistant cells. The detailed molecular analysis of these cells may help to identify key events of cancer cell dissemination, metastatic colony formation and systemic therapy escape. Methodology/Principal Findings: Using the Ampli1™ whole genome amplification (WGA) technology and high-resolution oligonucleotide aCGH microarrays we optimized conditions for the analysis of structural copy number changes. The protocol presented here enables reliable detection of numerical genomic alterations as small as 0.1 Mb in a single cell. Analysis of single cells from well-characterized cell lines and single normal cells confirmed the stringent quantitative nature of the amplification and hybridization protocol. Importantly, fixation and staining procedures used to detect DCCs showed no significant impact on the outcome of the analysis, proving the clinical usability of our method. In a proof-of-principle study we tracked the chromosomal changes of single DCCs over a full course of high-dose chemotherapy treatment by isolating and analyzing DCCs of an individual breast cancer patient at four different time points. Conclusions/Significance: The protocol enables detailed genome analysis of DCCs and thereby assessment of the clonal evolution during the natural course of the disease and under selection pressures. The results from an exemplary patient provide evidence that DCCs surviving selective therapeutic conditions may be recruited from a pool of genomically less advanced cells, which display a stable subset of specific genomic alterations. [ABSTRACT FROM AUTHOR]

Transcriptome sequencing using Illumina RNA-seq was performed on populations of black tiger shrimp from India. Samples were collected from (i) four landing centres around the east coastline (EC) of India, (ii) survivors of a severe WSSV infection during pond culture (SUR) and (iii) the Andaman Islands (AI) in the Bay of Bengal. Equal quantities of purified total RNA from homogenates of hepatopancreas, muscle, nervous tissue, intestinal tract, heart, gonad, gills, pleopod and lymphoid organs were combined to create AI, EC and SUR pools for RNA sequencing. De novo transcriptome assembly resulted in 136,223 contigs (minimum size 100 base pairs, bp) with a total length 61 Mb, an average length of 446 bp and an average coverage of 163× across all pools. Approximately 16% of contigs were annotated with BLAST hit information and gene ontology annotations. A total of 473,620 putative SNPs/indels were identified. An Illumina iSelect genotyping array containing 6,000 SNPs was developed and used to genotype 1024 offspring belonging to seven full-sibling families. A total of 3959 SNPs were mapped to 44 linkage groups. The linkage groups consisted of between 16–129 and 13–130 markers, of length between 139–10.8 and 109.1–10.5 cM and with intervals averaging between 1.2 and 0.9 cM for the female and male maps respectively. The female map was 28% longer than the male map (4060 and 2917 cM respectively) with a 1.6 higher recombination rate observed for female compared to male meioses. This approach has substantially increased expressed sequence and DNA marker resources for tiger shrimp and is a useful resource for QTL mapping and association studies for evolutionarily and commercially important traits. [ABSTRACT FROM AUTHOR]

Background: Prognostic biomarkers are needed for superficial gastroesophageal adenocarcinoma (EAC) to predict clinical outcomes and select therapy. Although recurrent mutations have been characterized in EAC, little is known about their clinical and prognostic significance. Aneuploidy is predictive of clinical outcome in many malignancies but has not been evaluated in superficial EAC. Methods: We quantified copy number changes in 41 superficial EAC using Affymetrix SNP 6.0 arrays. We identified recurrent chromosomal gains and losses and calculated the total copy number abnormality (CNA) count for each tumor as a measure of aneuploidy. We correlated CNA count with overall survival and time to first recurrence in univariate and multivariate analyses. Results: Recurrent segmental gains and losses involved multiple genes, including: HER2, EGFR, MET, CDK6, KRAS (recurrent gains); and FHIT, WWOX, CDKN2A/B, SMAD4, RUNX1 (recurrent losses). There was a 40-fold variation in CNA count across all cases. Tumors with the lowest and highest quartile CNA count had significantly better overall survival (p = 0.032) and time to first recurrence (p = 0.010) compared to those with intermediate CNA counts. These associations persisted when controlling for other prognostic variables. Significance: SNP arrays facilitate the assessment of recurrent chromosomal gain and loss and allow high resolution, quantitative assessment of segmental aneuploidy (total CNA count). The non-monotonic association of segmental aneuploidy with survival has been described in other tumors. The degree of aneuploidy is a promising prognostic biomarker in a potentially curable form of EAC. [ABSTRACT FROM AUTHOR]

Campylobacter jejuni is a major cause of bacterial gastroenteritis worldwide, and the capsular polysaccharide (CPS) of this organism is required for persistence and disease. C. jejuni produces over 47 different capsular structures, including a unique O-methyl phosphoramidate (MeOPN) modification present on most C. jejuni isolates. Although the MeOPN structure is rare in nature it has structural similarity to some synthetic pesticides. In this study, we have demonstrated, by whole genome comparisons and high resolution magic angle spinning NMR, that MeOPN modifications are common to several Campylobacter species. Using MeOPN biosynthesis and transferase mutants generated in C. jejuni strain 81–176, we observed that loss of MeOPN from the cell surface correlated with increased invasion of Caco-2 epithelial cells and reduced resistance to killing by human serum. In C. jejuni, the observed serum mediated killing was determined to result primarily from activation of the classical complement pathway. The C. jejuni MeOPN transferase mutant showed similar levels of colonization relative to the wild-type in chickens, but showed a five-fold drop in colonization when co-infected with the wild-type in piglets. In Galleria mellonella waxmoth larvae, the MeOPN transferase mutant was able to kill the insects at wild-type levels. Furthermore, injection of the larvae with MeOPN-linked monosaccharides or CPS purified from the wild-type strain did not result in larval killing, indicating that MeOPN does not have inherent insecticidal activity. [ABSTRACT FROM AUTHOR]

Systemic lupus erythematosus (SLE) is an autoimmune disease with known genetic, epigenetic, and environmental risk factors. To assess the role of DNA methylation in SLE, we collected CD4+ T-cells, CD19+ B-cells, and CD14+ monocytes from 49 SLE patients and 58 controls, and performed genome-wide DNA methylation analysis with Illumina Methylation450 microarrays. We identified 166 CpGs in B-cells, 97 CpGs in monocytes, and 1,033 CpGs in T-cells with highly significant changes in DNA methylation levels (p<1×10−8) among SLE patients. Common to all three cell-types were widespread and severe hypomethylation events near genes involved in interferon signaling (type I). These interferon-related changes were apparent in patients collected during active and quiescent stages of the disease, suggesting that epigenetically-mediated hypersensitivity to interferon persists beyond acute stages of the disease and is independent of circulating interferon levels. This interferon hypersensitivity was apparent in memory, naïve and regulatory T-cells, suggesting that this epigenetic state in lupus patients is established in progenitor cell populations. We also identified a widespread, but lower amplitude shift in methylation in CD4+ T-cells (>16,000 CpGs at FDR<1%) near genes involved in cell division and MAPK signaling. These cell type-specific effects are consistent with disease-specific changes in the composition of the CD4+ population and suggest that shifts in the proportion of CD4+ subtypes can be monitored at CpGs with subtype-specific DNA methylation patterns. [ABSTRACT FROM AUTHOR]

Next-generation sequencing projects continue to drive a vast accumulation of metagenomic sequence data. Given the growth rate of this data, automated approaches to functional annotation are indispensable and a cornerstone heuristic of many computational protocols is the concept of guilt by association. The guilt by association paradigm has been heavily exploited by genomic context methods that offer functional predictions that are complementary to homology-based annotations, thereby offering a means to extend functional annotation. In particular, operon methods that exploit co-directional intergenic distances can provide homology-free functional annotation through the transfer of functions among co-operonic genes, under the assumption that guilt by association is indeed applicable. Although guilt by association is a well-accepted annotative device, its applicability to metagenomic functional annotation has not been definitively demonstrated. Here a large-scale assessment of metagenomic guilt by association is undertaken where functional associations are predicted on the basis of co-directional intergenic distances. Specifically, functional annotations are compared within pairs of adjacent co-directional genes, as well as operons of various lengths (i.e. number of member genes), in order to reveal new information about annotative cohesion versus operon length. The results suggests that co-directional gene pairs offer reduced confidence for metagenomic guilt by association due to difficulty in resolving the existence of functional associations when intergenic distance is the sole predictor of pairwise gene interactions. However, metagenomic operons, particularly those with substantial lengths, appear to be capable of providing a superior basis for metagenomic guilt by association due to increased annotative stability. The need for improved recognition of metagenomic operons is discussed, as well as the limitations of the present work. [ABSTRACT FROM AUTHOR]

Campylobacter species.are phenotypically diverse in many aspects including host habitats and pathogenicities, which demands comprehensive characterization of the entire Campylobacter genus to study their underlying genetic diversification. Up to now, 34 Campylobacter strains have been sequenced and published in public databases, providing good opportunity to systemically analyze their genomic diversities. In this study, we first conducted genomic characterization, which includes genome-wide alignments, pan-genome analysis, and phylogenetic identification, to depict the genetic diversity of Campylobacter genus. Afterward, we improved the tetranucleotide usage pattern-based naïve Bayesian classifier to identify the abnormal composition fragments (ACFs, fragments with significantly different tetranucleotide frequency profiles from its genomic tetranucleotide frequency profiles) including horizontal gene transfers (HGTs) to explore the mechanisms for the genetic diversity of this organism. Finally, we analyzed the HGTs transferred via bacteriophage transductions. To our knowledge, this study is the first to use single nucleotide polymorphism information to construct liable microevolution phylogeny of 21 Campylobacter jejuni strains. Combined with the phylogeny of all the collected Campylobacter species based on genome-wide core gene information, comprehensive phylogenetic inference of all 34 Campylobacter organisms was determined. It was found that C. jejuni harbors a high fraction of ACFs possibly through intraspecies recombination, whereas other Campylobacter members possess numerous ACFs possibly via intragenus recombination. Furthermore, some Campylobacter strains have undergone significant ancient viral integration during their evolution process. The improved method is a powerful tool for bacterial genomic analysis. Moreover, the findings would provide useful information for future research on Campylobacter genus. [ABSTRACT FROM AUTHOR]

A workshop recently held at the École Polytechnique Fédérale de Lausanne ( EPFL, Switzerland) was dedicated to understanding the genetic basis of adaptive change, taking stock of the different approaches developed in theoretical population genetics and landscape genomics and bringing together knowledge accumulated in both research fields. Indeed, an important challenge in theoretical population genetics is to incorporate effects of demographic history and population structure. But important design problems (e.g. focus on populations as units, focus on hard selective sweeps, no hypothesis-based framework in the design of the statistical tests) reduce their capability of detecting adaptive genetic variation. In parallel, landscape genomics offers a solution to several of these problems and provides a number of advantages (e.g. fast computation, landscape heterogeneity integration). But the approach makes several implicit assumptions that should be carefully considered (e.g. selection has had enough time to create a functional relationship between the allele distribution and the environmental variable, or this functional relationship is assumed to be constant). To address the respective strengths and weaknesses mentioned above, the workshop brought together a panel of experts from both disciplines to present their work and discuss the relevance of combining these approaches, possibly resulting in a joint software solution in the future. [ABSTRACT FROM AUTHOR]

The early elongation checkpoint regulated by Positive Transcription Elongation Factor b (P-TEFb) is a critical control point for the expression of many genes. Spt5 interacts directly with RNA polymerase II and has an essential role in establishing this checkpoint, and also for further transcript elongation. Here we demonstrate that Drosophila Spt5 interacts both physically and genetically with the Polycomb Group (PcG) protein Pleiohomeotic (Pho), and the majority of Pho binding sites overlap with Spt5 binding sites across the genome in S2 cells. Our results indicate that Pho can interact with Spt5 to regulate transcription elongation in a gene specific manner. [ABSTRACT FROM AUTHOR]

Our bloodstream is considered to be an environment well separated from the outside world and the digestive tract. According to the standard paradigm large macromolecules consumed with food cannot pass directly to the circulatory system. During digestion proteins and DNA are thought to be degraded into small constituents, amino acids and nucleic acids, respectively, and then absorbed by a complex active process and distributed to various parts of the body through the circulation system. Here, based on the analysis of over 1000 human samples from four independent studies, we report evidence that meal-derived DNA fragments which are large enough to carry complete genes can avoid degradation and through an unknown mechanism enter the human circulation system. In one of the blood samples the relative concentration of plant DNA is higher than the human DNA. The plant DNA concentration shows a surprisingly precise log-normal distribution in the plasma samples while non-plasma (cord blood) control sample was found to be free of plant DNA. [ABSTRACT FROM AUTHOR]

Dermatofibrosarcoma protuberans (DFSP) is a very rare soft tissue sarcoma. DFSP often reveals a specific chromosome translocation, t(17;22)(q22;q13), which results in the fusion of collagen 1 alpha 1 (COL1A1) gene and platelet-derived growth factor-B (PDGFB) gene. The COL1A1-PDGFB fusion protein activates the PDGFB receptor and resultant constitutive activation of PDGFR receptor is essential in the pathogenesis of DFSP. Thus, blocking PDGFR receptor activation with imatinib has shown promising activity in the treatment of advanced and metastatic DFSP. Despite the success with targeted agents in cancers, acquired drug resistance eventually occurs. Here, we tried to identify potential drug resistance mechanisms against imatinib in a 46-year old female with DFSP who initially responded well to imatinib but suffered rapid disease progression. We performed whole-genome sequencing of both pre-treatment and post-treatment tumor tissue to identify the mutational events associated with imatinib resistance. No significant copy number alterations, insertion, and deletions were identified during imatinib treatment. Of note, we identified newly emerged 8 non-synonymous somatic mutations of the genes (ACAP2, CARD10, KIAA0556, PAAQR7, PPP1R39, SAFB2, STARD9, and ZFYVE9) in the imatinib-resistant tumor tissue. This study revealed diverse possible candidate mechanisms by which imatinib resistance to PDGFRB inhibition may arise in DFSP, and highlights the usefulness of whole-genome sequencing in identifying drug resistance mechanisms and in pursuing genome-directed, personalized anti-cancer therapy. [ABSTRACT FROM AUTHOR]

Copy number variations (CNVs) are important forms of genetic variation complementary to SNPs, and can be considered as promising markers for some phenotypic and economically important traits or diseases susceptibility in domestic animals. In the present study, we performed a genome-wide CNV identification in 14 individuals selected from diverse populations, including six types of Chinese indigenous breeds, one Asian wild boar population, as well as three modern commercial foreign breeds. We identified 63 CNVRs in total, which covered 9.98 Mb of polymorphic sequence and corresponded to 0.36% of the genome sequence. The length of these CNVRs ranged from 3.20 to 827.21 kb, with an average of 158.37 kb and a median of 97.85 kb. Functional annotation revealed these identified CNVR have important molecular function, and may play an important role in exploring the genetic basis of phenotypic variability and disease susceptibility among pigs. Additionally, to confirm these potential CNVRs, we performed qPCR for 12 randomly selected CNVRs and 8 of them (66.67%) were confirmed successfully. CNVs detected in diverse populations herein are essential complementary to the CNV map in the pig genome, which provide an important resource for studies of genomic variation and the association between various economically important traits and CNVs. [ABSTRACT FROM AUTHOR]

Vesicular trafficking plays a crucial role in protein localization and movement, signal transduction, and multiple developmental processes in eukaryotic cells. Vesicle fusion is the final and key step in vesicle-mediated trafficking and mainly relies on SNAREs (soluble N-ethylmaleimide-sensitive factor attachment protein receptors), the regulators including SM (Sec1/Munc18) family proteins, Rab GTPases and exocyst subunits. Verticillium dahliae is a widespread soil fungus that causes disruptive vascular diseases on a wide range of plants. To date, no genes involved in vesicular fusion process have been identified and characterized in V. dahliae. The recent publication of the draft genome sequence of V. dahliae allowed us to conduct a genome-wide identification, phylogeny and expression profile of genes encoding vesicular fusion components. Using compared genomics and phylogenetic methods, we identified 44 genes encoding vesicle fusion components in the V. dahliae genome. According to the structural features of their encoded proteins, the 44 V. dahliae genes were classified into 22 SNAREs (6 Qa-, 4 Qb-, 6 Qc-, 1 Qbc- and 5 R-types), 4 SM family proteins, 10 Rab GTPases and 8 exocyst proteins. Based on phylogeny and motif constitution analysis, orthologs of vesicle fusion component in filamentous fungi were generally clustered together into the same subclasses with well-supported bootstrap values. Analysis of the expression profiles of these genes indicated that many of them are significantly differentially expressed during vegetative growth and microsclerotia formation in V. dahliae. The analysis show that many components of vesicle fusion are well conserved in filamentous fungi and indicate that vesicle fusion plays a critical role in microsclerotia formation of smoke tree wilt fungus V. dahliae. The genome-wide identification and expression analysis of components involved in vesicle fusion should facilitate research in this gene family and give new insights toward elucidating their functions in growth, development and pathogenesis of V. dahliae. [ABSTRACT FROM AUTHOR]

Gene-based tests of association can increase the power of a genome-wide association study by aggregating multiple independent effects across a gene or locus into a single stronger signal. Recent gene-based tests have distinct approaches to selecting which variants to aggregate within a locus, modeling the effects of linkage disequilibrium, representing fractional allele counts from imputation, and managing permutation tests for p-values. Implementing these tests in a single, efficient framework has great practical value. Fast ASsociation Tests (Fast) addresses this need by implementing leading gene-based association tests together with conventional SNP-based univariate tests and providing a consolidated, easily interpreted report. Fast scales readily to genome-wide SNP data with millions of SNPs and tens of thousands of individuals, provides implementations that are orders of magnitude faster than original literature reports, and provides a unified framework for performing several gene based association tests concurrently and efficiently on the same data. Availability: https://bitbucket.org/baderlab/fast/downloads/FAST.tar.gz, with documentation at https://bitbucket.org/baderlab/fast/wiki/Home [ABSTRACT FROM AUTHOR]

Transcription Activators-Like Effectors (TALEs) belong to a family of virulence proteins from the Xanthomonas genus of bacterial plant pathogens that are translocated into the plant cell. In the nucleus, TALEs act as transcription factors inducing the expression of susceptibility genes. A code for TALE-DNA binding specificity and high-resolution three-dimensional structures of TALE-DNA complexes were recently reported. Accurate prediction of TAL Effector Binding Elements (EBEs) is essential to elucidate the biological functions of the many sequenced TALEs as well as for robust design of artificial TALE DNA-binding domains in biotechnological applications. In this work a program with improved EBE prediction performances was developed using an updated specificity matrix and a position weight correction function to account for the matching pattern observed in a validation set of TALE-DNA interactions. To gain a systems perspective on the large TALE repertoires from X. oryzae strains, this program was used to predict rice gene targets for 99 sequenced family members. Integrating predictions and available expression data in a TALE-gene network revealed multiple candidate transcriptional targets for many TALEs as well as several possible instances of functional convergence among TALEs. [ABSTRACT FROM AUTHOR]

Olfactory receptors (OR), responsible for detection of odor molecules, belong to the largest family of genes and are highly polymorphic in nature having distinct polymorphisms associated with specific regions around the globe. Since there are no reports on the presence of copy number variations in OR repertoire of Indian population, the present investigation in 43 Indians along with 270 HapMap and 31 Tibetan samples was undertaken to study genome variability and evolution. Analysis was performed using Affymetrix Genome-Wide Human SNP Array 6.0 chip, Affymterix CytoScan® High-Density array, HD-CNV, and MAFFT program. We observed a total of 1527 OR genes in 503 CNV events from 81.3% of the study group, which includes 67.6% duplications and 32.4% deletions encompassing more of genes than pseudogenes. We report human genotypic variation in functional OR repertoire size across populations and it was found that the combinatorial effect of both “orthologous obtained from closely related species” and “paralogous derived sequences” provide the complexity to the continuously occurring OR CNVs. [ABSTRACT FROM AUTHOR]

The early elongation checkpoint regulated by Positive Transcription Elongation Factor b (P-TEFb) is a critical control point for the expression of many genes. Spt5 interacts directly with RNA polymerase II and has an essential role in establishing this checkpoint, and also for further transcript elongation. Here we demonstrate that Drosophila Spt5 interacts both physically and genetically with the Polycomb Group (PcG) protein Pleiohomeotic (Pho), and the majority of Pho binding sites overlap with Spt5 binding sites across the genome in S2 cells. Our results indicate that Pho can interact with Spt5 to regulate transcription elongation in a gene specific manner. [ABSTRACT FROM AUTHOR]

Vesicular trafficking plays a crucial role in protein localization and movement, signal transduction, and multiple developmental processes in eukaryotic cells. Vesicle fusion is the final and key step in vesicle-mediated trafficking and mainly relies on SNAREs (soluble N-ethylmaleimide-sensitive factor attachment protein receptors), the regulators including SM (Sec1/Munc18) family proteins, Rab GTPases and exocyst subunits. Verticillium dahliae is a widespread soil fungus that causes disruptive vascular diseases on a wide range of plants. To date, no genes involved in vesicular fusion process have been identified and characterized in V. dahliae. The recent publication of the draft genome sequence of V. dahliae allowed us to conduct a genome-wide identification, phylogeny and expression profile of genes encoding vesicular fusion components. Using compared genomics and phylogenetic methods, we identified 44 genes encoding vesicle fusion components in the V. dahliae genome. According to the structural features of their encoded proteins, the 44 V. dahliae genes were classified into 22 SNAREs (6 Qa-, 4 Qb-, 6 Qc-, 1 Qbc- and 5 R-types), 4 SM family proteins, 10 Rab GTPases and 8 exocyst proteins. Based on phylogeny and motif constitution analysis, orthologs of vesicle fusion component in filamentous fungi were generally clustered together into the same subclasses with well-supported bootstrap values. Analysis of the expression profiles of these genes indicated that many of them are significantly differentially expressed during vegetative growth and microsclerotia formation in V. dahliae. The analysis show that many components of vesicle fusion are well conserved in filamentous fungi and indicate that vesicle fusion plays a critical role in microsclerotia formation of smoke tree wilt fungus V. dahliae. The genome-wide identification and expression analysis of components involved in vesicle fusion should facilitate research in this gene family and give new insights toward elucidating their functions in growth, development and pathogenesis of V. dahliae. [ABSTRACT FROM AUTHOR]

The cat (Felis silvestris catus) shows significant variation in pelage, morphological, and behavioral phenotypes amongst its over 40 domesticated breeds. The majority of the breed specific phenotypic presentations originated through artificial selection, especially on desired novel phenotypic characteristics that arose only a few hundred years ago. Variations in coat texture and color of hair often delineate breeds amongst domestic animals. Although the genetic basis of several feline coat colors and hair lengths are characterized, less is known about the genes influencing variation in coat growth and texture, especially rexoid – curly coated types. Cornish Rex is a cat breed defined by a fixed recessive curly coat trait. Genome-wide analyses for selection (di, Tajima’s D and nucleotide diversity) were performed in the Cornish Rex breed and in 11 phenotypically diverse breeds and two random bred populations. Approximately 63K SNPs were used in the analysis that aimed to localize the locus controlling the rexoid hair texture. A region with a strong signature of recent selective sweep was identified in the Cornish Rex breed on chromosome A1, as well as a consensus block of homozygosity that spans approximately 3 Mb. Inspection of the region for candidate genes led to the identification of the lysophosphatidic acid receptor 6 (LPAR6). A 4 bp deletion in exon 5, c.250_253_delTTTG, which induces a premature stop codon in the receptor, was identified via Sanger sequencing. The mutation is fixed in Cornish Rex, absent in all straight haired cats analyzed, and is also segregating in the German Rex breed. LPAR6 encodes a G protein-coupled receptor essential for maintaining the structural integrity of the hair shaft; and has mutations resulting in a wooly hair phenotype in humans. [ABSTRACT FROM AUTHOR]

Previous genome-wide association studies (GWAS) in multiple populations identified several genetic loci for coronary heart diseases (CHD). Here we utilized a 2-stage candidate gene association strategy in Chinese Han population to shed light on the putative association between several metabolic-related candidate genes and CHD. At the 1st stage, 190 patients with CHD and 190 controls were genotyped through the MassARRAY platform. At the 2nd stage, a larger sample including 400 patients and 392 controls was genotyped by the High Resolution Melt (HRM) method to confirm or rule out the associations with CHD. MLXIP expression level was quantified by the real time PCR in 65 peripheral blood samples. From the 21 studied single nucleotide polymorphisms (SNPs) of seven candidate genes: MLXIPL, MLXIP, MLX, ADIPOR1, VDR, SREBF1 and NR1H3, only one tag SNP rs4758685 (T→C) was found to be statistically associated with CHD (P-value = 0.02, Odds ratio (OR) of 0.83). After adjustment for the age, sex, lipid levels and diabetes, the association remained significant (P-value = 0.03). After adjustment for the hypertension, P-value became 0.20 although there was a significant difference in the allele distribution between the CHD patients with hypertension and the controls (P-value = 0.04, 406 vs 582). In conclusion, among the 21 tested SNPs, we identified a novel association between rs4758685 of MLXIP gene and CHD. The C allele of common variant rs4758685 interacted with hypertension, and was found to be protective against CHD in both allelic and genotypic models in Chinese Han population. [ABSTRACT FROM AUTHOR]

Along with computational approaches, NGS led technologies have caused a major impact upon the discoveries made in the area of miRNA biology, including novel miRNAs identification. However, to this date all microRNA discovery tools compulsorily depend upon the availability of reference or genomic sequences. Here, for the first time a novel approach, miReader, has been introduced which could discover novel miRNAs without any dependence upon genomic/reference sequences. The approach used NGS read data to build highly accurate miRNA models, molded through a Multi-boosting algorithm with Best-First Tree as its base classifier. It was comprehensively tested over large amount of experimental data from wide range of species including human, plants, nematode, zebrafish and fruit fly, performing consistently with >90% accuracy. Using the same tool over Illumina read data for Miscanthus, a plant whose genome is not sequenced; the study reported 21 novel mature miRNA duplex candidates. Considering the fact that miRNA discovery requires handling of high throughput data, the entire approach has been implemented in a standalone parallel architecture. This work is expected to cause a positive impact over the area of miRNA discovery in majority of species, where genomic sequence availability would not be a compulsion any more. [ABSTRACT FROM AUTHOR]

Waxy maize (Zea mays L. var. certaina Kulesh), with many excellent characters in terms of starch composition and economic value, has grown in China for a long history and its production has increased dramatically in recent decades. However, the evolution and origin of waxy maize still remains unclear. We studied the genetic diversity of Chinese waxy maize including typical landraces and inbred lines by SSR analysis and the results showed a wide genetic diversity in the Chinese waxy maize germplasm. We analyzed the origin and evolution of waxy maize by sequencing 108 samples, and downloading 52 sequences from GenBank for the waxy locus in a number of accessions from genus Zea. A sharp reduction of nucleotide diversity and significant neutrality tests (Tajima’s D and Fu and Li’s F*) were observed at the waxy locus in Chinese waxy maize but not in nonglutinous maize. Phylogenetic analysis indicated that Chinese waxy maize originated from the cultivated flint maize and most of the modern waxy maize inbred lines showed a distinct independent origin and evolution process compared with the germplasm from Southwest China. The results indicated that an agronomic trait can be quickly improved to meet production demand by selection. [ABSTRACT FROM AUTHOR]

Cetaceans are unique in being the only mammals completely adapted to an aquatic environment. This adaptation has required complex changes and sometimes a complete restructuring of physiology, behavior and morphology. Identifying genes that have been subjected to selection pressure during cetacean evolution would greatly enhance our knowledge of the ways in which genetic variation in this mammalian order has been shaped by natural selection. Here, we performed a genome-wide scan for positive selection in the dolphin lineage. We employed models of codon substitution that account for variation of selective pressure over branches on the tree and across sites in a sequence. We analyzed 7,859 nuclear-coding ortholog genes and using a series of likelihood ratio tests (LRTs), we identified 376 genes (4.8%) with molecular signatures of positive selection in the dolphin lineage. We used the cow as the sister group and compared estimates of selection in the cetacean genome to this using the same methods. This allowed us to define which genes have been exclusively under positive selection in the dolphin lineage. The enrichment analysis found that the identified positively selected genes are significantly over-represented for three exclusive functional categories only in the dolphin lineage: segment specification, mesoderm development and system development. Of particular interest for cetacean adaptation to an aquatic life are the following GeneOntology targets under positive selection: genes related to kidney, heart, lung, eye, ear and nervous system development. [ABSTRACT FROM AUTHOR]

The human retrotransposon with the highest copy number is the Alu element. The human genome contains over one million Alu elements that collectively account for over ten percent of our DNA. Full-length Alu elements are randomly distributed throughout the genome in both forward and reverse orientations. However, full-length widely spaced Alu pairs having two Alus in the same (direct) orientation are statistically more prevalent than Alu pairs having two Alus in the opposite (inverted) orientation. The cause of this phenomenon is unknown. It has been hypothesized that this imbalance is the consequence of anomalous inverted Alu pair interactions. One proposed mechanism suggests that inverted Alu pairs can ectopically interact, exposing both ends of each Alu element making up the pair to a potential double-strand break, or “hit”. This hypothesized “two-hit” (two double-strand breaks) potential per Alu element was used to develop a model for comparing the relative instabilities of human genes. The model incorporates both 1) the two-hit double-strand break potential of Alu elements and 2) the probability of exon-damaging deletions extending from these double-strand breaks. This model was used to compare the relative instabilities of 50 deletion-prone cancer genes and 50 randomly selected genes from the human genome. The output of the Alu element-based genomic instability model developed here is shown to coincide with the observed instability of deletion-prone cancer genes. The 50 cancer genes are collectively estimated to be 58% more unstable than the randomly chosen genes using this model. Seven of the deletion-prone cancer genes, ATM, BRCA1, FANCA, FANCD2, MSH2, NCOR1 and PBRM1, were among the most unstable 10% of the 100 genes analyzed. This algorithm may lay the foundation for comparing genetic risks posed by structural variations that are unique to specific individuals, families and people groups. [ABSTRACT FROM AUTHOR]

There is increasing evidence that heritable variation in gene expression underlies genetic variation in susceptibility to disease. Therefore, a comprehensive understanding of the similarity between relatives for transcript variation is warranted—in particular, dissection of phenotypic variation into additive and non-additive genetic factors and shared environmental effects. We conducted a gene expression study in blood samples of 862 individuals from 312 nuclear families containing MZ or DZ twin pairs using both pedigree and genotype information. From a pedigree analysis we show that the vast majority of genetic variation across 17,994 probes is additive, although non-additive genetic variation is identified for 960 transcripts. For 180 of the 960 transcripts with non-additive genetic variation, we identify expression quantitative trait loci (eQTL) with dominance effects in a sample of 339 unrelated individuals and replicate 31% of these associations in an independent sample of 139 unrelated individuals. Over-dominance was detected and replicated for a trans association between rs12313805 and ETV6, located 4MB apart on chromosome 12. Surprisingly, only 17 probes exhibit significant levels of common environmental effects, suggesting that environmental and lifestyle factors common to a family do not affect expression variation for most transcripts, at least those measured in blood. Consistent with the genetic architecture of common diseases, gene expression is predominantly additive, but a minority of transcripts display non-additive effects. [ABSTRACT FROM AUTHOR]

Fis, one of the most important nucleoid-associated proteins, functions as a global regulator of transcription in bacteria that has been comprehensively studied in Escherichia coli K12. Fis also influences the virulence of Salmonella enterica and pathogenic E. coli by regulating their virulence genes, however, the relevant mechanism is unclear. In this report, using combined RNA-seq and chromatin immunoprecipitation (ChIP)-seq technologies, we first identified 1646 Fis-regulated genes and 885 Fis-binding targets in the S. enterica serovar Typhimurium, and found a Fis regulon different from that in E. coli. Fis has been reported to contribute to the invasion ability of S. enterica. By using cell infection assays, we found it also enhances the intracellular replication ability of S. enterica within macrophage cell, which is of central importance for the pathogenesis of infections. Salmonella pathogenicity islands (SPI)-1 and SPI-2 are crucial for the invasion and survival of S. enterica in host cells. Using mutation and overexpression experiments, real-time PCR analysis, and electrophoretic mobility shift assays, we demonstrated that Fis regulates 63 of the 94 Salmonella pathogenicity island (SPI)-1 and SPI-2 genes, by three regulatory modes: i) binds to SPI regulators in the gene body or in upstream regions; ii) binds to SPI genes directly to mediate transcriptional activation of themselves and downstream genes; iii) binds to gene encoding OmpR which affects SPI gene expression by controlling SPI regulators SsrA and HilD. Our results provide new insights into the impact of Fis on SPI genes and the pathogenicity of S. enterica. [ABSTRACT FROM AUTHOR]

As the methodologies available for the detection of positive selection from genomic data vary in terms of assumptions and execution, weak correlations are expected among them. However, if there is any given signal that is consistently supported across different methodologies, it is strong evidence that the locus has been under past selection. In this paper, a straightforward frequentist approach based on the Stouffer Method to combine P-values across different tests for evidence of recent positive selection in common variations, as well as strategies for extracting biological information from the detected signals, were described and applied to high density single nucleotide polymorphism (SNP) data generated from dairy and beef cattle (taurine and indicine). The ancestral Bovinae allele state of over 440,000 SNP is also reported. Using this combination of methods, highly significant (P<3.17×10−7) population-specific sweeps pointing out to candidate genes and pathways that may be involved in beef and dairy production were identified. The most significant signal was found in the Cornichon homolog 3 gene (CNIH3) in Brown Swiss (P = 3.82×10−12), and may be involved in the regulation of pre-ovulatory luteinizing hormone surge. Other putative pathways under selection are the glucolysis/gluconeogenesis, transcription machinery and chemokine/cytokine activity in Angus; calpain-calpastatin system and ribosome biogenesis in Brown Swiss; and gangliosides deposition in milk fat globules in Gyr. The composite method, combined with the strategies applied to retrieve functional information, may be a useful tool for surveying genome-wide selective sweeps and providing insights in to the source of selection. [ABSTRACT FROM AUTHOR]

Background: It is thought that methylcytosine can be inherited through meiosis and mitosis, and that epigenetic variation may be under genetic control or correlation may be caused by neutral drift. However, DNA methylation also varies with tissue, developmental stage, and environmental factors. Eliminating these factors, we analyzed the levels and patterns, diversity and structure of genomic methylcytosine in the xylem of nine natural populations of Chinese white poplar. Principal Findings: On average, the relative total methylation and non-methylation levels were approximately 26.567% and 42.708% (P<0.001), respectively. Also, the relative CNG methylation level was higher than the relative CG methylation level. The relative methylation/non-methylation levels were significantly different among the nine natural populations. Epigenetic diversity ranged from 0.811 (Gansu) to 1.211 (Shaanxi), and the coefficients of epigenetic differentiation (GST = 0.159) were assessed by Shannon’s diversity index. Co-inertia analysis indicated that methylation-sensitive polymorphism (MSP) and genomic methylation pattern (CG-CNG) profiles gave similar distributions. Using a between-group eigen analysis, we found that the Hebei and Shanxi populations were independent of each other, but the Henan population intersected with the other populations, to some degree. Conclusions: Genome methylation in Populus tomentosa presented tissue-specific characteristics and the relative 5′-CCGG methylation level was higher in xylem than in leaves. Meanwhile, the genome methylation in the xylem shows great epigenetic variation and could be fixed and inherited though mitosis. Compared to genetic structure, data suggest that epigenetic and genetic variation do not completely match. [ABSTRACT FROM AUTHOR]

Tandem and segmental duplications significantly contribute to gene family expansion and genome evolution. Genome-wide identification of tandem and segmental genes has been analyzed before in several plant genomes. However, comparative studies in functional bias, expression divergence and their roles in species domestication are still lacking. We have carried out a genome-wide identification and comparative analysis of tandem and segmental genes in the rice genome. A total of 3,646 and 3,633 pairs of tandem and segmental genes, respectively, were identified in the genome. They made up around 30% of total annotated rice genes (excluding transposon-coding genes). Both tandem and segmental duplicates showed different physical locations and exhibited a biased subset of functions. These two types of duplicated genes were also under different functional constrains as shown by nonsynonymous substitutions per site (Ka) and synonymous substitutions per site (Ks) analysis. They are also differently regulated depending on the tissues and abiotic and biotic stresses based on transcriptomics data. The expression divergence might be related to promoter differentiation and DNA methylation status after tandem or segmental duplications. Both tandem and segmental duplications differ in their contribution to genetic novelty but evidence suggests that they play their role in species domestication and genome evolution. [ABSTRACT FROM AUTHOR]

Next-generation sequencing technologies expedited research to develop efficient computational tools for the identification of structural variants (SVs) and their use to study human diseases. As deeper data is obtained, the existence of higher complexity SVs in some genomes becomes more evident, but the detection and definition of most of these complex rearrangements is still in its infancy. The full characterization of SVs is a key aspect for discovering their biological implications. Here we present a pipeline (PeSV-Fisher) for the detection of deletions, gains, intra- and inter-chromosomal translocations, and inversions, at very reasonable computational costs. We further provide comprehensive information on co-localization of SVs in the genome, a crucial aspect for studying their biological consequences. The algorithm uses a combination of methods based on paired-reads and read-depth strategies. PeSV-Fisher has been designed with the aim to facilitate identification of somatic variation, and, as such, it is capable of analysing two or more samples simultaneously, producing a list of non-shared variants between samples. We tested PeSV-Fisher on available sequencing data, and compared its behaviour to that of frequently deployed tools (BreakDancer and VariationHunter). We have also tested this algorithm on our own sequencing data, obtained from a tumour and a normal blood sample of a patient with chronic lymphocytic leukaemia, on which we have also validated the results by targeted re-sequencing of different kinds of predictions. This allowed us to determine confidence parameters that influence the reliability of breakpoint predictions. Availability: PeSV-Fisher is available at http://gd.crg.eu/tools. [ABSTRACT FROM AUTHOR]

Expression of bacterial genes takes place under the control of RNA polymerase with exchangeable σ-subunits and multiple transcription factors. A typical promoter region contains one or several overlapping promoters. In the latter case promoters have the same or different σ-specificity and are often subjected to different regulatory stimuli. Genes, transcribed from multiple promoters, have on average higher expression levels. However, recently in the genome of Escherichia coli we found 78 regions with an extremely large number of potential transcription start points (promoter islands, PIs). It was shown that all PIs interact with RNA polymerase in vivo and are able to form transcriptionally competent open complexes both in vitro and in vivo but their transcriptional activity measured by oligonucleotide microarrays was very low, if any. Here we confirmed transcriptional defectiveness of PIs by analyzing the 5′-end specific RNA-seq data, but showed their ability to produce short oligos (9–14 bases). This combination of functional properties indicated a deliberate suppression of transcriptional activity within PIs. According to our data this suppression may be due to a specific conformation of the DNA double helix, which provides an ideal platform for interaction with both RNA polymerase and the histone-like nucleoid protein H-NS. The genomic DNA of E.coli contains therefore several dozen sites optimized by evolution for staying in a heterochromatin-like state. Since almost all promoter islands are associated with horizontally acquired genes, we offer them as specific components of bacterial evolution involved in acquisition of foreign genetic material by turning off the expression of toxic or useless aliens or by providing optimal promoter for beneficial genes. The putative molecular mechanism underlying the appearance of promoter islands within recipient genomes is discussed. [ABSTRACT FROM AUTHOR]

Keloid is benign fibroproliferative dermal tumors with unknown etiology. Recently, a genome-wide association study (GWAS) in Japanese population has identified 3 susceptibility loci (rs873549 at 1q41, rs940187 and rs1511412 at 3q22.3, rs8032158 at 15p21.3) for keloid. In order to examine whether these susceptibility loci are associated with keloid in the Chinese Han population, twelve previously reported SNPs were selected for replication in 714 cases and 2,944 controls by using Sequenom MassArray system. We found three SNPs in two regions showed significant association with keloid in the Chinese Han population: 1q41 (rs873549, P = 3.03×10−33, OR = 2.05, 95% CI: 1.82–2.31 and rs1442440, P = 9.85×10−18, OR = 0.56, 95% CI: 0.49–0.64, respectively) and 15q21.3 (rs2271289 located in NEDD4, P = 1.02×10−11, OR = 0.66, 95% CI: 0.58–0.74). We also detected one risk haplotype AG (P = 1.36×10−31, OR = 2.02) and two protective haplotypes of GA and AA (GA, P = 1.94×10−19, OR = 0.53, AA, P = 0.00043, OR = 0.78, respectively) from the two SNPs (rs873549 and rs1442440). Our study confirmed two previously reported loci 1q41 and 15q21.3 for keloid in the Chinese Han population, which suggested the common genetic factor predisposing to the development of keloid shared by the Chinese Han and Japanese populations. [ABSTRACT FROM AUTHOR]