Genome-wide evolution analysis reveals low CpG contents of fast-evolving genes and identifies antiviral microRNAs

Noncoding RNAs(ncRNAs) play important roles in many biological processes and provide materials for evolutionary adaptations beyond protein-coding genes, such as in the arms race between the host and pathogen. However, currently, a comprehensive high-resolution analysis of primate genomes that includes the latest annotated ncRNAs is not available. Here, we developed a computational pipeline to estimate the selections that act on noncoding regions based on comparisons with a large number of reference sequences in introns adjacent to the interested regions. Our method yields result comparable with those of the established codon-based method and phyloP method for coding genes; thus, it provides a holistic framework for estimating the selection on the entire genome. We further showed that fastevolving protein-coding genes and their corresponding 50 UTRs have a significantly lower frequency of the CpG dinucleotides than those evolving at an average pace, and these fast-evolving genes are enriched in the process of immunity and host defense. We also identified fast-evolving miRNAs with antiviral functions in cells. Our results provide a resource for high-resolution evolution analysis of the primate genomes.     

Sequence context at human single nucleotide polymorphisms: overrepresentation of CpG dinucleotide at polymorphic sites and suppression of variation in CpG islands

Human polymorphisms originate as mutations, and the influence of context on mutagenesis should be reflected in the distribution of sequences surrounding single nucleotide polymorphisms (SNPs). We have performed a computational survey of nearly two million human SNPs to determine if sequence-dependent hotspots for polymorphism exist in the human genome. Here we show that sequences containing CpG dinucleotides, which occur at low frequencies in the human genome, are 6.7-fold more abundant at polymorphic sites than expected. In contrast, polymorphisms in CpG sequences located within CpG islands, important regulatory regions that modulate gene expression, are 6.8-fold less prevalent than expected. The distribution of polymorphic alleles at CpGs in CpG islands is also significantly different from that in non-island regions. These data strongly support a role for 5-methylcytosine deamination in the generation of human variation, and suggest that variation at CpGs in islands is suppressed.     

奶牛产奶性状与乳房炎相关基因CpG含量及分布特征的比较分析

DNA甲基化是表观遗传学的重要组成部分,基因启动子区及第一外显子区的CpG甲基化通常抑制该基因的表达,而去甲基化则促进基因表达。已有的研究发现荷斯坦牛的乳房炎指标SCC(Somatic cell count)与产奶量呈较强负相关。文章分析并比较了这两类性状的相关基因的启动子区、第一外显子、下游2 000 bp序列中CpG含量及分布特征。结果表明,乳房炎相关基因的启动子、第一外显子中CpG含量显著低于产奶性状相关基因,而两类性状基因下游2 000 bp序列中CpG含量无显著性差异。另外,文中提出了两个量化基因序列中CpG特征的指标,一个是CpG平均距离,用来衡量序列中的CpG分布;另一个是条件概率p(G|C),用以量化序列中二核苷酸CpG随碱基C出现的可能性,并对两类基因的启动子和第一外显子区域的这两个指标做了统计检验。研究结果对产奶性状与乳房炎相关基因的DNA甲基化调控研究奠定了基础。     

Genome-wide identification of runs of homozygosity islands and associated genes in local dairy cattle breeds

Runs of homozygosity (ROH) are widely used as predictors of whole-genome inbreeding levels in cattle. They identify regions that have an unfavorable effect on a phenotype when homozygous, but also identify the genes associated with traits of economic interest present in these regions. Here, the distribution of ROH islands and enriched genes within these regions in four dairy cattle breeds were investigated. Cinisara (71), Modicana (72), Reggiana (168) and Italian Holstein (96) individuals were genotyped using the 50K v2 Illumina BeadChip. The genomic regions most commonly associated with ROHs were identified by selecting the top 1% of the single nucleotide polymorphisms (SNPs) most commonly observed in the ROH of each breed. In total, 11 genomic regions were identified in Cinisara and Italian Holstein, and eight in Modicana and Reggiana, indicating an increased ROH frequency level. Generally, ROH islands differed between breeds. The most homozygous region (>45% of individuals with ROH) was found in Modicana on chromosome 6 within a quantitative trail locus affecting milk fat and protein concentrations. We identified between 126 and 347 genes within ROH islands, which are involved in multiple signaling and signal transduction pathways in a wide variety of biological processes. The gene ontology enrichment provided information on possible molecular functions, biological processes and cellular components under selection related to milk production, reproduction, immune response and resistance/susceptibility to infection and diseases. Thus, scanning the genome for ROH could be an alternative strategy to detect genomic regions and genes related to important economic traits.     

Genome-wide analysis of DNA methylation status of CpG islands in embryoid bodies, teratomas, and fetuses

Differentiation of embryonic stem (ES) cells into embryoid bodies (EBs) provides an in vitro system for the study of early lineage determination during mammalian development. We have previously reported that there are 247 CpG islands that potentially have tissue-dependent and differentially methylated regions (T-DMRs). This provided evidence that the formation of DNA methylation patterns at CpG islands is a crucial epigenetic event underlying mammalian development. Here we present an analysis by the restriction landmark genomic scanning (RLGS) using NotI as a landmark enzyme of the genome-wide methylation status of CpG islands of ES cells and EBs and of teratomas produced from ES cells. These results are considered in relation to the methylation status of CpG islands of genomic DNA from normal fetus (10.5 dpc) and adult tissues. We have prepared a DNA methylation panel that consists of 259 T-DMRs and includes novel T-DMRs that are distinctly methylated or unmethylated in the teratomas. The DNA methylation pattern was complex and differed for the ES cells, EBs, and teratomas, providing evidence that differentiation of cells involves both de novo DNA methylation as well as demethylation. Comparison of the numbers of T-DMRs, that were differentially methylated or unmethylated among the cells and tissue types studied, revealed that the teratomas were the most epigenetically different from ES cells. Thus, analysis of the DNA methylation profiles prepared in this study provides new insights into the differentiation of ES cells and development of fetus, EB, teratoma, and somatic tissues.     

High‐resolution analysis of DNA synthesis start sites and nucleosome architecture at efficient mammalian replication origins

DNA replication origins are poorly characterized genomic regions that are essential to recruit and position the initiation complex to start DNA synthesis. Despite the lack of specific replicator sequences, initiation of replication does not occur at random sites in the mammalian genome. This has lead to the view that DNA accessibility could be a major determinant of mammalian origins. Here, we performed a high‐resolution analysis of nucleosome architecture and initiation sites along several origins of different genomic location and firing efficiencies. We found that mammalian origins are highly variable in nucleosome conformation and initiation patterns. Strikingly, initiation sites at efficient CpG island‐associated origins always occur at positions of high‐nucleosome occupancy. Origin recognition complex (ORC) binding sites, however, occur at adjacent but distinct positions marked by labile nucleosomes. We also found that initiation profiles mirror nucleosome architecture, both at endogenous origins and at a transgene in a heterologous system. Our studies provide a unique insight into the relationship between chromatin structure and initiation sites in the mammalian genome that has direct implications for how the replication programme can be accommodated to diverse epigenetic scenarios.     

比较分析草菇与其他真菌直系同源基因的关系

担子菌类的食用菌种类多、价值高、产量大,然而其产业的升级发展需要对食用菌生长发育相关生物学问题进行深入解析。目前多种食用菌完成了全基因组测序,然而作为非模式种其与模式丝状真菌间的直系同源基因目前尚缺乏全基因组水平的系统研究,在一定程度上限制了其分子生物学研究的深入。本研究以草菇为参照物种,将其与几种食用菌和模式丝状真菌进行两两直系同源基因分析,并对多物种间不同类型的直系同源基因进行功能富集。结果显示：一对一直系同源基因较多富集于基因复制、转录、翻译、修饰、加工等保守的基本功能类别;非一对一直系同源基因多属于基因家族,且包含了65%的转录因子,功能上富集在碳水化合物、脂质、氨基酸、次生代谢物及外源物质的代谢通路。无直系同源基因则较多富集在与基因重组、修复、信号转导相关的功能类别、特导性转录因子以及未知的预测基因。结果为食用菌分子生物学的深入研究提供有价值的参考。     

Genomic analysis of a migratory divide reveals candidate genes for migration and implicates selective sweeps in generating islands of differentiation

Differential gene flow, reductions in diversity following linked selection and/or features of the genome can structure patterns of genomic differentiation during the process of speciation. Possible sources of reproductive isolation are well studied between coastal and inland subspecies groups of Swainson's thrushes, with differences in seasonal migratory behaviour likely playing a key role in reducing hybrid fitness. We assembled and annotated a draft reference genome for this species and generated whole‐genome shotgun sequence data for populations adjacent to the hybrid zone between these groups. We documented substantial genomewide heterogeneity in relative estimates of genetic differentiation between the groups. Within population diversity was lower in areas of high relative differentiation, supporting a role for selective sweeps in generating this pattern. Absolute genetic differentiation was reduced in these areas, further suggesting that recurrent selective sweeps in the ancestral population and/or between divergent populations following secondary contact likely occurred. Relative genetic differentiation was also higher near centromeres and on the Z chromosome, suggesting that features of the genome also contribute to genomewide heterogeneity. Genes linked to migratory traits were concentrated in islands of differentiation, supporting previous suggestions that seasonal migration is under divergent selection between Swainson's thrushes. Differences in migratory behaviour likely play a central role in the speciation of many taxa; we developed the infrastructure here to permit future investigations into the role several candidate genes play in reducing gene flow between not only Swainson's thrushes but other species as well.     

Genomic analysis of Fas and FasL genes and absence of correlation with disease progression in AIDS

Apoptosis has been suggested as a major mechanism for the CD4⁺ T-lymphocyte depletion observed in patients infected with human immunodeficiency virus 1 (HIV-1). To evaluate the impact of genetic variations to apoptosis during progression of acquired immunodeficiency syndrome (AIDS), we have performed an extensive genetic analysis of Fas and Fas ligand (FasL) genes. The coding regions and promoters of these genes were resequenced in a cohort of 212 HIV-1-seropositive patients presenting extreme disease phenotypes and 155 healthy controls of Caucasian origin. Overall, 33 single nucleotide polymorphisms (SNPs) with an allele frequency >1% were identified and evaluated for their association with disease progression. Among them, 14 polymorphisms were newly characterized. We did not find any statistically significant association of Fas and FasL polymorphisms and haplotypes with AIDS progression.     

CpG island hypermethylation of multiple tumor suppressor genes associated with loss of their protein expression during rat lung carcinogenesis induced by 3-methylcholanthrene and diethylnitrosamine

The epigenetic mechanisms underlying the tumorigenesis caused by polycyclic aromatic hydrocarbons and nitrosamine compounds such as 3-methylcholanthrene (MCA) and diethylnitrosamine (DEN) are currently unknown. We reported previously that dynamic changes in DNA methylation occurred during MCA/DEN-induced rat lung carcinogenesis. Here, we used the same animal model to further study the evolution of methylation alterations in tumor suppressor genes (TSGs) DAPK1, FHIT, RASSF1A, and SOCS-3. We found that none of these genes were methylated in either normal or hyperplasia tissue. However, as the severity of the cancer progressed through squamous metaplasia and dysplasia to carcinoma in situ (CIS) and infiltrating carcinoma, so methylation became more prevalent. Particularly dramatic increases in the level of methylation, the average number of methylated genes, and the incidence of concurrent methylation in three genes were observed in CIS and infiltrating carcinoma. Similar but less profound changes were seen in squamous metaplasia and dysplasia. Furthermore, methylation status was closely correlated to loss of protein expression for these genes, with protein levels markedly declining along the continuum of carcinogenesis. These results suggest that progressive CpG island hypermethylation leading to inactivation of TSGs might be a vital molecular mechanism in the pathogenesis of MCA/DEN-induced multistep rat lung carcinogenesis.     

Altered tissue specificity of the revertant shrunken allele upon Dissociation (Ds) excision is associated with loss of expression and molecular rearrangement at the corresponding non-allelic isozyme locus in maize

Summary Transposable element Activator (Ac) induced wild-type stable revertants, derived from McClintock's Dissociation (Ds) insertion shrunken (sh) mutant sh-m5933, have been examined for sucrose synthases, SS1 and SS2, encoded by the revertant (Sh) locus and the non-allelic gene Sus (previously designated as Ss2), respectively. A structurally normal Sh locus has been previously described in these revertants. Immuno-blot (Western) and Southern hybridization analyses reported here identify one of the nine alleles, Sh-r5, as unique for several features. It showed altered tissue specificity, as the SS1 protein encoded by the Sh-r5 allele was readily detectable in the immature embryo which is otherwise characterized by the Sus expression only. The level of Sh-r5 expression at the protein and enzyme level was marked by endosperm specific SS1 abundance and a significant down-regulation in the embryo similar to the standard Sh and Sus loci in endosperm and embryo, respectively. We infer that tissue specific levels of gene expression among maize Ss genes is significantly determined by trans-regulatory factors present in these two tissues. The Sh-r5 strain also exhibited a complete loss of the Sus expression in all tissues tested in the plant. Lack of any detectable phenotypic abnormality in the Sh-r5 strain due to the loss of SS2 protein indicated that either the SS2 protein is nonessential or that the two SS isozymes are functionally compensatory. Genomic filter hybridizations with the Sus cDNA clone indicated that the Sus locus in the Sh-r5 strain was not deleted and was, in fact, unique among these revertants. Together, these data provide an unusual insight into the regulation and function of the two SS isozymes in the maize plant.     

Differentiation between isolates of Helicobacter pylori by PCR-RFLP analysis of urease A and B genes and comparison with ribosomal RNA gene patterns

Abstract Genetic diversity amongst 21 human gastric isolates of Helicobacter pylori was investigated by polymerase chain reaction amplification and Hae III digest (restriction fragment length polymorphism) analysis of an internal 2.4-kb segment of the urease A and urease B genes. H. pylori from 11 independent individuals yielded nine distinct restriction fragment patterns but only one pattern was common to H. pylori from two individuals. By contrast, multiple isolate sets of H. pylori from two patients each had common urease gene patterns. Most strains with the same urease gene patterns were distinguishable in their ribosomal RNA gene patterns. The study demonstrated diversity amongst H. pylori and established that PCR analysis of urease genes provided a novel method of identifying isolates. The profiles were reproducible and convenient to obtain and analyse, and were almost as discriminatory as Hae III ribopatterns.     

Haplotype analysis of key genes governing grain yield and quality traits across 3K RG panel reveals scope for the development of tailor‐made rice with enhanced genetic gains

Though several genes governing various major traits have been reported in rice, their superior haplotype combinations for developing ideal variety remains elusive. In this study, haplotype analysis of 120 previously functionally characterized genes, influencing grain yield (87 genes) and grain quality (33 genes) revealed significant variations in the 3K rice genome (RG) panel. For selected genes, meta‐expression analysis using already available datasets along with co‐expression network provided insights at systems level. Also, we conducted candidate gene based association study for the 120 genes and identified 21 strongly associated genes governing 10‐grain yield and quality traits. We report superior haplotypes upon phenotyping the subset of 3K RG panel, SD1‐H8 with haplotype frequency (HF) of 30.13% in 3K RG panel, MOC1‐H9 (HF: 23.08%), IPA1‐H14 (HF: 6.64%), DEP3‐H2 (HF: 5.59%), DEP1‐H2 (HF: 37.53%), SP1‐H3 (HF: 5.05%), LAX1‐H5 (HF: 1.56%), LP‐H13 (3.64%), OSH1‐H4 (5.52%), PHD1‐H14 (HF: 15.21%), AGO7‐H15 (HF: 3.33%), ROC5‐H2 (31.42%), RSR1‐H8 (HF: 4.20%) and OsNAS3‐H2 (HF: 1.00%). For heading date, Ghd7‐H8 (HF: 3.08%), TOB1‐H10 (HF: 4.60%) flowered early, Ghd7‐H14 (HF: 42.60%), TRX1‐H9 (HF: 27.97%), OsVIL3‐H14 (HF: 1.72%) for medium duration flowering, while Ghd7‐H6 (HF: 1.65%), SNB‐H9 (HF: 9.35%) were late flowering. GS5‐H4 (HF: 65.84%) attributed slender, GS5‐H5 (HF: 29.00%), GW2‐H2 (HF: 4.13%) were medium slender and GS5‐H9 (HF: 2.15%) for bold grains. Furthermore, haplotype analysis explained possible genetic basis for superiority of selected mega‐varieties. Overall, this study suggests the possibility for developing next‐generation tailor‐made rice with superior haplotype combinations of target genes suiting future food and nutritional demands via haplotype‐based breeding.     

The analysis of core and symbiotic genes of rhizobia nodulating Vicia from different continents reveals their common phylogenetic origin and suggests the distribution of Rhizobium leguminosarum strains together with Vicia seeds

In this work, we analysed the core and symbiotic genes of rhizobial strains isolated from Vicia sativa in three soils from the Northwest of Spain, and compared them with other Vicia endosymbionts isolated in other geographical locations. The analysis of rrs, recA and atpD genes and 16S–23S rRNA intergenic spacer showed that the Spanish strains nodulating V. sativa are phylogenetically close to those isolated from V. sativa and V. faba in different European, American and Asian countries forming a group related to Rhizobium leguminosarum. The analysis of the nodC gene of strains nodulating V. sativa and V. faba in different continents showed they belong to a phylogenetically compact group indicating that these legumes are restrictive hosts. The results of the nodC gene analysis allow the delineation of the biovar viciae showing a common phylogenetic origin of V. sativa and V. faba endosymbionts in several continents. Since these two legume species are indigenous from Europe, our results suggest a world distribution of strains from R. leguminosarum together with the V. sativa and V. faba seeds and a close coevolution among chromosome, symbiotic genes and legume host in this Rhizobium–Vicia symbiosis.