Molecular variability in the Celleporella hyalina (Bryozoa; Cheilostomata) species complex: evidence for cryptic speciation from complete mitochondrial genomes

The bryozoan Celleporella has been shown to be composed of multiple, often cryptic, lineages. We sequenced two complete mitochondrial (mt) genomes of the Celleporella hyalina species complex from Wales, UK and Norway (i) to determine genetic divergence at the complete mt genome level, and (ii) to design new molecular markers for examining the interrelationships amongst the major lineages. In addressing (i), we estimated genetic divergence at three levels: (a) nucleotide diversity (π), (b) genome size, and (c) gene order. Genes nad4L, nad6, and atp8 showed the highest levels of divergence, and rrnL, rrnS, and cox1 showed the lowest levels. Inter-genome nucleotide divergence of protein-coding and ribosomal RNA genes, measured as π, was 0.21. The two genomes differed substantially in size, with the Norwegian genome being 2,573 base pairs (bp) longer than the Welsh genome, 17,265 and 14,692 bp, respectively. This difference in size is attributable to long non-coding regions present in the Norwegian genome. Both genomes exhibit similar gene orders, except for the translocation of one transfer RNA (trnA). Considering the high nucleotide diversity, genome size difference and change in gene order, these mt genomes are considered sufficiently divergent to have originated from two distinct species. In addressing (ii) we designed PCR primers that flank the most conserved regions of the genome: 1,300 bp of cox1 and a contiguous 2,000 bp fragment of rrnL + rrnS. The primers have yielded products for tissue from Wales, Norway, New Zealand, Alaska and Chile and should provide useful tools in establishing species- and population-level diversity within the Celleporella complex.     

Low levels of nucleotide diversity at homoeologous Adh loci in allotetraploid cotton (Gossypium L.).

Levels of genetic diversity within and among populations and species are shaped by both external (population-level) and internal (genomic and genic) evolutionary forces. To address the effect of internal pressures, we estimated nucleotide diversity for a pair of homoeologous Adh loci in an allotetraploid species, Gossypium hirsutum. These data represent the first such estimates for a pair of homoeologous nuclear loci in plants. Estimates of nucleotide diversity for AdhA in Gossypium are lower than those for any plant nuclear gene yet described. This low diversity appears to reflect primarily a history of repeated, severe genetic bottlenecks associated with both speciation and recent domestication, supplemented by an unusually slow nucleotide substitution rate and an autogamous breeding system. While not statistically supportable, the sum of the observations also suggest differential evolutionary dynamics at each of the homoeologous loci.     

Individual haplotyping of whale sharks from seawater environmental DNA

Population genetic data can provide valuable information on the demography of a species. For rare and elusive marine megafauna, samples for generating the data are traditionally obtained from tissue biopsies, which can be logistically difficult and expensive to collect and require invasive sampling techniques. Analysis of environmental DNA (eDNA) offers an alternative, minimally invasive approach to provide important genetic information. Although eDNA approaches have been studied extensively for species detection and biodiversity monitoring in metabarcoding studies, the potential for the technique to address population-level questions remains largely unexplored. Here, we applied “eDNA haplotyping” to obtain estimates of the intraspecific genetic diversity of a whale shark (Rhincodon typus) aggregation at Ningaloo reef, Australia. Over 2 weeks, we collected seawater samples directly behind individual sharks prior to taking a tissue biopsy sample from the same animal. Our data showed a 100% match between mtDNA sequences recovered in the eDNA and tissue sample for all 28 individuals sampled. In the seawater samples, >97% of all reads were assigned to six dominant haplotypes, and a clear dominant signal (~99% of sample reads) was recovered in each sample. Our study demonstrates accurate individual-level haplotyping from seawater eDNA. When DNA from one individual clearly dominates each eDNA sample, it provides many of the same opportunities for population genetic analyses as a tissue sample, potentially removing the need for tissue sampling. Our results show that eDNA approaches for population-level analyses have the potential to supply critical demographic data for the conservation and management of marine megafauna.     

Estimating the contribution of mutation, recombination and gene conversion in the generation of haplotypic diversity

Recombination occurs through both homologous crossing over and homologous gene conversion during meiosis. The contribution of recombination relative to mutation is expected to be dramatically reduced in inbreeding organisms. We report coalescent-based estimates of the recombination parameter (rho) relative to estimates of the mutation parameter (theta) for 18 genes from the highly self-fertilizing grass, wild barley, Hordeum vulgare ssp. spontaneum. Estimates of rho/theta are much greater than expected, with a mean rho/theta approximately 1.5, similar to estimates from outcrossing species. We also estimate rho with and without the contribution of gene conversion. Genotyping errors can mimic the effect of gene conversion, upwardly biasing estimates of the role of conversion. Thus we report a novel method for identifying genotyping errors in nucleotide sequence data sets. We show that there is evidence for gene conversion in many large nucleotide sequence data sets including our data that have been purged of all detectable sequencing errors and in data sets from Drosophila melanogaster, D. simulans, and Zea mays. In total, 13 of 27 loci show evidence of gene conversion. For these loci, gene conversion is estimated to contribute an average of twice as much as crossing over to total recombination.     

Estimates of correlations among yield traits and somatic cell score with different models to adjust for bovine somatotropin effects on Holstein dairy cows

Records of Holstein cows from the Dairy Records Processing Center at Raleigh, NC were edited to obtain three data sets: 65,720 first, 50,694 second, and 65,445 later lactations. Correlations among yield traits and somatic cell score were estimated with three different models: 1) bovine somatotropin (bST) administration ignored, 2) bST administration as a fixed effect and 3) administration of bST as part of the contemporary group (herd-year-month-bST). Heritability estimates ranged from 0.13 to 0.17 for milk, 0.12 to 0.20 for fat, 0.14 to 0.16 for protein yields, and 0.08 to 0.09 for somatic cell score. Estimates were less for later than first lactations. Estimates of genetic correlations among yields ranged from 0.35 to 0.85 with no important differences between estimates with the 3 models. Estimates for lactation 2 agreed with estimates for lactation 1. Estimates of genetic correlations for later lactations were generally greater than for lactations 1 and 2 except between milk and protein yields. Estimates of genetic correlations between yields and somatic cell score were mostly negative or small (-0.45 to 0.11). Estimates of environmental correlations among yield traits were similar with all models (0.77 to 0.97). Estimates of environmental correlations between yields and somatic cell score were negative (-0.22 to -0.14). Estimates of phenotypic correlations among yield traits ranged from 0.70 to 0.95. Estimates of phenotypic correlations between yields and somatic cell score were small and negative. For all three data sets and all traits, no important differences in estimates of genetic parameters were found for the two models that adjusted for bST and the model that did not.     

Comparison of different nuclear DNA markers for estimating intraspecific genetic diversity in plants

A compilation was made of 307 studies using nuclear DNA markers for evaluating among- and within-population diversity in wild angiosperms and gymnosperms. Estimates derived by the dominantly inherited markers (RAPD, AFLP, ISSR) are very similar and may be directly comparable. STMS analysis yields almost three times higher values for within-population diversity whereas among-population diversity estimates are similar to those derived by the dominantly inherited markers. Number of sampled plants per population and number of scored microsatellite DNA alleles are correlated with some of the population genetics parameters. In addition, maximum geographical distance between sampled populations has a strong positive effect on among-population diversity. As previously verified with allozyme data, RAPD- and STMS-based analyses show that long-lived, outcrossing, late successional taxa retain most of their genetic variability within populations. By contrast, annual, selfing and/or early successional taxa allocate most of the genetic variability among populations. Estimates for among- and within-population diversity, respectively, were negatively correlated. The only major discrepancy between allozymes and STMS on the one hand, and RAPD on the other hand, concerns geographical range; within-population diversity was strongly affected when the former methods were used but not so in the RAPD-based studies. Direct comparisons between the different methods, when applied to the same plant material, indicate large similarities between the dominant markers and somewhat lower similarity with the STMS-based data, presumably due to insufficient number of analysed microsatellite DNA loci in many studies.     

Planar π-aromatic C3h B6H 3 + and π-antiaromatic C2h B8H2: boron hydride analogues of D3h C3H 3 + and D2h C4H4

Based upon extensive density functional theory and wave function theory calculations performed in this work, we predict the existence of the perfectly planar triangle C(3h) B(6)H(3)(+) (1, (1)A') and the double-chain stripe C(2h) B(8)H(2) (9, (1)A(g)) which are the ground states of the systems and the inorganic analogues of cyclopropene cation D(3h) C(3)H (3) (+) and cyclobutadiene D(2h) C(4)H(4), respectively. Detailed adaptive natural density partitioning (AdNDP) analyses indicate that C(3h) B(6)H (3) (+) is π plus σ doubly aromatic with two delocalized π-electrons and six delocalized σ-electrons formally conforming to the 4n + 2 aromatic rule, while C(2h) B(8)H(2) is π antiaromatic and σ aromatic with four delocalized π-electrons and ten delocalized σ-electrons. The perfectly planar C(2h) B(8)H(4) (5, (1)A(g)) also proves to be π antiaromatic analogous to D(2h) C(4)H(4), but it appears to be a local minimum about 50 kJ mol(-1) less stable than the three dimensional C(s) B(8)H(4)(6, (1)A'). AdNDP, nucleus independent chemical shifts (NICS) and electron localization function (ELF) analyses indicate that these boron hydride clusters form islands of both σ- and π-aromaticities and are overall aromatic in nature in ELF aromatic criteria.     

Combination of multiple microsatellite data sets to investigate genetic diversity and admixture of domestic cattle

Microsatellite markers are commonly used for population genetic analyses of livestock. However, up to now, combinations of microsatellite data sets or comparison of population genetic parameters from different studies and breeds has proven difficult. Often different genotyping methods have been employed, preventing standardization of microsatellite allele calling. In other cases different sets of markers have been genotyped, providing differing estimates of population genetic parameters. Here, we address these issues and illustrate a general two-step regression approach in cattle using three different sets of microsatellite data, to combine population genetics estimates of diversity and admixture. This regression-based method is independent of the loci genotyped but requires common breeds in the data sets. We show that combining microsatellite data sets can provide new insights on the origin and geographical distribution of genetic diversity and admixture in cattle, which will facilitate global management of this livestock species.     

Kinetic comparison of seven strains of 2,4-dichlorophenoxyacetic acid-degrading bacteria.

Seven strains of 2,4-dichlorophenoxyacetic acid-degrading bacteria, including Pseudomonas, Alcaligenes, and Bordetella spp., were compared on the basis of growth kinetics. Estimates of maximum growth rate (mu max, k1) and half-saturation growth constant (Ks, k3) were obtained by fitting substrate depletion curves to a four-parameter version of the integrated Monod equation. Estimates of Ks ranged from 2.2 micrograms/ml (10 microM) to 33.8 micrograms/ml (154 microM), and estimates of mu max ranged from 0.20 h-1 (Td = 3.5 h) to 0.32 h-1 (Td = 2.2 h). Estimates of mu max, but not Ks, were affected by changes in initial inoculum density. Maximum growth rates (mu max) were also estimated from turbidity measurements. They ranged from 0.10 h-1 (Td = 6.9 h) to 1.0 h-1 (Td = 0.7 h). There was no correlation between estimates of mu max derived from substrate depletion curves and those derived from turbidity measurements (P = 0.20).     

Kinetic comparison of seven strains of 2,4-dichlorophenoxyacetic acid-degrading bacteria.

Seven strains of 2,4-dichlorophenoxyacetic acid-degrading bacteria, including Pseudomonas, Alcaligenes, and Bordetella spp., were compared on the basis of growth kinetics. Estimates of maximum growth rate (mu max, k1) and half-saturation growth constant (Ks, k3) were obtained by fitting substrate depletion curves to a four-parameter version of the integrated Monod equation. Estimates of Ks ranged from 2.2 micrograms/ml (10 microM) to 33.8 micrograms/ml (154 microM), and estimates of mu max ranged from 0.20 h-1 (Td = 3.5 h) to 0.32 h-1 (Td = 2.2 h). Estimates of mu max, but not Ks, were affected by changes in initial inoculum density. Maximum growth rates (mu max) were also estimated from turbidity measurements. They ranged from 0.10 h-1 (Td = 6.9 h) to 1.0 h-1 (Td = 0.7 h). There was no correlation between estimates of mu max derived from substrate depletion curves and those derived from turbidity measurements (P = 0.20).     

Likely health outcomes for untreated acute febrile illness in the tropics in decision and economic models; a Delphi survey

Background

Modelling is widely used to inform decisions about management of malaria and acute febrile illnesses. Most models depend on estimates of the probability that untreated patients with malaria or bacterial illnesses will progress to severe disease or death. However, data on these key parameters are lacking and assumptions are frequently made based on expert opinion. Widely diverse opinions can lead to conflicting outcomes in models they inform.

Methods and Findings

A Delphi survey was conducted with malaria experts aiming to reach consensus on key parameters for public health and economic models, relating to the outcome of untreated febrile illnesses. Survey questions were stratified by malaria transmission intensity, patient age, and HIV prevalence. The impact of the variability in opinion on decision models is illustrated with a model previously used to assess the cost-effectiveness of malaria rapid diagnostic tests. Some consensus was reached around the probability that patients from higher transmission settings with untreated malaria would progress to severe disease (median 3%, inter-quartile range (IQR) 1–5%), and the probability that a non-malaria illness required antibiotics in areas of low HIV prevalence (median 20%). Children living in low transmission areas were considered to be at higher risk of progressing to severe malaria (median 30%, IQR 10–58%) than those from higher transmission areas (median 13%, IQR 7–30%). Estimates of the probability of dying from severe malaria were high in all settings (medians 60–73%). However, opinions varied widely for most parameters, and did not converge on resurveying.

Conclusions

This study highlights the uncertainty around potential consequences of untreated malaria and bacterial illnesses. The lack of consensus on most parameters, the wide range of estimates, and the impact of variability in estimates on model outputs, demonstrate the importance of sensitivity analysis for decision models employing expert opinion. Results of such models should be interpreted cautiously. The diversity of expert opinion should be recognised when policy options are debated.     

Genetic parameters for milkability from the first three lactations in Fleckvieh cows

Test-day records for average flow rate (AFR) from the routine dairy recording from Bavarian Fleckvieh cows were analysed. Two data sets with observations on approximately 20 000 cows each were sampled from the total data set. For the estimation of variance parameters, a two-step approach was applied. In a first step multiple-trait restricted maximum likelihood (REML) analyses were carried out. For each of the first three lactations, six time periods with up to 33 days were defined. An algorithm for iterative summing of expanded part matrices was applied in order to combine the estimates. In a second step covariance functions (CF) for additive-genetic variances and non-genetic animal variances were derived using second-order Legendre polynomials plus an exponential term. Estimates of test-day heritability for AFR ranged from 0.21 to 0.40, and were largest in lactation 1. For lactations 1 and 3, heritabilities decreased considerably towards the end of lactation. Genetic correlation estimates within lactation decreased as the distance between days in milk (DIM) increased. Genetic correlations between corresponding DIM in the three lactations were generally large, ranging from 0.80 to 0.99. The largest estimates were found between DIM from lactations 2 and 3. Results from this study suggest that including AFR data from second and third lactations in genetic evaluation systems could the improve accuracy of genetic selection.     

周公河齐口裂腹鱼种群的遗传多样性和遗传结构

为探讨周公河中齐口裂腹鱼(Schizothorax prenanti)的遗传多样性和遗传结构, 沿周公河连续设定6个采样点进行齐口裂腹鱼采集, 在进行基因组DNA提取后以线粒体控制区(D-loop)为分子标记进行种群遗传多样性和遗传结构的评估。结果从63尾齐口裂腹鱼中共检测到32个单倍型, 核苷酸多样性(π)为0.013, 单倍型多样性(h)为0.966。遗传多样性最高的为张家湾种群(PopF, π=0.018, h=1.000), 核苷酸多样性最低的为罗坝种群(PopC, π=0.010, h=0.970), 单倍型多样性最低的为瓦屋山大坝下游种群(PopB, π=0.015, h=0.867)。种群间共享13个单倍型, 多数齐口裂腹鱼种群间的遗传分化水平中等或较低, 仅Pop C和Pop F (F_st=0.195, P<0.01), 种群Pop A 和PopE (F_st=0.158, P<0.01)分化程度较高, 显示各种群间较密切的遗传关系。周公河不同河段齐口裂腹鱼群体的遗传多样性处于较高的水平, 各种群间具有较近的遗传关系。     

Rapid assessment of single-copy nuclear DNA variation in diverse species

We investigated the use of PCR primers designed to conserved exons within nuclear DNA to amplify potentially variable regions such as introns or hypervariable exons from a wide range of species. We then explored various approaches to assay population-level variation in these PCR products. Primers designed to amplify regions within the histone H2AF, myoglobin , MHC DQA , and aldolase (ALD) genes gave clean amplifications in diverse mammals (DQA) , and in birds, reptiles and mammals ( aldolase, H2AF, myoglobin ). The sequenced PCR products generally, but not always, confirmed that the correct locus had been amplified. Several primer sets produced smaller size fragments consistent with preferential amplification of intronless pseudogenes; this was confirmed by sequencing seal and reptile H2AF PCR products. Digestion with randomly selected four-base recognizing enzymes detected variation in some cases but not in others. In species/gene combinations with either low (e.g. seal H2AF, ALD-A ) or high (e.g. skink ALD-1 ) nucleotide diversity it was more efficient to sequence a small number of distantly related individuals (e.g. one per geographic population) and from these data to identify informative or potentially informative restriction enzymes for 'targeted' digestion. We conclude that for studies of population-level variation, the optimal approach is to use a battery of primers for initial PCR of both mtDNA and scnDNA loci, select those that give clean amplifications, and sequence one sample from each population to (i) confirm gene identity, (ii) estimate the amount of variation and, (iii) search for diagnostic restriction sites. This will allow determination of the most efficient approach for a large-scale study.     

Depression among Asian-American Adults in the Community: Systematic Review and Meta-Analysis

Objectives

In this systematic review, we provide an overview of the literature on depression among Asian-Americans and explore the possible variations in depression prevalence estimates by methodological and demographic factors.

Methods

Six databases were used to identify studies reporting a prevalence estimate for depression in Asian-American adults in non-clinical settings. Meta-analysis was used to calculate pooled estimates of rates of depression by assessment type. Statistical heterogeneity was assessed for subgroup analyses by gender, age, ethnicity, and other participant characteristics.

Results

A total of 58 studies met the review criteria (n = 21.731 Asian-American adults). Heterogeneity across the studies was considerably high. The prevalence of major depression assessed via standardized clinical interviews ranged between 4.5% and 11.3%. Meta-analyses revealed comparable estimated prevalence rates of depression as measured by the Center for Epidemiologic Studies Depression Scale (35.6%, 95% CI 27.6%–43.7%) and the Geriatric Depression Scale (33.1%, 95% CI 14.9%–51.3%). Estimates varied by Asian racial/ethnic group and other participant characteristics. Estimates of depression among special populations, which included maternity, caregivers, and homosexuals, were significantly higher than estimates obtained from other samples (58.8% vs 29.3%, p = .003). Estimates of depression among Korean and Filipino-Americans were similar (33.3%-34.4%); however, the estimates were twice as high as those for Chinese-Americans (15.7%; p = .012 for Korean, p = .049 for Filipino).

Conclusion

There appears to be wide variability in the prevalence rates of depression among Asian-Americans in the US. Practitioners and researchers who serve Asian-American adults need to be sensitive to the potential diversity of the expression of depression and treatment-seeking across Asian-American subgroups. Public health policies to increase Asian-American access to mental health care, including increased screening, are necessary. Further work is needed to determine whether strategies to reduce depression among specific Asian racial/ethnic groups is warranted.     

基于线粒体D-loop区的抚仙湖鱇浪白鱼遗传多样性分析

为了解鱇浪白鱼(Anabarilius grahami)种质资源和遗传多样性现状, 采集抚仙湖7个地理群体共计133尾鱇浪白鱼, 基于线粒体D-loop控制区全序列开展遗传多样性分析。结果表明: 鱇浪白鱼D-loop控制区序列的A+T含量(62.28%)高于G+C(37.72%); 共发现变异位点40个, 定义单倍型36个, 整体单倍型多样性指数(H_d)为0.791±0.036, 核酸多样性指数(π)为0.00254±0.00027, 呈“高H_d低π遗传分布”; 抚仙湖西岸各群体(H_d = 0.826—0.846, π=0.00236—0.0031)较东岸各群体(H_d =0.657—0.805, π=0.00204—0.00271)的平均遗传多样性呈现出“西高东低”的分布特征; 其中, 明星鱼洞(MX)群体的鱇浪白鱼单倍型多样性和核酸多样性指数最高, 下坝(XB)群体的鱇浪白鱼单倍型多样性和核酸多样性指数最低; 路岐(LQ)群体与海镜(HJ)群体遗传距离最远(0.00296), 但7个群体遗传分化不显著(F_st=0.01954, P＞0.05); 中性检验(Tajima’s D=–1.73617, P=0.03729; Fu’s F_s=–1.64259, P=0.23943)与核酸错配分布均表明抚仙湖鱇浪白鱼群体未经历种群扩张事件。综上所述, 抚仙湖鱇浪白鱼展现出高单倍型多样性和低核酸性的遗传多样性特征, 可能与线粒体进化速度较快有关; 而鱇浪白鱼扩散式生活史及人工增殖放流可能是造成其群体间遗传分化不显著的主要原因。研究利用线粒体DNA开展了抚仙湖鱇浪白鱼遗传多样性现状的初步评估, 为后续鱇浪白鱼种质资源保护及种业开发提供重要理论依据。     

S-allele diversity in Sorbus aucuparia and Crataegus monogyna (Rosaceae: Maloideae)

RT-PCR was used to obtain the first estimates from natural populations of allelic diversity at the RNase-based gametophytic self-incompatibility locus in the Rosaceae. A total of 20 alleles were retrieved from 20 Sorbus aucuparia individuals, whereas 17 alleles were found in 13 Crataegus monogyna samples. Estimates of population-level allele numbers fall within the range observed in the Solanaceae, the only other family with RNase-based incompatibility for which estimates are available. The nucleotide diversity of S-allele sequences was found to be much lower in the two Rosaceae species as compared with the Solanaceae. This was not due to a lower sequence divergence among most closely related alleles. Rather, it is the depth of the entire genealogy that differs markedly in the two families, with Rosaceae S-alleles exhibiting more recent apparent coalescence. We also investigated patterns of selection at the molecular level by comparing nucleotide diversity at synonymous and nonsynonymous sites. Stabilizing selection was inferred for the 5' region of the molecule, while evidence of diversifying selection was present elsewhere.     

Delimiting Species without Nuclear Monophyly in Madagascar's Mouse Lemurs

Background

Speciation begins when populations become genetically separated through a substantial reduction in gene flow, and it is at this point that a genetically cohesive set of populations attain the sole property of species: the independent evolution of a population-level lineage. The comprehensive delimitation of species within biodiversity hotspots, regardless of their level of divergence, is important for understanding the factors that drive the diversification of biota and for identifying them as targets for conservation. However, delimiting recently diverged species is challenging due to insufficient time for the differential evolution of characters—including morphological differences, reproductive isolation, and gene tree monophyly—that are typically used as evidence for separately evolving lineages.

Methodology

In this study, we assembled multiple lines of evidence from the analysis of mtDNA and nDNA sequence data for the delimitation of a high diversity of cryptically diverged population-level mouse lemur lineages across the island of Madagascar. Our study uses a multi-faceted approach that applies phylogenetic, population genetic, and genealogical analysis for recognizing lineage diversity and presents the most thoroughly sampled species delimitation of mouse lemur ever performed.

Conclusions

The resolution of a large number of geographically defined clades in the mtDNA gene tree provides strong initial evidence for recognizing a high diversity of population-level lineages in mouse lemurs. We find additional support for lineage recognition in the striking concordance between mtDNA clades and patterns of nuclear population structure. Lineages identified using these two sources of evidence also exhibit patterns of population divergence according to genealogical exclusivity estimates. Mouse lemur lineage diversity is reflected in both a geographically fine-scaled pattern of population divergence within established and geographically widespread taxa, as well as newly resolved patterns of micro-endemism revealed through expanded field sampling into previously poorly and well-sampled regions.     

Comparative phylogeography of Nearctic and Palearctic fishes

Combining phylogeographic data from mitochondrial DNA (mtDNA) of Nearctic and Palearctic freshwater and anadromous fishes, we used a comparative approach to assess the influence of historical events on evolutionary patterns and processes in regional fish faunas. Specifically, we (i) determined whether regional faunas differentially affected by Pleistocene glaciations show predictable differences in phylogeographic patterns; (ii) evaluated how processes of divergence and speciation have been influenced by such differential responses; and (iii) assessed the general contribution of phylogeographic studies to conservation issues. Comparisons among case studies revealed fundamental differences in phylogeographic patterns among regional faunas. Tree topologies were typically deeper for species from nonglaciated regions compared to northern species, whereas species with partially glaciated ranges were intermediate in their characteristics. Phylogeographic patterns were strikingly similar among southern species, whereas species in glaciated areas showed reduced concordance. The extent and locations of secondary contact among mtDNA lineages varied greatly among northern species, resulting in reduced intraspecific concordance of genetic markers for some northern species. Regression analysis of phylogeographic data for 42 species revealed significant latitudinal shifts in intraspecific genetic diversity. Both relative nucleotide diversity and estimates of evolutionary effective population size showed significant breakpoints matching the median latitude for the southern limit of the Pleistocene glaciations. Similarly, analysis of clade depth of phylogenetically distinct lineages vs. area occupied showed that evolutionary dispersal rates of species from glaciated and nonglaciated regions differed by two orders of magnitude. A negative relationship was also found between sequence divergence among sister species as a function of their median distributional latitude, indicating that recent bursts of speciation events have occurred in deglaciated habitats. Phylogeographic evidence for parallel evolution of sympatric northern species pairs in postglacial times suggested that differentiation of cospecific morphotypes may be driven by ecological release. Altogether, these results demonstrate that comparative phylogeography can be used to evaluate not only phylogeographic patterns but also evolutionary processes. As well as having significant implications for conservation programs, this approach enables new avenues of research for examining the regional, historical, and ecological factors involved in shaping intraspecific genetic diversity.