稻属C染色体组物种的群体遗传学和物种形成研究

植物近缘物种系统发育和物种形成过程一直以来都是植物进化生物学研究中最基本的问题之一，是人们理解自然界物种多样性产生和变化的前提。近缘物种间通常形态相似，遗传和分子水平的分化很小且常受诸如渐渗杂交、谱系分选、种内重组等微观进化事件的影响，导致植物近缘物种间系统发育和物种形成过程研究极为困难。在过去十多年中，生物技术的革新和理论方法的发展极大地推动了进化生物学研究，促进了人们对一些重要模式生物与其近缘种间的系统发育关系和物种分化过程的认识，如人、果蝇、线虫、拟南芥和玉米等。然而，迄今在植物中，许多重要类群及其近缘种的系统发育和物种形成过程研究仍不多见，包括在具有特殊重要性的栽培作物中，如稻属（Oryza L.）。由于属内包含有重要粮食作物水稻，稻属向来都是禾本科内备受关注的一个类群。本研究中，我们通过多基因序列的方法，探讨了稻属C染色体组三个近缘二倍体物种的系统发育和物种形成过程，主要研究结果如下。 为选择适宜的实验策略和保证序列数据的真实性，我们利用不同聚合酶扩增自交的栽培稻Oryza sativa ssp. japonica和异交的O. longistaminata不同类型的基因片段，采用克隆测序的方法，评估聚合酶链式扩增反应（PCR）中产生的非真实变异的状况。我们使用exTaq、exTaq和Pfu混和酶和PfuUltraTM酶三种不同聚合酶扩增了Adh1、GPA1和Waxy三个基因片段。在检测到的非真实变异中，PCR 错误的类型主要为单碱基变异和不同等位基因间重组，其中以单碱基变异为主，且突变类型以转换占绝大多数。比较不同酶扩增错误的结果表明，高保真PfuUltraTM酶对PCR反应错误有显著的改善，在扩增产物的单克隆中几乎检测不到PCR噪音，错误率仅为0.0001%，而exTaq和混和酶的错误率分别为0.096%和0.073%。从不同物种比较结果来看，exTaq酶和Pfu酶混用时在自交的O. sativa ssp. japonica内对PCR错误也表现明显的改善，但在异交的O. longistaminata中改善效应不太明显。在三个不同基因位点上，PCR扩增错误出现频率随扩增区域增长而变大。在PfuUltraTM酶的扩增产物中发现重组最少，exTaq和混和酶重组较多，且混和酶对重组改善效果不明显。基于不同聚合酶扩增错误对比研究结果，我们认为，由于能保证序列变异的真实性且不遗漏等位基因，核基因的克隆测序较宜于分离杂合个体中不同的等位基因。 基于随机挑取4个叶绿体和10个核基因位点，利用12份Oryza officinalis、8份O. eichingeri和4份O. rhizomatis材料，对稻属C染色体组三个近缘二倍体物种的系统发育关系作了深入分析。利用不同的系统发育分析方法对单基因位点序列和合并序列数据作了分析，结果表明，稻属C染色体组三个近缘种间呈多歧分支。因此，三个二倍体C染色体组物种可能经快速辐射分化形成。在不同核基因和叶绿体基因的系统发育树中各分支枝长均很短，也表明C染色体组的三个物种可能在较短时期内分化出。不同基因间拓扑结构不一致主要受谱系随机分选的影响。此外，C染色体组不同物种间的种间渐渗和不同等位基因重组对系统发育树的冲突也有影响。值得注意的是，C染色体组三个物种的辐射分化不排除由于相邻两次物种形成事件间隔时间太近、目前数据量不够而无法分辨的可能。在本研究中，我们发现，对于系统发育重建困难的类群，利用等位基因构建物种谱系树有助于挖掘不同基因间结果不一致的因素。 基于10个随机选取的核基因序列数据，利用物种水平的取样方式和群体遗传学分析方法，我们研究了稻属C染色体组三个近缘物种O. officinalis、O. eichingeri和O. rhizomatis的核苷酸多态性，并根据多态性水平和式样，推测了三个近缘种分化的历史。结果表明，在C染色体组的三个近缘种中，仅分布于斯里兰卡的O. rhizomatis的核苷酸多态性水平相对最低（θsil = 0.0038），而间断分布于非洲和斯里兰卡的O. eichingeri最高（θsil = 0.0057）。与被子植物其他类群相比，稻属C染色体组三个物种的核苷酸多态性水平显得较低，O. eichingeri的核苷酸多态性仅约为玉米及其近缘野生种的23-46%和拟南芥的35%。C染色体组内三个野生种相对较低的核苷酸多态性水平可能起因于其较小的祖先有效群体。物种形成模型分析表明，O. officinalis和其近缘种从最近共同祖先分化开后可能经历了居群缩减的历史，且自最近共同祖先分开后，三个物种彼此间并无显著的基因交流。基于分子钟粗略估算了C染色体三个物种分化时间，结果表明，三个物种彼此在很短的时期内分开，约0.63-0.68 Myr。同时，O. eichingeri在非洲和斯里兰卡两个地理宗的分异时间约为0.37 Myr，且推测斯里兰卡的O. eichingeri主要由西非经长距离扩散传播到斯里兰卡。

Complete genome sequence of lymphocystis disease virus isolated from China

Lymphocystis diseases in fish throughout the world have been extensively described. Here we report the complete genome sequence of lymphocystis disease virus isolated in China (LCDV-C), an LCDV isolated from cultured flounder (Paralichthys olivaceus) with lymphocystis disease in China. The LCDV-C genome is 186,250 bp, with a base composition of 27.25% G+C. Computer-assisted analysis revealed 240 potential open reading frames (ORFs) and 176 nonoverlapping putative viral genes, which encode polypeptides ranging from 40 to 1,193 amino acids. The percent coding density is 67%, and the average length of each ORF is 702 bp. A search of the GenBank database using the 176 individual putative genes revealed 103 homologues to the corresponding ORFs of LCDV-1 and 73 potential genes that were not found in LCDV-1 and other iridoviruses. Among the 73 genes, there are 8 genes that contain conserved domains of cellular genes and 65 novel genes that do not show any significant homology with the sequences in public databases. Although a certain extent of similarity between putative gene products of LCDV-C and corresponding proteins of LCDV-1 was revealed, no colinearity was detected when their ORF arrangements and coding strategies were compared to each other, suggesting that a high degree of genetic rearrangements between them has occurred. And a large number of tandem and overlapping repeated sequences were observed in the LCDV-C genome. The deduced amino acid sequence of the major capsid protein (MCP) presents the highest identity to those of LCDV-1 and other iridoviruses among the LCDV-C gene products. Furthermore, a phylogenetic tree was constructed based on the multiple alignments of nine MCP amino acid sequences. Interestingly, LCDV-C and LCDV-1 were clustered together, but their amino acid identity is much less than that in other clusters. The unexpected levels of divergence between their genomes in size, gene organization, and gene product identity suggest that LCDV-C and LCDV-1 shouldn't belong to a same species and that LCDV-C should be considered a species different from LCDV-1.

The downy mildew resistance locus Pp523 is located on chromosome C8 of Brassica oleracea L.

We have previously constructed a genetic map of Brassica oleracea L. containing the Pp523 locus that confers downy mildew resistance to adult plants. In this work, 44 SSR markers of reference for the Brassica C genome chromosomes were added to the map, allowing the nine major linkage groups to be assigned to the nine chromosomes of B. oleracea. Locus Pp523 was located on chromosome C8, and a locus determining flower colour was mapped to chromosome C3. In comparison with the first version of the map, the new map is denser and more compact. The available genomic information on B. oleracea was enriched with the chromosome location of two phenotypic traits and 421 DNA markers (RAPD, ISSR, AFLP, SCAR, BAC-end derived STS, SSR and other PCR markers). Conversely, the genomic information on B. oleracea chromosome C8 is being used as an additional tool for the map-based cloning of Pp523, the first gene for adult plant resistance to downy mildew precisely located to a specific chromosome of this crop species.

Phylogenetic insights of Avena genus based on genomic analysis of 45S rDNA molecular organization

Tese de mestrado. Biologia (Biologia Molecular e Genética). Universidade de Lisboa, Faculdade de Ciências, 2014

Simple sequence repeats reveal uneven distribution of genetic diversity in chloroplast genomes of Brassica oleracea L. and (n=9) wild relatives

Diversity in the chloroplast genome of 171 accessions representing the Brassica 'C' (n = 9) genome, including domesticated and wild B. oleracea and nine inter-fertile related wild species, was investigated using six chloroplast SSR (microsatellite) markers. The lack of diversity detected among 105 cultivated and wild accessions of B. oleracea contrasted starkly with that found within its wild relatives. The vast majority of B. oleracea accessions shared a single haplotype, whereas as many as six haplotypes were detected in two wild species, B. villosa Biv. and B. cretica Lam.. The SSRs proved to be highly polymorphic across haplotypes, with calculated genetic diversity values (H) of 0.23-0.87. In total, 23 different haplotypes were detected in C genome species, with an additional five haplotypes detected in B. rapa L. (A genome n = 10) and another in B. nigra L. (B genome, n = 8). The low chloroplast diversity of B. oleracea is not suggestive of multiple domestication events. The predominant B. oleracea haplotype was also common in B. incana Ten. and present in low frequencies in B. villosa, B. macrocarpa Guss, B. rupestris Raf. and B. cretica. The chloroplast SSRs reveal a wealth of diversity within wild Brassica species that will facilitate further evolutionary and phylogeographic studies of this important crop genus.

Shoot calcium and magnesium concentrations differ between subtaxa, are highly heritable, and associate with potentially pleiotropic loci in Brassica oleracea

Calcium (Ca) and magnesium (Mg) are the most abundant group II elements in both plants and animals. Genetic variation in shoot Ca and shoot Mg concentration (shoot Ca and Mg) in plants can be exploited to biofortify food crops and thereby increase dietary Ca and Mg intake for humans and livestock. We present a comprehensive analysis of within-species genetic variation for shoot Ca and Mg, demonstrating that shoot mineral concentration differs significantly between subtaxa (varietas). We established a structured diversity foundation set of 376 accessions to capture a high proportion of species-wide allelic diversity within domesticated Brassica oleracea, including representation of wild relatives (C genome, 1n = 9) from natural populations. These accessions and 74 modern F-1 hybrid cultivars were grown in glasshouse and field environments. Shoot Ca and Mg varied 2- and 2.3-fold, respectively, and was typically not inversely correlated with shoot biomass, within most subtaxa. The closely related capitata (cabbage) and sabauda (Savoy cabbage) subtaxa consistently had the highest mean shoot Ca and Mg. Shoot Ca and Mg in glasshouse-grown plants was highly correlated with data from the field. To understand and dissect the genetic basis of variation in shoot Ca and Mg, we studied homozygous lines from a segregating B. oleracea mapping population. Shoot Ca and Mg was highly heritable (up to 40). Quantitative trait loci (QTL) for shoot Ca and Mg were detected on chromosomes C2, C6, C7, C8, and, in particular, C9, where QTL accounted for 14 to 55 of the total genetic variance. The presence of QTL on C9 was substantiated by scoring recurrent backcross substitution lines, derived from the same parents. This also greatly increased the map resolution, with strong evidence that a 4-cM region on C9 influences shoot Ca. This region corresponds to a 0.41-Mb region on Arabidopsis (Arabidopsis thaliana) chromosome 5 that includes 106 genes. There is also evidence that pleiotropic loci on C8 and C9 affect shoot Ca and Mg. Map-based cloning of these loci will reveal how shoot-level phenotypes relate to Ca 21 and Mg 21 uptake and homeostasis at the molecular level.

Lack of evidence to support the association of a single IL28B genotype SNP rs12979860 with the HTLV-1 clinical outcomes and proviral load

Background: The Interleukin 28B (IL28B) rs12979860 polymorphisms was recently reported to be associated with the human T-cell leukemia virus type 1 (HTLV-1) proviral load (PvL) and the development of the HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP). Methods: In an attempt to examine this hypothesis, we assessed the association of the rs12979860 genotypes with HTLV-1 PvL levels and clinical status in 112 unrelated Brazilian subjects (81 HTLV-1 asymptomatic carriers, 24 individuals with HAM/TSP and 7 with Adult T cell Leukemia/Lymphoma (ATLL)). Results: All 112 samples were successfully genotyped and their PvLs compared. Neither the homozygote TT nor the heterozygote CT mutations nor the combination genotypes (TT/CT) were associated with a greater PvL. We also observed no significant difference in allele distribution between asymptomatic carriers and patients with HTLV-1 associated HAM/TSP. Conclusions: Our study failed to support the previously reported positive association between the IL28B rs12979860 polymorphisms and an increased risk of developing HAM/TSP in the Brazilian population.

Discrimination of the closely related A and D genomes of the hexaploid oat Avena sativa L.

A satellite DNA sequence, As120a, specific to the A-genome chromosomes in the hexaploid oat, Avena sativa L., was isolated by subcloning a fragment with internal tandem repeats from a plasmid, pAs120, that had been obtained from an Avena strigosa (As genome) genomic library. Southern and in situ hybridization showed that sequences with homology to sequences within pAs120 were dispersed throughout the genome of diploid (A and C genomes), tetraploid (AC genomes), and hexaploid (ACD genomes) Avena species. In contrast, sequences homologous to As120a were found in two A-genome species (A. strigosa and Avena longiglumis) and in the hexaploid A. sativa whereas this sequence was little amplified in the tetraploid Avena murphyi and was absent in the remaining A- and C-genome diploid species. In situ hybridization of pAs120a to hexaploid oat species revealed the distribution of elements of the As120a repeated family over both arms of 14 of 42 chromosomes of this species. By using double in situ hybridization with pAs120a and a C genome-specific probe, three sets of 14 chromosomes were revealed corresponding to the A, C, and D genomes of the hexaploid species. Simultaneous in situ hybridizations with pAs120a and ribosomal probes were used to assign the SAT chromosomes of hexaploid species to their correct genomes. This work reports a sequence able to distinguish between the closely related A and D genomes of hexaploid oats. This sequence offers new opportunities to analyze the relationships of Avena species and to explore the possible evolution of various polyploid oat species.

A genome-wide asociation study identifies new psoriasis susceptibility loci and an interaction betwEn HLA-C and ERAP1

To identify new susceptibility loci for psoriasis, we undertOk a genome-wide asociation study of 594,224 SNPs in 2,622 individuals with psoriasis and 5,667 controls. We identified asociations at eight previously unreported genomic loci. Seven loci harbored genes with recognized iMune functions (IL28RA, REL, IFIH1, ERAP1, TRAF3IP2, NFKBIA and TYK2). These asociations were replicated in 9,079 European samples (six loci with a combined P < 5-10 -8 and two loci with a combined P < 5-10-7). We also report compeLing evidence for an interaction betwEn the HLA-C and ERAP1 loci (combined P = 6.95-10-6). ERAP1 plays an important role in MHC claS I peptide proceSing. ERAP1 variants only influenced psoriasis susceptibility in individuals carrying the HLA-C risk aLele. Our findings implicate pathways that integrate epidermal barrier dysfunction with iNate and adaptive iMune dysregulation in psoriasis pathogenesis.

Elucidation of hepatitis C virus transmission and early diversification by single genome sequencing.

A precise molecular identification of transmitted hepatitis C virus (HCV) genomes could illuminate key aspects of transmission biology, immunopathogenesis and natural history. We used single genome sequencing of 2,922 half or quarter genomes from plasma viral RNA to identify transmitted/founder (T/F) viruses in 17 subjects with acute community-acquired HCV infection. Sequences from 13 of 17 acute subjects, but none of 14 chronic controls, exhibited one or more discrete low diversity viral lineages. Sequences within each lineage generally revealed a star-like phylogeny of mutations that coalesced to unambiguous T/F viral genomes. Numbers of transmitted viruses leading to productive clinical infection were estimated to range from 1 to 37 or more (median = 4). Four acutely infected subjects showed a distinctly different pattern of virus diversity that deviated from a star-like phylogeny. In these cases, empirical analysis and mathematical modeling suggested high multiplicity virus transmission from individuals who themselves were acutely infected or had experienced a virus population bottleneck due to antiviral drug therapy. These results provide new quantitative and qualitative insights into HCV transmission, revealing for the first time virus-host interactions that successful vaccines or treatment interventions will need to overcome. Our findings further suggest a novel experimental strategy for identifying full-length T/F genomes for proteome-wide analyses of HCV biology and adaptation to antiviral drug or immune pressures.

Genome-wide Association Study Implicates HLA-C*01:02 as a Risk Factor at the Major Histocompatibility Complex Locus in Schizophrenia

BACKGROUND: We performed a genome-wide association study (GWAS) to identify common risk variants for schizophrenia. METHODS: The discovery scan included 1606 patients and 1794 controls from Ireland, using 6,212,339 directly genotyped or imputed single nucleotide polymorphisms (SNPs). A subset of this sample (270 cases and 860 controls) was subsequently included in the Psychiatric GWAS Consortium-schizophrenia GWAS meta-analysis. RESULTS: One hundred eight SNPs were taken forward for replication in an independent sample of 13,195 cases and 31,021 control subjects. The most significant associations in discovery, corrected for genomic inflation, were (rs204999, p combined = 1.34 × 10(-9) and in combined samples (rs2523722 p combined = 2.88 × 10(-16)) mapped to the major histocompatibility complex (MHC) region. We imputed classical human leukocyte antigen (HLA) alleles at the locus; the most significant finding was with HLA-C*01:02. This association was distinct from the top SNP signal. The HLA alleles DRB1*03:01 and B*08:01 were protective, replicating a previous study. CONCLUSIONS: This study provides further support for involvement of MHC class I molecules in schizophrenia. We found evidence of association with previously reported risk alleles at the TCF4, VRK2, and ZNF804A loci.

Genetic variation in IL28B is associated with chronic hepatitis C and treatment failure: a genome-wide association study.

BACKGROUND & AIMS: Hepatitis C virus (HCV) induces chronic infection in 50% to 80% of infected persons; approximately 50% of these do not respond to therapy. We performed a genome-wide association study to screen for host genetic determinants of HCV persistence and response to therapy. METHODS: The analysis included 1362 individuals: 1015 with chronic hepatitis C and 347 who spontaneously cleared the virus (448 were coinfected with human immunodeficiency virus [HIV]). Responses to pegylated interferon alfa and ribavirin were assessed in 465 individuals. Associations between more than 500,000 single nucleotide polymorphisms (SNPs) and outcomes were assessed by multivariate logistic regression. RESULTS: Chronic hepatitis C was associated with SNPs in the IL28B locus, which encodes the antiviral cytokine interferon lambda. The rs8099917 minor allele was associated with progression to chronic HCV infection (odds ratio [OR], 2.31; 95% confidence interval [CI], 1.74-3.06; P = 6.07 x 10(-9)). The association was observed in HCV mono-infected (OR, 2.49; 95% CI, 1.64-3.79; P = 1.96 x 10(-5)) and HCV/HIV coinfected individuals (OR, 2.16; 95% CI, 1.47-3.18; P = 8.24 x 10(-5)). rs8099917 was also associated with failure to respond to therapy (OR, 5.19; 95% CI, 2.90-9.30; P = 3.11 x 10(-8)), with the strongest effects in patients with HCV genotype 1 or 4. This risk allele was identified in 24% of individuals with spontaneous HCV clearance, 32% of chronically infected patients who responded to therapy, and 58% who did not respond (P = 3.2 x 10(-10)). Resequencing of IL28B identified distinct haplotypes that were associated with the clinical phenotype. CONCLUSIONS: The association of the IL28B locus with natural and treatment-associated control of HCV indicates the importance of innate immunity and interferon lambda in the pathogenesis of HCV infection.

Enhanced hepatitis C virus genome replication and lipid accumulation mediated by inhibition of AMP-activated protein kinase

Hepatitis C virus (HCV) infection is associated with dysregulation of both lipid and glucose metabolism. As well as contributing to viral replication, these perturbations influence the pathogenesis associated with the virus, including steatosis, insulin resistance, and type 2 diabetes. AMP-activated protein kinase (AMPK) plays a key role in regulation of both lipid and glucose metabolism. We show here that, in cells either infected with HCV or harboring an HCV subgenomic replicon, phosphorylation of AMPK at threonine 172 and concomitant AMPK activity are dramatically reduced. We demonstrate that this effect is mediated by activation of the serine/threonine kinase, protein kinase B, which inhibits AMPK by phosphorylating serine 485. The physiological significance of this inhibition is demonstrated by the observation that pharmacological restoration of AMPK activity not only abrogates the lipid accumulation observed in virus-infected and subgenomic replicon-harboring cells but also efficiently inhibits viral replication. These data demonstrate that inhibition of AMPK is required for HCV replication and that the restoration of AMPK activity may present a target for much needed anti-HCV therapies.

A Novel Diagnostic Target in the Hepatitis C Virus Genome

BackgroundDetection and quantification of hepatitis C virus (HCV) RNA is integral to diagnostic and therapeutic regimens. All molecular assays target the viral 5'-noncoding region (59-NCR), and all show genotype-dependent variation of sensitivities and viral load results. Non-western HCV genotypes have been under-represented in evaluation studies. An alternative diagnostic target region within the HCV genome could facilitate a new generation of assays.Methods and FindingsIn this study we determined by de novo sequencing that the 3'-X-tail element, characterized significantly later than the rest of the genome, is highly conserved across genotypes. To prove its clinical utility as a molecular diagnostic target, a prototype qualitative and quantitative test was developed and evaluated multicentrically on a large and complete panel of 725 clinical plasma samples, covering HCV genotypes 1-6, from four continents (Germany, UK, Brazil, South Africa, Singapore). To our knowledge, this is the most diversified and comprehensive panel of clinical and genotype specimens used in HCV nucleic acid testing (NAT) validation to date. The lower limit of detection (LOD) was 18.4 IU/ml (95% confidence interval, 15.3-24.1 IU/ml), suggesting applicability in donor blood screening. The upper LOD exceeded 10(-9) IU/ml, facilitating viral load monitoring within a wide dynamic range. In 598 genotyped samples, quantified by Bayer VERSANT 3.0 branched DNA (bDNA), X-tail-based viral loads were highly concordant with bDNA for all genotypes. Correlation coefficients between bDNA and X-tail NAT, for genotypes 1-6, were: 0.92, 0.85, 0.95, 0.91, 0.95, and 0.96, respectively; X-tail-based viral loads deviated by more than 0.5 log10 from 5'-NCR-based viral loads in only 12% of samples (maximum deviation, 0.85 log10). The successful introduction of X-tail NAT in a Brazilian laboratory confirmed the practical stability and robustness of the X-tail-based protocol. The assay was implemented at low reaction costs (US$8.70 per sample), short turnover times (2.5 h for up to 96 samples), and without technical difficulties.ConclusionThis study indicates a way to fundamentally improve HCV viral load monitoring and infection screening. Our prototype assay can serve as a template for a new generation of viral load assays. Additionally, to our knowledge this study provides the first open protocol to permit industry-grade HCV detection and quantification in resource-limited settings.

Evaluation of canonical siRNA and dicer substrate RNA for inhibition of hepatitis C virus genome replication: a comparative study

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)