Genetic determinants of complement activation in the general population.

Complement is a fundamental innate immune response component. Its alterations are associated with severe systemic diseases. To illuminate the complement's genetic underpinnings, we conduct genome-wide association studies of the functional activity of the classical (CP), lectin (LP), and alternative (AP) complement pathways in the Cooperative Health Research in South Tyrol study (n = 4,990). We identify seven loci, encompassing 13 independent, pathway-specific variants located in or near complement genes (CFHR4, C7, C2, MBL2) and non-complement genes (PDE3A, TNXB, ABO), explaining up to 74% of complement pathways' genetic heritability and implicating long-range haplotypes associated with LP at MBL2. Two-sample Mendelian randomization analyses, supported by transcriptome- and proteome-wide colocalization, confirm known causal pathways, establish within-complement feedback loops, and implicate causality of ABO on LP and of CFHR2 and C7 on AP. LP causally influences collectin-11 and KAAG1 levels and the risk of mouth ulcers. These results build a comprehensive resource to investigate the role of complement in human health.

Relative Telomere Length Is Associated With Age-Related Macular Degeneration in Women.

Purpose: Relative telomere length (RTL) is a biomarker for physiological aging. Premature shortening of telomeres is associated with oxidative stress, which is one possible pathway that might contribute to age-related macular degeneration (AMD). We therefore aimed to investigate the association between RTL and AMD in a well-characterized group of elderly individuals. Methods: We measured RTL in participants of the AugUR study using a multiplex quantitative PCR-based assay determining the ratio between the telomere product and a single-copy gene product (T/S ratio). AMD was assessed by manual grading of color fundus images using the Three Continent AMD Consortium Severity Scale. Results: Among the 2262 individuals 70 to 95 years old (627 with AMD and 1635 without AMD), RTL was significantly shorter in individuals with AMD compared to AMD-free participants. In age- and sex-adjusted logistic regression analyses, we observed an 8% higher odds for AMD per 0.1 unit shorter RTL (odds ratio [OR] = 1.08; 95% confidence interval [CI], 1.02-1.14; P = 0.005). The estimates remained stable when adjusted for smoking, high-density lipoprotein cholesterol, cardiovascular disease, diabetes, and hypertension. Interestingly, this association was only present in women (OR = 1.14; 95% CI, 1.06-1.23; P < 0.001), but not in men (OR = 1.01; 95% CI, 0.93-1.10; P = 0.76). A significant sex-by-RTL interaction on AMD was detected (P = 0.043). Conclusions: Our results show an association of RTL with AMD that was restricted to women. This is in line with altered reactive oxygen species levels and higher telomerase activity in women and provides an indication for a sex-differential pathway for oxidative stress and AMD.

An in-depth analysis of the mitochondrial phylogenetic landscape of Cambodia.

Cambodia harbours a variety of human aboriginal populations that have scarcely been studied in terms of genetic diversity of entire mitochondrial genomes. Here we present the matrilineal gene pool of 299 Cambodian refugees from three different ethnic groups (Cham, Khmer, and Khmer Loeu) deriving from 16 Cambodian districts. After establishing a DNA-saving high-throughput strategy for mitochondrial whole-genome Sanger sequencing, a HaploGrep based workflow was used for quality control, haplogroup classification and phylogenetic reconstruction. The application of diverse phylogenetic algorithms revealed an exciting picture of the genetic diversity of Cambodia, especially in relation to populations from Southeast Asia and from the whole world. A total of 224 unique haplotypes were identified, which were mostly classified under haplogroups B5a1, F1a1, or categorized as newly defined basal haplogroups or basal sub-branches of R, N and M clades. The presence of autochthonous maternal lineages could be confirmed as reported in previous studies. The exceptional homogeneity observed between and within the three investigated Cambodian ethnic groups indicates genetic isolation of the whole population. Between ethnicities, genetic barriers were not detected. The mtDNA data presented here increases the phylogenetic resolution in Cambodia significantly, thereby highlighting the need for an update of the current human mtDNA phylogeny.

Autotrophic carbon fixation strategies used by nitrifying prokaryotes in freshwater lakes.

Niche specialization of nitrifying prokaryotes is usually studied with tools targeting molecules involved in the oxidation of ammonia and nitrite. The ecological significance of diverse CO2 fixation strategies used by nitrifiers is, however, mostly unexplored. By analyzing autotrophy-related genes in combination with amoA marker genes based on droplet digitial PCR and CARD-FISH counts targeting rRNA, we quantified the distribution of nitrifiers in eight stratified lakes. Ammonia oxidizing (AO) Thaumarchaeota using the 3-hydroxypropionate/4-hydroxybutyrate pathway dominated deep and oligotrophic lakes, whereas Nitrosomonas-related taxa employing the Calvin cycle were important AO bacteria in smaller lakes. The occurrence of nitrite oxidizing Nitrospira, assimilating CO2 with the reductive TCA cycle, was strongly correlated with the distribution of Thaumarchaeota. Recently discovered complete ammonia-oxidizing bacteria (comammox) belonging to Nitrospira accounted only for a very small fraction of ammonia oxidizers (AOs) present at the study sites. Altogether, this study gives a first insight on how physicochemical characteristics in lakes are associated to the distribution of nitrifying prokaryotes with different CO2 fixation strategies. Our investigations also evaluate the suitability of functional genes associated with individual CO2 assimilation pathways to study niche preferences of different guilds of nitrifying microorganisms based on an autotrophic perspective.

Bioaccumulation of ultraviolet sunscreen compounds (mycosporine-like amino acids) by the heterotrophic freshwater ciliate Bursaridium living in alpine lakes.

Ciliates in shallow alpine lakes are exposed to high levels of incident solar ultraviolet radiation (UVR). We observed the presence of specific sunscreen compounds, the mycosporine-like amino acids (MAAs), in several populations of Bursaridium, a relatively large ciliate species found in such lakes. The populations from 3 highly UV transparent lakes revealed the presence of 7 MAAs (MG, SH, PR, PI, AS, US, and PE) in total concentrations of 3.6-52.4 10-5 µg µg-1 dry weight (DW) per individual, whereas in one glacially turbid and less UV transparent lake, no MAAs were detected in the Bursaridium population. The MAAs in the ciliates generally reflected the composition and relative amounts of the lakes' seston MAAs, assuming that the ciliates fed on MAA-rich plankton. We experimentally found that naturally acquired MAAs prevented ciliate mortality under simulated UVR and photosynthetically active radiation (PAR) conditions. We further tested the dietary regulation of the MAAs-content in the ciliates under artificial UVR and PAR exposure and found an increase in MAAs concentrations in all treatments. Our assumption was that several stress factors other than irradiation were involved in the synthesis or up-regulation of MAAs.

CO2 assimilation strategies in stratified lakes: Diversity and distribution patterns of chemolithoautotrophs.

While mechanisms of different carbon dioxide (CO2 ) assimilation pathways in chemolithoautotrohic prokaryotes are well understood for many isolates under laboratory conditions, the ecological significance of diverse CO2 fixation strategies in the environment is mostly unexplored. Six stratified freshwater lakes were chosen to study the distribution and diversity of the Calvin-Benson-Bassham (CBB) cycle, the reductive tricarboxylic acid (rTCA) cycle, and the recently discovered archaeal 3-hydroxypropionate/4-hydroxybutyrate (HP/HB) pathway. Eleven primer sets were used to amplify and sequence genes coding for selected key enzymes in the three pathways. Whereas the CBB pathway with different forms of RubisCO (IA, IC and II) was ubiquitous and related to diverse bacterial taxa, encompassing a wide range of potential physiologies, the rTCA cycle in Epsilonproteobacteria and Chloribi was exclusively detected in anoxic water layers. Nitrifiying Nitrosospira and Thaumarchaeota, using the rTCA and HP/HB cycle respectively, are important residents in the aphotic and (micro-)oxic zone of deep lakes. Both taxa were of minor importance in surface waters and in smaller lakes characterized by an anoxic hypolimnion. Overall, this study provides a first insight on how different CO2 fixation strategies and chemical gradients in lakes are associated to the distribution of chemoautotrophic prokaryotes with different functional traits.

Heme oxygenase-1 gene promoter microsatellite polymorphism is associated with progressive atherosclerosis and incident cardiovascular disease.

OBJECTIVE: The enzyme heme oxygenase-1 (HO-1) exerts cytoprotective effects in response to various cellular stressors. A variable number tandem repeat polymorphism in the HO-1 gene promoter region has previously been linked to cardiovascular disease. We examined this association prospectively in the general population. APPROACH AND RESULTS: Incidence of stroke, myocardial infarction, or vascular death was registered between 1995 and 2010 in 812 participants of the Bruneck Study aged 45 to 84 years (49.4% males). Carotid atherosclerosis progression was quantified by high-resolution ultrasound. HO-1 variable number tandem repeat length was determined by polymerase chain reaction. Subjects with ≥32 tandem repeats on both HO-1 alleles compared with the rest of the population (recessive trait) featured substantially increased cardiovascular disease risk (hazard ratio [95% confidence interval], 5.45 [2.39, 12.42]; P<0.0001), enhanced atherosclerosis progression (median difference in atherosclerosis score [interquartile range], 2.1 [0.8, 5.6] versus 0.0 [0.0, 2.2] mm; P=0.0012), and a trend toward higher levels of oxidized phospholipids on apolipoprotein B-100 (median oxidized phospholipids/apolipoprotein B level [interquartile range], 11364 [4160, 18330] versus 4844 [3174, 12284] relative light units; P=0.0554). Increased cardiovascular disease risk in those homozygous for ≥32 repeats was also detected in a pooled analysis of 7848 participants of the Bruneck, SAPHIR, and KORA prospective studies (hazard ratio [95% confidence interval], 3.26 [1.50, 7.33]; P=0.0043). CONCLUSIONS: This study found a strong association between the HO-1 variable number tandem repeat polymorphism and cardiovascular disease risk confined to subjects with a high number of repeats on both HO-1 alleles and provides evidence for accelerated atherogenesis and decreased antioxidant defense in this vascular high-risk group.

Large-scale mitochondrial DNA analysis in Southeast Asia reveals evolutionary effects of cultural isolation in the multi-ethnic population of Myanmar.

BACKGROUND: Myanmar is the largest country in mainland Southeast Asia with a population of 55 million people subdivided into more than 100 ethnic groups. Ruled by changing kingdoms and dynasties and lying on the trade route between India and China, Myanmar was influenced by numerous cultures. Since its independence from British occupation, tensions between the ruling Bamar and ethnic minorities increased. RESULTS: Our aim was to search for genetic footprints of Myanmar's geographic, historic and sociocultural characteristics and to contribute to the picture of human colonization by describing and dating of new mitochondrial DNA (mtDNA) haplogroups. Therefore, we sequenced the mtDNA control region of 327 unrelated donors and the complete mitochondrial genome of 44 selected individuals according to highest quality standards. CONCLUSION: Phylogenetic analyses of the entire mtDNA genomes uncovered eight new haplogroups and three unclassified basal M-lineages. The multi-ethnic population and the complex history of Myanmar were reflected in its mtDNA heterogeneity. Population genetic analyses of Burmese control region sequences combined with population data from neighboring countries revealed that the Myanmar haplogroup distribution showed a typical Southeast Asian pattern, but also Northeast Asian and Indian influences. The population structure of the extraordinarily diverse Bamar differed from that of the Karen people who displayed signs of genetic isolation. Migration analyses indicated a considerable genetic exchange with an overall positive migration balance from Myanmar to neighboring countries. Age estimates of the newly described haplogroups point to the existence of evolutionary windows where climatic and cultural changes gave rise to mitochondrial haplogroup diversification in Asia.

Blockade of receptor activator of nuclear factor-κB (RANKL) signaling improves hepatic insulin resistance and prevents development of diabetes mellitus.

Hepatic insulin resistance is a driving force in the pathogenesis of type 2 diabetes mellitus (T2DM) and is tightly coupled with excessive storage of fat and the ensuing inflammation within the liver. There is compelling evidence that activation of the transcription factor nuclear factor-κB (NF-κB) and downstream inflammatory signaling pathways systemically and in the liver are key events in the etiology of hepatic insulin resistance and ß-cell dysfunction, although the molecular mechanisms involved are incompletely understood. We here test the hypothesis that receptor activator of NF-κB ligand (RANKL), a prototypic activator of NF-κB, contributes to this process using both an epidemiological and experimental approach. In the prospective population-based Bruneck Study, a high serum concentration of soluble RANKL emerged as a significant (P<0.001) and independent risk predictor of T2DM manifestation. In close agreement, systemic or hepatic blockage of RANKL signaling in genetic and nutritional mouse models of T2DM resulted in a marked improvement of hepatic insulin sensitivity and amelioration or even normalization of plasma glucose concentrations and glucose tolerance. Overall, this study provides evidence for a role of RANKL signaling in the pathogenesis of T2DM. If so, translation to the clinic may be feasible given current pharmacological strategies to lower RANKL activity to treat osteoporosis.

Genetic determinants of the ankle-brachial index: a meta-analysis of a cardiovascular candidate gene 50K SNP panel in the candidate gene association resource (CARe) consortium.

BACKGROUND: Candidate gene association studies for peripheral artery disease (PAD), including subclinical disease assessed with the ankle-brachial index (ABI), have been limited by the modest number of genes examined. We conducted a two stage meta-analysis of ∼50,000 SNPs across ∼2100 candidate genes to identify genetic variants for ABI. METHODS AND RESULTS: We studied subjects of European ancestry from 8 studies (n=21,547, 55% women, mean age 44-73 years) and African American ancestry from 5 studies (n=7267, 60% women, mean age 41-73 years) involved in the candidate gene association resource (CARe) consortium. In each ethnic group, additive genetic models were used (with each additional copy of the minor allele corresponding to the given beta) to test each SNP for association with continuous ABI (excluding ABI>1.40) and PAD (defined as ABI<0.90) using linear or logistic regression with adjustment for known PAD risk factors and population stratification. We then conducted a fixed-effects inverse-variance weighted meta-analyses considering a p<2×10(-6) to denote statistical significance. RESULTS: In the European ancestry discovery meta-analyses, rs2171209 in SYTL3 (ß=-0.007, p=6.02×10(-7)) and rs290481 in TCF7L2 (ß=-0.008, p=7.01×10(-7)) were significantly associated with ABI. None of the SNP associations for PAD were significant, though a SNP in CYP2B6 (p=4.99×10(-5)) was among the strongest associations. These 3 genes are linked to key PAD risk factors (lipoprotein(a), type 2 diabetes, and smoking behavior, respectively). We sought replication in 6 population-based and 3 clinical samples (n=15,440) for rs290481 and rs2171209. However, in the replication stage (rs2171209, p=0.75; rs290481, p=0.19) and in the combined discovery and replication analysis the SNP-ABI associations were no longer significant (rs2171209, p=1.14×10(-3); rs290481, p=8.88×10(-5)). In African Americans, none of the SNP associations for ABI or PAD achieved an experiment-wide level of significance. CONCLUSIONS: Genetic determinants of ABI and PAD remain elusive. Follow-up of these preliminary findings may uncover important biology given the known gene-risk factor associations. New and more powerful approaches to PAD gene discovery are warranted.

Factors involved in the distribution pattern of ciliates in the water column of a transparent alpine lake.

The recurrent depth preference of three ciliate species (two prostomatids and one haptorid) in a transparent alpine lake indicates the existence of niche partitioning among them involving potential factors such as avoidance of high ultraviolet radiation levels and zooplankton predation, as well as competition for food resources.

Expression of the IGF axis is decreased in local prostate cancer but enhanced after benign prostate epithelial differentiation and TGF-ß treatment.

The insulin-like growth factor (IGF) axis is a molecular pathway intensively investigated in cancer research. Clinical trials targeting the IGF1 receptor (IGF1R) in different tumors, including prostate cancer, are under way. Although studies on the IGF axis in prostate cancer have already entered into clinical trials, the expression and functional role of the IGF axis in benign prostate and in prostate cancer needs to be better defined. We determined mRNA expression levels of the IGF axis in microdissected tissue specimens of local prostate cancer using quantitative PCR. All members of the IGF axis, including IGF1, IGF2, IGF binding proteins 1 through 6, and insulin receptor, were measured in both the stromal and epithelial compartments of the prostate. IGF1, IGF2, IGF1R, and insulin receptor were down-regulated in local prostate cancer tissue compared with matched benign tissue, suggesting that the IGF axis is not induced during prostate cancer development. Using a new prostate epithelial differentiation model, we demonstrate that the expression of the IGF axis is enhanced during normal prostate epithelial differentiation and regulated by tumor growth factor (TGF)-ß. Our data reveal a functional role of the IGF axis in prostate differentiation, underscoring the importance of the IGF axis in normal development and emphasizing the importance of accurate target validation before moving to advanced clinical trials.

Meta analysis of candidate gene variants outside the LPA locus with Lp(a) plasma levels in 14,500 participants of six White European cohorts.

BACKGROUND: Both genome-wide association studies and candidate gene studies have reported that the major determinant of plasma levels of the Lipoprotein (a) [Lp(a)] reside within the LPA locus on chromosome 6. We have used data from the HumanCVD BeadChip to explore the contribution of other candidate genes determining Lp(a) levels. METHODS: 48,032 single nucleotide polymorphisms (SNPs) from the Illumina HumanCVD BeadChip were genotyped in 5059 participants of the Whitehall II study (WHII) of randomly ascertained healthy men and women. SNPs showing association with Lp(a) levels of p<10(-4) outside the LPA locus were selected for replication in a total of an additional 9463 participants of five European based studies (EAS, EPIC-Norfolk, NPHSII, PROCARDIS, and SAPHIR). RESULTS: In Whitehall II, apart from the LPA locus (where p values for several SNPs were <10(-30)) there was significant association at four loci GALNT2, FABP1, PPARGC1A and TNFRSFF11A. However, a meta-analysis of the six studies did not confirm any of these findings. CONCLUSION: Results from this meta analysis of 14,522 participants revealed no candidate genes from the HumanCVD BeadChip outside the LPA locus to have an effect on Lp(a) levels. Further studies with genome-wide and denser SNP coverage are required to confirm or refute this finding.

Are freshwater mixotrophic ciliates less sensitive to solar ultraviolet radiation than heterotrophic ones?

We tested whether mixotrophic ciliates are more resistant to solar ultraviolet radiation (UVR) than heterotrophic ones because symbiotic algae can provide self-shading by cell matter absorption and eventually by direct UV screening from mycosporine-like amino acids (MAAs). Sensitivity of a natural assemblage to solar radiation was tested in experiments in the original lake and in a more UV transparent alpine lake after transplantation of the ciliates. In both lakes, the assemblage was exposed either to full sunlight, to photosynthetically active radiation only, or kept in the dark. In each lake, exposure was for 5 h at the surface and at the depth corresponding to the 10% attenuation depth at 320 nm. Overall, when the assemblage was exposed to surface UVR, only one out of four dominant mixotrophic ciliates, Vorticella chlorellata, was more resistant than heterotrophic species. The higher UV resistance in V. chlorellata was related to the presence of MAAs and the high percentage of ciliate volume occupied by algal symbionts. Our results indicate that effects of UVR were species-specific and depended on efficient screening of these wavelengths, but also on the depth preference of the ciliates and thus, on their previous exposure history to UVR.

Somatic mutations throughout the entire mitochondrial genome are associated with elevated PSA levels in prostate cancer patients.

The genetic etiology of prostate cancer, the most common form of male cancer in western countries, is complex and the interplay of disease genes with environmental factors is far from being understood. Studies on somatic mitochondrial DNA (mtDNA) mutations have become an important aspect of cancer research because these mutations might have functional consequences and/or might serve as biosensors for tumor detection and progression. We sequenced the entire mitochondrial genome (16,569 bp) from 30 prospectively collected pairs of macrodissected cancerous and benign cells from prostate cancer patients and compared their genetic variability. Given recent concerns regarding the authenticity of newly discovered mtDNA mutations, we implemented a high-quality procedure for mtDNA whole-genome sequencing. In addition, the mitochondrial genes MT-CO2, MT-CO3, MT-ATP6, and MT-ND6 were sequenced in further 35 paired samples from prostate cancer patients. We identified a total of 41 somatic mutations in 22 out of 30 patients: the majority of these mutations have not previously been observed in the human phylogeny. The presence of somatic mutations in transfer RNAs (tRNAs) was found to be associated with elevated PSA levels (14.25 ± 5.44 versus 7.15 ± 4.32 ng/ml; p = 0.004). The level and degree of heteroplasmy increased with increasing tumor activity. In summary, somatic mutations in the mitochondrial genome are frequent events in prostate cancer. Mutations mapping to mitochondrial tRNAs, ribosomal RNAs, and protein coding genes might impair processes that occur within the mitochondrial compartment (e.g., transcription, RNA processing, and translation) and might finally affect oxidative phosphorylation.

Association of STR polymorphisms in CMA1 and IL-4 with asthma and atopy: the SAPALDIA cohort.

Asthma is a chronic pulmonary disorder that is characterized by airway inflammation and bronchial hyperreactivity. Several genetic loci have been associated with asthma, and some of these associations have been replicated in independent studies. However, larger population-based replication studies for the association of short tandem repeat (STR) polymorphisms with asthma are limited. In this study, we investigated the association of STR polymorphisms in genes encoding mast cell chymase (CMA1), uteroglobin (UGB), tumor necrosis factor-α (TNF-α) and interleukin-4 (IL-4) with asthma and atopic phenotypes in the large population-based Swiss Cohort Study SAPALDIA. Our results show that the STR polymorphism in the CMA1 gene is associated with asthma and that this association is even stronger with atopic asthma. Similarly, we observed a weak association of the IL-4 2-allele with asthma that tended to be stronger for atopic asthma than for nonatopic asthma. This minor IL-4 2-allele was also associated with higher IgE levels, with a higher risk for a positive skin prick test and with a trend for a higher risk for bronchial hyperresponsiveness. These results support previous findings suggesting a role for CMA1 and IL-4 in atopic asthma and for IL-4 in atopy in general.

Sex and age interaction with genetic association of atherogenic uric acid concentrations.

BACKGROUND: High serum uric acid levels are associated with gout, atherosclerosis and cardiovascular disease. Three genes (SLC2A9, ABCG2, and SLC17A3) were reported to be involved in the regulation of uric acid levels. DESIGN AND METHODS: SNPs rs2231142 (ABCG2) and rs1165205 (SLC17A3) were genotyped in three cohorts (n=4492) and combined with previously genotyped SNPs within SLC2A9 (rs6855911, rs7442295, rs6449213, rs12510549). RESULTS: Each copy of the minor allele decreased uric acid levels by 0.30-0.38 mg/dL for SLC2A9 (p values: 10(-20)-10(-36)) and increased levels by 0.34 mg/dL for ABCG2 (p=1.1x10(-16)). SLC17A3 influenced uric acid levels only modestly. Together the SNPs showed graded associations with uric acid levels of 0.111 mg/dL per risk allele (p=3.8x10(-42)). In addition, we observed a sex-specific interaction of age with the association of SLC2A9 SNPs with uric acid levels, where increasing age strengthened the association of SNPs in women and decreased the association in men. CONCLUSIONS: Genetic variants within SLC2A9,ABCG2 and SLC17A3 show highly significant associations with uric acid levels, and for SNPs within SLC2A9 this association is strongly modified by age and sex.

Symbiotic ciliates receive protection against UV damage from their algae: a test with Paramecium bursaria and Chlorella.

We assessed the photoprotective role of symbiotic Chlorella in the ciliate Paramecium bursaria by comparing their sensitivity to UV radiation (UVR) with Chlorella-reduced and Chlorella-free (aposymbiotic) cell lines of the same species. Aposymbiotic P. bursaria had significantly higher mortality than the symbiotic cell lines when exposed to UVR. To elucidate the protection mechanism, we assessed the algal distribution within the ciliate using thin-sections and transmission electron microscopy and estimated the screening factor by Chlorella based on an optical model. These analyses evidenced a substantial screening factor ranging, from 59.2% to 93.2% (320nm) for regular algal distribution. This screening efficiency reached up to approximately 100% when Chlorella algae were dislocated to the posterior region of the ciliate. The dislocation was observed in symbiotic ciliates only under exposure to UV plus photosynthetically active radiation (PAR) or to high PAR levels. Moreover, under exposure to UVB radiation and high PAR, symbiotic P. bursaria aggregated into dense spots. This behavior could represent an efficient avoidance strategy not yet described for ciliates. Analyses of the intact symbiosis and their algal symbionts for UV-screening compounds (mycosporine-like amino acids and sporopollenin) proved negative. Overall, our results show that photoprotection in this ciliate symbiosis represents an additional advantage to the hitherto postulated nutritional benefits.

CILIATE-SYMBIONT SPECIFICITY OF FRESHWATER ENDOSYMBIOTIC CHLORELLA (TREBOUXIOPHYCEAE, CHLOROPHYTA)(1).

The nature of Chlorella symbioses in invertebrates and protists has attracted much interest, but the uncertain taxonomy of the algal partner has constrained a deeper ecological understanding of this symbiosis. We sequenced parts of the nuclear 18S rDNA, the internal transcribed spacer (ITS)-1 region, and the chloroplast 16S rDNA of several Chlorella isolated from pelagic ciliate species of different lakes, Paramecium bursaria symbionts, and free-living Chlorella to elucidate phylogenetic relationships of Chlorella-like algae and to assess their host specificity. Sequence analyses resulted in well-resolved phylogenetic trees providing strong statistical support for a homogenous 'zoochlorellae' group of different ciliate species from one lake, but clearly different Chlorella in one of those ciliate species occurring in another lake. The two Chlorella strains isolated from the same ciliate species, but from lakes having a 10-fold difference in underwater UV transparency, also presented a distinct physiological trait, such as the ability to synthesize UV-absorbing substances known as mycosporine-like amino acids (MAAs). Algal symbionts of all P. bursaria strains of different origin resolved in one clade apart from the other ciliate symbionts but split into two distinct lineages, suggesting the existence of a biogeographic pattern. Overall, our results suggest a high degree of species specificity but also hint at the importance of physiological adaptation in symbiotic Chlorella.

An experimental test of the symbiosis specificity between the ciliate Paramecium bursaria and strains of the unicellular green alga Chlorella.

The ciliate Paramecium bursaria living in mutualistic relationship with the unicellular green alga Chlorella is known to be easily infected by various potential symbionts/parasites such as bacteria, yeasts and other algae. Permanent symbiosis, however, seems to be restricted to Chlorella taxa. To test the specificity of this association, we designed infection experiments with two aposymbiotic P. bursaria strains and Chlorella symbionts isolated from four Paramecium strains, seven other ciliate hosts and two Hydra strains, as well as three free-living Chlorella species. Paramecium bursaria established stable symbioses with all tested Chlorella symbionts of ciliates, but never with symbiotic Chlorella of Hydra viridissima or with free-living Chlorella. Furthermore, we tested the infection specificity of P. bursaria with a 1:1:1 mixture of three compatible Chlorella strains, including the native symbiont, and then identified the strain of the newly established symbiosis by sequencing the internal transcribed spacer region 1 of the 18S rRNA gene. The results indicated that P. bursaria established symbiosis with its native symbiont. We conclude that despite clear preferences for their native Chlorella, the host-symbiont relationship in P. bursaria is flexible.