Frequency Of Nucleotide Substitutions Research Articles

Evolution of retrocopies in the context of HUSH silencing

Retrotransposition is one of the main factors responsible for gene duplication and thus genome evolution. However, the sequences that undergo this process are not only an excellent source of biological diversity, but in certain cases also pose a threat to the integrity of the DNA. One of the mechanisms that protects against the incorporation of mobile elements is the HUSH complex, which is responsible for silencing long, intronless, transcriptionally active transposed sequences that are rich in adenine on the sense strand. In this study, broad sets of human and porcine retrocopies were analysed with respect to the above factors, taking into account evolution of these molecules. Analysis of expression pattern, genomic structure, transcript length, and nucleotide substitution frequency showed the strong relationship between the expression level and exon length as well as the protective nature of introns. The results of the studies also showed that there is no direct correlation between the expression level and adenine content. However, protein-coding retrocopies, which have a lower adenine content, have a significantly higher expression level than the adenine-rich non-coding but expressed retrocopies. Therefore, although the mechanism of HUSH silencing may be an important part of the regulation of retrocopy expression, it is one component of a more complex molecular network that remains to be elucidated.

Characterization of the Spectrum of Mitochondrial DNA Nucleotide Substitutions in Human Populations in High Altitude Environments

Using phylogenetic analysis of the nucleotide sequences of whole mitochondrial genomes (mtDNA), the spectra of germinal nucleotide substitutions (on the L-chain of mtDNA) were reconstructed in highland populations of the Pamirs and Tibet in comparison with regional indigenous groups of West Asia, Northeast Siberia, and South Siberia. No differences were found in the distribution of nucleotide substitution frequencies in the mtDNA spectra depending on the population distribution by altitude. Pyrimidine transitions dominate in all mtDNA spectra, and T → C substitutions are the most frequent among them. Next in frequency in most regional groups are A → G substitutions, but in the Pamir and northeast Asian groups G → A substitutions are prevalent. Of the transversions in all populations studied C → A replacements were found to be predominant, except for the Tibetan one, where A → C substitutions are more frequent. The lack of differences in the distribution of mtDNA mutations in high-altitude and non-highland populations indicates that the structure of mtDNA spectra in human populations is independent of the oxidative stress intensity in mitochondria.

Conserved G-Quadruplex-Forming Sequences in Mammalian TERT Promoters and Their Effect on Mutation Frequency.

Somatic mutations in the promoter region of the human telomerase reverse transcriptase (hTERT) gene have been identified in many types of cancer. The hTERT promoter is known to be enriched with sequences that enable the formation of G-quadruplex (G4) structures, whose presence is associated with elevated mutagenicity and genome instability. Here, we used a bioinformatics tool (QGRS mapper) to search for G4-forming sequences (G4 motifs) in the 1000 bp TERT promoter regions of 141 mammalian species belonging to 20 orders, 5 of which, including primates and predators, contain more than 10 species. Groups of conserved G4 motifs and single-nucleotide variants within these groups were discovered using a block alignment approach (based on the Nucleotide PanGenome explorer). It has been shown that: (i) G4 motifs are predominantly located in the region proximal to the transcription start site (up to 400 bp) and are over-represented on the non-coding strand of the TERT promoters, (ii) 11 to 22% of the G4 motifs found are evolutionarily conserved across the related organisms, and (iii) a statistically significant higher frequency of nucleotide substitutions in the conserved G4 motifs compared to the surrounding regions was confirmed only for the order Primates. These data support the assumption that G4s can interfere with the DNA repair process and affect the evolutionary adaptation of organisms and species.

Characterization of Potential Molecular Markers in Lac Insect Kerria lacca (Kerr) Responsible for Lac Production.

Simple SummaryLac insect Kerria lacca (Hemiptera, Kerriidae) is widely used to produce lac in tropical and subtropical areas. This species is genomically understudied and there is no information present on single-nucleotide polymorphisms and insertions and deletions of lac insects. This study involves distribution and characterization of molecular markers in K. lacca. Nucleotide substitution frequency and transition to transversion ratio were highest at the late adult stage. The maximum numbers of SNPs and InDels were distributed on chromosome 8 and 1, respectively. Most variants at low frequencies were detected in InDels compared to SNPs. Gene function analysis showed several gene variations, and candidate SNPs that were involved in lac synthesis. This work paves the way to determining further mechanisms involved in increased lac production.Kerria lacca (Kerr) is an important lac insect extensively used in industrial products in the form of resin, wax and dye. The scarce knowledge on molecular markers for K. lacca is a barrier in elucidating genetic information. Our study identified a total of 16,921 single-nucleotide polymorphisms (SNPs), and 6231 insertions and deletions (InDels)—of which, intergenic variation accounted for 41.22% and 56.30%, and exonic variation accounted for 39.10% and 17.46%, of SNPs and InDels, respectively. Observation of SNPs suggested that nucleotide substitution frequency and transition to transversion (Ts/Tv) ratio were highest at the late adult stage, 3.97, compared to at the other stages, with a genome-wide Ts/Tv ratio of 2.95. The maximum number of SNPs, 2853 (16.86%), was identified in chromosome 8, while the lowest, 1126 (6.65%), was identified in chromosome 7. The maximum and minimum numbers of InDels were located on chromosome 1 and 7, with 834 (13.38%) and 519 (8.33%), respectively. Annotation showed that highest numbers of exonic and intergenic SNPs were present at the late adult stage, whereas the maximum number of InDels was found at the larval stage. On the basis of gene function, 47 gene variations were screened and 23 candidate genes were identified in associations with lac production. Concluding work will enhance knowledge on molecular markers to facilitate an increase in lac production in K. lacca as well as other lac insects.

Genetic diversity of the cultivated Salvia miltiorrhiza populations revealed by four intergenic spacers.

For better understanding the genetic diversity and phylogeny of the cultivated Salvia miltiorrhiza populations, four intergenic spacer sequences, ETS, psbA-trnH, trnL-trnF, and ycf1-rps15 of the 40 populations collected from China were Polymerase Chain Reaction (PCR) amplified, analyzed both individually and in combination. Haplotype diversity analysis showed that the cultivated S. miltiorrhiza populations had a very rich genetic diversity and an excellent capacity to resist environmental pressure. The best-fit nucleotide substitution models for ETS, psbA-trnH, trnL-trnF, ycf1-rps15, and their combined sequences were HKY+I, T92, T92, T92+G, and T92+G, respectively; the nucleotide conversion frequency in the combined sequences was lower than the transversion, and the relatively high nucleotide substitution frequencies suggests its high genetic variability. Neutral tests showed that the spacer sequences of the populations conform with the neutral evolution model, and there has been no current expansion events occurred. Phylogeny analyses based on both the individual and the combined sequences showed that the 40 populations were clustered in two clades with a very similar topological structure. The discrimination rate of the combined sequence marker is significantly increased to 52.5% (21 populations) over the highest 35% (13 populations) by the single marker of ETS, though still inadequate but a big step forward. Further exploration of more DNA markers is needed. This study for the first time revealed the rich genetic diversity and phylogeny of the currently cultivated S. miltiorrhiza populations in China and provides novel alternative molecular markers for the genetic identification and resources evaluation of the cultivated S. miltiorrhiza populations.

Non-B DNA: a major contributor to small- and large-scale variation in nucleotide substitution frequencies across the genome.

Approximately 13% of the human genome can fold into non-canonical (non-B) DNA structures (e.g. G-quadruplexes, Z-DNA, etc.), which have been implicated in vital cellular processes. Non-B DNA also hinders replication, increasing errors and facilitating mutagenesis, yet its contribution to genome-wide variation in mutation rates remains unexplored. Here, we conducted a comprehensive analysis of nucleotide substitution frequencies at non-B DNA loci within noncoding, non-repetitive genome regions, their ±2 kb flanking regions, and 1-Megabase windows, using human-orangutan divergence and human single-nucleotide polymorphisms. Functional data analysis at single-base resolution demonstrated that substitution frequencies are usually elevated at non-B DNA, with patterns specific to each non-B DNA type. Mirror, direct and inverted repeats have higher substitution frequencies in spacers than in repeat arms, whereas G-quadruplexes, particularly stable ones, have higher substitution frequencies in loops than in stems. Several non-B DNA types also affect substitution frequencies in their flanking regions. Finally, non-B DNA explains more variation than any other predictor in multiple regression models for diversity or divergence at 1-Megabase scale. Thus, non-B DNA substantially contributes to variation in substitution frequencies at small and large scales. Our results highlight the role of non-B DNA in germline mutagenesis with implications to evolution and genetic diseases.

Comparative insight into the genomic landscape of SARS-CoV-2 and identification of mutations associated with the origin of infection and diversity.

The analyses of 2325 severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) genomes revealed 107, 162, and 65 nucleotide substitutions in the coding region of SARS-CoV-2 from the three continents America, Europe, and Asia, respectively. Of these nucleotide substitutions 58, 94, and 37 were nonsynonymous types mostly present in the Nsp2, Nsp3, Spike, and ORF9. A continent-specific phylogram analyses clustered the SARS-CoV-2 in the different group based on the frequency of nucleotide substitutions. Detailed analyses about the continent-specific amino acid changes and their effectiveness by SNAP2 software was investigated. We found 11 common nonsynonymous mutations;amongthem, two novel effective mutations were identified in ORF9 (S194L and S202N). Intriguingly, ORF9 encodes nucleocapsid phosphoprotein possessing many effective mutations across continents and could be a potential candidate after the spike protein for studying the role of mutation in viral assembly and pathogenesis. Among the two forms of certain frequent mutation, one form is more prevalent in Europe continents (Nsp12:L314, Nsp13:P504, Nsp13:Y541, Spike:G614, and ORF8:L84) while other forms aremore prevalent in American (Nsp12:P314, Nsp13:L504, Nsp13:C541, Spike:D614, and ORF8:L84) and Asian continents (Spike:D614), indicating the spatial and temporal dynamics of SARS-CoV-2. We identified highly conserved 38 regions and among these regions, 11 siRNAs were predicted on stringent criteria that can be used to suppress the expression of viral genes and the corresponding reduction of human viral infections. The present investigation provides information ondifferent mutations and will pave the way for differentiating strains based on virulence and their use in the development of better antiviral therapy.

Characterization of mutational processes in ER+ metastatic breast cancer.

1019 Background: Diverse mutational processes are active at different stages of cancer development and cause the mutational burden in cancer cells. These processes leave distinctive signatures, characterized by varying frequencies of nucleotide substitutions and their flanking base-pairs. Mutational signatures found in early stage breast cancer have been described, but less is known about the mutational processes in metastatic breast cancer (MBC). Methods: We performed whole exome sequencing (WES) on 209 metastatic tumor biopsies from patients with ER+ MBC. In a subset of patients, we obtained and performed WES on the matched primary biopsy (n = 62). Decomposition of mutational signatures was done by determining the linear combination of 30 previously reported mutational signatures (Alexandrov et al., 2013) that most accurately reconstructs the mutational profile of each tumor sample. For comparison, signature decomposition was also performed on 695 treatment-naïve primary ER+ breast cancers form The Cancer Genome Atlas (TCGA). Results: Cohort-level comparison of signature activity between our cohort of MBC biopsies versus the primary ER+ breast cancers from TCGA demonstrated enrichment of Signature 13 in MBC samples (q = 0.002), attributed to activity of the APOBEC family of cytidine deaminases and implicated in resistance to endocrine therapies. Activity of Signature 22, associated to activity of the nucleotide-excision repair mechanism, was depleted in MBC tumors relative to TCGA ER+ breast tumors (q = 0.02). Decomposition of mutational signatures at the clone level in the subset of patients with both metastatic and matched primary samples allowed us to infer signatures acquired in the metastatic setting. Among mutations in acquired clones, we found enrichment of APOBEC signatures (Sig. 2, q = 0.018; Sig. 13, q = 0.003). Conversely, activity of Signature 1, attributed to the normal aging process, was enriched in mutations shared between primary and MBC samples (q = 0.004). This suggests that mutations developed early in tumor evolution may be driven primarily by natural processes, while later mutations may reflect selective pressures imposed by therapies and cellular defects acquired during disease progression. Conclusions: The profile of somatic mutations in MBC differs from early stage breast cancer. APOBEC-related mutations arising late during tumor evolution are one possible mechanism for this difference. Further studying the associations between clinical features and therapies to mutational processes could inform design of clinical trials and management of MBC.

The non-random patterns of genetic variation induced by asymmetric somatic hybridization in wheat

BackgroundAsymmetric somatic hybridization is an efficient crop breeding approach by introducing several exogenous chromatin fragments, which leads to genomic shock and therefore induces genome-wide genetic variation. However, the fundamental question concerning the genetic variation such as whether it occurs randomly and suffers from selection pressure remains unknown.ResultsHere, we explored this issue by comparing expressed sequence tags of a common wheat cultivar and its asymmetric somatic hybrid line. Both nucleotide substitutions and indels (insertions and deletions) had lower frequencies in coding sequences than in un-translated regions. The frequencies of nucleotide substitutions and indels were both comparable between chromosomes with and without introgressed fragments. Nucleotide substitutions distributed unevenly and were preferential to indel-flanking sequences, and the frequency of nucleotide substitutions at 5′-flanking sequences of indels was obviously higher in chromosomes with introgressed fragments than in those without exogenous fragment. Nucleotide substitutions and indels both had various frequencies among seven groups of allelic chromosomes, and the frequencies of nucleotide substitutions were strongly negatively correlative to those of indels. Among three sets of genomes, the frequencies of nucleotide substitutions and indels were both heterogeneous, and the frequencies of nucleotide substitutions exhibited drastically positive correlation to those of indels.ConclusionsOur work demonstrates that the genetic variation induced by asymmetric somatic hybridization is attributed to both whole genomic shock and local chromosomal shock, which is a predetermined and non-random genetic event being closely associated with selection pressure. Asymmetric somatic hybrids provide a worthwhile model to further investigate the nature of genomic shock induced genetic variation.

Flanking monomer repeats determine decreased context complexity of single nucleotide polymorphism sites in the human genome

The study of the dependence of the mutation frequency in human genome was conducted by the example of a set of documented single nucleotide polymorphisms (SNPs) from the 1000 genomes project. The tasks of the development of new computer methods for the statistical analysis of genetic texts based on estimations of sequences complexity were considered. The application of the complexity profiles in a sliding window to the analysis of sites containing single nucleotide polymorphisms in the human genome was demonstrated. A local decrease in the text complexity near SNPs was established. Based on the analysis of the complexity profiles in the regions containing SNPs, it was demonstrated that the flanking monomer repeats determine the decreased context complexity of single nucleotide polymorphism sites in human genome. The effect of local decrease in the text complexity level for sequences flanking SNP sites was confirmed for the data on polymorphisms in the rat and mouse genomes. Differences in the context organization for coding and regulatory sequences (that are reflected in the text complexity of nucleotide sequences containing SNPs) were determined. Changes in the point mutation frequencies were previously demonstrated for the sequences containing microsatellites. Using more general mathematical apparatus and more complete data, a saturation of local genomic surroundings containing SNPs with polytracts and simple repeated sequences was demonstrated in this work. Oligonucleotides with increased frequency in the genomic SNP surroundings in human were determined; their association with polytracts was demonstrated. The presence of polytracts can indicate a greater probability of a break in the double DNA strand at this point (resulting in an increased frequency of nucleotide substitutions). The obtained estimations were determined by a previously developed complex of computer programs, which allows us to efficiently determine the frequency spectrum of oligonucleotides with a fixed length, to compare nucleotide frequencies in a larger sample (in addition to estimating the complexity of the phased sequences).

DNA barcoding reflects the diversity and variety of brooding traits of fish species in the family Syngnathidae along China’s coast

DNA barcoding offers a rapid and accurate assessment for species labelling and identification. Here, we investigated the performance of DNA barcodes in a sample of Syngnathidae, a unique group of fish with male pregnancy. A total of 1002 DNA barcodes using 649 base-pair fragments of the cytochrome c oxidase subunit I (COI) gene were generated. All species were associated with distinct DNA barcode and could be readily distinguished; seven of the COI barcode clusters represented the first species records submitted to the Barcode of Life Data Systems (BOLD) and GenBank databases. In the Neighbor-joining tree of COI sequences, two major clusters (Gastrophori and Urophori) were displayed, which could also be classified by their brood pouches. However, the tail-brooding Hippichthys cyanospilus and trunk-brooding Syngnathoides biaculeatus were reverse-clustered together with the Gastrophori and the Urophori, respectively. Both seahorse and pipefish sequences showed high frequencies of nucleotide substitutions. The probability of nucleotide substitutions, those in pipefish were higher than that of all seahorses. However, we identified no signal for positive selection based on the COI gene in any of the data sets. Our results supported DNA barcoding as an efficient molecular tool for achieving better monitoring, conservation, and management of fisheries.

Genetic differentiation of sockeye salmon Oncorhynchus nerka from Kamchatka River basin and the lake–river systems of the west coast of the bering sea as inferred from data on single nucleotide polymorphism

The variability of 45 single nucleotide polymorphism loci (SNP) was studied in sockeye salmon from the Kamchatka River basin and four lake-river systems of the west coast of the Bering Sea. Based on the genetic differentiation estimates for the largest sockeye salmon populations of Eastern Kamchatka and Chukotka, the examined samples were combined into two regional groups represented by the population of the Kamchatka River drainage, which included numerous local subpopulations and seasonal races, and the northern population grouping from the rivers of Olutorsko-Navarinsky raion, wherein the sockeye salmon from Maynypilginskaya Lake-River system was relatively isolated. Considerable divergence was observed between the island (Sarannoe Lake, Bering Island) and continental populations. Genetic heterogeneity was revealed and groups of early- and late-maturing individuals were isolated in the sample of late-run sockeye salmon from Kamchatka River. In Apuka River, subdivision of the spawning run into two genetically distinct spatial and temporal groupings was also observed. The results suggest that the differentiation of sockeye salmon samples by single nucleotide substitution frequencies was largely due to differences in the direction and strength of local selection at some loci in the population complexes and intrapopulation groupings from the examined river basins of Eastern Kamchatka, Chukotka, and Commander Islands.

The complete nucleotide sequence of the mitochondrial genome of Bactrocera minax (Diptera: Tephritidae).

The complete 16,043 bp mitochondrial genome (mitogenome) of Bactrocera minax (Diptera: Tephritidae) has been sequenced. The genome encodes 37 genes usually found in insect mitogenomes. The mitogenome information for B. minax was compared to the homologous sequences of Bactrocera oleae, Bactrocera tryoni, Bactrocera philippinensis, Bactrocera carambolae, Bactrocera papayae, Bactrocera dorsalis, Bactrocera correcta, Bactrocera cucurbitae and Ceratitis capitata. The analysis indicated the structure and organization are typical of, and similar to, the nine closely related species mentioned above, although it contains the lowest genome-wide A+T content (67.3%). Four short intergenic spacers with a high degree of conservation among the nine tephritid species mentioned above and B. minax were observed, which also have clear counterparts in the control regions (CRs). Correlation analysis among these ten tephritid species revealed close positive correlation between the A+T content of zero-fold degenerate sites (P0FD), the ratio of nucleotide substitution frequency at P0FD sites to all degenerate sites (zero-fold degenerate sites, two-fold degenerate sites and four-fold degenerate sites) and amino acid sequence distance (ASD) were found. Further, significant positive correlation was observed between the A+T content of four-fold degenerate sites (P4FD) and the ratio of nucleotide substitution frequency at P4FD sites to all degenerate sites; however, we found significant negative correlation between ASD and the A+T content of P4FD, and the ratio of nucleotide substitution frequency at P4FD sites to all degenerate sites. A higher nucleotide substitution frequency at non-synonymous sites compared to synonymous sites was observed in nad4, the first time that has been observed in an insect mitogenome. A poly(T) stretch at the 5′ end of the CR followed by a [TA(A)]n-like stretch was also found. In addition, a highly conserved G+A-rich sequence block was observed in front of the poly(T) stretch among the ten tephritid species and two tandem repeats were present in the CR.

Comparative analysis of the genomes of Clostera anastomosis (L.) granulovirus and Clostera anachoreta granulovirus

Clostera anastomosis (L.) granulovirus (CaLGV) and Clostera anachoreta granulovirus (ClanGV) are both capable of infecting each other's native host insects. Despite this, we have little information on their genetic relationship. The complete nucleotide sequence of CaLGV was determined and compared with that of the genome of ClanGV. The circular, double-stranded DNA CaLGV genome (GenBank accession no. KC179784) had a G+C content of 46.7% and was 101,818 bp in size (331 bp larger than the ClanGV genome). Overall, the CaLGV nucleotide sequence was found to be 90% identical to that of ClanGV. It contained a total of 123 ORFs, 119 of which had ClanGV homologues, with an identical transcription direction and ORF organization. Seventy-five of the 119 ORFs showed 90% or greater identity to their ClanGV homologues. CaLGV contained only a single identifiable homologous region (hrs)/repeat region (similar to ClanGV hr4). The mean frequency of nucleotide substitutions in the CaLGV/ClanGV coding regions was 8.33%. CaLGV contained four unique ORFs (CaL23, CaL39, CaL48 and CaL92). Eight ORFs found in both CaLGV and ClanGV have no homologues in other baculoviruses. Intergenic regions of CaLGV and ClanGV occupied 6.6% and 7% of their respective genomes. CaLGV appears closer phylogenetically to ClanGV than to any other baculoviruses.

Ethyl Methane Sulphonate (EMS) Induced Mutagenesis in Malaysian Rice (cv. MR219) for Lethal Dose Determination

Chemical and physical mutagenesis has been used to increase genetic variability in crop plants. More than 430 new varieties have been derived as mutants of rice (Oryza sativa L.) via the application of different mutagenic agents. Chemical mutagens such as ethyl methane sulphonate (EMS), diepoxybutane-derived (DEB), sodium azide and irradiation (Gamma rays, X-rays and fast neutrons) have been widely used to induce a large number of functional variations in rice and others crops. Among chemical mutagens, the alkylating agent, ethyl methane sulfonate (EMS) is the most commonly used in plants as it causes a high frequency of nucleotide substitutions, as detected in different genomes. In this study, seeds of potential genotype of the popular variety, (Oryza sativa L. spp. Indica cv. MR219) were treated with EMS at concentrations of 0.25%, 0.50%, 0.75%, 1%, 1.25%, 1.5% and 2%. Sensitivity to EMS was determined by various measurements on the M1 generation. As concentration of applied EMS increased, will decrease in germination, seedling height, root length and emergence under field conditions was observed in M1 generation as compared to the non-treatment control. Plant height and root length also decreased with increases in EMS mutagenesis in an approximately linear fashion. The LD25 and LD50 values were observed based on growth reduction of seedlings after EMS treatment with 0.25% and 0.50% on the rice variety (Oryza sativa L. spp. Indica cv. MR219).

Genome evolution in diploid and tetraploid Coffea species as revealed by comparative analysis of orthologous genome segments

Sequence comparison of orthologous regions enables estimation of the divergence between genomes, analysis of their evolution and detection of particular features of the genomes, such as sequence rearrangements and transposable elements. Despite the economic importance of Coffea species, little genomic information is currently available. Coffea is a relatively young genus that includes more than one hundred diploid species and a single tetraploid species. Three Coffea orthologous regions of 470-900 kb were analyzed and compared: both subgenomes of allotetraploid Coffea arabica (contributed by the diploid species Coffea eugenioides and Coffea canephora) and the genome of diploid C. canephora. Sequence divergence was calculated on global alignments or on coding and non-coding sequences separately. A search for transposable elements detected 43 retrotransposons and 198 transposons in the sequences analyzed. Comparative insertion analysis made it possible to locate 165 TE insertions in the phylogenetic tree of the three genomes/subgenomes. In the tetraploid C. arabica, a homoeologous non-reciprocal transposition (HNRT) was detected and characterized: a 50 kb region of the C. eugenioides derived subgenome replaced the C. canephora derived counterpart. Comparative sequence analysis on three Coffea genomes/subgenomes revealed almost perfect gene synteny, low sequence divergence and a high number of shared transposable elements. Compared to the results of similar analysis in other genera (Aegilops/Triticum and Oryza), Coffea genomes/subgenomes appeared to be dramatically less diverged, which is consistent with the relatively recent radiation of the Coffea genus. Based on nucleotide substitution frequency, the HNRT was dated at 10,000-50,000 years BP, which is also the most recent estimation of the origin of C. arabica.

Characterization of the gilthead seabream ( Sparus aurata L.) transferrin gene: Genomic structure, constitutive expression and SNP variation

Transferrin (Tf) is a multi-function protein with a central role in iron metabolism, and it is this function that is associated with a role in the innate immune system response. The clear link between Tf and immune defense mechanism lead to propose Tf as a candidate gene for disease resistance. In this study, genomic and cDNA sequences of Tf gene in gilthead seabream ( Sparus aurata L.) ( SaTf gene), were identified and characterized. SaTf gene structure consists of a coding region of 2076 nucleotides divided into 17 exons and a no-coding region that includes 16 introns and spans 5495 nucleotides. The deduced Tf protein for gilthead seabream is composed of 691 amino acids and consists of an initial peptide and two lobes (N- and C-lobes). This gene structure is similar to that of previously described Tf genes in other fish species. RT-PCR analyses carried out in different tissues and two developmental stages showed tissue-and stage-specific Tf expression in gilthead seabream. Finally, by sequencing the transferrin genomic sequences of 20 unrelated seabreams, 31 SNPs were identified. These data allowed the estimation of the frequency of nucleotide substitution in the SaTf gene as 1SNP per 253 bp. SNPs were detected in different regions of the genomic sequence but they were mainly localized in non-coding regions, specifically, SNP frequency in non-coding regions was fifteen-fold higher than within coding regions.

Pentatricopeptide Repeat Proteins Constrain Genome Evolution in Chloroplasts

Higher plants encode hundreds of pentatricopeptide repeat proteins (PPRs) that are involved in several types of RNA processing reactions. Most PPR genes are predicted to be targeted to chloroplasts or mitochondria, and many are known to affect organellar gene expression. In some cases, RNA binding has been directly demonstrated, and the sequences of the cis-elements are known. In this work, we demonstrate that RNA cis-elements recognized by PPRs are constrained in chloroplast genome evolution. Cis-elements for two PPR genes and several RNA editing sites were analyzed for sequence changes by pairwise nucleotide substitution frequency, pairwise indel frequency, and maximum likelihood (ML) phylogenetic distances. All three of these analyses demonstrated that sequences within the cis-element are highly conserved compared with surrounding sequences. In addition, we have compared sequences around chloroplast editing sites and homologous sequences in species that lack an editing site due to the presence of a genomic T. Cis-elements for RNA editing sites are highly conserved in angiosperms; by contrast, comparable sequences around a genomically encoded T exhibit higher rates of nucleotide substitution, higher frequencies of indels, and greater ML distances. The loss in requirement for editing to create the ndhD start codon has resulted in the conversion of the PPR gene responsible for editing that site to a pseudogene. We show that organellar dependence on nuclear-encoded PPR proteins for gene expression has constrained the evolution of cis-elements that are required at the level of RNA processing. Thus, the expansion of the PPR gene family in plants has had a dramatic effect on the evolution of plant organelle genomes.

Selective pressure for allelic diversity in SeM of Streptococcus equi does not affect immunoreactive proteins SzPSe or Se18.9

Streptococcus equi, a clone or biovar of an ancestral Streptococcus zooepidemicus of Lancefield group C causes equine strangles, a highly contagious tonsillitis and lymphadenitis of the head and neck. At least 74 alleles based on N-terminal amino acid sequence of the anti-phagocytic SeM have been observed among isolates of S. equi from N. America, Europe and Japan. A d N/ d S ratio of 5.93 for the 5′ region of sem is indicative of positive selective pressure. The aim of this study was to determine whether variations in SeM were accompanied by variations in the surface exposed SzPSe and secreted Se18.9, both of which bind to equine tonsillar epithelium and, along with SeM, elicit strong nasopharyngeal IgA responses during convalescence. Sequences of genes for these proteins from 25 S. equi expressing 19 different SeM alleles isolated over 40 years in different countries were compared. No variation was observed in szpse, except for an Australian isolate with a deletion of a single repeat in the 3′ end of the gene. Interestingly, only two SNP loci were detected in se18.9 compared to 93 and 55 in sem and szpse, respectively. The high frequency of nucleotide substitutions in s zpse may be related to its mosaic structure since this gene in S. zooepidemicus exists in a variety of combinations of sequence segments and has a central hypervariable region that includes exogenous DNA sequence based on an atypical G–C percentage. In summary, the results of this study document very different responses of streptococcal genes for 3 immunoreactive proteins to selection pressure of the nasopharyngeal mucosal immune response.

Single-nucleotide polymorphism, linkage disequilibrium and geographic structure in the malaria parasite Plasmodium vivax: prospects for genome-wide association studies

BackgroundThe ideal malaria parasite populations for initial mapping of genomic regions contributing to phenotypes such as drug resistance and virulence, through genome-wide association studies, are those with high genetic diversity, allowing for numerous informative markers, and rare meiotic recombination, allowing for strong linkage disequilibrium (LD) between markers and phenotype-determining loci. However, levels of genetic diversity and LD in field populations of the major human malaria parasite P. vivax remain little characterized.ResultsWe examined single-nucleotide polymorphisms (SNPs) and LD patterns across a 100-kb chromosome segment of P. vivax in 238 field isolates from areas of low to moderate malaria endemicity in South America and Asia, where LD tends to be more extensive than in holoendemic populations, and in two monkey-adapted strains (Salvador-I, from El Salvador, and Belem, from Brazil). We found varying levels of SNP diversity and LD across populations, with the highest diversity and strongest LD in the area of lowest malaria transmission. We found several clusters of contiguous markers with rare meiotic recombination and characterized a relatively conserved haplotype structure among populations, suggesting the existence of recombination hotspots in the genome region analyzed. Both silent and nonsynonymous SNPs revealed substantial between-population differentiation, which accounted for ~40% of the overall genetic diversity observed. Although parasites clustered according to their continental origin, we found evidence for substructure within the Brazilian population of P. vivax. We also explored between-population differentiation patterns revealed by loci putatively affected by natural selection and found marked geographic variation in frequencies of nucleotide substitutions at the pvmdr-1 locus, putatively associated with drug resistance.ConclusionThese findings support the feasibility of genome-wide association studies in carefully selected populations of P. vivax, using relatively low densities of markers, but underscore the risk of false positives caused by population structure at both local and regional levels.See commentary: http://www.biomedcentral.com/1741-7007/8/90