高山松的群体历史及其适应性进化研究

物种形成一直以来是进化生物学的一个中心议题。因为杂交能够快速创造出新的遗传变异和提高基因组重组率，加速物种形成和适应性进化，所以杂交物种形成是其重要的组成部分。但是直到最近几年，人们才逐渐意识到同倍体杂种形成的重要性和普遍性。此外，人工杂交种群结果证实了杂种的适应性起源于亲本基因的超亲分离现象，这些基因之间通过加和效应或上位效应相互连锁，并在杂种后代中逐渐固定，因此，探讨维持这种适应机制的杂种的连锁不平衡模式就显得尤为重要。到目前为止，关于自然杂种的核苷酸多态性和连锁不平衡结构的报道非常稀见，对杂种的进化历史知之更少。高山松能生长于亲本种不能正常繁殖生长的青藏高原上，这种极端环境为杂种的适应性进化和与亲本的生态隔离造就了绝佳的条件，其生态和生殖的稳定性极大地方便了我们研究植物物种形成和适应性进化的遗传机理，对高山松开展的生态学、遗传学的系列分析已使高山松成为同倍体杂种形成的经典范例。 该论文通过随机挑选7个核基因位点对5个高山松群体、3个油松群体和3个云南松群体的164个大配子体单倍体基因组进行了核苷酸多态性分析。研究发现：所有基因座位的单倍型组成和基因谱系的拓扑结构都支持高山松的多次杂交起源。高山松的平均核苷酸多态性与油松相当，θW达到0.0107;比云南松高出一倍;高于已报道的裸子植物类群遗传变异水平。高山松如此高的遗传多态性与其杂交特性和有效群体大有关。我们运用分子钟和共祖模拟分析，推测出高山松的有效群体为7.32 x 105。等位基因共祖时间模拟分析发现，杂交过程早于青藏高原的隆升，也即是说，在杂种稳定成种以前，两个亲本种之间存在广泛的渐渗杂交。此外，高山松杂种的群体间分化严重，并且在不同的群体中，我们找到了多个偏离中性的基因座位。由此说明，高山松复杂的进化历史和适应的策略区域化。 在显著偏离中性的基因座位上，选取群体历史清晰的高山松群体，对其全基因序列多态性和连锁不平衡结构进行了深入调查，以期找到选择作用位点或区域。Ara-like和Dhn1两个基因在高山松与亲本之间都没有发现固定变异。在高山松中，Ara-like基因和Dhn1基因核苷酸多态性式样和连锁不平衡结构存在明显差异，这说明在两个基因上选择作用的程度和方向不尽相同。 我们对Dhn1基因的多态性分布、LD结构和中性偏离水平进行了分析，结果表明，该基因可能在PdNX和PdLZ群体中受到平衡选择的作用。从Dhn1基因谱系结构可以看出有两种来源的等位基因造成这种平衡多态性，一种属于祖先类型，另一种是从云南松继承衍生而来。这两种等位基因之间的分化很大，极有可能造成编码蛋白的亚功能化。事实上，Dhn1编码的脱水蛋白在植物对环境的抗逆过程中发挥着重要作用。由此我们推测高山松在高海拔极端环境中的适应性进化与Dhn1基因的平衡作用有关。 对Ara-like和Dhn1基因进行HKA检测，结果表明，与Dhn1基因相比，在Ara-like基因上，双维管束亚属与单维管束亚属之间存在显著分化，沉默突变位点的分化Ksil达到了0.1392，远远高于平均水平0.0508;在双维管束亚属共发现了43个固定突变位点，其中有6个能导致氨基酸突变，它们有可能导致了Ara-like基因功能的分化。结合裸子植物近缘物种间共享多态性的普遍性，我们推测Ara-like基因可能在单维管束亚属中的进化速率加快，暗示其在单维管束亚属中的适应性进化。 最后，基因内连锁不平衡分析结果显示，随机筛选的基因座位之间不存在连锁不平衡。高山松平均的基因内LD程度非常低，仅有18%的信息位点之间显著连锁。平均基因内LD在油松和高山松中的衰减速率很快，尤其在杂种中下降最快，在不到200bp以内就降到0.1以下。 LD的结果印证了高山松的有效群体大和多次起源特性。另外，我们也怀疑杂交物种形成过程中染色体组的重组和重排频繁发生，也是造成自然杂种现有群体的LD水平低的一个主要原因。 通过LD衰减曲线估计高山松的平均单位重组率,比亲本油松高17倍，比云南松高45倍。这个结果表明，杂种共适应的基因间要维持超亲分离，需要强烈的自然选择压力，才能保持它们之间的连锁不平衡。

Mancha areolada da seringueira (Hevea brasiliensis) na Amazônia: evolução do patógeno (Thanatephorus cucumeris/Rhizoctonia solani grupo de anastomose AG 2-2 Hb) num patossistema tropical

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

Ancient Dispersal of the Human Fungal Pathogen Cryptococcus gattii from the Amazon Rainforest

Over the past two decades, several fungal outbreaks have occurred, including the high-profile 'Vancouver Island' and 'Pacific Northwest' outbreaks, caused by Cryptococcus gattii, which has affected hundreds of otherwise healthy humans and animals. Over the same time period, C. gattii was the cause of several additional case clusters at localities outside of the tropical and subtropical climate zones where the species normally occurs. In every case, the causative agent belongs to a previously rare genotype of C. gattii called AFLP6/VGII, but the origin of the outbreak clades remains enigmatic. Here we used phylogenetic and recombination analyses, based on AFLP and multiple MLST datasets, and coalescence gene genealogy to demonstrate that these outbreaks have arisen from a highly-recombining C. gattii population in the native rainforest of Northern Brazil. Thus the modern virulent C. gattii AFLP6/VGII outbreak lineages derived from mating events in South America and then dispersed to temperate regions where they cause serious infections in humans and animals.

Genetic consequences of habitat fragmentation during a range expansion

We investigate the effect of habitat fragmentation on the genetic diversity of a species experiencing a range expansion. These two evolutionary processes have not been studied yet, at the same time, owing to the difficulties of deriving analytic results for non-equilibrium models. Here we provide a description of their interaction by using extensive spatial and temporal coalescent simulations and we suggest guidelines for a proper genetic sampling to detect fragmentation. To model habitat fragmentation, we simulated a two-dimensional lattice of demes partitioned into groups (patches) by adding barriers to dispersal. After letting a population expand on this grid, we sampled lineages from the lattice at several scales and studied their coalescent history. We find that in order to detect fragmentation, one needs to extensively sample at a local level rather than at a landscape level. This is because the gene genealogy of a scattered sample is less sensitive to the presence of genetic barriers. Considering the effect of temporal changes of fragmentation intensities, we find that at least 10, but often >100, generations are needed to affect local genetic diversity and population structure. This result explains why recent habitat fragmentation does not always lead to detectable signatures in the genetic structure of populations. Finally, as expected, long-distance dispersal increases local genetic diversity and decreases levels of population differentiation, efficiently counteracting the effects of fragmentation.

The phylogeny of closely related species as revealed by the genealogy of a speciation gene, Odysseus

Molecular differentiation between races or closely related species is often incongruent with the reproductive divergence of the taxa of interest. Shared ancient polymorphism and/or introgression during secondary contact may be responsible for the incongruence. At loci contributing to speciation, these two complications should be minimized (1, 2); hence, their variation may more faithfully reflect the history of the species' reproductive differentiation. In this study, we analyzed DNA polymorphism at the Odysseus (OdsH) locus of hybrid sterility between Drosophila mauritiana and Drosophila simulans and were able to verify such a prediction. Interestingly, DNA variation only a short distance away (1.8 kb) appears not to be influenced by the forces that shape the recent evolution of the OdsH coding region. This locus thus may represent a test case of inferring phylogeny of very closely related species.

Variants in the human potassium channel gene (KCNN3) are associated with migraine in a high risk genetic isolate

The calcium-activated potassium ion channel gene (KCNN3) is located in the vicinity of the familial hemiplegic migraine type 2 locus on chromosome 1q21.3. This gene is expressed in the central nervous system and plays a role in neural excitability. Previous association studies have provided some, although not conclusive, evidence for involvement of this gene in migraine susceptibility. To elucidate KCNN3 involvement in migraine, we performed gene-wide SNP genotyping in a high-risk genetic isolate from Norfolk Island, a population descended from a small number of eighteenth century Isle of Man ‘Bounty Mutineer’ and Tahitian founders. Phenotype information was available for 377 individuals who are related through the single, well-defined Norfolk pedigree (96 were affected: 64 MA, 32 MO). A total of 85 SNPs spanning the KCNN3 gene were genotyped in a sub-sample of 285 related individuals (76 affected), all core members of the extensive Norfolk Island ‘Bounty Mutineer’ genealogy. All genotyping was performed using the Illumina BeadArray platform. The analysis was performed using the statistical program SOLAR v4.0.6 assuming an additive model of allelic effect adjusted for the effects of age and sex. Haplotype analysis was undertaken using the program HAPLOVIEW v4.0. A total of four intronic SNPs in the KCNN3 gene displayed significant association (P < 0.05) with migraine. Two SNPs, rs73532286 and rs6426929, separated by approximately 0.1 kb, displayed complete LD (r 2 = 1.00, D′ = 1.00, D′ 95% CI = 0.96–1.00). In all cases, the minor allele led to a decrease in migraine risk (beta coefficient = 0.286–0.315), suggesting that common gene variants confer an increased risk of migraine in the Norfolk pedigree. This effect may be explained by founder effect in this genetic isolate. This study provides evidence for association of variants in the KCNN3 ion channel gene with migraine susceptibility in the Norfolk genetic isolate with the rarer allelic variants conferring a possible protective role. This the first comprehensive analysis of this potential candidate gene in migraine and also the first study that has utilised the unique Norfolk Island large pedigree isolate to implicate a specific migraine gene. Studies of additional variants in KCNN3 in the Norfolk pedigree are now required (e.g. polyglutamine variants) and further analyses in other population data sets are required to clarify the association of the KCNN3 gene and migraine risk in the general outbred population.

Gene flow from Corymbia hybrids in northern New South Wales

Some of the most productive taxa for forestry are interspecific F1 hybrids grown as exotics in the tropics and subtropics. Attributes of resilience, adaptability and vigour which engender the hybrids for wood production, may also exacerbate the risk they present from gene flow to native species gene pools or to local ecologies as weeds. To determine the biological and genetic factors that influence the extent of hybridisation, we examine the distribution and genealogy of wildlings surrounding plantings of locally-exotic Corymbia torelliana (Section Cadageria) near native C. henryi (Section Maculatae) in northern New South Wales. Our study showed pre-mating and pre- and post-zygotic barriers were incomplete, with in situ generation and natural establishment of both F1 hybrids (n = 3) and advanced generation hybrids under the disturbed conditions bordering native forest. As hybrids were located on alluvial flats exposed to frost, they also likely have an extended ecological range relative to native C. henryi. Despite the likely generation of large viable seed crops on F1 trees at the site over many years, establishment success and survival of advanced generation hybrids may be low, as only 5 immature and no mature advanced generation hybrids were identified. Propagation and genetic analysis of a seed crop from one F1 wildling showed early survival and vigour of seedlings in cultivation was high, and that at least for some F1 in some seasons, backcrossing to the recurrent native C. henryi parent is favoured (60%), whereas selfing (10%) and crossing with other F1 (30%) was less frequent. Transport of seed by stingless bees probably accounted for long distance dispersal from C. torelliana, but this mechanism does not appear to supplement gravity-dispersal of seed from the F1. Coupled with other evidence from studies of bee behaviour, controlled pollination in Corymbia sp., and long-term fitness in second generation eucalypt hybrids, we anticipate gene flow via pollen rather than seed will be the greater challenge for managing the risk of introgression of C. torelliana ancestry into native species from the planted F1 hybrid. If large sources of F1 pollen become available to compete with native pollen, gene flow will probably be frequent and hybrids may establish in disturbed conditions and in habitats beyond the ecological range of their native parent. Further study is needed to determine the degree to which outbreeding depression and poor survival inhibits on-going gene flow.

Progress toward a global genealogy of common carp (Cyprinus carpio L.) strains using mitochondrial nucleotide sequences data

As part of a study of genetic variation in the Vietnamese strains of the common carp (Cyprinus carpio L.) using direct DNA sequencing of mitochondrial control and ATPase6/8 gene regions, samples from a number of other countries were analyzed for comparison. Results show that the levels of sequence divergence in common carp is low on a global scale, with the Asian carp having the highest diversity while Koi and European carp are invariant. A genealogical analysis supports a close relationship among Vietnamese, Koi, Chinese Color and, to a lesser extent, European carp. Koi carp appear to have originated from a strain of Chinese red carp. There is considerable scope to extend this research through the analysis of additional samples of carp from around the world, especially from China, in order to generate a comprehensive global genealogy of common carp strains.

DNA variation in a 13-Mb region including the F9 gene: inferring the genealogical history and causal role of a hemophilia B mutation (IVS 5+13 A -> G)

About 5.5% of all UK hemophilia B patients have the base substitution IVS 5+13 A-->G as the only change in their factor (F)IX gene (F9). This generates a novel donor splice site which fits the consensus better than the normal intron 5 donor splice. Use of the novel splice site should result in a missense mutation followed by the abnormal addition of four amino acids to the patients' FIX. In order to explain the prevalence of this mutation, its genealogical history is examined. Analysis of restriction fragment length polymorphism in the 21 reference UK individuals (from different families) with the above mutation showed identical haplotypes in 19 while two differed from the rest and from each other. In order to investigate the history of the mutation and to verify that it had occurred independently more than once, the sequence variation in 1.5-kb segments scattered over a 13-Mb region including F9 was examined in 18 patients and 15 controls. This variation was then analyzed with a recently developed Bayesian approach that reconstructs the genealogy of the gene investigated while providing evidence of independent mutations that contribute disconnected branches to the genealogical tree. The method also provides minimum estimates of the age of the mutation inherited by the members of coherent trees. This revealed that 17 or 18 mutant genes descend from a founder who probably lived 450 years ago, while one patient carries an independent mutation. The independent recurrence of the IVS5+13 A-->G mutation strongly supports the conclusion that it is the cause of these patients' mild hemophilia.

Growth hormone 1 gene (GH1) polymorphisms as possible markers of the production potential of beef cattle using the Brazilian Canchim breed as a model

The growth hormone 1 gene (GH1) is a candidate gene for body weight and weight gain in cattle since it plays a fundamental role in growth regulation. We investigated the GH1 gene AluI and DdeI restriction enzyme polymorphisms, located 149 bp apart in the cattle genome, as possible markers of the production potential of Canchim crossbreed cattle, a 5/8 Charolais (Bos taurus) and 3/8 Nelore (Bos indicus) breed developed in Brazil, by evaluating the birth weight, weaning weight, yearling weight and plasma insulin-like growth factor-1 (IGF-1) concentration of 7 month to 10 months old Canchim calves (n = 204) of known genealogy and which had been genotyped for the AluI and DdeI markers. Our results showed significant effect (p < 0.05) between the homozygous DdeI+/DdeI+ polymorphism and the estimated breeding value for weaning weight (ESB-WW), while the AluI leucine homozygous (L/L) and leucine/valine (L/V) heterozygous polymorphisms showed no significant effect on the traits studied. The restriction sites of the two enzymes led to the formation of haplotypes which also exerted a significant effect (p < 0.05) on the ESB-WW, with the largest difference being 8.5 kg in favor of the homozygous L plus DdeI+/L plus DdeI+ genotype over the heterozygous L plus DdeI-/V plus DdeI+ genotype.

Cryptic speciation and recombination in the fungus Paracoccidioides brasiliensis as revealed by gene genealogies

Paracoccidioides brasiliensis is the etiologic agent of paracoccidioidomycosis, a disease confined to Latin America and of marked importance in the endemic areas due to its frequency and severity. This species is considered to be clonal according to mycological criteria and has been shown to vary in virulence. To characterize natural genetic variation and reproductive mode in this fungus, we analyzed P. brasiliensis phylogenetically in search of cryptic species and possible recombination using concordance and nondiscordance of gene genealogies with respect to phylogenies of eight regions in five nuclear loci. Our data indicate that this fungus consists of at least three distinct, previously unrecognized species: S1 (species 1 with 38 isolates), PS2 (phylogenetic species 2 with six isolates), and PS3 (phylogenetic species 3 with 21 isolates). Genealogies of four of the regions studied strongly supported the PS2 clade, composed of five Brazilian and one Venezuelan isolate. The second clade, PS3, composed solely of 21 Colombian isolates, was strongly supported by the alpha-tubulin genealogy. The remaining 38 individuals formed S1. Two of the three lineages of P. brasiliensis, S1 and PS2, are sympatric across their range, suggesting barriers to gene flow other than geographic isolation. Our study provides the first evidence for possible sexual reproduction in P. brasiliensis S1, but does not rule it out in the other two species.

Horizontal gene transfer confers fermentative metabolism in the respiratory-deficient plant trypanosomatid Phytomonas serpens

Among trypanosomatids, the genus Phytomonas is the only one specifically adapted to infect plants. These hosts provide a particular habitat with a plentiful supply of carbohydrates. Phytomonas sp. lacks a cytochrome-mediated respiratory chain and Krebs cycle, and ATP production relies predominantly on glycolysis. We have characterised the complete gene encoding a putative pyruvate/indolepyruvate decarboxylase (PDC/IPDC) (548 amino acids) of P. serpens, that displays high amino acid sequence similarity with phytobacteria and Leishmania enzymes. No orthologous PDC/IPDC genes were found in Trypanosoma cruzi or T. brucei. Conservation of the PDC/IPDC gene sequence was verified in 14 Phytomonas isolates. A phylogenetic analysis shows that Phytomonas protein is robustly monophyletic with Leishmania spp. and C. fasciculata enzymes. In the trees this clade appears as a sister group of indolepyruvate decarboxylases of gamma-proteobacteria. This supports the proposition that a horizontal gene transfer event from a donor phytobacteria to a recipient ancestral trypanosome has occurred prior to the separation between Phytomonas. Leishmania and Crithidia. We have measured the PDC activity in P. serpens cell extracts. The enzyme has a Km value for pyruvate of 1.4 mM. The acquisition of a PDC, a key enzyme in alcoholic fermentation, explains earlier observations that ethanol is one of the major end-products of glucose catabolism under aerobic and anaerobic conditions. This represents an alternative and necessary route to reoxidise part of the NADH produced in the highly demanding glycolytic pathway and highlights the importance of this type of event in metabolic adaptation. (C) 2012 Elsevier B.V. All rights reserved.

Phylogeny of genes for secretion NTPases: Identification of the widespread tadA subfamily and development of a diagnostic key for gene classification

Macromolecular transport systems in bacteria currently are classified by function and sequence comparisons into five basic types. In this classification system, type II and type IV secretion systems both possess members of a superfamily of genes for putative NTP hydrolase (NTPase) proteins that are strikingly similar in structure, function, and sequence. These include VirB11, TrbB, TraG, GspE, PilB, PilT, and ComG1. The predicted protein product of tadA, a recently discovered gene required for tenacious adherence of Actinobacillus actinomycetemcomitans, also has significant sequence similarity to members of this superfamily and to several unclassified and uncharacterized gene products of both Archaea and Bacteria. To understand the relationship of tadA and tadA-like genes to those encoding the putative NTPases of type II/IV secretion, we used a phylogenetic approach to obtain a genealogy of 148 NTPase genes and reconstruct a scenario of gene superfamily evolution. In this phylogeny, clear distinctions can be made between type II and type IV families and their constituent subfamilies. In addition, the subgroup containing tadA constitutes a novel and extremely widespread subfamily of the family encompassing all putative NTPases of type IV secretion systems. We report diagnostic amino acid residue positions for each major monophyletic family and subfamily in the phylogenetic tree, and we propose an easy method for precisely classifying and naming putative NTPase genes based on phylogeny. This molecular key-based method can be applied to other gene superfamilies and represents a valuable tool for genome analysis.

No evidence for association of ataxia-telangiectasia mutated gene T2119C and C3161G amino acid substitution variants with risk of breast cancer

Background There is evidence that certain mutations in the double-strand break repair pathway ataxia-telangiectasia mutated gene act in a dominant-negative manner to increase the risk of breast cancer. There are also some reports to suggest that the amino acid substitution variants T2119C Ser707Pro and C3161G Pro1054Arg may be associated with breast cancer risk. We investigate the breast cancer risk associated with these two nonconservative amino acid substitution variants using a large Australian population-based case–control study. Methods The polymorphisms were genotyped in more than 1300 cases and 600 controls using 5' exonuclease assays. Case–control analyses and genotype distributions were compared by logistic regression. Results The 2119C variant was rare, occurring at frequencies of 1.4 and 1.3% in cases and controls, respectively (P = 0.8). There was no difference in genotype distribution between cases and controls (P = 0.8), and the TC genotype was not associated with increased risk of breast cancer (adjusted odds ratio = 1.08, 95% confidence interval = 0.59–1.97, P = 0.8). Similarly, the 3161G variant was no more common in cases than in controls (2.9% versus 2.2%, P = 0.2), there was no difference in genotype distribution between cases and controls (P = 0.1), and the CG genotype was not associated with an increased risk of breast cancer (adjusted odds ratio = 1.30, 95% confidence interval = 0.85–1.98, P = 0.2). This lack of evidence for an association persisted within groups defined by the family history of breast cancer or by age. Conclusion The 2119C and 3161G amino acid substitution variants are not associated with moderate or high risks of breast cancer in Australian women.

The Expanded Human Kallikrein (KLK) gene family : genomic organisation, tissue specific expression and potential functions

The tissue kallikreins are serine proteases encoded by highly conserved multigene families. The rodent kallikrein (KLK) families are particularly large, consisting of 13 26 genes clustered in one chromosomal locus. It has been recently recognised that the human KLK gene family is of a similar size (15 genes) with the identification of another 12 related genes (KLK4-KLK15) within and adjacent to the original human KLK locus (KLK1-3) on chromosome 19q13.4. The structural organisation and size of these new genes is similar to that of other KLK genes except for additional exons encoding 5 or 3 untranslated regions. Moreover, many of these genes have multiple mRNA transcripts, a trait not observed with rodent genes. Unlike all other kallikreins, the KLK4-KLK15 encoded proteases are less related (25–44%) and do not contain a conventional kallikrein loop. Clusters of genes exhibit high prostatic (KLK2-4, KLK15) or pancreatic (KLK6-13) expression, suggesting evolutionary conservation of elements conferring tissue specificity. These genes are also expressed, to varying degrees, in a wider range of tissues suggesting a functional involvement of these newer human kallikrein proteases in a diverse range of physiological processes.