Bioinformatic Amplicon Read Processing Strategies Strongly Affect Eukaryotic Diversity and the Taxonomic Composition of Communities.

Amplicon read sequencing has revolutionized the field of microbial diversity studies. The technique has been developed for bacterial assemblages and has undergone rigorous testing with mock communities. However, due to the great complexity of eukaryotes and the numbers of different rDNA copies, analyzing eukaryotic diversity is more demanding than analyzing bacterial or mock communities, so studies are needed that test the methods of analyses on taxonomically diverse natural communities. In this study, we used 20 samples collected from the Baltic Sea ice, slush and under-ice water to investigate three program packages (UPARSE, mothur and QIIME) and 18 different bioinformatic strategies implemented in them. Our aim was to assess the impact of the initial steps of bioinformatic strategies on the results when analyzing natural eukaryotic communities. We found significant differences among the strategies in resulting read length, number of OTUs and estimates of diversity as well as clear differences in the taxonomic composition of communities. The differences arose mainly because of the variable number of chimeric reads that passed the pre-processing steps. Singleton removal and denoising substantially lowered the number of errors. Our study showed that the initial steps of the bioinformatic amplicon read processing strategies require careful consideration before applying them to eukaryotic communities.

Full likelihood analysis of genetic risk with variable age at onset disease--combining population-based registry data and demographic information.

BACKGROUND: In genetic studies of rare complex diseases it is common to ascertain familial data from population based registries through all incident cases diagnosed during a pre-defined enrollment period. Such an ascertainment procedure is typically taken into account in the statistical analysis of the familial data by constructing either a retrospective or prospective likelihood expression, which conditions on the ascertainment event. Both of these approaches lead to a substantial loss of valuable data. METHODOLOGY AND FINDINGS: Here we consider instead the possibilities provided by a Bayesian approach to risk analysis, which also incorporates the ascertainment procedure and reference information concerning the genetic composition of the target population to the considered statistical model. Furthermore, the proposed Bayesian hierarchical survival model does not require the considered genotype or haplotype effects be expressed as functions of corresponding allelic effects. Our modeling strategy is illustrated by a risk analysis of type 1 diabetes mellitus (T1D) in the Finnish population-based on the HLA-A, HLA-B and DRB1 human leucocyte antigen (HLA) information available for both ascertained sibships and a large number of unrelated individuals from the Finnish bone marrow donor registry. The heterozygous genotype DR3/DR4 at the DRB1 locus was associated with the lowest predictive probability of T1D free survival to the age of 15, the estimate being 0.936 (0.926; 0.945 95% credible interval) compared to the average population T1D free survival probability of 0.995. SIGNIFICANCE: The proposed statistical method can be modified to other population-based family data ascertained from a disease registry provided that the ascertainment process is well documented, and that external information concerning the sizes of birth cohorts and a suitable reference sample are available. We confirm the earlier findings from the same data concerning the HLA-DR3/4 related risks for T1D, and also provide here estimated predictive probabilities of disease free survival as a function of age.

Evolution of the capsular regulatory genes in Streptococcus pneumoniae.

The major pneumococcal virulence determinant is its capsule, and pneumococcal epidemiology is based on 91 capsular serotypes, each corresponding to the structure of the capsular polysaccharide determined by the type-specific capsular genome. Here, we provide the beginnings of an approach to intertwine serotype epidemiology, capsular regulatory gene characteristics on the basis of existing sequence information, and the reanalysis of published epidemiological data. We present an approach to explain epidemiological characteristics of serotypes on the basis of genetic differences in their capsular regulatory genes. The part of the capsular genome that regulates capsular expression falls into 2 highly divergent sequence clans: the ancestral pneumococcal capsular regulatory gene sequences (present in 49 serotypes) and laterally transferred sequences (present in 32 serotypes). Our survey of epidemiological data showed a tendency of the ancestral type of the capsular regulatory genome to be associated with carriage and the laterally transferred sequences to be associated with invasive disease isolates. The regulatory gene region shows mosaic structures that have signatures of recent recombination events, reminiscent of structures known from antibiotic resistance genes.

Phylogeny of Mysis (Crustacea, Mysida): history of continental invasions inferred from molecular and morphological data.

We studied the phylogenetic history of opossum shrimps of the genus Mysis Latreille, 1802 (Crustacea: Mysida) using parsimony analyses of morphological characters, DNA sequence data from mitochondrial (16S, COI and CytB) and nuclear genes (ITS2, 18S), and eight allozyme loci. With these data we aimed to resolve a long-debated question of the origin of the non-marine (continental) taxa in the genus, i.e., "glacial relicts" in circumpolar postglacial lakes and "arctic immigrants" in the Caspian Sea. A simultaneous analysis of the data sets gave a single tree supporting monophyly of all continental species, as well as monophyly of the taxa from circumpolar lakes and from the Caspian Sea. A clade of three circumarctic marine species was sister group to the continental taxa, whereas Atlantic species had more distant relationships to the others. Small molecular differentiation among the morphologically diverse endemic species from the Caspian Sea suggested their recent speciation, while the phenotypically more uniform "glacial relict" species from circumpolar lakes (Mysis relicta group) showed deep molecular divergences. For the length-variable ITS2 region both direct optimization and a priori alignment procedures gave similar topologies, although the former approach provided a better overall resolution. In terms of partitioned Bremer support (PBS), mitochondrial protein coding genes provided the largest contribution (83%) to the total tree resolution. This estimate however, appears to be partly spurious, due to the concerted inheritance of mitochondrial characters and probable cases of introgression or ancestral polymorphism.