Genome-Wide Comparative Analysis of Microsatellites in Pineapple

Pineapple (Ananas comosus (L.) Merrill) is the second most important tropical fruit in term of international trade. The availability of whole genomic sequences and expressed sequence tags (ESTs) offers an opportunity to identify and characterize microsatellite or simple sequence repeat (SSR) markers in pineapple. A total of 278,245 SSRs and 41,962 SSRs with an overall density of 728.57 SSRs/Mb and 619.37 SSRs/Mb were mined from genomic and ESTs sequences, respectively. 5′-untranslated regions (5′-UTRs) had the greatest amount of SSRs, 3.6–5.2 fold higher SSR density than other regions. For repeat length, 12 bp was the predominant repeat length in both assembled genome and ESTs. Class I SSRs were underrepresented compared with class II SSRs. For motif length, dinucleotide repeats were the most abundant in genomic sequences, whereas trinucleotides were the most common motif in ESTs. Tri- and hexanucleotides of total SSRs were more prevalent in ESTs than in the whole genome. The SSR frequency decreased dramatically as repeat times increased. AT was the most frequent single motif across the entire genome while AG was the most abundant motif in ESTs. Across six examined plant species, the pineapple genome displayed the highest density, substantially more than the second-place cucumber. Annotation and expression analyses were also conducted for genes containing SSRs. This thorough analysis of SSR markers in pineapple provided valuable information on the frequency and distribution of SSRs in the pineapple genome. This genomic resource will expedite genomic research and pineapple improvement.     

Genome-Wide Analysis of Repetitive Elements in Papaya

Papaya (Carica papaya L.) is an important fruit crop cultivated in tropical and subtropical regions worldwide. A first draft of its genome sequence has been recently released. Together with Arabidopsis, rice, poplar, grapevine and other genomes in the pipeline, it represents a good opportunity to gain insight into the organization of plant genomes. Here we report a detailed analysis of repetitive elements in the papaya genome, including transposable elements (TEs), tandemly-arrayed sequences, and high copy number genes. These repetitive sequences account for ～56% of the papaya genome with TEs being the most abundant at 52%, tandem repeats at 1.3% and high copy number genes at 3%. Most common types of TEs are represented in the papaya genome with retrotransposons being the dominant class, accounting for 40% of the genome. The most prevalent retrotransposons are Ty3-gypsy (27.8%) and Ty1-copia (5.5%). Among the tandem repeats, microsatellites are the most abundant in number, but represent only 0.19% of the genome. Minisatellites and satellites are less abundant, but represent 0.68% and 0.43% of the genome, respectively, due to greater repeat length. Despite an overall smaller gene repertoire in papaya than many other angiosperms, a significant fraction of genes (>2%) are present in large gene families with copy number greater than 20. This repeat database clarified a major part of the papaya genome organization and partly explained the lower gene repertoire in papaya than in Arabidopsis.     

转基因番木瓜的抗病性及分子鉴定

对T1代转番木瓜环斑病毒(PRV)复制酶(RP)突变体基因的两个番木瓜株系,进行了抗病性和分子生物学分析。结果表明,转基因番木瓜对PRV抗性达到高抗或免疫,目的基因RP遗传至转基因后代并在RNA水平表达,PCR可检测到CaMV35S启动子序列、标记基因NPTII。 Abstract:Virus resistance in field and molecular biological characterizations of the transgenes were analyzed for two lines of T1 generation of transgenic papaya with the replicase mutant gene from papaya ringspot virus (PRV).The transgenic plants showed highly resistant or immune against PRV.Results indicated that the transgenes inherited to and expressed at RNA level in the progenies.     

Genome-Wide Diet-Gene Interaction Analyses for Risk of Colorectal Cancer

Dietary factors, including meat, fruits, vegetables and fiber, are associated with colorectal cancer; however, there is limited information as to whether these dietary factors interact with genetic variants to modify risk of colorectal cancer. We tested interactions between these dietary factors and approximately 2.7 million genetic variants for colorectal cancer risk among 9,287 cases and 9,117 controls from ten studies. We used logistic regression to investigate multiplicative gene-diet interactions, as well as our recently developed Cocktail method that involves a screening step based on marginal associations and gene-diet correlations and a testing step for multiplicative interactions, while correcting for multiple testing using weighted hypothesis testing. Per quartile increment in the intake of red and processed meat were associated with statistically significant increased risks of colorectal cancer and vegetable, fruit and fiber intake with lower risks. From the case-control analysis, we detected a significant interaction between rs4143094 (10p14/near GATA3) and processed meat consumption (OR = 1.17; p = 8.7E-09), which was consistently observed across studies (p heterogeneity = 0.78). The risk of colorectal cancer associated with processed meat was increased among individuals with the rs4143094-TG and -TT genotypes (OR = 1.20 and OR = 1.39, respectively) and null among those with the GG genotype (OR = 1.03). Our results identify a novel gene-diet interaction with processed meat for colorectal cancer, highlighting that diet may modify the effect of genetic variants on disease risk, which may have important implications for prevention.     

Genome-Wide Analyses of Gene Expression during Mouse Endochondral Ossification

Background

Endochondral ossification is a complex process involving a series of events that are initiated by the establishment of a chondrogenic template and culminate in its replacement through the coordinated activity of osteoblasts, osteoclasts and endothelial cells. Comprehensive analyses of in vivo gene expression profiles during these processes are essential to obtain a complete understanding of the regulatory mechanisms involved.

Methodology/Principal Findings

To address these issues, we completed a microarray screen of three zones derived from manually segmented embryonic mouse tibiae. Classification of genes differentially expressed between each respective zone, functional categorization as well as characterization of gene expression patterns, cytogenetic loci, signaling pathways and functional motifs both confirmed reported data and provided novel insights into endochondral ossification. Parallel comparisons of the microdissected tibiae data set with our previously completed micromass culture screen further corroborated the suitability of micromass cultures for modeling gene expression in chondrocyte development. The micromass culture system demonstrated striking similarities to the in vivo microdissected tibiae screen; however, the micromass system was unable to accurately distinguish gene expression differences in the hypertrophic and mineralized zones of the tibia.

Conclusions/Significance

These studies allow us to better understand gene expression patterns in the growth plate and endochondral bones and provide an important technical resource for comparison of gene expression in diseased or experimentally-manipulated cartilages. Ultimately, this work will help to define the genomic context in which genes are expressed in long bones and to understand physiological and pathological ossification.     

Insight in Genome-Wide Association of Metabolite Quantitative Traits by Exome Sequence Analyses

Metabolite quantitative traits carry great promise for epidemiological studies, and their genetic background has been addressed using Genome-Wide Association Studies (GWAS). Thus far, the role of less common variants has not been exhaustively studied. Here, we set out a GWAS for metabolite quantitative traits in serum, followed by exome sequence analysis to zoom in on putative causal variants in the associated genes. ¹H Nuclear Magnetic Resonance (¹H-NMR) spectroscopy experiments yielded successful quantification of 42 unique metabolites in 2,482 individuals from The Erasmus Rucphen Family (ERF) study. Heritability of metabolites were estimated by SOLAR. GWAS was performed by linear mixed models, using HapMap imputations. Based on physical vicinity and pathway analyses, candidate genes were screened for coding region variation using exome sequence data. Heritability estimates for metabolites ranged between 10% and 52%. GWAS replicated three known loci in the metabolome wide significance: CPS1 with glycine (P-value  = 1.27×10⁻³²), PRODH with proline (P-value  = 1.11×10⁻¹⁹), SLC16A9 with carnitine level (P-value  = 4.81×10⁻¹⁴) and uncovered a novel association between DMGDH and dimethyl-glycine (P-value  = 1.65×10⁻¹⁹) level. In addition, we found three novel, suggestively significant loci: TNP1 with pyruvate (P-value  = 1.26×10⁻⁸), KCNJ16 with 3-hydroxybutyrate (P-value  = 1.65×10⁻⁸) and 2p12 locus with valine (P-value  = 3.49×10⁻⁸). Exome sequence analysis identified potentially causal coding and regulatory variants located in the genes CPS1, KCNJ2 and PRODH, and revealed allelic heterogeneity for CPS1 and PRODH. Combined GWAS and exome analyses of metabolites detected by high-resolution ¹H-NMR is a robust approach to uncover metabolite quantitative trait loci (mQTL), and the likely causative variants in these loci. It is anticipated that insight in the genetics of intermediate phenotypes will provide additional insight into the genetics of complex traits.     

HCV Genome-Wide Genetic Analyses in Context of Disease Progression and Hepatocellular Carcinoma

Hepatitis C virus (HCV) is a major cause of hepatitis and hepatocellular carcinoma (HCC) world-wide. Most HCV patients have relatively stable disease, but approximately 25% have progressive disease that often terminates in liver failure or HCC. HCV is highly variable genetically, with seven genotypes and multiple subtypes per genotype. This variation affects HCV’s sensitivity to antiviral therapy and has been implicated to contribute to differences in disease. We sequenced the complete viral coding capacity for 107 HCV genotype 1 isolates to determine whether genetic variation between independent HCV isolates is associated with the rate of disease progression or development of HCC. Consensus sequences were determined by sequencing RT-PCR products from serum or plasma. Positions of amino acid conservation, amino acid diversity patterns, selection pressures, and genome-wide patterns of amino acid covariance were assessed in context of the clinical phenotypes. A few positions were found where the amino acid distributions or degree of positive selection differed between in the HCC and cirrhotic sequences. All other assessments of viral genetic variation and HCC failed to yield significant associations. Sequences from patients with slow disease progression were under a greater degree of positive selection than sequences from rapid progressors, but all other analyses comparing HCV from rapid and slow disease progressors were statistically insignificant. The failure to observe distinct sequence differences associated with disease progression or HCC employing methods that previously revealed strong associations with the outcome of interferon α-based therapy implies that variable ability of HCV to modulate interferon responses is not a dominant cause for differential pathology among HCV patients. This lack of significant associations also implies that host and/or environmental factors are the major causes of differential disease presentation in HCV patients.     

Genome-Wide Analyses Suggest Mechanisms Involving Early B-Cell Development in Canine IgA Deficiency

Immunoglobulin A deficiency (IgAD) is the most common primary immune deficiency disorder in both humans and dogs, characterized by recurrent mucosal tract infections and a predisposition for allergic and other immune mediated diseases. In several dog breeds, low IgA levels have been observed at a high frequency and with a clinical resemblance to human IgAD. In this study, we used genome-wide association studies (GWAS) to identify genomic regions associated with low IgA levels in dogs as a comparative model for human IgAD. We used a novel percentile groups-approach to establish breed-specific cut-offs and to perform analyses in a close to continuous manner. GWAS performed in four breeds prone to low IgA levels (German shepherd, Golden retriever, Labrador retriever and Shar-Pei) identified 35 genomic loci suggestively associated (p <0.0005) to IgA levels. In German shepherd, three genomic regions (candidate genes include KIRREL3 and SERPINA9) were genome-wide significantly associated (p <0.0002) with IgA levels. A ~20kb long haplotype on CFA28, significantly associated (p = 0.0005) to IgA levels in Shar-Pei, was positioned within the first intron of the gene SLIT1. Both KIRREL3 and SLIT1 are highly expressed in the central nervous system and in bone marrow and are potentially important during B-cell development. SERPINA9 expression is restricted to B-cells and peaks at the time-point when B-cells proliferate into antibody-producing plasma cells. The suggestively associated regions were enriched for genes in Gene Ontology gene sets involving inflammation and early immune cell development.     

利用全测序基因组数据对被子植物纤维素合成酶基因的进化分析

利用已经分离的植物纤维素合成酶基因的Cellulose_synt结构域为检索序列,从NCBI和其他数据库中调取已完成测序的物种的纤维素合成酶的氨基酸序列,共涉及10个物种的171个基因,基于以上氨基酸序列,应用MEGA4.0生成系统进化树。结果表明:CesA基因和Csl基因直向的相似度远大于平行的相似度,且它们的分化可能在单子叶和真双子叶植物分化之前,单子叶和真双子叶植物的最近共同祖先中至少存在7个CesA基因,综合已知的模式植物CesA基因的功能(初生壁或次生壁形成特异性),可为推测其他物种中该基因的功能提供帮助。     

Genome-Wide Identification and Comprehensive Analyses of the Kinomes in Four Pathogenic Microsporidia Species

Microsporidia have attracted considerable attention because they infect a wide range of hosts, from invertebrates to vertebrates, and cause serious human diseases and major economic losses in the livestock industry. There are no prospective drugs to counteract this pathogen. Eukaryotic protein kinases (ePKs) play a central role in regulating many essential cellular processes and are therefore potential drug targets. In this study, a comprehensive summary and comparative analysis of the protein kinases in four microsporidia–Enterocytozoon bieneusi, Encephalitozoon cuniculi, Nosema bombycis and Nosema ceranae–was performed. The results show that there are 34 ePKs and 4 atypical protein kinases (aPKs) in E. bieneusi, 29 ePKs and 6 aPKs in E. cuniculi, 41 ePKs and 5 aPKs in N. bombycis, and 27 ePKs and 4 aPKs in N. ceranae. These data support the previous conclusion that the microsporidian kinome is the smallest eukaryotic kinome. Microsporidian kinomes contain only serine-threonine kinases and do not contain receptor-like and tyrosine kinases. Many of the kinases related to nutrient and energy signaling and the stress response have been lost in microsporidian kinomes. However, cell cycle-, development- and growth-related kinases, which are important to parasites, are well conserved. This reduction of the microsporidian kinome is in good agreement with genome compaction, but kinome density is negatively correlated with proteome size. Furthermore, the protein kinases in each microsporidian genome are under strong purifying selection pressure. No remarkable differences in kinase family classification, domain features, gain and/or loss, and selective pressure were observed in these four species. Although microsporidia adapt to different host types, the coevolution of microsporidia and their hosts was not clearly reflected in the protein kinases. Overall, this study enriches and updates the microsporidian protein kinase database and may provide valuable information and candidate targets for the design of treatments for pathogenic diseases.     

Genome-Wide Analyses of Recombination Suggest That Giardia intestinalis Assemblages Represent Different Species

Giardia intestinalis is a major cause of waterborne enteric disease in humans. The species is divided into eight assemblages suggested to represent separate Giardia species based on host specificities and the genetic divergence of marker genes. We have investigated whether genome-wide recombination occurs between assemblages using the three available G. intestinalis genomes. First, the relative nonsynonymous substitution rates of the homologs were compared for 4,009 positional homologs. The vast majority of these comparisons indicate genetic isolation without interassemblage recombinations. Only a region of 6 kbp suggests genetic exchange between assemblages A and E, followed by gene conversion events. Second, recombination-detecting software fails to identify within-gene recombination between the different assemblages for most of the homologs. Our results indicate very low frequency of recombination between the syntenic core genes, suggesting that G. intestinalis assemblages are genetically isolated lineages and thus should be viewed as separated Giardia species.     

有尾噬菌体目中微卫星和复合型微卫星的发生及分析

简单重复序列亦称微卫星,被成功应用于许多真核生物、原核生物和病毒的基因组和进化研究,但是噬菌体中的微卫星目前很少被研究。因此对60条尾病毒目基因组中的微卫星和和复合型微卫星(由两个或两个以上直接相邻的微卫星组成)做综合性分析,在这60个基因组中总共观察到11 874个微卫星和449个复合型微卫星。相关性分析表明微卫星个数与基因组大小成正线性相关(ρ=0.899, P<0.01)。参考序列中的微卫星个数少于对应的随机序列中微卫星个数,这种反常现象主要是因为参考序列含有较少的单核苷酸和二核苷酸重复。A/T和AT/TA重复是单核苷酸和二核苷酸重复中最主要的类型,因此单核苷酸重复中的GC含量明显低于相应的序列中的GC含量;相比之下,微卫星中的二核苷酸和三核苷酸重复的GC含量与对应的参考序列的GC含量无明显区别。尾病毒目基因组中的这些结果与其它生物体基因组存在一定的差别。有助于了解尾病毒目中微卫星的分布、进化和生物学功能。     

Genome-Wide Comparative Analysis of pogo-Like Transposable Elements in Different Fusarium Species

The recent availability of genome sequences of four different Fusarium species offers the opportunity to perform extensive comparative analyses, in particular of repeated sequences. In a recent work, the overall content of such sequences in the genomes of three phylogenetically related Fusarium species, F. graminearum, F. verticillioides, and F. oxysporum f. sp. lycopersici has been estimated. In this study, we present an exhaustive characterization of pogo-like elements, named Fots, in four Fusarium genomes. Overall 10 Fot and two Fot-related miniature inverted-repeat transposable element families were identified, revealing a diversification of multiple lineages of pogo-like elements, some of which accompanied by a gain of introns. This analysis also showed that such elements are present in an unusual high proportion in the genomes of F. oxysporum f. sp. lycopersici and Nectria haematococca (anamorph F. solani f. sp. pisi) in contrast with most other fungal genomes in which retroelements are the most represented. Interestingly, our analysis showed that the most numerous Fot families all contain potentially active or mobilisable copies, thus conferring a mutagenic potential of these transposable elements and consequently a role in strain adaptation and genome evolution. This role is strongly reinforced when examining their genomic distribution which is clearly biased with a high proportion (more than 80%) located on strain- or species-specific regions enriched in genes involved in pathogenicity and/or adaptation. Finally, the different reproductive characteristics of the four Fusarium species allowed us to investigate the impact of the process of repeat-induced point mutations on the expansion and diversification of Fot elements.     

Comparative Genome-Wide Analysis of the Malate Dehydrogenase Gene Families in Cotton

Malate dehydrogenases (MDHs) play crucial roles in the physiological processes of plant growth and development. In this study, 13 and 25 MDH genes were identified from Gossypium raimondii and Gossypium hirsutum, respectively. Using these and 13 previously reported Gossypium arboretum MDH genes, a comparative molecular analysis between identified MDH genes from G. raimondii, G. hirsutum, and G. arboretum was performed. Based on multiple sequence alignments, cotton MDHs were divided into five subgroups: mitochondrial MDH, peroxisomal MDH, plastidial MDH, chloroplastic MDH and cytoplasmic MDH. Almost all of the MDHs within the same subgroup shared similar gene structure, amino acid sequence, and conserved motifs in their functional domains. An analysis of chromosomal localization suggested that segmental duplication played a major role in the expansion of cotton MDH gene families. Additionally, a selective pressure analysis indicated that purifying selection acted as a vital force in the evolution of MDH gene families in cotton. Meanwhile, an expression analysis showed the distinct expression profiles of GhMDHs in different vegetative tissues and at different fiber developmental stages, suggesting the functional diversification of these genes in cotton growth and fiber development. Finally, a promoter analysis indicated redundant but typical cis-regulatory elements for the potential functions and stress activity of many MDH genes. This study provides fundamental information for a better understanding of cotton MDH gene families and aids in functional analyses of the MDH genes in cotton fiber development.     

Comparative Corticotype Analyses in Tetrahymena*

SYNOPSIS. Strains of Tetrahymena pyriformis, including amicronucleate strain GL and representatives of 9 syngens, have been examined to determine the patterns whereby cortical features vary with numbers of ciliary meridians. The characteristics scored were the positions of the contractile vacuole pores (CVP's), the extent of the area within which CVP's develop, the incidence of supernumerary CVP's, and the number of postoral meridians. Intrasyngenic comparisons were possible in 6 syngens and permitted an assessment of intrasyngenic variation for these characteristics. Only the CVP positions appear to be reasonably constant within syngens in the strains examined. On the basis of this criterion the syngens can be arrayed in an approximate order of 1, 3, 7, 6, 8, 9, 2, 5, GL and 4; the angle formed between the central axis, the stomatogenic meridian and the CVP's is most acute in syngen 1.