Decoupled differentiation of gene expression and coding sequence among Drosophila populations

Owing to the relevance to evolutionary theories of genotypic and phenotypic evolution, the correspondence of differentiation among natural populations in complex phenotypic traits and genetic markers has been studied extensively, and generally found to be poor. In contrast, the correspondence of differentiation among natural populations in gene expression, now often considered a genomic era proxy for the phenotype, and genetic markers, remains largely unexplored. Here, an analysis of expression and nucleotide sequence polymorphism of 106 genes in Drosophila melanogaster strains of the Cosmopolitan (M) and Zimbabwe, Africa (Z) mating races showed that differentiation of gene expression and of coding sequences, measured as QST and GST, respectively, were uncorrelated and, generally, QST > GST. However, an exploratory analysis showed that GST of the 5 prime sequences of genes was correlated with QST calculated from expression data, while GST of the coding sequences remained uncorrelated with QST. This scenario is consistent with the population differentiation at cis-regulatory regions that is decoupled from differentiation of the coding regions. However, despite evidence for selection on global levels of gene expression (deduced from QST > GST), 5 prime sequence polymorphisms generally were compatible with selective neutrality, suggesting differentiation in cis-regulated gene expression for these genes has been promoted by drift or selection too weak or too long ago to be detected, or higher organizational levels underlying the genetic architecture of expression are targets of selection. In all, this raises the question how selection on the expression changes (i.e. the phenotype) can be so obvious yet elusive at the level of the nucleotide sequence. Our contrasts between genetic differentiation of populations in expression and sequences revealed that even when genotype and phenotype can be connected the sources of variation that are the target of selection remain to be identified.     

Complex patterns of cis‐regulatory polymorphisms in ebony underlie standing pigmentation variation in Drosophila melanogaster

Pigmentation traits in adult Drosophila melanogaster were used in this study to investigate how phenotypic variations in continuous ecological traits can be maintained in a natural population. First, pigmentation variation in the adult female was measured at seven different body positions in 20 strains from the Drosophila melanogaster Genetic Reference Panel (DGRP) originating from a natural population in North Carolina. Next, to assess the contributions of cis‐regulatory polymorphisms of the genes involved in the melanin biosynthesis pathway, allele‐specific expression levels of four genes were quantified by amplicon sequencing using a 454 GS Junior. Among those genes, ebony was significantly associated with pigmentation intensity of the thoracic segment. Detailed sequence analysis of the gene regulatory regions of this gene indicated that many different functional cis‐regulatory alleles are segregating in the population and that variations outside the core enhancer element could potentially play important roles in the regulation of gene expression. In addition, a slight enrichment of distantly associated SNP pairs was observed in the ~10 kb cis‐regulatory region of ebony, which suggested the presence of interacting elements scattered across the region. In contrast, sequence analysis in the core cis‐regulatory region of tan indicated that SNPs within the region are significantly associated with allele‐specific expression level of this gene. Collectively, the data suggest that the underlying genetic differences in the cis‐regulatory regions that control intraspecific pigmentation variation can be more complex than those of interspecific pigmentation trait differences, where causal genetic changes are typically confined to modular enhancer elements.     

Adaptive protein evolution and regulatory divergence in Drosophila

Two recent studies demonstrated a positive correlation between divergence in gene expression and protein sequence in Drosophila. This correlation could be driven by positive selection or variation in functional constraint. To distinguish between these alternatives, we compared patterns of molecular evolution for 1,862 genes with two previously reported estimates of expression divergence in Drosophila. We found a slight negative trend (nonsignificant) between positive selection on protein sequence and divergence in expression levels between Drosophila melanogaster and Drosophila simulans. Conversely, shifts in expression patterns during Drosophila development showed a positive association with adaptive protein evolution, though as before the relationship was weak and not significant. Overall, we found no strong evidence for an increase in the incidence of positive selection on protein-coding regions in genes with divergent expression in Drosophila, suggesting that the previously reported positive association between protein and regulatory divergence primarily reflects variation in functional constraint.     

Advanced intercross line mapping of Eae5 reveals Ncf-1 and CLDN4 as candidate genes for experimental autoimmune encephalomyelitis

Eae5 in rats was originally identified in two F(2) intercrosses, (DA x BN) and (E3 x DA), displaying linkage to CNS inflammation and disease severity in experimental autoimmune encephalomyelitis (EAE), respectively. This region overlaps with an arthritis locus, Pia4, which was also identified in the (E3 x DA) cross. Two congenic strains, BN.DA-Eae5 and BN.DA-Eae5.R1, encompassing the previously described Eae5 and Pia4, were established. DA alleles within the chromosome 12 fragment conferred an increase in disease susceptibility as well as increased inflammation and demyelination in the CNS as compared with BN alleles. To enable a more precise fine mapping of EAE regulatory genes, we used a rat advanced intercross line between the EAE-susceptible DA strain and the EAE-resistant PVG.1AV1 strain. Linkage analysis performed in the advanced intercross line considerably narrowed down the myelin oligodendrocyte glycoprotein-EAE regulatory locus (Eae5) to a approximately 1.3-megabase region with a defined number of candidate genes. In this study we demonstrate a regulatory effect of Eae5 on MOG-EAE by using both congenic strains as well as fine mapping these effects to a region containing Ncf-1, a gene associated with arthritis. In addition to structural polymorphisms in Ncf-1, both sequence polymorphisms and expression differences were identified in CLDN4. CLDN4 is a tight junction protein involved in blood-brain barrier integrity. In conclusion, our data strongly suggests Ncf-1 to be a gene shared between two organ-specific inflammatory diseases with a possible contribution by CLDN4 in encephalomyelitis.     

生防菌Bs-916突变菌株M49芽胞形成和感受态形成能力

摘要：生防枯草芽胞杆菌Bs-916（Bacillus subtilis）在水稻纹枯病的防治上效果显著。应用离子注入突变对Bs-916进行了突变,获得了一系列的突变菌株。其中突变菌株M49,其表面活性素Surfactin分泌量比出发菌株Bs-916大大降低并导致其防效降低。【目的】为了确认影响该菌株防效降低的影响因子,对其表型和相关基因表达水平进行了研究。【方法】应用生孢培养基,通过芽胞形成能力评测方法比较该菌株和野生菌株Bs-916的芽胞形成能力;通过转化质粒的实验评测突变菌株M49和野生型Bs-916的     

Wolbachia interactions that determine Drosophila melanogaster survival

Abstract.— We have recently described a mutualistic symbiosis in which Wolbachia bacteria were shown to improve the fitness of some Drosophila melanogaster stocks. Wolbachia did not extend longevity in all Drosophila genotypes, even though 16s rDNA sequences indicated that our Drosophila stocks were infected with the same Wolbachia strain. Here, we use reciprocal hybrid crosses between two Drosophila strains, one that lived longer with Wolbachia (Z53) and one that did not (Z2), to investigate the inheritance of the survival phenotype and its dependence on the host genotype, sex, and mating conditions. Wolbachia's positive effects were more apparent in hybrid flies than in parental flies, ruling out exclusive maternal inheritance or the dependence of the survival phenotype on Wolbachia strain differences. The Wolbachia survival effects were more apparent in single-sex cages, where courtship and mating were not permitted. In these cages, nearly all flies with Wolbachia lived longer than uninfected flies, even though strain Z2 showed no Wolbachia effect in mixed-sex mating cages. We used comparisons between single- and mixed-sex cages to estimate the cost of reproduction for both sexes. Our data suggest that Wolbachia infection may increase the inferred cost of reproduction, particularly in males. Wolbachia can even produce a positive survival effect almost as large as the negative survival effect associated with reproduction. We discuss the implications of our experiments for the study of insect symbioses.     

Sequence polymorphisms within the pMGA genes and pMGA antigenic variants in Mycoplasma gallisepticum

Antigenic variants of Mycoplasma gallisepticum major surface lipoprotein, pMGA, are encoded by a large gene family. In this study sequence analyses of the PCR-amplified pMGA genes showed two types of sequences similar to the pMGA1.2 gene in M. gallisepticum strains. They differed in the sequence encoding a proline-rich region (PRR) at the N-terminus of the pMGA protein. The type A genes had sequences similar to the published pMGA1.2 gene sequence of strain S6, whereas the type B genes lacked the second repetitive segment encoding PTPN sequence within PRR and were similar to the published sequence of PG31 strain. Low in vitro passages of M. gallisepticum strains isolated recently in Slovenia from four avian species showed very different expression patterns of pMGA1.2 and pMGA1.9 genes. Among isogenic populations of S6(B) and IHB1 strains a high frequency of pMGA antigenic variants lacking an epitope for monoclonal antibody (mAb) 71 was found. Strain IHB1 clones, which synthesized pMGA recognized by mAb 71, transcribed pMGA genes whose partial sequence encoded the amino acid sequence (262)TNGDEPRSVS of the mAb 71 epitope. Other IHB1 clones synthesized pMGA variants with different isoelectric points, lacking the epitope for mAb 71, but expressing downstream epitopes for other mAbs. Our study suggests that a molecular basis for pMGA antigenic variation lies in the corresponding changes at the DNA level.     

Tissue-Specific Expression Phenotypes of Hawaiian Drosophila Adh Genes in Drosophila Melanogaster Transformants

Interspecific differences in the tissue-specific patterns of expression displayed by the alcohol dehydrogenase (Adh) genes within the Hawaiian picture-winged Drosophila represent a rich source of evolutionary variation in gene regulation. Study of the cis-acting elements responsible for regulatory differences between Adh genes from various species is greatly facilitated by analyzing the behavior of the different Adh genes in a homogeneous background. Accordingly, the Adh gene from Drosophila grimshawi was introduced into the germ line of Drosophila melanogaster by means of P element-mediated transformation, and transformants carrying this gene were compared to transformants carrying the Adh genes from Drosophila affinidisjuncta and Drosophila hawaiiensis. The results indicate that the D. affinidisjuncta and D. grimshawi genes have relatively higher levels of expression and broader tissue distribution of expression than the D. hawaiiensis gene in larvae. All three genes are expressed at similar overall levels in adults, with differences in tissue distribution of enzyme activity corresponding to the pattern in the donor species. However, certain systematic differences between Adh gene expression in transformants and in the Hawaiian Drosophila are noted along with tissue-specific position effects in some cases. The implications of these findings for the understanding of evolved regulatory variation are discussed.     

Integration of QTL and bioinformatic tools to identify candidate genes for triglycerides in mice

To identify genetic loci influencing lipid levels, we performed quantitative trait loci (QTL) analysis between inbred mouse strains MRL/MpJ and SM/J, measuring triglyceride levels at 8 weeks of age in F2 mice fed a chow diet. We identified one significant QTL on chromosome (Chr) 15 and three suggestive QTL on Chrs 2, 7, and 17. We also carried out microarray analysis on the livers of parental strains of 282 F2 mice and used these data to find cis-regulated expression QTL. We then narrowed the list of candidate genes under significant QTL using a "toolbox" of bioinformatic resources, including haplotype analysis; parental strain comparison for gene expression differences and nonsynonymous coding single nucleotide polymorphisms (SNP); cis-regulated eQTL in livers of F2 mice; correlation between gene expression and phenotype; and conditioning of expression on the phenotype. We suggest Slc25a7 as a candidate gene for the Chr 7 QTL and, based on expression differences, five genes (Polr3 h, Cyp2d22, Cyp2d26, Tspo, and Ttll12) as candidate genes for Chr 15 QTL. This study shows how bioinformatics can be used effectively to reduce candidate gene lists for QTL related to complex traits.     

Deletion of the Schizophyllum Commune Aα Locus: The Roles of Aα Y and Z Mating-Type Genes

The Aα locus is one of four master regulatory loci that determine mating type and regulate sexual development in Schizophyllum commune. We have made a plasmid containing a URA1 gene disruption of the Aα Y1 gene. Y1 is the sole Aα gene in Aα1 strains. We used the plasmid construction to produce an Aα null (i.e., AαΔ) strain by replacing the genomic Y1 gene with URA1 in an Aα1 strain. To characterize the role of the Aα genes in the regulation of sexual development, we transformed various Aα Y and Z alleles into AαΔ strains and examined the acquired mating types and mating abilities of the transformants. These experiments demonstrate that the Aα Y gene is not essential for fungal viability and growth, that a solitary Z Aα mating-type gene does not itself activate development, that Aβ proteins are sufficient to activate the A developmental pathway in the absence of Aα proteins and confirm that Y and Z genes are the sole determinants of Aα mating type. The data from these experiments support and refine our model of the regulation of A-pathway events by Y and Z proteins.     

Structural analysis of the Drosophila rpA1 gene, a member of the eucaryotic ''A'' type ribosomal protein family.

The expression of ribosomal protein (r-protein) genes is uniquely regulated at the translational level during early development of Drosophila. Here we report results of a detailed analysis of the r-protein rpA1 gene. A cloned DNA sequence coding for rpA1 has been identified by hybrid-selected translation and amino acid composition analysis. The rpA1 gene was localized to polytene chromosome band 53CD. The nucleotide sequence of the rpA1 gene and its cDNA have been determined. rpA1 is a single copy gene and sequence comparison between the gene and its cDNA indicates that this r-protein gene is intronless. Allelic restriction site polymorphisms outside of the gene were observed, while the coding sequence is well conserved between two Drosophila strains. The protein has unusual domains rich in Ala and charged residues. The rpA1 is homologous to the "A" family of eucaryotic acidic r-proteins which are known to play a key role in the initiation and elongation steps of protein synthesis.     

DNA-based classification and sequence heterogeneities in the 16S rRNA genes of Lactobacillus casei/paracasei and related species

The sequence differences within the 16S rRNA genes of Lactobacillus casei/paracasei and related species, Lactobacillus zeae and Lactobacillus rhamnosus, were investigated. Thirty-seven strains of mostly human or cheese origin were grouped by restriction endonuclease analysis (REA) of the total chromosomal DNA and by temporal temperature gradient gel electrophoresis (TTGE) of PCR-amplified 16S rRNA gene fragments. REA verified that all strains were genomically unique and singled out three major clusters, one L. rhamnosus-cluster and two clusters containing L. paracasei strains. The groups obtained by TTGE corresponded with one exception to the REA-clusters. In the TTGE clustering all L. paracasei strains formed one general group with one TTGE-band in common, and this group was sub-divided into five subgroups due to the presence of more than one TTGE-band in four of the subgroups. The occurrence of multiple TTGE-bands was investigated by amplifying and cloning of the 16S rRNA genes from the strains showing this phenomenon, thereby 12 clones from each strain were sequenced, demonstrating polymorphisms in almost all the cases. Subjecting the clones displaying sequence variations to TTGE as well as sequencing of 16S rDNA revealed by ribotyping of the strains, verified the presence of polymorphisms within the 16S rRNA genes. The migration characteristic of amplified DNA from a single clone corresponded to a specific band in the TTGE-pattern of the strain from which the clone originated. Southern blot hybridisation with a 16S rDNA probe demonstrated the presence of at least five 16S rRNA genes in L. casei/paracasei. A higher degree of variable positions than previously reported was observed in the 16S rRNA gene fragments of the members in the complex. Sequence comparison between the 16S rRNA gene copies of L. casei (CCUG 21451T) and L. zeae (CCUG 35515T) demonstrated that the two species shared almost the same sequence in some copies while the others were more different. Our results provide one explanation for the difficulties in reaching clear-cut taxa within the L. casei/paracasei complex.     

格特隐球菌核基因、交配型位点和线粒体基因的多位点序列分型及重组分析

研究格特隐球菌VGI和VGII基因型菌株间线粒体基因重组与核基因、交配型位点重组间的关系。采用分子鉴定区分受试格特隐球菌的血清型、基因型和交配型;选取包含核基因和线粒体基因10个位点的多位点序列分型（ MLST）进行基因型和系统发育分析;采用同质性检验评价基因系谱间的一致性。多位点的序列和系统发育分析结合同质性检验表明,受试位点中,仅ATP6位点发现VGI和VGII菌株间基因的杂交和重组,该基因系谱与其余位点的基因系谱呈现不一致性。结果显示,受试的部分VGI菌株中的ATP6位点含有VGII菌株的基因序列,表明VGI和VGII菌株间线粒体基因的重组;没有发现核基因及交配型位点中两者间的重组,提示VGI和VGII菌株间线粒体基因的重组现象与交配行为无关。     

Mating Genes of the Trichophyton mentagrophytes Complex

The mating type (-)-specific gene of the alpha-box and the mating type (+)-specific gene of the high-mobility group (HMG) DNA-binding domain were confirmed in zoophilic dematophytes of Arthroderma simii and A. vanbreuseghemii. The sequence of the alpha-box gene was 1,375 bp, containing 2 exons (from 172 to 463 bp and from 513 to 1,375 bp) in the A. simii (-) mating type strain and 1,380 bp, containing 2 exons (from 177 to 468 bp and from 518 to 1,380 bp) in the A. vanbreuseghemii (-) mating type strain. The sequence of the HMG gene was 1,871 bp, containing 2 exons (from 181 to 362 bp and from 426 to 1,440 bp, coding a protein of 398 amino acids) in the A. simii (+) mating type strain and 1,811 bp containing 2 exons (from 158 to 339 bp and from 403 to 1,381 bp, coding a protein of 386 amino acids) in the A. vanbreuseghemii (+) mating type strain. Of 15 animal isolates and 72 human isolates examined, the alpha-box gene was detected in five of the animal isolates and in none of the human isolates, while the HMG gene was detected in the other 10 of the animal isolates and in all of the human isolates. Phylogenetic analysis of the alpha-box and HMG genes of Trichophyton mentagrophytes complex strains and the Microsporum gypseum strain revealed that these strains were divided into 4 clusters; the first cluster consisting of A. vanbreuseghemii and the isolates from animals and humans, the second cluster consisting of A. simii, the third cluster consisting of A. benhamiae and the fourth cluster consisting of M. gypseum. These results indicate that anthropophilic T. mentagrophytes evolved from the A. vanbreuseghemii (+) mating strain.     

纳他霉素高产菌株Streptomyces gilvosporeus F607基因组及其生物合成基因簇分析

【背景】纳他霉素(Natamycin)是一种天然、广谱、高效的多烯大环内酯类抗真菌剂,褐黄孢链霉菌(Streptomyces gilvosporeus)是一种重要的纳他霉素产生菌。目前S. gilvosporeus基因组序列分析还未有报道,限制了该菌中纳他霉素及其他次级代谢产物合成及调控的研究。【目的】解析纳他霉素高产菌株S. gilvosporeus F607的基因组序列信息,挖掘其次级代谢产物基因资源,为深入研究该菌株的纳他霉素高产机理及生物合成调控机制奠定基础。【方法】利用相关软件对F607菌株的基因组序列进行基因预测、功能注释、进化分析和共线性分析,并预测次级代谢产物合成基因簇;对纳他霉素生物合成基因簇进行注释分析,比较分析不同菌种中纳他霉素生物合成基因簇的差异;分析预测S.gilvosporeusF607中纳他霉素生物合成途径。【结果】F607菌株基因组总长度为8482298bp,(G+C)mol%为70.95%,分别在COG、GO、KEGG数据库提取到5 062、4 428、5063个基因的注释信息。同时,antiSMASH软件预测得到29个次级代谢产物合成基因簇,其中纳他霉素基因簇与S.natalensis、S. chattanoogensis等菌株的纳他霉素基因簇相似性分别为81%和77%。除2个参与调控的sngT和sgnH基因和9个未知功能的orf基因有差异外,S. gilvosporeus F607基因簇中其他纳他霉素生物合成基因及其排列顺序与已知的纳他霉素基因簇高度一致。【结论】分析了S. gilvosporeus全基因组信息,预测了S. gilvosporeus F607中纳他霉素生物合成的途径,为从基因组层面上解析S. gilvosporeus F607菌株高产纳他霉素的内在原因提供了基础数据,为揭示纳他霉素高产的机理及工业化生产和未来新药的发现奠定了良好的基础。