Genome-Wide Variation Patterns Uncover the Origin and Selection in Cultivated Ginseng (Panax ginseng Meyer)

Chinese ginseng (Panax ginseng Meyer) is a medicinally important herb and plays crucial roles in traditional Chinese medicine. Pharmacological analyses identified diverse bioactive components from Chinese ginseng. However, basic biological attributes including domestication and selection of the ginseng plant remain under-investigated. Here, we presented a genome-wide view of the domestication and selection of cultivated ginseng based on the whole genome data. A total of 8,660 protein-coding genes were selected for genome-wide scanning of the 30 wild and cultivated ginseng accessions. In complement, the 45s rDNA, chloroplast and mitochondrial genomes were included to perform phylogenetic and population genetic analyses. The observed spatial genetic structure between northern cultivated ginseng (NCG) and southern cultivated ginseng (SCG) accessions suggested multiple independent origins of cultivated ginseng. Genome-wide scanning further demonstrated that NCG and SCG have undergone distinct selection pressures during the domestication process, with more genes identified in the NCG (97 genes) than in the SCG group (5 genes). Functional analyses revealed that these genes are involved in diverse pathways, including DNA methylation, lignin biosynthesis, and cell differentiation. These findings suggested that the SCG and NCG groups have distinct demographic histories. Candidate genes identified are useful for future molecular breeding of cultivated ginseng.

Diversity and Genome Analysis of Australian and Global Oilseed Brassica napus L. Germplasm Using Transcriptomics and Whole Genome Re-sequencing

Intensive breeding of Brassica napus has resulted in relatively low diversity, such that B. napus would benefit from germplasm improvement schemes that sustain diversity. As such, samples representative of global germplasm pools need to be assessed for existing population structure, diversity and linkage disequilibrium (LD). Complexity reduction genotyping-by-sequencing (GBS) methods, including GBS-transcriptomics (GBS-t), enable cost-effective screening of a large number of samples, while whole genome re-sequencing (WGR) delivers the ability to generate large numbers of unbiased genomic single nucleotide polymorphisms (SNPs), and identify structural variants (SVs). Furthermore, the development of genomic tools based on whole genomes representative of global oilseed diversity and orientated by the reference genome has substantial industry relevance and will be highly beneficial for canola breeding. As recent studies have focused on European and Chinese varieties, a global diversity panel as well as a substantial number of Australian spring types were included in this study. Focusing on industry relevance, 633 varieties were initially genotyped using GBS-t to examine population structure using 61,037 SNPs. Subsequently, 149 samples representative of global diversity were selected for WGR and both data sets used for a side-by-side evaluation of diversity and LD. The WGR data was further used to develop genomic resources consisting of a list of 4,029,750 high-confidence SNPs annotated using SnpEff, and SVs in the form of 10,976 deletions and 2,556 insertions. These resources form the basis of a reliable and repeatable system allowing greater integration between canola genomics studies, with a strong focus on breeding germplasm and industry applicability.

Selective constraint and adaptive potential of West Nile virus within and among naturally infected avian hosts and mosquito vectors

Arthropod-borne viruses are among the most genetically constrained RNA viruses, yet they have a remarkable propensity to adapt and emerge. We studied wild birds and mosquitoes naturally infected with West Nile virus (WNV) in a 'hot spot' of virus transmission in Chicago, IL, USA. We generated full coding WNV genome sequences from spatiotemporally matched bird and mosquito samples using high-throughput sequencing, allowing a molecular evolutionary assessment with deep coverage. Mean FST among samples was 0.66 (±0.02 SE) and was bimodal, with mean nucleotide diversity being higher between samples (interhost πN = 0.001; πS = 0.024) than within them (intrahost πN < 0.0001; πS < 0.001). Eight genomic sites with FST > 1.01 (in the PrM, NS2a, NS3, NS4b, and 5'-noncoding genomic regions) showed bird versus mosquito variant frequency differences of >30 per cent and/or polymorphisms fixed in ≥5 host or vector individuals, suggesting host tropism for these variants. However, phylogenetic analyses demonstrated a lack of grouping by bird or mosquito, most inter-sample differences were synonymous (mean interhost πN/πS = 0.04), and there was no significant difference between hosts and vectors in either their nucleotide diversities or levels of purifying selection (mean intrahost πN/πS = 0.28 in birds and πN/πS = 0.21 in mosquitoes). This finding contrasts with the 'trade-off' and 'selective sieve' hypotheses that have been proposed and tested in the laboratory, which predict strong host versus vector effects on WNV genetic variation, with heightened selective constraint in birds alternating with heightened viral diversity in mosquitoes. Overall, our data show WNV to be highly selectively constrained within and between both hosts and vectors but still able to vary at a limited number of sites across the genome. Such site-specific plasticity in the face of overall selective constraint may offer a mechanism whereby highly constrained viruses such as WNV and its relatives can still adapt and emerge.

Chromosome arm-specific patterns of polymorphism associated with chromosomal inversions in the major African malaria vector, Anopheles funestus

Chromosomal inversions facilitate local adaptation of beneficial mutations and modulate genetic polymorphism, but the extent of their effects within the genome is still insufficiently understood. The genome of Anopheles funestus, a malaria mosquito endemic to sub-Saharan Africa, contains an impressive number of paracentric polymorphic inversions, which are unevenly distributed among chromosomes and provide an excellent framework for investigating the genomic impacts of chromosomal rearrangements. Here, we present results of a fine-scale analysis of genetic variation within the genome of two weakly differentiated populations of Anopheles funestus inhabiting contrasting moisture conditions in Cameroon. Using population genomic analyses, we found that genetic divergence between the two populations is centred on regions of the genome corresponding to three inversions, which are characterized by high values of FST , absolute sequence divergence and fixed differences. Importantly, in contrast to the 2L chromosome arm, which is collinear, nucleotide diversity is significantly reduced along the entire length of three autosome arms bearing multiple overlapping chromosomal rearrangements. These findings support the idea that interactions between reduced recombination and natural selection within inversions contribute to sculpt nucleotide polymorphism across chromosomes in An. funestus.

Population transcriptomic characterization of the genetic and expression variation of a candidate progenitor of Miscanthus energy crops

The use of transcriptome data in the study of the population genetics of a species can capture faint signals of both genetic variation and expression variation and can provide a broad picture of a species' genomic response to environmental conditions. In this study, we characterized the genetic and expression diversity of Miscanthus lutarioriparius by comparing more than 16,225 transcripts obtained from 78 individuals, belonging to 10 populations distributed across the species' entire geographic range. We only observed a low level of nucleotide diversity (π = 0.000434) among the transcriptome data of these populations, which is consistent with highly conserved sequences of functional elements and protein-coding genes captured with this method. Tests of population divergence using the transcriptome data were consistent with previous microsatellite data but proved to be more sensitive, particularly if gene expression variation was considered as well. For example, the analysis of expression data showed that genes involved in photosynthetic processes and responses to temperature or reactive oxygen species stimuli were significantly enriched in certain populations. This differential gene expression was primarily observed among populations and not within populations. Interestingly, nucleotide diversity was significantly negatively correlated with expression diversity within populations, while this correlation was positive among populations. This suggests that genetic and expression variation play separate roles in adaptation and population persistence. Combining analyses of genetic and gene expression variation represents a promising approach for studying the population genetics of wild species and may uncover both adaptive and nonadaptive processes.

The amount of DNA polymorphism when population size changes linearly

Population size is one of the main factors that determine the amount of DNA polymorphism. We examined a model under which the population size changed linearly. Because of the simplicity of this model, we could analytically obtain the expectation of nucleotide diversity, E(π), and the expectation of the amount of DNA polymorphism, E(θ), based on the number of segregating sites. The results suggest that E(π) is larger than E(θ) when the population size decreased and that E(π) is smaller than E(θ) when the population size increased. The expected time to the most recent common ancestor could also be obtained under this model.

Unexpected high genetic diversity in small populations suggests maintenance by associative overdominance

The effective population size (N_e ) is a central factor in determining maintenance of genetic variation. The neutral theory predicts that loss of variation depends on N_e , with less genetic drift in larger populations. We monitored genetic drift in 42 Drosophila melanogaster populations of different adult census population sizes (10, 50 or 500) using pooled RAD sequencing. In small populations, variation was lost at a substantially lower rate than expected. This observation was consistent across two ecological relevant thermal regimes, one stable and one with a stressful increase in temperature across generations. Estimated ratios between N_e and adult census size were consistently higher in small than in larger populations. The finding provides evidence for a slower than expected loss of genetic diversity and consequently a higher than expected long-term evolutionary potential in small fragmented populations. More genetic diversity was retained in areas of low recombination, suggesting that associative overdominance, driven by disfavoured homozygosity of recessive deleterious alleles, is responsible for the maintenance of genetic diversity in smaller populations. Consistent with this hypothesis, the X-chromosome, which is largely free of recessive deleterious alleles due to hemizygosity in males, fits neutral expectations even in small populations. Our experiments provide experimental answers to a range of unexpected patterns in natural populations, ranging from variable diversity on X-chromosomes and autosomes to surprisingly high levels of nucleotide diversity in small populations.

Hill-Robertson Interference Reduces Genetic Diversity on a Young Plant Y-Chromosome

X and Y chromosomes differ in effective population size (N_e ), rates of recombination, and exposure to natural selection, all of which can affect patterns of genetic diversity. On Y chromosomes with suppressed recombination, natural selection is expected to eliminate linked neutral variation, and lower the N_e of Y compared to X chromosomes or autosomes. However, female-biased sex ratios and high variance in male reproductive success can also reduce Y-linked N_e , making it difficult to infer the causes of low Y-diversity. Here, we investigate the factors affecting levels of polymorphism during sex chromosome evolution in the dioecious plant Rumexhastatulus (Polygonaceae). Strikingly, we find that neutral diversity for genes on the Y chromosome is, on average, 2.1% of the value for their X-linked homologs, corresponding to a chromosome-wide reduction of 93% compared to the standard neutral expectation. We demonstrate that the magnitude of this diversity loss is inconsistent with reduced male N_e caused by neutral processes. Instead, using forward simulations and estimates of the distribution of deleterious fitness effects, we show that Y chromosome diversity loss can be explained by purifying selection acting in aggregate over a large number of genetically linked sites. Simulations also suggest that our observed level of Y-diversity is consistent with the joint action of purifying and positive selection, but only for models in which there were fewer constrained sites than we empirically estimated. Given the relatively recent origin of R. hastatulus sex chromosomes, our results imply that Y-chromosome degeneration in the early stages may be largely driven by selective interference rather than by neutral genetic drift of silenced Y-linked genes.

Moderate nucleotide diversity in the Atlantic herring is associated with a low mutation rate

The Atlantic herring is one of the most abundant vertebrates on earth but its nucleotide diversity is moderate (π = 0.3%), only three-fold higher than in human. Here, we present a pedigree-based estimation of the mutation rate in this species. Based on whole-genome sequencing of four parents and 12 offspring, the estimated mutation rate is 2.0 × 10-9 per base per generation. We observed a high degree of parental mosaicism indicating that a large fraction of these de novo mutations occurred during early germ cell development. The estimated mutation rate - the lowest among vertebrates analyzed to date - partially explains the discrepancy between the rather low nucleotide diversity in herring and its huge census population size. But a species like the herring will never reach its expected nucleotide diversity because of fluctuations in population size over the millions of years it takes to build up high nucleotide diversity.

Cis- and Trans-regulatory Effects on Gene Expression in a Natural Population of Drosophila melanogaster

Cis- and trans-regulatory mutations are important contributors to transcriptome evolution. Quantifying their relative contributions to intraspecific variation in gene expression is essential for understanding the population genetic processes that underlie evolutionary changes in gene expression. Here, we have examined this issue by quantifying genome-wide, allele-specific expression (ASE) variation using a crossing scheme that produces F₁ hybrids between 18 different Drosophila melanogaster strains sampled from the Drosophila Genetic Reference Panel and a reference strain from another population. Head and body samples from F₁ adult females were subjected to RNA sequencing and the subsequent ASE quantification. Cis- and trans-regulatory effects on expression variation were estimated from these data. A higher proportion of genes showed significant cis-regulatory variation (∼28%) than those that showed significant trans-regulatory variation (∼9%). The sizes of cis-regulatory effects on expression variation were 1.98 and 1.88 times larger than trans-regulatory effects in heads and bodies, respectively. A generalized linear model analysis revealed that both cis- and trans-regulated expression variation was strongly associated with nonsynonymous nucleotide diversity and tissue specificity. Interestingly, trans-regulated variation showed a negative correlation with local recombination rate. Also, our analysis on proximal transposable element (TE) insertions suggested that they affect transcription levels of ovary-expressed genes more pronouncedly than genes not expressed in the ovary, possibly due to defense mechanisms against TE mobility in the germline. Collectively, our detailed quantification of ASE variations from a natural population has revealed a number of new relationships between genomic factors and the effects of cis- and trans-regulatory factors on expression variation.

Patterns of Inter-Chromosomal Gene Conversion on the Male-Specific Region of the Human Y Chromosome

The male-specific region of the human Y chromosome (MSY) is characterized by the lack of meiotic recombination and it has long been considered an evolutionary independent region of the human genome. In recent years, however, the idea that human MSY did not have an independent evolutionary history begun to emerge with the discovery that inter-chromosomal gene conversion (ICGC) can modulate the genetic diversity of some portions of this genomic region. Despite the study of the dynamics of this molecular mechanism in humans is still in its infancy, some peculiar features and consequences of it can be summarized. The main effect of ICGC is to increase the allelic diversity of MSY by generating a significant excess of clustered single nucleotide polymorphisms (SNPs) (defined as groups of two or more SNPs occurring in close proximity and on the same branch of the Y phylogeny). On the human MSY, 13 inter-chromosomal gene conversion hotspots (GCHs) have been identified so far, involving donor sequences mainly from the X-chromosome and, to a lesser extent, from autosomes. Most of the GCHs are evolutionary conserved and overlap with regions involved in aberrant X–Y crossing-over. This review mainly focuses on the dynamics and the current knowledge concerning the recombinational landscape of the human MSY in the form of ICGC, on how this molecular mechanism may influence the evolution of the MSY, and on how it could affect the information enclosed within a genomic region which, until recently, appeared to be an evolutionary independent unit.

[Divergence of paralogous growth-hormone-encoding genes and their promoters in Salmonidae]

In many fish species, including salmonids, the growth-hormone is encoded by two duplicated paralogous genes, gh1 and gh2. Both genes were already in place at the time of divergence of species in this group. A comparison of the entire sequence of these genes of salmonids has shown that their conserved regions are associated with exons, while their most variable regions correspond to introns. Introns C and D include putative regulatory elements (sites Pit-1, CRE, and ERE), that are also conserved. In chars, the degree of polymorphism of gh2 gene is 2-3 times as large as that in gh1 gene. However, a comparison across all Salmonidae species would not extent this observation to other species. In both these chars' genes, the promoters are conserved mainly because they correspond to putative regulatory sequences (TATA box, binding sites for the pituitary transcription factor Pit-1 (F1-F4), CRE, GRE and RAR/RXR elements). The promoter of gh2 gene has a greater degree of polymorphism compared with gh1 gene promoter in all investigated species of salmonids. The observed differences in the rates of accumulation of changes in growth hormone encoding paralogs could be explained by differences in the intensity of selection.

Covariation in levels of nucleotide diversity in homologous regions of the avian genome long after completion of lineage sorting

Closely related species may show similar levels of genetic diversity in homologous regions of the genome owing to shared ancestral variation still segregating in the extant species. However, after completion of lineage sorting, such covariation is not necessarily expected. On the other hand, if the processes that govern genetic diversity are conserved, diversity may potentially covary even among distantly related species. We mapped regions of conserved synteny between the genomes of two divergent bird species-collared flycatcher and hooded crow-and identified more than 600 Mb of homologous regions (66% of the genome). From analyses of whole-genome resequencing data in large population samples of both species we found nucleotide diversity in 200 kb windows to be well correlated (Spearman's ρ = 0.407). The correlation remained highly similar after excluding coding sequences. To explain this covariation, we suggest that a stable avian karyotype and a conserved landscape of recombination rate variation render the diversity-reducing effects of linked selection similar in divergent bird lineages. Principal component regression analysis of several potential explanatory variables driving heterogeneity in flycatcher diversity levels revealed the strongest effects from recombination rate variation and density of coding sequence targets for selection, consistent with linked selection. It is also possible that a stable karyotype is associated with a conserved genomic mutation environment contributing to covariation in diversity levels between lineages. Our observations imply that genetic diversity is to some extent predictable.

Sequence diversity patterns suggesting balancing selection in partially sex-linked genes of the plant Silene latifolia are not generated by demographic history or gene flow

DNA sequence diversity in genes in the partially sex-linked pseudoautosomal region (PAR) of the sex chromosomes of the plant Silene latifolia is higher than expected from within-species diversity of other genes. This could be the footprint of sexually antagonistic (SA) alleles that are maintained by balancing selection in a PAR gene (or genes) and affect polymorphism in linked genome regions. SA selection is predicted to occur during sex chromosome evolution, but it is important to test whether the unexpectedly high sequence polymorphism could be explained without it, purely by the combined effects of partial linkage with the sex-determining region and the population's demographic history, including possible introgression from Silene dioica. To test this, we applied approximate Bayesian computation-based model choice to autosomal sequence diversity data, to find the most plausible scenario for the recent history of S. latifolia and then to estimate the posterior density of the most relevant parameters. We then used these densities to simulate variation to be expected at PAR genes. We conclude that an excess of variants at high frequencies at PAR genes should arise in S. latifolia populations only for genes with strong associations with fully sex-linked genes, which requires closer linkage with the fully sex-linked region than that estimated for the PAR genes where apparent deviations from neutrality were observed. These results support the need to invoke selection to explain the S. latifolia PAR gene diversity, and encourage further work to test the possibility of balancing selection due to sexual antagonism.

Distinct patterns of natural selection in Na+/H+ antiporter genes in Populus euphratica and Populus pruinosa

Salt tolerance genes constitute an important class of loci in plant genomes. Little is known about the extent to which natural selection in saline environments has acted upon these loci, and what types of nucleotide diversity such selection has given rise to. Here, we surveyed genetic diversity in three types of Na⁺/H⁺ antiporter gene (SOS, NhaD, and NHX, belonging to the cation/proton antiporter 1 family), which have well‐characterized essential roles in plant salt tolerance. Ten Na⁺/H⁺ antiporter genes and 16 neutral loci randomly selected as controls were sequenced from 17 accessions of two closely related members of the genus Populus, Populus euphratica and Populus pruinosa, section Turanga, which are native to northwest China. The results show that salt tolerance genes are common targets of natural selection in P. euphratica and P. pruinosa. Moreover, the patterns of nucleotide variation across the three types of Na⁺/H⁺ antiporter gene are distinctly different in these two closely related Populus species, and gene flow from P. pruinosa to P. euphratica is highly restricted. Our results suggest that natural selection played an important role in shaping the current distinct patterns of Na⁺/H⁺ antiporter genes, resulting in adaptive evolution in P. euphratica and P. pruinosa.

The Impact of Linked Selection in Chimpanzees: A Comparative Study

Levels of nucleotide diversity vary greatly across the genomes of most species owing to multiple factors. These include variation in the underlying mutation rates, as well as the effects of both direct and linked selection. Fundamental to interpreting the relative importance of these forces is the common observation of a strong positive correlation between nucleotide diversity and recombination rate. While indeed observed in humans, the interpretation of this pattern has been difficult in the absence of high-quality polymorphism data and recombination maps in closely related species. Here, we characterize genetic features driving nucleotide diversity in Western chimpanzees using a recently generated whole genome polymorphism data set. Our results suggest that recombination rate is the primary predictor of nucleotide variation with a strongly positive correlation. In addition, telomeric distance, regional GC-content, and regional CpG-island content are strongly negatively correlated with variation. These results are compared with humans, with both similarities and differences interpreted in the light of the estimated effective population sizes of the two species as well as their strongly differing recent demographic histories.

Loss of genetic diversity as a consequence of selection in response to high pCO2

Standing genetic variation may allow for rapid evolutionary response to the geologically unprecedented changes in global conditions. However, there is little known about the consequences of such rapid evolutionary change. Here, we measure genetic responses to experimental low and high pCO 2 levels in purple sea urchin larvae, Strongylocentrotus purpuratus. We found greater loss of nucleotide diversity in high pCO 2 levels (18.61%; 900 μatm) compared to low pCO 2 levels (10.12%; 400 μatm). In the wild, this loss could limit the evolutionary capacity of future generations. In contrast, we found minimal evidence that purple sea urchin larvae physiologically respond to high pCO 2 through alternative splicing of transcripts (11 genes), despite a strong signal of alternative splicing between different developmental stages (1193 genes). However, in response to high pCO 2, four of the 11 alternatively spliced transcripts encoded ribosomal proteins, suggesting the regulation of translation as a potential response mechanism. The results of this study indicate that while the purple urchin presently may have enough standing genetic variation in response to rapid environmental change, this reservoir of resilience is a finite resource and could quickly diminish.

Genetic Diversity and Molecular Evolution of a Violaxanthin De-epoxidase Gene in Maize

Violaxanthin de-epoxidase (VDE) has a critical role in the carotenoid biosynthesis pathway, which is involved in protecting the photosynthesis apparatus from damage caused by excessive light. Here, a VDE gene in maize, ZmVDE1, was cloned and shown to have functional domains in common with the gramineous VDE protein. Candidate gene association analysis indicated that no polymorphic sites in ZmVDE1 were significant association with any of the examined carotenoid-related traits at P = 0.05 in an association panel containing 155 maize inbred lines. Nucleotide diversity analysis of VDE1 in maize and teosinte indicated that its exon had less genetic variation, consistent with the conserved function of VDE1 in plants. In addition, dramatically reduced nucleotide diversity, fewer haplotypes and a significantly negative parameter deviation for Tajima’s D test of ZmVDE1 in maize and teosinte suggested that a potential selective force had acted across the ZmVDE1 locus. We further identified a 4.2 Mb selective sweep with low recombination surrounding the ZmVDE1 locus that resulted in severely reduced nucleotide diversity on chromosome 2. Collectively, natural selection and the conserved domains of ZmVDE1 might show an important role in the xanthophyll cycle of the carotenoid biosynthesis pathway.

Population Genomic Analysis Reveals Differential Evolutionary Histories and Patterns of Diversity across Subgenomes and Subpopulations of Brassica napus L

The allotetraploid species Brassica napus L. is a global crop of major economic importance, providing canola oil (seed) and vegetables for human consumption and fodder and meal for livestock feed. Characterizing the genetic diversity present in the extant germplasm pool of B. napus is fundamental to better conserve, manage and utilize the genetic resources of this species. We used sequence-based genotyping to identify and genotype 30,881 SNPs in a diversity panel of 782 B. napus accessions, representing samples of winter and spring growth habits originating from 33 countries across Europe, Asia, and America. We detected strong population structure broadly concordant with growth habit and geography, and identified three major genetic groups: spring (SP), winter Europe (WE), and winter Asia (WA). Subpopulation-specific polymorphism patterns suggest enriched genetic diversity within the WA group and a smaller effective breeding population for the SP group compared to WE. Interestingly, the two subgenomes of B. napus appear to have different geographic origins, with phylogenetic analysis placing WE and WA as basal clades for the other subpopulations in the C and A subgenomes, respectively. Finally, we identified 16 genomic regions where the patterns of diversity differed markedly from the genome-wide average, several of which are suggestive of genomic inversions. The results obtained in this study constitute a valuable resource for worldwide breeding efforts and the genetic dissection and prediction of complex B. napus traits.

Genetic Diversity on the Human X Chromosome Does Not Support a Strict Pseudoautosomal Boundary

Unlike the autosomes, recombination between the X chromosome and the Y chromosome is often thought to be constrained to two small pseudoautosomal regions (PARs) at the tips of each sex chromosome. PAR1 spans the first 2.7 Mb of the proximal arm of the human sex chromosomes, whereas the much smaller PAR2 encompasses the distal 320 kb of the long arm of each sex chromosome. In addition to PAR1 and PAR2, there is a human-specific X-transposed region that was duplicated from the X to the Y chromosome. The X-transposed region is often not excluded from X-specific analyses, unlike the PARs, because it is not thought to routinely recombine. Genetic diversity is expected to be higher in recombining regions than in nonrecombining regions because recombination reduces the effect of linked selection. In this study, we investigated patterns of genetic diversity in noncoding regions across the entire X chromosome of a global sample of 26 unrelated genetic females. We found that genetic diversity in PAR1 is significantly greater than in the nonrecombining regions (nonPARs). However, rather than an abrupt drop in diversity at the pseudoautosomal boundary, there is a gradual reduction in diversity from the recombining through the nonrecombining regions, suggesting that recombination between the human sex chromosomes spans across the currently defined pseudoautosomal boundary. A consequence of recombination spanning this boundary potentially includes increasing the rate of sex-linked disorders (e.g., de la Chapelle) and sex chromosome aneuploidies. In contrast, diversity in PAR2 is not significantly elevated compared to the nonPARs, suggesting that recombination is not obligatory in PAR2. Finally, diversity in the X-transposed region is higher than in the surrounding nonPARs, providing evidence that recombination may occur with some frequency between the X and Y chromosomes in the X-transposed region.

Duplication of OsHAP family genes and their association with heading date in rice

Heterotrimeric Heme Activator Protein (HAP) family genes are involved in the regulation of flowering in plants. It is not clear how many HAP genes regulate heading date in rice. In this study, we identified 35 HAP genes, including seven newly identified genes, and performed gene duplication and candidate gene-based association analyses. Analyses showed that segmental duplication and tandem duplication are the main mechanisms of HAP gene duplication. Expression profiling and functional identification indicated that duplication probably diversifies the functions of HAP genes. A nucleotide diversity analysis revealed that 13 HAP genes underwent selection. A candidate gene-based association analysis detected four HAP genes related to heading date. An investigation of transgenic plants or mutants of 23 HAP genes confirmed that overexpression of at least four genes delayed heading date under long-day conditions, including the previously cloned Ghd8/OsHAP3H. Our results indicate that the large number of HAP genes in rice was mainly produced by gene duplication, and a few HAP genes function to regulate heading date. Selection of HAP genes is probably caused by their diverse functions rather than regulation of heading.

The impacts of polyploidy, geographic and ecological isolations on the diversification of Panax (Araliaceae)

BACKGROUND

Panax L. is a medicinally important genus within family Araliaceae, where almost all species are of cultural significance for traditional Chinese medicine. Previous studies suggested two independent origins of the East Asia and North America disjunct distribution of this genus and multiple rounds of whole genome duplications (WGDs) might have occurred during the evolutionary process.

RESULTS

We employed multiple chloroplast and nuclear markers to investigate the evolution and diversification of Panax. Our phylogenetic analyses confirmed previous observations of the independent origins of disjunct distribution and both ancient and recent WGDs have occurred within Panax. The estimations of divergence time implied that the ancient WGD might have occurred before the establishment of Panax. Thereafter, at least two independent recent WGD events have occurred within Panax, one of which has led to the formation of three geographically isolated tetraploid species P. ginseng, P. japonicus and P. quinquefolius. Population genetic analyses showed that the diploid species P. notoginseng harbored significantly lower nucleotide diversity than those of the two tetraploid species P. ginseng and P. quinquefolius and the three species showed distinct nucleotide variation patterns at exon regions.

CONCLUSION

Our findings based on the phylogenetic and population genetic analyses, coupled with the species distribution patterns of Panax, suggested that the two rounds of WGD along with the geographic and ecological isolations might have together contributed to the evolution and diversification of this genus.

A simulation study of sample size for DNA barcoding

For some groups of organisms, DNA barcoding can provide a useful tool in taxonomy, evolutionary biology, and biodiversity assessment. However, the efficacy of DNA barcoding depends on the degree of sampling per species, because a large enough sample size is needed to provide a reliable estimate of genetic polymorphism and for delimiting species. We used a simulation approach to examine the effects of sample size on four estimators of genetic polymorphism related to DNA barcoding: mismatch distribution, nucleotide diversity, the number of haplotypes, and maximum pairwise distance. Our results showed that mismatch distributions derived from subsamples of ≥20 individuals usually bore a close resemblance to that of the full dataset. Estimates of nucleotide diversity from subsamples of ≥20 individuals tended to be bell‐shaped around that of the full dataset, whereas estimates from smaller subsamples were not. As expected, greater sampling generally led to an increase in the number of haplotypes. We also found that subsamples of ≥20 individuals allowed a good estimate of the maximum pairwise distance of the full dataset, while smaller ones were associated with a high probability of underestimation. Overall, our study confirms the expectation that larger samples are beneficial for the efficacy of DNA barcoding and suggests that a minimum sample size of 20 individuals is needed in practice for each population.

Accurate Genetic Detection of Hepatitis C Virus Transmissions in Outbreak Settings

Hepatitis C is a major public health problem in the United States and worldwide. Outbreaks of hepatitis C virus (HCV) infections are associated with unsafe injection practices, drug diversion, and other exposures to blood and are difficult to detect and investigate. Here, we developed and validated a simple approach for molecular detection of HCV transmissions in outbreak settings. We obtained sequences from the HCV hypervariable region 1 (HVR1), using end-point limiting-dilution (EPLD) technique, from 127 cases involved in 32 epidemiologically defined HCV outbreaks and 193 individuals with unrelated HCV strains. We compared several types of genetic distances and calculated a threshold, using minimal Hamming distances, that identifies transmission clusters in all tested outbreaks with 100% accuracy. The approach was also validated on sequences obtained using next-generation sequencing from HCV strains recovered from 239 individuals, and findings showed the same accuracy as that for EPLD. On average, the nucleotide diversity of the intrahost population was 6.2 times greater in the source case than in any incident case, allowing the correct detection of transmission direction in 8 outbreaks for which source cases were known. A simple and accurate distance-based approach developed here for detecting HCV transmissions streamlines molecular investigation of outbreaks, thus improving the public health capacity for rapid and effective control of hepatitis C.

Genome-wide patterns of recombination, linkage disequilibrium and nucleotide diversity from pooled resequencing and single nucleotide polymorphism genotyping unlock the evolutionary history of Eucalyptus grandis

We used high-density single nucleotide polymorphism (SNP) data and whole-genome pooled resequencing to examine the landscape of population recombination (ρ) and nucleotide diversity (ϴw ), assess the extent of linkage disequilibrium (r(2) ) and build the highest density linkage maps for Eucalyptus. At the genome-wide level, linkage disequilibrium (LD) decayed within c. 4-6 kb, slower than previously reported from candidate gene studies, but showing considerable variation from absence to complete LD up to 50 kb. A sharp decrease in the estimate of ρ was seen when going from short to genome-wide inter-SNP distances, highlighting the dependence of this parameter on the scale of observation adopted. Recombination was correlated with nucleotide diversity, gene density and distance from the centromere, with hotspots of recombination enriched for genes involved in chemical reactions and pathways of the normal metabolic processes. The high nucleotide diversity (ϴw = 0.022) of E. grandis revealed that mutation is more important than recombination in shaping its genomic diversity (ρ/ϴw = 0.645). Chromosome-wide ancestral recombination graphs allowed us to date the split of E. grandis (1.7-4.8 million yr ago) and identify a scenario for the recent demographic history of the species. Our results have considerable practical importance to Genome Wide Association Studies (GWAS), while indicating bright prospects for genomic prediction of complex phenotypes in eucalypt breeding.