A transcriptome resource for the koala (Phascolarctos cinereus) : insights into koala retrovirus transcription and sequence diversity

Background The koala, Phascolarctos cinereus, is a biologically unique and evolutionarily distinct Australian arboreal marsupial. The goal of this study was to sequence the transcriptome from several tissues of two geographically separate koalas, and to create the first comprehensive catalog of annotated transcripts for this species, enabling detailed analysis of the unique attributes of this threatened native marsupial, including infection by the koala retrovirus. Results RNA-Seq data was generated from a range of tissues from one male and one female koala and assembled de novo into transcripts using Velvet-Oases. Transcript abundance in each tissue was estimated. Transcripts were searched for likely protein-coding regions and a non-redundant set of 117,563 putative protein sequences was produced. In similarity searches there were 84,907 (72%) sequences that aligned to at least one sequence in the NCBI nr protein database. The best alignments were to sequences from other marsupials. After applying a reciprocal best hit requirement of koala sequences to those from tammar wallaby, Tasmanian devil and the gray short-tailed opossum, we estimate that our transcriptome dataset represents approximately 15,000 koala genes. The marsupial alignment information was used to look for potential gene duplications and we report evidence for copy number expansion of the alpha amylase gene, and of an aldehyde reductase gene. Koala retrovirus (KoRV) transcripts were detected in the transcriptomes. These were analysed in detail and the structure of the spliced envelope gene transcript was determined. There was appreciable sequence diversity within KoRV, with 233 sites in the KoRV genome showing small insertions/deletions or single nucleotide polymorphisms. Both koalas had sequences from the KoRV-A subtype, but the male koala transcriptome has, in addition, sequences more closely related to the KoRV-B subtype. This is the first report of a KoRV-B-like sequence in a wild population. Conclusions This transcriptomic dataset is a useful resource for molecular genetic studies of the koala, for evolutionary genetic studies of marsupials, for validation and annotation of the koala genome sequence, and for investigation of koala retrovirus. Annotated transcripts can be browsed and queried at http://koalagenome.org

Sequence diversity and differential expression of major phenylpropanoid-flavonoid biosynthetic genes among three mango varieties

Background: Mango fruits contain a broad spectrum of phenolic compounds which impart potential health benefits; their biosynthesis is catalysed by enzymes in the phenylpropanoid-flavonoid (PF) pathway. The aim of this study was to reveal the variability in genes involved in the PF pathway in three different mango varieties Mangifera indica L., a member of the family Anacardiaceae: Kensington Pride (KP), Irwin (IW) and Nam Doc Mai (NDM) and to determine associations with gene expression and mango flavonoid profiles. Results: A close evolutionary relationship between mango genes and those from the woody species poplar of the Salicaceae family (Populus trichocarpa) and grape of the Vitaceae family (Vitis vinifera), was revealed through phylogenetic analysis of PF pathway genes. We discovered 145 SNPs in total within coding sequences with an average frequency of one SNP every 316bp. Variety IW had the highest SNP frequency (one SNP every 258bp) while KP and NDM had similar frequencies (one SNP every 369bp and 360bp, respectively). The position in the PF pathway appeared to influence the extent of genetic diversity of the encoded enzymes. The entry point enzymes phenylalanine lyase (PAL), cinnamate 4-mono-oxygenase (C4H) and chalcone synthase (CHS) had low levels of SNP diversity in their coding sequences, whereas anthocyanidin reductase (ANR) showed the highest SNP frequency followed by flavonoid 3'-hydroxylase (F3'H). Quantitative PCR revealed characteristic patterns of gene expression that differed between mango peel and flesh, and between varieties. Conclusions: The combination of mango expressed sequence tags and availability of well-established reference PF biosynthetic genes from other plant species allowed the identification of coding sequences of genes that may lead to the formation of important flavonoid compounds in mango fruits and facilitated characterisation of single nucleotide polymorphisms between varieties. We discovered an association between the extent of sequence variation and position in the pathway for up-stream genes. The high expression of PAL, C4H and CHS genes in mango peel compared to flesh is associated with high amounts of total phenolic contents in peels, which suggest that these genes have an influence on total flavonoid levels in mango fruit peel and flesh. In addition, the particularly high expression levels of ANR in KP and NDM peels compared to IW peel and the significant accumulation of its product epicatechin gallate (ECG) in those extracts reflects the rate-limiting role of ANR on ECG biosynthesis in mango. © 2015 Hoang et al.

Demographic history of India and mtDNA-sequence diversity

The demographic history of India was examined by comparing mtDNA sequences obtained from members of three culturally divergent Indian subpopulations (endogamous caste groups). While an inferred tree revealed some clustering according to caste affiliation, there was no clear separation into three genetically distinct groups along caste lines. Comparison of pairwise nucleotide difference distributions, however, did indicate a difference in growth patterns between two of the castes. The Brahmin population appears to have undergone either a rapid expansion or steady growth. The low-ranking Mukri caste, however, may have either maintained a roughly constant population size or undergone multiple bottlenecks during that period. Comparison of the Indian sequences to those obtained from other populations, using a tree, revealed that the Indian sequences, along with ah other non-African samples, form a starlike cluster. This cluster may represent a major expansion, possibly originating in southern Asia, taking place at some point after modern humans initially left Africa.

Mitochondrial cytochrome b gene sequence diversity in the Korean hare, Lepus coreanus Thomas (Mammalia, Lagomorpha)

Partial sequences of the mitochondrial cytochrome b gene of the Korean hare (Lepus coreanus) were analyzed to determine the degree of genetic diversity. Nine haPlotyes were observed, and the maximum Tamura-Nei nucleotide distance among them was 2.8%, indicating that genetic diversity of L. coreanus is moderate. In order to clarify the Korean hare's taxonomic status and relationship with the Manchurian hare (L. mandshuricus) and the Chinese hare (L. sinensis), these nine haplotypes of the Korean hare were compared with 13 haplotypes from five other species of eastern Asian Lepus including L. mandshwicus and L. sinensis. The Korean hare was distinct in its cytochrome b gene, and it is confirmed that L. coreanus is a valid species, as noted by Jones and Johnson (1965, Univ. Kansas Publ. (Mus. Nat. Hist.) 16:357). Further analyses of mtDNA cytochrome b gene with additional specimens of L. coreanus from North Korea and other species of Lepus from eastern Asia are needed to clarify the taxonomic status of the divergent mtDNA clades of L. mandshuricus and L. sinensis.

Genetic relationship of Chinese ethnic populations revealed by mtDNA sequence diversity

The origin and demographic history of the ethnic populations of China have not been clearly resolved. In this study, we examined the hypervariable segment I sequences (HVSI) of the mitochondrial DNA control region in 372 individuals from nine Chinese popu

Mitochondrial DNA sequence diversity and origin of Chinese domestic yak

In order to clarify the origin and genetic diversity of yak in China, we analysed mitochondrial DNA (mtDNA) control region sequences (similar to 891 bp) in 52 individuals from four domestic yak (Poephagus grunniens) breeds, as well as from a hybrid betwee

Relation of pretreatment sequence diversity in NS5A region of HCV genotype 1 with immune response between pegylated-INF/ribavirin therapy outcomes

Hepatitis C virus (HCV) infection frequently persists despite substantial virus-specific immune responses and the combination of pegylated interferon (INF)-alpha and ribavirin therapy. Major histocompatibility complex class I restricted CD8+ T cells are responsible for the control of viraemia in HCV infection, and several studies suggest protection against viral infection associated with specific HLAs. The reason for low rates of sustained viral response (SVR) in HCV patients remains unknown. Escape mutations in response to cytotoxic T lymphocyte are widely described; however, its influence in the treatment outcome is ill understood. Here, we investigate the differences in CD8 epitopes frequencies from the Los Alamos database between groups of patients that showed distinct response to pegylated alpha-INF with ribavirin therapy and test evidence of natural selection on the virus in those who failed treatment, using five maximum likelihood evolutionary models from PAML package. The group of sustained virological responders showed three epitopes with frequencies higher than Non-responders group, all had statistical support, and we observed evidence of selection pressure in the last group. No escape mutation was observed. Interestingly, the epitope VLSDFKTWL was 100% conserved in SVR group. These results suggest that the response to treatment can be explained by the increase in immune pressure, induced by interferon therapy, and the presence of those epitopes may represent an important factor in determining the outcome of therapy.

Sequence diversity and evolutionary dynamics of the dimorphic antigen merozoite surface protein-6 and other Msp genes of Plasmodium falciparum

Immune evasion by Plasmodium falciparum is favored by extensive allelic diversity of surface antigens. Some of them, most notably the vaccine-candidate merozoite surface protein (MSP)-1, exhibit a poorly understood pattern of allelic dimorphism, in which all observed alleles group into two highly diverged allelic families with few or no inter-family recombinants. Here we describe contrasting levels and patterns of sequence diversity in genes encoding three MSP-1-associated surface antigens of P. falciparum, ranging from an ancient allelic dimorphism in the Msp-6 gene to a near lack of allelic divergence in Msp-9 to a more classical multi-allele polymorphism in Msp-7 Other members of the Msp-7 gene family exhibit very little polymorphism in non-repetitive regions. A comparison of P. falciparum Msp-6 sequences to an orthologous sequence from P. reichenowi provided evidence for distinct evolutionary histories of the 5` and 3` segments of the dimorphic region in PfMsp-6, consistent with one dimorphic lineage having arisen from recombination between now-extinct ancestral alleles. In addition. we uncovered two surprising patterns of evolution in repetitive sequence. Firsts in Msp-6, large deletions are associated with (nearly) identical sequence motifs at their borders. Second, a comparison of PfMsp-9 with the P. reichenowi ortholog indicated retention of a significant inter-unit diversity within an 18-base pair repeat within the coding region of P. falciparum, but homogenization in P. reichenowi. (C) 2009 Elsevier B.V. All rights reserved.

Plasmodium falciparum: sequence diversity and antibody recognition of the Merozoite surface protein-2 (MSP-2) in Brazilian Amazonia

The merozoite surface protein-2 (MSP-2) of Plasmodium falciparum comprises repeats flanked by dimorphic domains defining the allelic families FC27 and IC1. Here, we examined sequence diversity at the msp-2 locus in Brazil and its impact on MSP-2 antibody recognition by local patients. Only 25 unique partial sequences of msp-2 were found in 61 isolates examined. The finding of identical msp-2 sequences in unrelated parasites, collected 6-13 years apart, suggests that no major directional selection is exerted by variant-specific immunity in this malaria-endemic area. To examine antibody cross-reactivity, recombinant polypeptides derived from locally prevalent and foreign MSP-2 variants were used in ELISA. Foreign IC1-type variants, such as 3D7 (currently tested for human vaccination), were less frequently recognized than FC27-type and local IC1-type variants. Antibodies discriminated between local and foreign IC1-type variants, but cross-recognized structurally different local IC1-type variants. The use of evolutionary models of MSP-2 is suggested to design vaccines that minimize differences between local parasites and vaccine antigens. (C) 2004 Elsevier B.V. All rights reserved.

Relation of pretreatment sequence diversity in NS5A region of HCV genotype 1 with immune response between pegylated-INF/ribavirin therapy outcomes

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Origins of sequence diversity in the malaria vaccine candidate merozoite surface protein-2 (MSP-2) in Amazonian isolates of Plasmodium falciparum

The recent evolution of Plasmodium falciparum is at odds with the extensive polymorphism found in most genes coding for antigens. Here, we examined the patterns and putative mechanisms of sequence diversification in the merozoite surface protein-2 (MSP-2), a major malarial repetitive surface antigen. We compared the msp-2 gene sequences from closely related clones derived from sympatric parasite isolates from Brazilian Amazonia and used microsatellite typing to examine, in these same clones, the haplotype background of chromosome 2, where msp-2 is located. We found examples of msp-2 sequence rearrangements putatively created by nonreciprocal recombinational events, such as replication slippage and gene conversion, while maintaining the chromosome haplotype. We conclude that these nonreciprocal recombination events may represent a major source of antigenic diversity in MSP-2 in P falciparum populations with low rates of classical meiotic recombination. (c) 2006 Elsevier B.V. All rights reserved.

A genome-to-genome analysis of associations between human genetic variation, HIV-1 sequence diversity, and viral control

HIV-1 sequence diversity is affected by selection pressures arising from host genomic factors. Using paired human and viral data from 1071 individuals, we ran >3000 genome-wide scans, testing for associations between host DNA polymorphisms, HIV-1 sequence variation and plasma viral load (VL), while considering human and viral population structure. We observed significant human SNP associations to a total of 48 HIV-1 amino acid variants (p<2.4 × 10−12). All associated SNPs mapped to the HLA class I region. Clinical relevance of host and pathogen variation was assessed using VL results. We identified two critical advantages to the use of viral variation for identifying host factors: (1) association signals are much stronger for HIV-1 sequence variants than VL, reflecting the ‘intermediate phenotype’ nature of viral variation; (2) association testing can be run without any clinical data. The proposed genome-to-genome approach highlights sites of genomic conflict and is a strategy generally applicable to studies of host–pathogen interaction.

Unexpected beta2-microglobulin sequence diversity in individual rainbow trout.

For mammals beta2-microglobulin (beta2m), the light chain of major histocompatibility complex (MHC) class I molecules, is invariant (or highly conserved) and is encoded by a single gene unlinked to the MHC. We find that beta2m of a salmonid fish, the rainbow trout (Oncorhynchus mykiss), does not conform to the mammalian paradigm. Ten of 12 randomly selected beta2m cDNA clones from an individual fish have different nucleotide sequences. A complex restriction fragment length polymorphism pattern is observed with rainbow trout, suggesting multiple beta2m genes in the genome, in excess of the two genes expected from the ancestral salmonid tetraploidy. Additional duplication and diversification of the beta2m genes might have occurred subsequently. Variation in the beta2m cDNA sequences is mainly at sites that do not perturb the structure of the mature beta2m protein, showing that the observed diversity of the trout beta2m genes is not primarily a result of pathogen selection.