Overcoming taxonomic challenges in DNA barcoding for improvement of identification and preservation of clariid catfish species

DNA barcoding without assessing reliability and validity causes taxonomic errors of species identification, which is responsible for disruptions of their conservation and aquaculture industry. Although DNA barcoding facilitates molecular identification and phylogenetic analysis of species, its availability in clariid catfish lineage remains uncertain. In this study, DNA barcoding was developed and validated for clariid catfish. 2,970 barcode sequences from mitochondrial cytochrome c oxidase I (COI) and cytochrome b (Cytb) genes and D-loop sequences were analyzed for 37 clariid catfish species. The highest intraspecific nearest neighbor distances were 85.47%, 98.03%, and 89.10% for COI, Cytb, and D-loop sequences, respectively. This suggests that the Cytb gene is the most appropriate for identifying clariid catfish and can serve as a standard region for DNA barcoding. A positive barcoding gap between interspecific and intraspecific sequence divergence was observed in the Cytb dataset but not in the COI and D-loop datasets. Intraspecific variation was typically less than 4.4%, whereas interspecific variation was generally more than 66.9%. However, a species complex was detected in walking catfish and significant intraspecific sequence divergence was observed in North African catfish. These findings suggest the need to focus on developing a DNA barcoding system for classifying clariid catfish properly and to validate its efficacy for a wider range of clariid catfish. With an enriched database of multiple sequences from a target species and its genus, species identification can be more accurate and biodiversity assessment of the species can be facilitated.

The effects of inversion polymorphisms on patterns of neutral genetic diversity

The strong reduction in the frequency of recombination in heterozygotes for an inversion and a standard gene arrangement causes the arrangements to become partially isolated genetically, resulting in sequence divergence between them and changes in the levels of neutral variability at nucleotide sites within each arrangement class. Previous theoretical studies on the effects of inversions on neutral variability have assumed either that the population is panmictic or that it is divided into 2 populations subject to divergent selection. Here, the theory is extended to a model of an arbitrary number of demes connected by migration, using a finite island model with the inversion present at the same frequency in all demes. Recursion relations for mean pairwise coalescent times are used to obtain simple approximate expressions for diversity and divergence statistics for an inversion polymorphism at equilibrium under recombination and drift, and for the approach to equilibrium following the sweep of an inversion to a stable intermediate frequency. The effects of an inversion polymorphism on patterns of linkage disequilibrium are also examined. The reduction in effective recombination rate caused by population subdivision can have significant effects on these statistics. The theoretical results are discussed in relation to population genomic data on inversion polymorphisms, with an emphasis on Drosophila melanogaster. Methods are proposed for testing whether or not inversions are close to recombination-drift equilibrium, and for estimating the rate of recombinational exchange in heterozygotes for inversions; difficulties involved in estimating the ages of inversions are also discussed.

Plastome evolution and phylogeny of the tribe Ruteae (Rutaceae)

Rutaceae is a large family, and the genus-level classification in the subfamilies or tribes of this family is not unified based on different taxonomic treatments. Until now, phylogenetic relationships of some genera in traditional tribe Ruteae have not been clearly resolved. In this study, seven new complete plastomes of this tribe were sequenced, and a comparative analysis was performed to investigate their plastome characteristics and evolution. In addition, we inferred the phylogenetic relationships of Ruteae based on complete plastome and nuclear ITS data. All plastomes exhibited a typical quadripartite structure and were relatively conserved in their structure and gene arrangement. Their genome sizes ranged from 154,656 bp to 160,677 bp, and the size variation was found to be associated with differences in IR expansion and gene loss. A total of 112 to 114 genes were identified in the genomes, including 78 to 79 protein-coding genes, 30 tRNA genes, 4 rRNA genes, and 2 pseudogenes. Sequence divergence analysis indicated that non-coding regions exhibited a higher percentage of variable characters, and nine non-coding and six coding regions were identified as divergent hotspots. Phylogenetic results based on different datasets showed that this tribe was divided into three reciprocally exclusive groups. The phylogenetic analyses between plastome and nuclear ITS data were partly incongruent with each other. This study provides new insights into plastome evolution of Ruteae as well as Rutaceae. The availability of these plastomes provides useful genomic resources for molecular DNA barcodes and phylogenetically informative markers and deepens our understanding of the phylogeny in Ruteae.

Complete Chloroplast Genomes of 14 Subspecies of D. glomerata: Phylogenetic and Comparative Genomic Analyses

Orchardgrass (Dactylis glomerata L.) is a species in the Gramineae family that is highly important economically and valued for its role in ecology. However, the phylogeny and taxonomy of D. glomerata are still controversial based on current morphological and molecular evidence. The study of chloroplast (cp) genomes has developed into a powerful tool to develop molecular markers for related species and reveal the relationships between plant evolution and phylogenetics. In this study, we conducted comparative genomic analyses and phylogenetic inferences on 14 cp genomes of D. glomerata originating from the Mediterranean and Eurasia. The genome size ranged from 134,375 bp to 134,993 bp and exhibited synteny of gene organization and order. A total of 129–131 genes were identified, including 85–87 protein coding genes, 38 tRNA genes and 8 rRNA genes. The cp sequences were highly conserved, and key sequence variations were detected at the junctions of inverted repeats (IRs)/small single–copy (SSC) regions. Moreover, nine highly variable regions were identified among the subspecies based on a sequence divergence analysis. A total of 285 RNA editing sites were detected that were relevant to 52 genes, where rpoB exhibited the most abundant RNA editing sites. The phylogenetic analysis revealed that all Dactylis subspecies clustered into a monophyletic group and most branches provided a high support bootstrap. The main divergence time of D. glomerata was dated to the Miocene era, and this could have been due to changes in the climate. These findings will provide useful insights for further studies on phylogeny, the identification of subspecies and the development of hypotheses for the evolutionary history of the genus Dactylis and of the Gramineae family.

Insights Into Comparative Analyses and Phylogenomic Implications of Acer (Sapindaceae) Inferred From Complete Chloroplast Genomes

Acer L. (Sapindaceae) is one of the most diverse and widespread plant genera in the Northern Hemisphere. It comprises 124-156 recognized species, with approximately half being native to Asia. Owing to its numerous morphological features and hybridization, this genus is taxonomically and phylogenetically ranked as one of the most challenging plant taxa. Here, we report the complete chloroplast genome sequences of five Acer species and compare them with those of 43 published Acer species. The chloroplast genomes were 149,103-158,458 bp in length. We conducted a sliding window analysis to find three relatively highly variable regions (psbN-rps14, rpl32-trnL, and ycf1) with a high potential for developing practical genetic markers. A total of 76-103 SSR loci were identified in 48 Acer species. The positive selection analysis of Acer species chloroplast genes showed that two genes (psaI and psbK) were positively selected, implying that light level is a selection pressure for Acer species. Using Bayes empirical Bayes methods, we also identified that 20 cp gene sites have undergone positive selection, which might result from adaptation to specific ecological niches. In phylogenetic analysis, we have reconfirmed that Acer pictum subsp. mono and A. truncatum as sister species. Our results strongly support the sister relationships between sections Platanoidea and Macrantha and between sections Trifoliata and Pentaphylla. Moreover, series Glabra and Arguta are proposed to promote to the section level. The chloroplast genomic resources provided in this study assist taxonomic and phylogenomic resolution within Acer and the Sapindaceae family.

Quantitative modelling of fine-scale variations in the Arabidopsis thaliana crossover landscape

In, essentially, all species where meiotic crossovers (COs) have been studied, they occur preferentially in open chromatin, typically near gene promoters and to a lesser extent, at the end of genes. Here, in the case of Arabidopsis thaliana, we unveil further trends arising when one considers contextual information, namely summarised epigenetic status, gene or intergenic region size, and degree of divergence between homologs. For instance, we find that intergenic recombination rate is reduced if those regions are less than 1.5 kb in size. Furthermore, we propose that the presence of single nucleotide polymorphisms enhances the rate of CO formation compared to when homologous sequences are identical, in agreement with previous works comparing rates in adjacent homozygous and heterozygous blocks. Lastly, by integrating these different effects, we produce a quantitative and predictive model of the recombination landscape that reproduces much of the experimental variation.

Large X-Linked Palindromes Undergo Arm-to-Arm Gene Conversion across Mus Lineages

Large (>10 kb), nearly identical (>99% nucleotide identity), palindromic sequences are enriched on mammalian sex chromosomes. Primate Y-palindromes undergo high rates of arm-to-arm gene conversion, a proposed mechanism for maintaining their sequence integrity in the absence of X-Y recombination. It is unclear whether X-palindromes, which can freely recombine in females, undergo arm-to-arm gene conversion and, if so, at what rate. We generated high-quality sequence assemblies of Mus molossinus and M. spretus X-palindromic regions and compared them with orthologous M. musculus X-palindromes. Our evolutionary sequence comparisons find evidence of X-palindrome arm-to-arm gene conversion at rates comparable to autosomal allelic gene conversion rates in mice. Mus X-palindromes also carry more derived than ancestral variants between species, suggesting that their sequence is rapidly diverging. We speculate that in addition to maintaining genes' sequence integrity via sequence homogenization, palindrome arm-to-arm gene conversion may also facilitate rapid sequence divergence.

Resolving the Early Divergence Pattern of Teleost Fish Using Genome-Scale Data

Regarding the phylogenetic relationship of the three primary groups of teleost fishes, Osteoglossomorpha (bonytongues and others), Elopomorpha (eels and relatives), Clupeocephala (the remaining teleost fish), early morphological studies hypothesized the first divergence of Osteoglossomorpha, whereas the recent prevailing view is the first divergence of Elopomorpha. Molecular studies supported all the possible relationships of the three primary groups. This study analyzed genome-scale data from four previous studies: 1) 412 genes from 12 species, 2) 772 genes from 15 species, 3) 1,062 genes from 30 species, and 4) 491 UCE loci from 27 species. The effects of the species, loci, and models used on the constructed tree topologies were investigated. In the analyses of the data sets (1)–(3), although the first divergence of Clupeocephala that left the other two groups in a sister relationship was supported by concatenated sequences and gene trees of all the species and genes, the first divergence of Elopomorpha among the three groups was supported using species and/or genes with low divergence of sequence and amino-acid frequencies. This result corresponded to that of the UCE data set (4), whose sequence divergence was low, which supported the first divergence of Elopomorpha with high statistical significance. The increase in accuracy of the phylogenetic construction by using species and genes with low sequence divergence was predicted by a phylogenetic informativeness approach and confirmed by computer simulation. These results supported that Elopomorpha was the first basal group of teleost fish to have diverged, consistent with the prevailing view of recent morphological studies.

Functional Conservation of Divergent p63-Bound cis-Regulatory Elements

The transcription factor p63 is an essential regulator of vertebrate ectoderm development, including epidermis, limbs, and craniofacial tissues. Here, we have investigated the evolutionary conservation of p63 binding sites (BSs) between zebrafish and human. First, we have analyzed sequence conservation of p63 BSs by comparing ChIP-seq data from human keratinocytes and zebrafish embryos, observing a very poor conservation. Next, we compared the gene regulatory network orchestrated by p63 in both species and found a high overlap between them, suggesting a high degree of functional conservation during evolution despite sequence divergence and the large evolutionary distance. Finally, we used transgenic reporter assays in zebrafish embryos to functionally validate a set of equivalent p63 BSs from zebrafish and human located close to genes involved in epidermal development. Reporter expression was driven by human and zebrafish BSs to many common tissues related to p63 expression domains. Therefore, we conclude that the gene regulatory network controlled by p63 is highly conserved across vertebrates despite the fact that p63-bound regulatory elements show high divergence.

How new genes are born

Analysis of yeast, fly and human genomes suggests that sequence divergence is not the main source of orphan genes.

Synteny-based analyses indicate that sequence divergence is not the main source of orphan genes

The origin of 'orphan' genes, species-specific sequences that lack detectable homologues, has remained mysterious since the dawn of the genomic era. There are two dominant explanations for orphan genes: complete sequence divergence from ancestral genes, such that homologues are not readily detectable; and de novo emergence from ancestral non-genic sequences, such that homologues genuinely do not exist. The relative contribution of the two processes remains unknown. Here, we harness the special circumstance of conserved synteny to estimate the contribution of complete divergence to the pool of orphan genes. By separately comparing yeast, fly and human genes to related taxa using conservative criteria, we find that complete divergence accounts, on average, for at most a third of eukaryotic orphan and taxonomically restricted genes. We observe that complete divergence occurs at a stable rate within a phylum but at different rates between phyla, and is frequently associated with gene shortening akin to pseudogenization.

The Complete Chloroplast Genomes of Echinacanthus Species (Acanthaceae): Phylogenetic Relationships, Adaptive Evolution, and Screening of Molecular Markers

Among the four species of Echinacanthus (Acanthaceae), one distributed in the West Himalayan region and three restricted to the Sino-Vietnamese karst region. Because of its ecological significance, molecular markers are necessary for proper assessment of its genetic diversity and phylogenetic relationships. Herein, the complete chloroplast genomes of four Echinacanthus species were determined for the first time. The results indicated that all the chloroplast genomes were mapped as a circular structure and each genomes included 113 unique genes, of which 80 were protein-coding, 29 were tRNAs, and 4 were rRNAs. However, the four cp genomes ranged from 151,333 to 152,672 bp in length. Comparison of the four cp genomes showed that the divergence level was greater between geographic groups. We also analyzed IR expansion or contraction in the four cp genomes and the fifth type of the large single copy/inverted repeat region in Lamiales was suggested. Furthermore, based on the analyses of comparison and nucleotide variability, six most divergent sequences (rrn16, ycf1, ndhA, rps16-trnQ-UUG, trnS-GCU-trnG-UCC, and psaA-ycf3) were identified. A total of 37-45 simple sequence repeats were discovered in the four species and 22 SSRs were identified as candidate effective molecular markers for detecting interspecies polymorphisms. These SSRs and hotspot regions could be used as potential molecular markers for future study. Phylogenetic analysis based on Bayesian and parsimony methods did not support the monophyly of Echinacanthus. The phylogenetic relationships among the four species were clearly resolved and the results supported the recognition of the Sino-Vietnamese Echinacanthus species as a new genus. Based on the protein sequence evolution analysis, 12 genes (rpl14, rpl16, rps4, rps15, rps18, rps19, psbK, psbN, ndhC, ndhJ, rpoB, and infA) were detected under positive selection in branch of Sino-Vietnamese Echinacanthus species. These genes will lead to understanding the adaptation of Echinacanthus species to karst environment. The study will help to resolve the phylogenetic relationship and understand the adaptive evolution of Echinacanthus. It will also provide genomic resources and potential markers suitable for future species identification and speciation studies of the genus.

Species dilemma of musk deer (Moschus spp) in India: molecular data on cytochrome c oxidase I suggests distinct genetic lineage in Uttarakhand compared to other Moschus species

Musk deer are of high conservation priority owing to poaching pressure because of its musk pod. Representation of musk deer status using genetics is poorly documented in India, and it is not confirmed as to how many species of musk deer are present. We characterize for the first time, the genetic diversity of musk deer from Uttarakhand using Cytochrome Oxidase sub-unit (COI) gene (486 bp) and compared with the data available for other species. Results revealed the presence of six haplotypes in the Uttarakhand population amongst 17 sequences. Of these, 12 sequences shared the single haplotype. The intra-species sequences divergence was 0.003-0.017, whereas divergence with other species of musk deer was 0.071-0.081. Bayesian phylogenetic tree revealed that samples from Uttarakhand formed a separate clade with respect to other species of musk deer, whereas three species distributed in China clustered in the same clade and showed low sequences divergence, i.e., 0.002-0.061. Because of different ecomorph reported, we suggest using the barcoding based approach for inter and intra-species distinction and delineating species boundaries across the range for effective conservation. Besides, systematic classification, DNA barcoding would also help in dealing wildlife offence cases for disposal of the legal report in court.

Atomistic simulations indicate the functional loop-to-coiled-coil transition in influenza hemagglutinin is not downhill

Influenza hemagglutinin (HA) mediates viral entry into host cells through a large-scale conformational rearrangement at low pH that leads to fusion of the viral and endosomal membranes. Crystallographic and biochemical data suggest that a loop-to-coiled-coil transition of the B-loop region of HA is important for driving this structural rearrangement. However, the microscopic picture for this proposed "spring-loaded" movement is missing. In this study, we focus on understanding the transition of the B loop and perform a set of all-atom molecular dynamics simulations of the full B-loop trimeric structure with the CHARMM36 force field. The free-energy profile constructed from our simulations describes a B loop that stably folds half of the postfusion coiled coil in tens of microseconds, but the full coiled coil is unfavorable. A buried hydrophilic residue, Thr59, is implicated in destabilizing the coiled coil. Interestingly, this conserved threonine is the only residue in the B loop that strictly differentiates between the group 1 and 2 HA molecules. Microsecond-scale constant temperature simulations revealed that kinetic traps in the structural switch of the B loop can be caused by nonnative, intramonomer, or intermonomer β-sheets. The addition of the A helix stabilized the postfusion state of the B loop, but introduced the possibility for further β-sheet structures. Overall, our results do not support a description of the B loop in group 2 HAs as a stiff spring, but, rather, it allows for more structural heterogeneity in the placement of the fusion peptides during the fusion process.

DNA barcoding revealed a new species of Neolissochilus Rainboth, 1985 from the Kaladan River of Mizoram, North East India

Neolissochilus kaladanensis sp. nov., a new cyprinid species, is described from the Kaladan River drainage of Mizoram. It differs from all other valid Neolissochilus species in having higher number of gill rakers on the lower arm of the first gill arch (13-14 vs. 12 or below in all the species). The analysis of mitochondrial gene cytochrome c oxidase subunit I (COI) sequences separated N. kaladanensis sp. nov. from all other Neolissochilus and Tor species with an average genetic distance of 6.0%. It is further separated from the morphologically most similar species N. hendersoni and N. soroides by a genetic distance of 6.7% and 6.8%, respectively. Based on the lowest BIC and AICc scores, best fit model for COI dataset was TN93 + G + I, out of 24 different nucleotide substitution models tested. The maximum-likelihood (ML) phylogenetic tree was constructed using the COI sequences of representative Neolissochilus and Tor species. The anomalies observed among the GenBank sequences of the genera Tor and Neolissochilus are also discussed.

The Complete Chloroplast Genome of Heimia myrtifolia and Comparative Analysis within Myrtales

Heimia myrtifolia is an important medicinal plant with several pharmacologically active alkaloids and is also used as an ornamental landscape plant. The purpose of this study is to complete and characterize the chloroplast (cp) genome of H. myrtifolia and compare genomic features to other Myrtales species' cp genomes. The analysis showed that H. myrtifolia has a total length of 159,219 bp with a typical quadripartite structure containing two identical inverted repeats (IRs) of 25,643 bp isolated by one large single copy (LSC) of 88,571 bp and one small single copy (SSC) of 18,822 bp. The H. myrtifolia cp genome contains 129 genes with eight ribosomal RNAs, 30 transfer RNAs, and 78 protein coding genes, in which 17 genes are duplicated in two IR regions. The genome organization including gene type and number and guanine-cytosine (GC) content is analyzed among the 12 cp genomes in this study. Approximately 255 simple sequence repeats (SSRs) and 16 forward, two reverses, and two palindromic repeats were identified in the H. myrtifolia cp genome. By comparing the whole H. myrtifolia cp genome with 11 other Myrtales species, the results showed that the sequence similarity was high between coding regions while sequence divergence was high between intergenic regions. By employing the full cp genomes for phylogenetic analysis, structural and sequence differences were characterized between H. myrtifolia and 11 Myrtales species illustrating what patterns are common in the evolution of cp genomes within the Myrtales. The first entire cp genome in the genus Heimia provides a valuable resource for further studies in these medicinally and ornamentally important taxa.

Comparative Plastid Genomes of Primula Species: Sequence Divergence and Phylogenetic Relationships

Compared to traditional DNA markers, genome-scale datasets can provide mass information to effectively address historically difficult phylogenies. Primula is the largest genus in the family Primulaceae, with members distributed mainly throughout temperate and arctic areas of the Northern Hemisphere. The phylogenetic relationships among Primula taxa still maintain unresolved, mainly due to intra- and interspecific morphological variation, which was caused by frequent hybridization and introgression. In this study, we sequenced and assembled four complete plastid genomes (Primula handeliana, Primula woodwardii, Primula knuthiana, and Androsace laxa) by Illumina paired-end sequencing. A total of 10 Primula species (including 7 published plastid genomes) were analyzed to investigate the plastid genome sequence divergence and their inferences for the phylogeny of Primula. The 10 Primula plastid genomes were similar in terms of their gene content and order, GC content, and codon usage, but slightly different in the number of the repeat. Moderate sequence divergence was observed among Primula plastid genomes. Phylogenetic analysis strongly supported that Primula was monophyletic and more closely related to Androsace in the Primulaceae family. The phylogenetic relationships among the 10 Primula species showed that the placement of P. knuthiana–P. veris clade was uncertain in the phylogenetic tree. This study indicated that plastid genome data were highly effective to investigate the phylogeny.

Different Natural Selection Pressures on the atpF Gene in Evergreen Sclerophyllous and Deciduous Oak Species: Evidence from Comparative Analysis of the Complete Chloroplast Genome of Quercus aquifolioides with Other Oak Species

Quercus is an economically important and phylogenetically complex genus in the family Fagaceae. Due to extensive hybridization and introgression, it is considered to be one of the most challenging plant taxa, both taxonomically and phylogenetically. Quercus aquifolioides is an evergreen sclerophyllous oak species that is endemic to, but widely distributed across, the Hengduanshan Biodiversity Hotspot in the Eastern Himalayas. Here, we compared the fully assembled chloroplast (cp) genome of Q. aquifolioides with those of three closely related species. The analysis revealed a cp genome ranging in size from 160,415 to 161,304 bp and with a typical quadripartite structure, composed of two inverted repeats (IRs) separated by a small single copy (SSC) and a large single copy (LSC) region. The genome organization, gene number, gene order, and GC content of these four Quercus cp genomes are similar to those of many angiosperm cp genomes. We also analyzed the Q. aquifolioides repeats and microsatellites. Investigating the effects of selection events on shared protein-coding genes using the Ka/Ks ratio showed that significant positive selection had acted on the atpF gene of Q. aquifolioides compared to two deciduous oak species, and that there had been significant purifying selection on the atpF gene in the chloroplast of evergreen sclerophyllous oak trees. In addition, site-specific selection analysis identified positively selected sites in 12 genes. Phylogenetic analysis based on shared protein-coding genes from 14 species defined Q. aquifolioides as belonging to sect. Heterobalanus and being closely related to Q. rubra and Q. aliena. Our findings provide valuable genetic information for use in accurately identifying species, resolving taxonomy, and reconstructing the phylogeny of the genus Quercus.

Saurogobio punctatus sp. nov., a new cyprinid gudgeon (Teleostei: Cypriniformes) from the Yangtze River, based on both morphological and molecular data

A new cyprinid gudgeon, Saurogobio punctatus sp. nov., is described based on specimens collected from the Yangtze River, China. The new species can be distinguished from its congeners by differences in both morphology and the cytochrome b (cytb) gene sequence. Numerous minute blackish spots are scattered on dorsal and caudal fins in S. punctatus sp. nov. v. absent in the other seven valid Saurogobio species. The new species can be further distinguished from its congeners by the following unique combination of characters: a dorsal fin with eight branched rays; absence of scales in chest area before pectoral origin; upper and lower lips thick, covered with papillae; and a papillose mental pad approximately triangular. Morphologically, the new species most resembles the Chinese lizard gudgeon Saurogobio dabryi, but the new species lays yellowish adhesive eggs v. white pelagic eggs in S. dabryi. A phylogenetic analysis of all Saurogobio species based on cytb gene sequences indicated that S. punctatus sp. nov was distinctly separated from its congeners, with mean sequence divergence ranging from 12·6 to 21·0%. Therefore, molecular data further supported the distinctiveness of the new species.

Development of Chloroplast Genomic Resources for Oryza Species Discrimination

Rice is the most important crop in the world as the staple food for over half of the population. The wild species of Oryza represent an enormous gene pool for genetic improvement of rice cultivars. Accurate and rapid identification of these species is critical for effective utilization of the wild rice germplasm. In this study, we developed valuable chloroplast molecular markers by comparing the chloroplast genomes for species identification. Four chloroplast genomes of Oryza were newly sequenced on the Illumina HiSeq platform and other 14 Oryza species chloroplast genomes from Genbank were simultaneously taken into consideration for comparative analyses. Among 18 Oryza chloroplast genomes, five variable regions (rps16-trnQ, trnTEYD, psbE-petL, rpoC2 and rbcL-accD) were detected for DNA barcodes, in addition to differences in simple sequence repeats (SSR) and repeat sequences. The highest species resolution (72.22%) was provided by rpoC2 and rbcL-accD with distance-based methods. Three-marker combinations (rps16-trnQ + trnTEYD + rbcL-accD, rps16-trnQ + trnTEYD + rpoC2 and rpoC2 + trnTEYD + psbE-petL) showed the best species resolution (100%). Phylogenetic analysis based on the chloroplast genome provided the best resolution of Oryza. In the comparison of chloroplast genomes in this study, identification of the most variable regions and assessment of the focal regions of divergence were efficient in developing species-specific DNA barcodes. Based on evaluation of the chloroplast genomic resources, we conclude that chloroplast genome sequences are a reliable and valuable molecular marker for exploring the wild rice genetic resource in rice improvement.

The Complete Chloroplast Genome of Wild Rice (Oryza minuta) and Its Comparison to Related Species

Oryza minuta, a tetraploid wild relative of cultivated rice (family Poaceae), possesses a BBCC genome and contains genes that confer resistance to bacterial blight (BB) and white-backed (WBPH) and brown (BPH) plant hoppers. Based on the importance of this wild species, this study aimed to understand the phylogenetic relationships of O. minuta with other Oryza species through an in-depth analysis of the composition and diversity of the chloroplast (cp) genome. The analysis revealed a cp genome size of 135,094 bp with a typical quadripartite structure and consisting of a pair of inverted repeats separated by small and large single copies, 139 representative genes, and 419 randomly distributed microsatellites. The genomic organization, gene order, GC content and codon usage are similar to those of typical angiosperm cp genomes. Approximately 30 forward, 28 tandem and 20 palindromic repeats were detected in the O. minuta cp genome. Comparison of the complete O. minuta cp genome with another eleven Oryza species showed a high degree of sequence similarity and relatively high divergence of intergenic spacers. Phylogenetic analyses were conducted based on the complete genome sequence, 65 shared genes and matK gene showed same topologies and O. minuta forms a single clade with parental O. punctata. Thus, the complete O. minuta cp genome provides interesting insights and valuable information that can be used to identify related species and reconstruct its phylogeny.

Comparative Analysis of Six Lagerstroemia Complete Chloroplast Genomes

Crape myrtles are economically important ornamental trees of the genus Lagerstroemia L. (Lythraceae), with a distribution from tropical to northern temperate zones. They are positioned phylogenetically to a large subclade of rosids (in the eudicots) which contain more than 25% of all the angiosperms. They commonly bloom from summer till fall and are of significant value in city landscape and environmental protection. Morphological traits are shared inter-specifically among plants of Lagerstroemia to certain extent and are also influenced by environmental conditions and different developmental stages. Thus, classification of plants in Lagerstroemia at species and cultivar levels is still a challenging task. Chloroplast (cp) genome sequences have been proven to be an informative and valuable source of cp DNA markers for genetic diversity evaluation. In this study, the complete cp genomes of three Lagerstroemia species were newly sequenced, and three other published cp genome sequences of Lagerstroemia were retrieved for comparative analyses in order to obtain an upgraded understanding of the application value of genetic information from the cp genomes. The six cp genomes ranged from 152,049 bp (L. subcostata) to 152,526 bp (L. speciosa) in length. We analyzed nucleotide substitutions, insertions/deletions, and simple sequence repeats in the cp genomes, and discovered 12 relatively highly variable regions that will potentially provide plastid markers for further taxonomic, phylogenetic, and population genetics studies in Lagerstroemia. The phylogenetic relationships of the Lagerstroemia taxa inferred from the datasets from the cp genomes obtained high support, indicating that cp genome data may be useful in resolving relationships in this genus.

Molecular systematics and DNA barcoding of Altai osmans, oreoleuciscus (pisces, cyprinidae, and leuciscinae), and their nearest relatives, inferred from sequences of cytochrome b (Cyt-b), cytochrome oxidase c (Co-1), and complete mitochondrial genome

Mitochondrial DNA (mtDNA) at the protein-coding Cyt-b gene along with data retrieved from GenBank for Co-1 gene fragments and complete mitochondrial genome (mitogenome) of Altai osmans and the nearest relatives of Leuciscinae fish species were compared for the estimation of variability and phylogenetic tree building. Phylogenetic trees were built by four techniques: Bayesian (BA), maximum likelihood (ML), maximum parsimony (MP), and neighbor-joining (NJ). Resolution of Cyt-b trees for species of two genera (Oreoleuciscus and Phoxinus) was quite distinct at all the approaches. For Tribolodon, the single gene trees were not well resolved; however, the mitogenome tree was resolved. Species identification on per individual basis (DNA barcoding) was high for both Cyt-b and Co-1 genes. The trees built using the data for 13 protein mitochondrial genes revealed a complicated phylogenetic pattern within the subfamily Leuciscinae. Scores of the average p-distances at three taxonomic levels were considerably different: (1) 1.16 ± 0.96, (2) 8.21 ± 1.01, and (3) 16.41 ± 0.85 for Cyt-b and (1) 1.04 ± 0.78, (2) 8.30 ± 0.92, and (3) 10.74 ± 0.79 for 13 protein genes of mitogenome, where (1) is intraspecies, (2) is intragenus, and (3) is intrasubfamily levels. Data on mitogenome distances were summarized for the taxonomic hierarchy for the first time. A concordant increase in distance score with growth of the rank of taxa (having the minimum score at the intraspecies level), both for a single gene and the whole mitogenome, substantiates the concept that speciation in the subfamily Leuciscinae in most cases follows the geographic mode. The distinct clustering of Altai osmans, Oreoleuciscus potanini and O. humilis, in the Cyt-b and Co-1 gene trees with small overall genetic distances, obtained for both genes, allows us to consider these taxa as separate but genetically sister species.

Complete Chloroplast Genome of Nicotiana otophora and its Comparison with Related Species

Nicotiana otophora is a wild parental species of Nicotiana tabacum, an interspecific hybrid of Nicotiana tomentosiformis and Nicotiana sylvestris. However, N. otophora is least understood as an alternative paternal donor. Here, we compared the fully assembled chloroplast (cp) genome of N. otophora and with those of closely related species. The analysis showed a cp genome size of 156,073 bp and exhibited a typical quadripartite structure, which contains a pair of inverted repeats separated by small and large single copies, containing 163 representative genes, with 165 microsatellites distributed unevenly throughout the genome. Comparative analysis of a gene with known function across Nicotiana species revealed 76 protein-coding sequences, 20 tRNA sequences, and 3 rRNA sequence shared between the cp genomes. The analysis revealed that N. otophora is a sister species to N. tomentosiformis within the Nicotiana genus, and Atropha belladonna and Datura stramonium are their closest relatives. These findings provide a valuable analysis of the complete N. otophora cp genome, which can identify species, elucidate taxonomy, and reconstruct the phylogeny of genus Nicotiana.

Sequence variability is correlated with weak immunogenicity in Streptococcus pyogenes M protein

The M protein of Streptococcus pyogenes, a major bacterial virulence factor, has an amino-terminal hypervariable region (HVR) that is a target for type-specific protective antibodies. Intriguingly, the HVR elicits a weak antibody response, indicating that it escapes host immunity by two mechanisms, sequence variability and weak immunogenicity. However, the properties influencing the immunogenicity of regions in an M protein remain poorly understood. Here, we studied the antibody response to different regions of the classical M1 and M5 proteins, in which not only the HVR but also the adjacent fibrinogen-binding B repeat region exhibits extensive sequence divergence. Analysis of antisera from S. pyogenes-infected patients, infected mice, and immunized mice showed that both the HVR and the B repeat region elicited weak antibody responses, while the conserved carboxy-terminal part was immunodominant. Thus, we identified a correlation between sequence variability and weak immunogenicity for M protein regions. A potential explanation for the weak immunogenicity was provided by the demonstration that protease digestion selectively eliminated the HVR-B part from whole M protein-expressing bacteria. These data support a coherent model, in which the entire variable HVR-B part evades antibody attack, not only by sequence variability but also by weak immunogenicity resulting from protease attack.

Evidence of Statistical Inconsistency of Phylogenetic Methods in the Presence of Multiple Sequence Alignment Uncertainty

Evolutionary studies usually use a two-step process to investigate sequence data. Step one estimates a multiple sequence alignment (MSA) and step two applies phylogenetic methods to ask evolutionary questions of that MSA. Modern phylogenetic methods infer evolutionary parameters using maximum likelihood or Bayesian inference, mediated by a probabilistic substitution model that describes sequence change over a tree. The statistical properties of these methods mean that more data directly translates to an increased confidence in downstream results, providing the substitution model is adequate and the MSA is correct. Many studies have investigated the robustness of phylogenetic methods in the presence of substitution model misspecification, but few have examined the statistical properties of those methods when the MSA is unknown. This simulation study examines the statistical properties of the complete two-step process when inferring sequence divergence and the phylogenetic tree topology. Both nucleotide and amino acid analyses are negatively affected by the alignment step, both through inaccurate guide tree estimates and through overfitting to that guide tree. For many alignment tools these effects become more pronounced when additional sequences are added to the analysis. Nucleotide sequences are particularly susceptible, with MSA errors leading to statistical support for long-branch attraction artifacts, which are usually associated with gross substitution model misspecification. Amino acid MSAs are more robust, but do tend to arbitrarily resolve multifurcations in favor of the guide tree. No inference strategies produce consistently accurate estimates of divergence between sequences, although amino acid MSAs are again more accurate than their nucleotide counterparts. We conclude with some practical suggestions about how to limit the effect of MSA uncertainty on evolutionary inference.

Role of AtMSH7 in UV-B-induced DNA damage recognition and recombination

The mismatch repair (MMR) system maintains genome integrity by correcting replication-associated errors and inhibiting recombination between divergent DNA sequences. The basic features of the pathway have been highly conserved throughout evolution, although the nature and number of the proteins involved in this DNA repair system vary among organisms. Plants have an extra mismatch recognition protein, MutSγ, which is a heterodimer: MSH2-MSH7. To further understand the role of MSH7 in vivo, we present data from this protein in Arabidopsis thaliana. First, we generated transgenic plants that express β-glucuronidase (GUS) under the control of the MSH7 promoter. Histochemical staining of the transgenic plants indicated that MSH7 is preferentially expressed in proliferating tissues. Then, we identified msh7 T-DNA insertion mutants. Plants deficient in MSH7 show increased levels of UV-B-induced cyclobutane pyrimidine dimers relative to wild-type (WT) plants. Consistent with the patterns of MSH7 expression, we next analysed the role of the protein during somatic and meiotic recombination. The frequency of somatic recombination between homologous or homeologous repeats (divergence level of 1.6%) was monitored using a previously described GUS recombination reporter assay. Disruption of MSH7 has no effect on the rates of somatic homologous or homeologous recombination under control conditions or after UV-B exposure. However, the rate of meiotic recombination between two genetically linked seed-specific fluorescent markers was 97% higher in msh7 than in WT plants. Taken together, these results suggest that MSH7 is involved in UV-B-induced DNA damage recognition and in controlling meiotic recombination.

Genome-wide analysis reveals diverged patterns of codon bias, gene expression, and rates of sequence evolution in picea gene families

The recent sequencing of several gymnosperm genomes has greatly facilitated studying the evolution of their genes and gene families. In this study, we examine the evidence for expression-mediated selection in the first two fully sequenced representatives of the gymnosperm plant clade (Picea abies and Picea glauca). We use genome-wide estimates of gene expression (>50,000 expressed genes) to study the relationship between gene expression, codon bias, rates of sequence divergence, protein length, and gene duplication. We found that gene expression is correlated with rates of sequence divergence and codon bias, suggesting that natural selection is acting on Picea protein-coding genes for translational efficiency. Gene expression, rates of sequence divergence, and codon bias are correlated with the size of gene families, with large multicopy gene families having, on average, a lower expression level and breadth, lower codon bias, and higher rates of sequence divergence than single-copy gene families. Tissue-specific patterns of gene expression were more common in large gene families with large gene expression divergence than in single-copy families. Recent family expansions combined with large gene expression variation in paralogs and increased rates of sequence evolution suggest that some Picea gene families are rapidly evolving to cope with biotic and abiotic stress. Our study highlights the importance of gene expression and natural selection in shaping the evolution of protein-coding genes in Picea species, and sets the ground for further studies investigating the evolution of individual gene families in gymnosperms.

Molecular and morphological differentiation between Aphis gossypii Glover (Hemiptera, Aphididae) and related species, with particular reference to the North American Midwest

The cotton aphid, Aphis gossypii, is one of the most biologically diverse species of aphids; a polyphagous species in a family where most are host specialists. It is economically important and belongs to a group of closely related species that has challenged aphid taxonomy. The research presented here seeks to clarify the taxonomic relationships and status of species within the Aphid gossypii group in the North American Midwest. Sequences of the mitochondrial cytochrome oxidase 1 (COI), nuclear elongation factor 1-α (EF1-α), and nuclear sodium channel para-type (SCP) genes were used to differentiate between Aphid gossypii and related species. Aphis monardae, previously synonymised with Aphid gossypii, is re-established as a valid species. Phylogenetic analyses support the close relationship of members of the Aphid gossypii group native to North America (Aphid forbesi, Aphid monardae, Aphid oestlundi, Aphid rubifolii, and Aphid rubicola), Europe (Aphid nasturtii, Aphid urticata and Aphid sedi), and Asia (Aphid agrimoniae, Aphid clerodendri, Aphid glycines, Aphid gossypii, Aphid hypericiphaga, Aphid ichigicola, Aphid ichigo, Aphid sanguisorbicola, Aphid sumire and Aphid taraxicicola). The North American species most closely related to Aphid gossypii are Aphid monardae and Aphid oestlundi. The cosmopolitan Aphid gossypii and Aphid sedi identified in the USA are genetically very similar using COI and EF1-α sequences, but the SCP gene shows greater genetic distance between them. We present a discussion of the biological and morphological differentiation of these species.

Phylogeny based on two mtDNA genes (Co-1, Cyt-B) among Sculpins (Scorpaeniformes, Cottidae) and some other scorpionfish in the Russian Far East

The average scores of p-distances at Co-1 and Cyt-b genes revealed a pattern of increased nucleotide diversity for four categories: (1) intraspecies, (2) intragenus, (3) intrafamily, and (4) intraorder. Scores for Co-1 genes were following: (1) 0.24 ± 0.19% (n = 66), (2) 1.54 ± 0.27% (n = 33), (3) 11.57 ± 0.24% (n = 42), and (4) 12.91 ± 0.16% (n = 102); while for Cyt-b genes, those values were (1) 1.92 ± 0.55% (n = 24), (2) 8.19 ± 0.09% (n = 812), (3) 15.42 ± 0.07% (n = 1573), and (4) 23.13 ± 0.19% (n = 211). A concordant increase in the distance score, as the rank of taxa grows, with the minimum at intraspecies level substantiates the concept that speciation in the order Scorpaeniformes mostly follows a geographic mode. Phylogenetic trees for 40 sequences of scorpionfish and four other fishes belonging to Actinopterigii were developed using Co-1 gene and four different analytical approaches: Bayesian, ML, MP, and NJ. Resolution of the Co-1 trees for the three genera that included different species branches was quite distinct in all approaches. Intrafamily branching for Cyt-b in Cottidae was also supported by all four kinds of trees but only for less numerous genera in the research.

The first complete chloroplast genome of the Genistoid legume Lupinus luteus: evidence for a novel major lineage-specific rearrangement and new insights regarding plastome evolution in the legume family

BACKGROUND AND AIMS

To date chloroplast genomes are available only for members of the non-protein amino acid-accumulating clade (NPAAA) Papilionoid lineages in the legume family (i.e. Millettioids, Robinoids and the 'inverted repeat-lacking clade', IRLC). It is thus very important to sequence plastomes from other lineages in order to better understand the unusual evolution observed in this model flowering plant family. To this end, the plastome of a lupine species, Lupinus luteus, was sequenced to represent the Genistoid lineage, a noteworthy but poorly studied legume group.

METHODS

The plastome of L. luteus was reconstructed using Roche-454 and Illumina next-generation sequencing. Its structure, repetitive sequences, gene content and sequence divergence were compared with those of other Fabaceae plastomes. PCR screening and sequencing were performed in other allied legumes in order to determine the origin of a large inversion identified in L. luteus.

KEY RESULTS

The first sequenced Genistoid plastome (L. luteus: 155 894 bp) resulted in the discovery of a 36-kb inversion, embedded within the already known 50-kb inversion in the large single-copy (LSC) region of the Papilionoideae. This inversion occurs at the base or soon after the Genistoid emergence, and most probably resulted from a flip-flop recombination between identical 29-bp inverted repeats within two trnS genes. Comparative analyses of the chloroplast gene content of L. luteus vs. Fabaceae and extra-Fabales plastomes revealed the loss of the plastid rpl22 gene, and its functional relocation to the nucleus was verified using lupine transcriptomic data. An investigation into the evolutionary rate of coding and non-coding sequences among legume plastomes resulted in the identification of remarkably variable regions.

CONCLUSIONS

This study resulted in the discovery of a novel, major 36-kb inversion, specific to the Genistoids. Chloroplast mutational hotspots were also identified, which contain novel and potentially informative regions for molecular evolutionary studies at various taxonomic levels in the legumes. Taken together, the results provide new insights into the evolutionary landscape of the legume plastome.

Heterosis is prevalent among domesticated but not wild strains of Saccharomyces cerevisiae

Crosses between inbred but unrelated individuals often result in an increased fitness of the progeny. This phenomenon is known as heterosis and has been reported for wild and domesticated populations of plants and animals. Analysis of heterosis is often hindered by the fact that the genetic relatedness between analyzed organisms is only approximately known. We studied a collection of Saccharomyces cerevisiae isolates from wild and human-created habitats whose genomes were sequenced and thus their relatedness was fully known. We reasoned that if these strains accumulated different deleterious mutations at an approximately constant rate, then heterosis should be most visible in F1 heterozygotes from the least related parents. We found that heterosis was substantial and positively correlated with sequence divergence, but only in domesticated strains. More than 80% of the heterozygous hybrids were more fit than expected from the mean of their homozygous parents, and approximately three-quarters of those exceeded even the fittest parent. Our results support the notion that domestication brings about relaxation of selection and accumulation of deleterious mutations. However, other factors may have contributed as well. In particular, the observed build-up of genetic load might be facilitated by a decrease, and not increase, in the rate of inbreeding.

Four cases of Taenia saginata infection with an analysis of COX1 gene

Human taeniases had been not uncommon in the Republic of Korea (=Korea) until the 1980s. The prevalence decreased and a national survey in 2004 revealed no Taenia egg positive cases. However, a subsequent national survey in 2012 showed 0.04% (10 cases) prevalence of Taenia spp. eggs suggesting its resurgence in Korea. We recently encountered 4 cases of Taenia saginata infection who had symptoms of taeniasis that included discharge of proglottids. We obtained several proglottids from each case. Because the morphological features of T. saginata are almost indistinguishable from those of Taenia asiatica, molecular analyses using the PCR-RFLP and DNA sequencing of the cytochrome c oxidase subunit 1 (cox1) were performed to identify the species. The PCR-RFLP patterns of all of the 4 specimens were consistent with T. saginata, and the cox1 gene sequence showed 99.8-100% identity with that of T. saginata reported previously from Korea, Japan, China, and Cambodia. All of the 4 patients had the history of travel abroad but its relation with contracting taeniasis was unclear. Our findings may suggest resurgence of T. saginata infection among people in Korea.

Genome differentiation of Drosophila melanogaster from a microclimate contrast in Evolution Canyon, Israel

The opposite slopes of "Evolution Canyon" in Israel have served as a natural model system of adaptation to a microclimate contrast. Long-term studies of Drosophila melanogaster populations inhabiting the canyon have exhibited significant interslope divergence in thermal and drought stress resistance, candidate genes, mobile elements, habitat choice, mating discrimination, and wing-shape variation, all despite close physical proximity of the contrasting habitats, as well as substantial interslope migration. To examine patterns of genetic differentiation at the genome-wide level, we used high coverage sequencing of the flies' genomes. A total of 572 genes were significantly different in allele frequency between the slopes, 106 out of which were associated with 74 significantly overrepresented gene ontology (GO) terms, particularly so with response to stimulus and developmental and reproductive processes, thus corroborating previous observations of interslope divergence in stress response, life history, and mating functions. There were at least 37 chromosomal "islands" of interslope divergence and low sequence polymorphism, plausible signatures of selective sweeps, more abundant in flies derived from one (north-facing) of the slopes. Positive correlation between local recombination rate and the level of nucleotide polymorphism was also found.

Rapid degeneration of noncoding DNA regions surrounding SlAP3X/Y after recombination suppression in the dioecious plant Silene latifolia

Silene latifolia is a dioecious plant with heteromorphic XY sex chromosomes. Previous studies of sex chromosome–linked genes have suggested a gradual divergence between the X-linked and the Y-linked genes in proportion to the distance from the pseudoautosomal region. However, such a comparison has yet to be made for the noncoding regions. To better characterize the nonrecombining region of the X and Y chromosomes, we sequenced bacterial artificial chromosome clones containing the sex chromosome–linked paralogs SlAP3X and SlAP3Y, including 115 kb and 73 kb of sequences, respectively, flanking these genes. The synonymous nucleotide divergence between SlAP3X and SlAP3Y indicated that recombination stopped approximately 3.4 million years ago. Sequence homology analysis revealed the presence of six long terminal repeat retrotransposon-like elements. Using the nucleotide divergence calculated between left and right long terminal repeat sequences, insertion dates were estimated to be 0.083–1.6 million years ago, implying that all elements detected were inserted after recombination stopped. A reciprocal sequence homology search facilitated the identification of four homologous noncoding DNA regions between the X and Y chromosomes, spanning 6.7% and 10.6% of the X chromosome–derived and Y chromosome–derived sequences, respectively, investigated. Genomic Southern blotting and fluorescence in situ hybridization showed that the noncoding DNA flanking SlAP3X/Y has homology to many regions throughout the genome, regardless of whether they were homologous between the X and Y chromosomes. This finding suggests that most noncoding DNA regions rapidly lose their counterparts because of the introduction of transposable elements and indels (insertion–deletions) after recombination has stopped.

Divergent allele advantage at MHC-DRB through direct and maternal genotypic effects and its consequences for allele pool composition and mating

It is still debated whether main individual fitness differences in natural populations can be attributed to genome-wide effects or to particular loci of outstanding functional importance such as the major histocompatibility complex (MHC). In a long-term monitoring project on Galápagos sea lions (Zalophus wollebaeki), we collected comprehensive fitness and mating data for a total of 506 individuals. Controlling for genome-wide inbreeding, we find strong associations between the MHC locus and nearly all fitness traits. The effect was mainly attributable to MHC sequence divergence and could be decomposed into contributions of own and maternal genotypes. In consequence, the population seems to have evolved a pool of highly divergent alleles conveying near-optimal MHC divergence even by random mating. Our results demonstrate that a single locus can significantly contribute to fitness in the wild and provide conclusive evidence for the 'divergent allele advantage' hypothesis, a special form of balancing selection with interesting evolutionary implications.

Gene divergence of homeologous regions associated with a major seed protein content QTL in soybean

Understanding several modes of duplication contributing on the present genome structure is getting an attention because it could be related to numerous agronomically important traits. Since soybean serves as a rich protein source for animal feeds and human consumption, breeding efforts in soybean have been directed toward enhancing seed protein content. The publicly available soybean sequences and its genomically featured elements facilitate comprehending of quantitative trait loci (QTL) for seed protein content in concordance with homeologous regions in soybean genome. Although parts of chromosome (Chr) 20 and Chr 10 showed synteny, QTLs for seed protein content present only on Chr 20. Using comparative analysis of gene contents in recently duplicated genomic regions harboring QTL for protein/oil content on Chrs 20 and 10, a total of 27 genes are present in duplicated regions of both Chrs. Notably, 4 tandem duplicates of the putative homeobox protein 22 (HB22) are present only on Chr 20 and this Medicago truncatula homolog expressed in endosperm at seed filling stage. These tandem duplicates could contribute on the protein/oil QTL of Chr 20. Our study suggests that non-shared gene contents within the duplicated genomic regions might lead to absence/presence of QTL related to protein/oil content.

Social discrimination by quantitative assessment of immunogenetic similarity

Genes of the major histocompatibility complex (MHC) that underlie the adaptive immune system may allow vertebrates to recognize their kin. True kin-recognition genes should produce signalling products to which organisms can respond. Allelic variation in the peptide-binding region (PBR) of MHC molecules determines the pool of peptides that can be presented to trigger an immune response. To examine whether these MHC peptides also might underlie assessments of genetic similarity, we tested whether Xenopus laevis tadpoles socially discriminate between pairs of siblings with which they differed in PBR amino acid sequences. We found that tadpoles (four sibships, n = 854) associated preferentially with siblings with which they were more similar in PBR amino acid sequence. Moreover, the strength of their preference for a conspecific was directly proportional to the sequence similarity between them. Discrimination was graded, and correlated more closely with functional sequence differences encoded by MHC class I and class II alleles than with numbers of shared haplotypes. Our results thus suggest that haplotype analyses may fail to reveal fine-scale behavioural responses to divergence in functionally expressed sequences. We conclude that MHC-PBR gene products mediate quantitative social assessment of immunogenetic similarity that may facilitate kin recognition in vertebrates.

Coadaptation of isoacceptor tRNA genes and codon usage bias for translation efficiency in Aedes aegypti and Anopheles gambiae

The transfer RNAs (tRNAs) are essential components of translational machinery. We determined that tRNA isoacceptors (tRNAs with different anticodons but incorporating the same amino acid in protein synthesis) show differential copy number abundance, genomic distribution patterns and sequence evolution between Aedes aegypti and Anopheles gambiae mosquitoes. The tRNA-Ala genes are present in unusually high copy number in the Ae. aegypti genome but not in An. gambiae. Many of the tRNA-Ala genes of Ae. aegypti are flanked by a highly conserved sequence that is not observed in An. gambiae. The relative abundance of tRNA isoacceptor genes is correlated with preferred (or optimal) and nonpreferred (or rare) codons for ∼2-4% of the predicted protein coding genes in both species. The majority (∼74-85%) of these genes are related to pathways involved with translation, energy metabolism and carbohydrate metabolism. Our results suggest that these genes and the related pathways may be under translational selection in these mosquitoes.

USING POPULATION GENOMICS TO DETECT SELECTION IN NATURAL POPULATIONS: KEY CONCEPTS AND METHODOLOGICAL CONSIDERATIONS

Natural selection shapes patterns of genetic variation among individuals, populations, and species, and it does so differentially across genomes. The field of population genomics provides a comprehensive genome-scale view of the action of selection, even beyond traditional model organisms. However, even with nearly complete genomic sequence information, our ability to detect the signature of selection on specific genomic regions depends on choosing experimental and analytical tools appropriate to the biological situation. For example, processes that occur at different timescales, such as sorting of standing genetic variation, mutation-selection balance, or fixed interspecific divergence, have different consequences for genomic patterns of variation. Inappropriate experimental or analytical approaches may fail to detect even strong selection or falsely identify a signature of selection. Here we outline the conceptual framework of population genomics, relate genomic patterns of variation to evolutionary processes, and identify major biological factors to be considered in studies of selection. As data-gathering technology continues to advance, our ability to understand selection in natural populations will be limited more by conceptual and analytical weaknesses than by the amount of molecular data. Our aim is to bring critical biological considerations to the fore in population genomics research and to spur the development and application of analytical tools appropriate to diverse biological systems.

Comparative analysis of nuclear tRNA genes of Nasonia vitripennis and other arthropods, and relationships to codon usage bias

Using bioinformatics methods, we identified a total of 221 and 199 tRNA genes in the nuclear genomes of Nasonia vitripennis and honey bee (Apis mellifera), respectively. We performed comparative analyses of Nasonia tRNA genes with honey bee and other selected insects to understand genomic distribution, sequence evolution and relationship of tRNA copy number with codon usage patterns. Many tRNA genes are located physically close to each other in the form of small clusters in the Nasonia genome. However, the number of clusters and the tRNA genes that form such clusters vary from species to species. In particular, the Ala-, Pro-, Tyr- and His-tRNA genes tend to accumulate in clusters in Nasonia but not in honey bee, whereas the bee contains a long cluster of 15 tRNA genes (of which 13 are Gln-tRNAs) that is absent in Nasonia. Though tRNA genes are highly conserved, contrasting patterns of nucleotide diversity are observed among the arm and loop regions of tRNAs between Nasonia and honey bee. Also, the sequence convergence between the reconstructed ancestral tRNAs and the present day tRNAs suggests a common ancestral origin of Nasonia and honey bee tRNAs. Furthermore, we also present evidence that the copy number of isoacceptor tRNAs (those having a different anticodon but charge the same amino acid) is correlated with codon usage patterns of highly expressed genes in Nasonia.

Restriction of ectopic recombination by interhomolog interactions during Saccharomyces cerevisiae meiosis

In Saccharomyces cerevisiae meiosis, recombination occurs frequently between sequences at the same location on homologs (allelic recombination) and can take place between dispersed homologous sequences (ectopic recombination). Ectopic recombination occurs less often than does allelic, especially when homologous sequences are on heterologous chromosomes. To account for this, it has been suggested that homolog pairing (homolog colocalization and alignment) either promotes allelic recombination or restricts ectopic recombination. The latter suggestion was tested by examining ectopic recombination in two cases where normal interhomolog relationships are disrupted. In the first case, one member of a homolog pair was replaced by a homologous (related but not identical) chromosome that has diverged sufficiently to prevent allelic recombination. In the second case, ndj1 mutants were used to delay homolog pairing and synapsis. Both circumstances resulted in a substantial increase in the frequency of ectopic recombination between arg4-containing plasmid inserts located on heterologous chromosomes. These findings suggest that, during normal yeast meiosis, progressive homolog colocalization, alignment, synapsis, and allelic recombination restrict the ability of ectopically located sequences to find each other and recombine. In the absence of such restrictions, the meiotic homology search may encompass the entire genome.