The selection landscape and genetic legacy of ancient Eurasians.

The Holocene (beginning around 12,000 years ago) encompassed some of the most significant changes in human evolution, with far-reaching consequences for the dietary, physical and mental health of present-day populations. Using a dataset of more than 1,600 imputed ancient genomes1, we modelled the selection landscape during the transition from hunting and gathering, to farming and pastoralism across West Eurasia. We identify key selection signals related to metabolism, including that selection at the FADS cluster began earlier than previously reported and that selection near the LCT locus predates the emergence of the lactase persistence allele by thousands of years. We also find strong selection in the HLA region, possibly due to increased exposure to pathogens during the Bronze Age. Using ancient individuals to infer local ancestry tracts in over 400,000 samples from the UK Biobank, we identify widespread differences in the distribution of Mesolithic, Neolithic and Bronze Age ancestries across Eurasia. By calculating ancestry-specific polygenic risk scores, we show that height differences between Northern and Southern Europe are associated with differential Steppe ancestry, rather than selection, and that risk alleles for mood-related phenotypes are enriched for Neolithic farmer ancestry, whereas risk alleles for diabetes and Alzheimer's disease are enriched for Western hunter-gatherer ancestry. Our results indicate that ancient selection and migration were large contributors to the distribution of phenotypic diversity in present-day Europeans.

Genomic-based microsatellite development for Ternstroemia (Pentaphylacaceae) and transferability to other Ericales.

BACKGROUND: The genus Ternstroemia is associated with the vulnerable tropical montane cloud forest in Mexico and with other relevant vegetation types worldwide. It contains threatened and pharmacologically important species and has taxonomic issues regarding its species limits. This study describes 38 microsatellite markers generated using a genomic-based approach. METHODS AND RESULTS: We tested 23 of these markers in a natural population of Ternstroemia lineata. These markers are highly polymorphic (all loci polymorphic with 3-14 alleles per locus and expected heterozygosity between 0.202 and 0.908), most of them (19 out of 23) are in Hardy-Weinberg Equilibrium and free of null alleles (18 out of 23). Also we found no evidence of linkage among them. Finally, we tested the transferability to six other American species of Ternstroemia, two other Pentaphylacaceae species, and four species from different families within the order Ericales. CONCLUSIONS: These molecular resources are promising tools to investigate genetic diversity loss and as barcodes for ethnopharmacological applications and species delimitation in the family Pentaphylacaceae and some Ericales, among other applications.

The genetic architecture of the human thalamus and its overlap with ten common brain disorders.

The thalamus is a vital communication hub in the center of the brain and consists of distinct nuclei critical for consciousness and higher-order cortical functions. Structural and functional thalamic alterations are involved in the pathogenesis of common brain disorders, yet the genetic architecture of the thalamus remains largely unknown. Here, using brain scans and genotype data from 30,114 individuals, we identify 55 lead single nucleotide polymorphisms (SNPs) within 42 genetic loci and 391 genes associated with volumes of the thalamus and its nuclei. In an independent validation sample (n = 5173) 53 out of the 55 lead SNPs of the discovery sample show the same effect direction (sign test, P = 8.6e-14). We map the genetic relationship between thalamic nuclei and 180 cerebral cortical areas and find overlapping genetic architectures consistent with thalamocortical connectivity. Pleiotropy analyses between thalamic volumes and ten psychiatric and neurological disorders reveal shared variants for all disorders. Together, these analyses identify genetic loci linked to thalamic nuclei and substantiate the emerging view of the thalamus having central roles in cortical functioning and common brain disorders.

Preclinical Nanomedicines for Polyglutamine-Based Neurodegenerative Diseases.

Polyglutamine (polyQ) diseases, such as Huntington's disease and several types of spinocerebellar ataxias, are dominantly inherited progressive neurodegenerative disorders and characterized by the presence of expanded CAG trinucleotide repeats in the respective disease locus of the patient genomes. Patients with polyQ diseases currently need to rely on symptom-relieving treatments because disease-modifying therapeutic interventions remain scarce. Many disease-modifying therapeutic agents are now under clinical testing for treating polyQ diseases, but their delivery to the brain is often too invasive (e.g., intracranial injection) or inefficient, owing to in vivo degradation and clearance by physiological barriers (e.g., oral and intravenous administration). Nanoparticles provide a feasible solution for improving drug delivery to the brain, as evidenced by an increasing number of preclinical studies that document the efficacy of nanomedicines for polyQ diseases over the past 5-6 years. In this review, we present the pathogenic mechanisms of polyQ diseases, the common animal models of polyQ diseases for evaluating the efficacy of nanomedicines, and the common administration routes for delivering nanoparticles to the brain. Next, we summarize the recent preclinical applications of nanomedicines for treating polyQ diseases and improving neurological conditions in vivo, placing emphasis on antisense oligonucleotides, small peptide inhibitors, and small molecules as the disease-modifying agents. We conclude with our perspectives of the burgeoning field of "nanomedicines for polyQ diseases", including the use of inorganic nanoparticles and potential drugs as next-generation nanomedicines, development of higher-order animal models of polyQ diseases, and importance of "brain-nano" interactions.

A functional molecular marker for detecting blister blight disease resistance in tea (Camellia sinensis L.).

KEY MESSAGE: Identification of an EST-SSR molecular marker associated with Blister blight, a common fungal disease of tea, facilitating marker-assisted selection, marking a milestone in tea molecular breeding. lister blight (BB) leaf disease of tea, caused by the fungus Exobasidium vexans, results in 25-30% crop loss annually. BB is presently controlled by Cu based fungicides, but genetic resistance is the most viable option in disease management. Tea is a naturally out-crossing, woody perennial necessitating a long time for completion of a breeding programme. Marker-assisted selection (MAS) is vital to expedite breeding programmes and also for better accuracy in gene identification. The aim of the current research was to derive marker-trait associations using an F1 population segregating for BB. The population was genotyped at 11 expressed sequence tag simple sequence repeat loci followed by detecting the alleles by fragment analysis. The genotypic and phenotypic data were subjected to single-marker analysis resulting in the identification of EST-SSR073 as a diagnostic marker amplifying three alleles of the sizes, 168, 170 and 190 bp in F1. Of them, alleles 190 and 168 bp were confirmed to concur BB resistance and susceptibility, respectively. The alleles were validated in a panel of 64 tea cultivars, resulting in the amplification of 12 alleles at EST-SSR073. The EST-SSR073 allele sequences matched with Camellia sinensis photosystem-I reaction center subunit-II. The marker EST-SSR073 can be effectively used in breeding tea against BB, recording a milestone in MAS in tea.

The molecular basis for allelic differences suggests Restorer-of-fertility 1 is a complex locus in sugar beet (Beta vulgaris L.).

BACKGROUND: Cytoplasmic male sterility (CMS) is a widely used trait for hybrid seed production in many crops. Sugar beet CMS is associated with a unique mitochondrial protein named preSATP6 that forms a 250-kDa complex. Restorer-of-fertility 1 (Rf1) is a nuclear gene that suppresses CMS and is, hence, one of the targets of sugar beet breeding. Rf1 has dominant, semi-dominant and recessive alleles, suggesting that it may be a multi-allelic locus; however, the molecular basis for differences in genetic action is obscure. Molecular cloning of Rf1 revealed a gene (orf20) whose protein products produced in transgenics can bind with preSATP6 to generate a novel 200-kDa complex. The complex is also detected in fertility-restored anthers concomitant with a decrease in the amount of the 250-kDa complex. Molecular diversity of the Rf1 locus involves organizational diversity of a gene cluster composed of orf20-like genes (RF-Oma1s). We examined the possibility that members of the clustered RF-Oma1 in this locus could be associated with fertility restoration. RESULTS: Six yet uncharacterized RF-Oma1s from dominant and recessive alleles were examined to determine whether they could generate the 200-kDa complex. Analyses of transgenic calli revealed that three RF-Oma1s from a dominant allele could generate the 200-kDa complex, suggesting that clustered RF-Oma1s in the dominant allele can participate in fertility restoration. None of the three copies from two recessive alleles was 200-kDa generative. The absence of this ability was confirmed by analyzing mitochondrial complexes in anthers of plants having these recessive alleles. Together with our previous data, we designed a set of PCR primers specific to the 200-kDa generative RF-Oma1s. The amount of mRNA measured by this primer set inversely correlated with the amount of the 250-kDa complex in anthers and positively correlated with the strength of the Rf1 alleles. CONCLUSIONS: Fertility restoration by sugar beet Rf1 can involve multiple RF-Oma1s clustered in the locus, implying that stacking 200-kDa generative copies in the locus strengthens the efficacy, whereas the absence of 200-kDa generative copies in the locus makes the allele recessive irrespective of the copy number. We propose that sugar beet Rf1 is a complex locus.

A Genome-Wide Association Study Identifies the Association between the 12q24 Locus and Black Tea Consumption in Japanese Populations.

Several genome-wide association studies (GWASs) have reported the association between genetic variants and the habitual consumption of foods and drinks; however, no association data are available regarding the consumption of black tea. The present study aimed to identify genetic variants associated with black tea consumption in 12,258 Japanese participants. Data on black tea consumption were collected by a self-administered questionnaire, and genotype data were obtained from a single nucleotide polymorphism array. In the discovery GWAS, two loci met suggestive significance (p < 1.0 × 10-6). Three genetic variants (rs2074356, rs144504271, and rs12231737) at 12q24 locus were also significantly associated with black tea consumption in the replication stage (p < 0.05) and during the meta-analysis (p < 5.0 × 10-8). The association of rs2074356 with black tea consumption was slightly attenuated by the additional adjustment for alcohol drinking frequency. In conclusion, genetic variants at the 12q24 locus were associated with black tea consumption in Japanese populations, and the association is at least partly mediated by alcohol drinking frequency.

Genetic variants and underlying mechanisms influencing variance heterogeneity in maize.

Traditional genetic studies focus on identifying genetic variants associated with the mean difference in a quantitative trait. Because genetic variants also influence phenotypic variation via heterogeneity, we conducted a variance-heterogeneity genome-wide association study to examine the contribution of variance heterogeneity to oil-related quantitative traits. We identified 79 unique variance-controlling single nucleotide polymorphisms (vSNPs) from the sequences of 77 candidate variance-heterogeneity genes for 21 oil-related traits using the Levene test (P < 1.0 × 10-5 ). About 30% of the candidate genes encode enzymes that work in lipid metabolic pathways, most of which define clear expression variance quantitative trait loci. Of the vSNPs specifically associated with the genetic variance heterogeneity of oil concentration, 89% can be explained by additional linked mean-effects genetic variants. Furthermore, we demonstrated that gene × gene interactions play important roles in the formation of variance heterogeneity for fatty acid compositional traits. The interaction pattern was validated for one gene pair (GRMZM2G035341 and GRMZM2G152328) using yeast two-hybrid and bimolecular fluorescent complementation analyses. Our findings have implications for uncovering the genetic basis of hidden additive genetic effects and epistatic interaction effects, and we indicate opportunities to stabilize efficient breeding and selection of high-oil maize (Zea mays L.).

A high density SLAF-seq SNP genetic map and QTL for seed size, oil and protein content in upland cotton.

BACKGROUND: Cotton is a leading natural fiber crop. Beyond its fiber, cottonseed is a valuable source of plant protein and oil. Due to the much higher value of cotton fiber, there is less consideration of cottonseed quality despite its potential value. Though some QTL controlling cottonseed quality have been identified, few of them that warrant further study are known. Identifying stable QTL controlling seed size, oil and protein content is necessary for improvement of cottonseed quality. RESULTS: In this study, a recombinant inbred line (RIL) population was developed from a cross between upland cotton cultivars/lines Yumian 1 and M11. Specific locus amplified fragment sequencing (SLAF-seq) technology was used to construct a genetic map that covered 3353.15 cM with an average distance between consecutive markers of 0.48 cM. The seed index, together with kernel size, oil and protein content were further used to identify QTL. In total, 58 QTL associated with six traits were detected, including 13 stable QTL detected in all three environments and 11 in two environments. CONCLUSION: A high resolution genetic map including 7033 SNP loci was constructed through specific locus amplified fragment sequencing technology. A total of 13 stable QTL associated with six cottonseed quality traits were detected. These stable QTL have the potential for fine mapping, identifying candidate genes, elaborating molecular mechanisms of cottonseed development, and application in cotton breeding programs.

Identification of a gene encoding polygalacturonase expressed specifically in short styles in distylous common buckwheat (Fagopyrum esculentum).

Common buckwheat (Fagopyrum esculentum) is a heteromorphic self-incompatible (SI) species with two types of floral architecture: thrum (short style) and pin (long style). The floral morphology and intra-morph incompatibility are controlled by a single genetic locus, S. However, the molecular mechanisms underlying the heteromorphic self-incompatibility of common buckwheat remain unclear. To identify these mechanisms, we performed proteomic, quantitative reverse-transcription PCR, and linkage analyses. Comparison of protein profiles between the long and short styles revealed a protein unique to the short style. Amino-acid sequencing revealed that it was a truncated form of polygalacturonase (PG); we designated the gene encoding this protein FePG1. Phylogenetic analysis classified FePG1 into the same clade as PGs that function in pollen development and floral morphology. FePG1 expression was significantly higher in short styles than in long styles. It was expressed in flowers of a short-homostyle line but not in flowers of a long-homostyle line. Linkage analysis indicated that FePG1 was not linked to the S locus; it could be a factor downstream of this locus. Our finding of a gene putatively working under the regulation of the S locus provides useful information for elucidation of the mechanism of heteromorphic self-incompatibility.

Metabolic and skeletal homeostasis are maintained in full locus GPRC6A knockout mice.

The G protein-coupled receptor class C, group 6, subtype A (GPRC6A) is suggested to have a physiological function in glucose and bone metabolism, although the precise role lacks consensus due to varying findings in different knockout (KO) mouse models and inconsistent findings on the role of osteocalcin, a proposed GPRC6A agonist. We have further characterized a full locus GPRC6A KO model with respect to energy metabolism, including a long-term high-dose glucocorticoid metabolic challenge. Additionally, we analyzed the microarchitecture of tibiae from young, middle-aged and aged GPRC6A KO mice and wildtype (WT) littermates. Compared to WT, vehicle-treated KO mice presented with normal body composition, unaltered insulin sensitivity and basal serum insulin and glucose levels. Corticosterone (CS) treatment resulted in insulin resistance, abnormal fat accrual, loss of lean mass and suppression of serum osteocalcin levels in both genotypes. Interestingly, serum osteocalcin and skeletal osteocalcin mRNA levels were significantly lower in vehicle-treated GPRC6A KO mice compared to WT animals. However, WT and KO age groups did not differ in long bone mass and structure assessed by micro-computed tomography. We conclude that GPRC6A is not involved in glucose metabolism under normal physiological conditions, nor does it mediate glucocorticoid-induced dysmetabolism in mice. Moreover, GPRC6A does not appear to possess a direct, non-compensable role in long bone microarchitecture under standard conditions.

Detection of major loci associated with the variation of 18 important agronomic traits between Solanum pimpinellifolium and cultivated tomatoes.

Wild species can be used to improve various agronomic traits in cultivars; however, a limited understanding of the genetic basis underlying the morphological differences between wild and cultivated species hinders the integration of beneficial traits from wild species. In the present study, we generated and sequenced recombinant inbred lines (RILs, 201 F10 lines) derived from a cross between Solanum pimpinellifolium and Solanum lycopersicum tomatoes. Based on a high-resolution recombination bin map to uncover major loci determining the phenotypic variance between wild and cultivated tomatoes, 104 significantly associated loci were identified for 18 agronomic traits. On average, these loci explained ~39% of the phenotypic variance of the RILs. We further generated near-isogenic lines (NILs) for four identified loci, and the lines exhibited significant differences for the associated traits. We found that two loci could improve the flower number and inflorescence architecture in the cultivar following introgression of the wild-species alleles. These findings allowed us to construct a trait-locus network to help explain the correlations among different traits based on the pleiotropic or linked loci. Our results provide insights into the morphological changes between wild and cultivated tomatoes, and will help to identify key genes governing important agronomic traits for the molecular selection of elite tomato varieties.

Endogenous Locus Reporter Assays.

Reporter gene assays are widely used in high-throughput screening (HTS) to identify compounds that modulate gene expression. Traditionally a reporter gene assay is built by cloning an endogenous promoter sequence or synthetic response elements in the regulatory region of a reporter gene to monitor transcriptional activity of a specific biological process (exogenous reporter assay). In contrast, an endogenous locus reporter has a reporter gene inserted in the endogenous gene locus that allows the reporter gene to be expressed under the control of the same regulatory elements as the endogenous gene, thus more accurately reflecting the changes seen in the regulation of the actual gene. In this chapter, we introduce some of the considerations behind building a reporter gene assay for high-throughput compound screening and describe the methods we have utilized to establish 1536-well format endogenous locus reporter and exogenous reporter assays for the screening of compounds that modulate Myc pathway activity.

Dissection of the genetic architecture of three seed-quality traits and consequences for breeding in Brassica napus.

Genome-wide association studies (GWASs) combining high-throughput genome resequencing and phenotyping can accelerate the dissection of genetic architecture and identification of genes for plant complex traits. In this study, we developed a rapeseed genomic variation map consisting of 4 542 011 SNPs and 628 666 INDELs. GWAS was performed for three seed-quality traits, including erucic acid content (EAC), glucosinolate content (GSC) and seed oil content (SOC) using 3.82 million polymorphisms in an association panel. Six, 49 and 17 loci were detected to be associated with EAC, GSC and SOC in multiple environments, respectively. The mean total contribution of these loci in each environment was 94.1% for EAC and 87.9% for GSC, notably higher than that for SOC (40.1%). A high correlation was observed between phenotypic variance and number of favourable alleles for associated loci, which will contribute to breeding improvement by pyramiding these loci. Furthermore, candidate genes were detected underlying associated loci, based on functional polymorphisms in gene regions where sequence variation was found to correlate with phenotypic variation. Our approach was validated by detection of well-characterized FAE1 genes at each of two major loci for EAC on chromosomes A8 and C3, along with MYB28 genes at each of three major loci for GSC on chromosomes A9, C2 and C9. Four novel candidate genes were detected by correlation between GSC and SOC and observed sequence variation, respectively. This study provides insights into the genetic architecture of three seed-quality traits, which would be useful for genetic improvement of B. napus.

Synteny analysis of genes and distribution of loci controlling oil content and fatty acid profile based on QTL alignment map in Brassica napus.

BACKGROUND: Deciphering the genetic architecture of a species is a good way to understand its evolutionary history, but also to tailor its profile for breeding elite cultivars with desirable traits. Aligning QTLs from diverse population in one map and utilizing it for comparison, but also as a basis for multiple analyses assure a stronger evidence to understand the genetic system related to a given phenotype. RESULTS: In this study, 439 genes involved in fatty acid (FA) and triacylglycerol (TAG) biosyntheses were identified in Brassica napus. B. napus genome showed mixed gene loss and insertion compared to B. rapa and B. oleracea, and C genome had more inserted genes. Identified QTLs for oil (OC-QTLs) and fatty acids (FA-QTLs) from nine reported populations were projected on the physical map of the reference genome "Darmor-bzh" to generate a map. Thus, 335 FA-QTLs and OC-QTLs could be highlighted and 82 QTLs were overlapping. Chromosome C3 contained 22 overlapping QTLs with all trait studied except for C18:3. In total, 218 candidate genes which were potentially involved in FA and TAG were identified in 162 QTLs confidence intervals and some of them might affect many traits. Also, 76 among these candidate genes were found inside 57 overlapping QTLs, and candidate genes for oil content were in majority (61/76 genes). Then, sixteen genes were found in overlapping QTLs involving three populations, and the remaining 60 genes were found in overlapping QTLs of two populations. Interaction network and pathway analysis of these candidate genes indicated ten genes that might have strong influence over the other genes that control fatty acids and oil formation. CONCLUSION: The present results provided new information for genetic basis of FA and TAG formation in B. napus. A map including QTLs from numerous populations was built, which could serve as reference to study the genome profile of B. napus, and new potential genes emerged which might affect seed oil. New useful tracks were showed for the selection of population or/and selection of interesting genes for breeding improvement purpose.

A disulfide bond A-like oxidoreductase is a strong candidate gene for self-incompatibility in apricot (Prunus armeniaca) pollen.

S-RNase based gametophytic self-incompatibility (SI) is a widespread prezygotic reproductive barrier in flowering plants. In the Solanaceae, Plantaginaceae and Rosaceae gametophytic SI is controlled by the pistil-specific S-RNases and the pollen S-locus F-box proteins but non-S-specific factors, namely modifiers, are also required. In apricot, Prunus armeniaca (Rosaceae), we previously mapped two pollen-part mutations that confer self-compatibility in cultivars Canino and Katy at the distal end of chromosome 3 (M-locus) unlinked to the S-locus. Here, we used high-resolution mapping to identify the M-locus with an ~134 kb segment containing ParM-1-16 genes. Gene expression analysis identified four genes preferentially expressed in anthers as modifier gene candidates, ParM-6, -7, -9 and -14. Variant calling of WGS Illumina data from Canino, Katy, and 10 self-incompatible cultivars detected a 358 bp miniature inverted-repeat transposable element (MITE) insertion in ParM-7 shared only by self-compatible apricots, supporting ParM-7 as strong candidate gene required for SI. ParM-7 encodes a disulfide bond A-like oxidoreductase protein, which we named ParMDO. The MITE insertion truncates the ParMDO ORF and produces a loss of SI function, suggesting that pollen rejection in Prunus is dependent on redox regulation. Based on phylogentic analyses we also suggest that ParMDO may have originated from a tandem duplication followed by subfunctionalization and pollen-specific expression.

Use of genotyping-by-sequencing to determine the genetic structure in the medicinal plant chamomile, and to identify flowering time and alpha-bisabolol associated SNP-loci by genome-wide association mapping.

BACKGROUND: Chamomile (Matricaria recutita L.) has a long history of use in herbal medicine with various applications, and the flower heads contain numerous secondary metabolites which are medicinally active. In the major crop plants, next generation sequencing (NGS) approaches are intensely applied to exploit genetic resources, to develop genomic resources and to enhance breeding. Here, genotyping-by-sequencing (GBS) has been used in the non-model medicinal plant chamomile to evaluate the genetic structure of the cultivated varieties/populations, and to perform genome wide association study (GWAS) focusing on genes with large effect on flowering time and the medicinally important alpha-bisabolol content. RESULTS: GBS analysis allowed the identification of 6495 high-quality SNP-markers in our panel of 91 M. recutita plants from 33 origins (2-4 genotypes each) and 4 M. discoidea plants as outgroup, grown in the greenhouse in Gatersleben, Germany. M. recutita proved to be clearly distinct from the outgroup, as was demonstrated by different cluster and principal coordinate analyses using the SNP-markers. Chamomile genotypes from the same origin were mostly genetically similar. Model-based cluster analysis revealed one large group of tetraploid genotypes with low genetic differentiation including 39 plants from 14 origins. Tetraploids tended to display lower genetic diversity than diploids, probably reflecting their origin by artificial polyploidisation from only a limited set of genetic backgrounds. Analyses of flowering time demonstrated that diploids generally flowered earlier than tetraploids, and the analysis of alpha-bisabolol identified several tetraploid genotypes with a high content. GWAS identified highly significant (P < 0.01) SNPs for flowering time (9) and alpha-bisabolol (71). One sequence harbouring SNPs associated with flowering time was described to play a role in self-pollination in Arabidopsis thaliana, whereas four sequences harbouring SNPs associated with alpha-bisabolol were identified to be involved in plant biotic and abiotic stress response in various plants species. CONCLUSIONS: The first genomic resource for future applications to enhance breeding in chamomile was created, andanalyses of diversity will facilitate the exploitation of these genetic resources. The GWAS data pave the way for future research towards the genetics underlying important traits in chamomile, the identification of marker-trait associations, and development of reliable markers for practical breeding.

Intricacies in arrangement of SNP haplotypes suggest "Great Admixture" that created modern humans.

BACKGROUND: Inferring history from genomic sequences is challenging and problematic because chromosomes are mosaics of thousands of small Identicalby-descent (IBD) fragments, each of them having their own unique story. However, the main events in recent evolution might be deciphered from comparative analysis of numerous loci. A paradox of why humans, whose effective population size is only 104, have nearly three million frequent SNPs is formulated and examined. RESULTS: We studied 5398 loci evenly covering all human autosomes. Common haplotypes built from frequent SNPs that are present in people from various populations have been examined. We demonstrated highly non-random arrangement of alleles in common haplotypes. Abundance of mutually exclusive pairs of common haplotypes that have different alleles at every polymorphic position (so-called Yin/Yang haplotypes) was found in 56% of loci. A novel widely spread category of common haplotypes named Mosaic has been described. Mosaic consists of numerous pieces of Yin/Yang haplotypes and represents an ancestral stage of one of them. Scenarios of possible appearance of large number of frequent human SNPs and their habitual arrangement in Yin/Yang common haplotypes have been evaluated with an advanced genomic simulation algorithm. CONCLUSIONS: Computer modeling demonstrated that the observed arrangement of 2.9 million frequent SNPs could not originate from a sole stand-alone population. A "Great Admixture" event has been proposed that can explain peculiarities with frequent SNP distributions. This Great Admixture presumably occurred 100-300 thousand years ago between two ancestral populations that had been separated from each other about a million years ago. Our programs and algorithms can be applied to other species to perform evolutionary and comparative genomics.

Molecular characterization of aspartylglucosaminidase, a lysosomal hydrolase upregulated during strobilation in the moon jellyfish, Aurelia aurita.

The life cycle of the moon jellyfish, Aurelia aurita, alternates between a benthic asexual polyp stage and a planktonic sexual medusa (jellyfish) stage. Transition from polyp to medusa is called strobilation. To investigate the molecular mechanisms of strobilation, we screened for genes that are upregulated during strobilation using the differential display method and we identified aspartylglucosaminidase (AGA), which encodes a lysosomal hydrolase. Similar to AGAs from other species, Aurelia AGA possessed an N-terminal signal peptide and potential N-glycosylation sites. The genomic region of Aurelia AGA was approximately 9.8 kb in length and contained 12 exons and 11 introns. Quantitative RT-PCR analysis revealed that AGA expression increased during strobilation, and was then decreased in medusae. To inhibit AGA function, we administered the lysosomal acidification inhibitors, chloroquine or bafilomycin A1, to animals during strobilation. Both inhibitors disturbed medusa morphogenesis at the oral end, suggesting involvement of lysosomal hydrolases in strobilation.

Isolation and characterization of 14 tetranucleotide microsatellite loci for the cannonball jellyfish (Stomolophus sp.) by next generation sequencing.

The Cannonball jellyfish (Stomolophus sp.) is a species of jellyfish with high relevance in artisanal fishing. Studies of their populations do not extend beyond the morphological descriptions knowing that presents a great morphological variability. However, there are no genetic studies to determine the number of independent populations, so microsatellite markers become a suitable option. Since there are no species-specific microsatellite loci, in this paper, 14 new microsatellite loci are characterized. Microsatellite loci were isolated de novo through next generation sequencing, by two runs on Illumina MiSeq. A total of 506,771,269 base pair were obtained, from which 142,616 were microsatellite loci, and 1546 of them could design primers. We tested 14 primer pairs on 32 individuals from Bahía de La Paz, Gulf of California. We observed low genetic variation among loci (mean number of alleles per locus = 4.33, mean observed heterozygosity 0.381, mean expected heterozygosity 0.501). These loci are the first ones described for the species and will be helpful to carry out genetic diversity and population genetics studies.