Genetic Diversity and Population Structure in Türkiye Bread Wheat Genotypes Revealed by Simple Sequence Repeats (SSR) Markers

Wheat genotypes should be improved through available germplasm genetic diversity to ensure food security. This study investigated the molecular diversity and population structure of a set of Türkiye bread wheat genotypes using 120 microsatellite markers. Based on the results, 651 polymorphic alleles were evaluated to determine genetic diversity and population structure. The number of alleles ranged from 2 to 19, with an average of 5.44 alleles per locus. Polymorphic information content (PIC) ranged from 0.031 to 0.915 with a mean of 0.43. In addition, the gene diversity index ranged from 0.03 to 0.92 with an average of 0.46. The expected heterozygosity ranged from 0.00 to 0.359 with a mean of 0.124. The unbiased expected heterozygosity ranged from 0.00 to 0.319 with an average of 0.112. The mean values of the number of effective alleles (Ne), genetic diversity of Nei (H) and Shannon's information index (I) were estimated at 1.190, 1.049 and 0.168, respectively. The highest genetic diversity (GD) was estimated between genotypes G1 and G27. In the UPGMA dendrogram, the 63 genotypes were grouped into three clusters. The three main coordinates were able to explain 12.64, 6.38 and 4.90% of genetic diversity, respectively. AMOVA revealed diversity within populations at 78% and between populations at 22%. The current populations were found to be highly structured. Model-based cluster analyses classified the 63 genotypes studied into three subpopulations. The values of F-statistic (Fst) for the identified subpopulations were 0.253, 0.330 and 0.244, respectively. In addition, the expected values of heterozygosity (He) for these sub-populations were recorded as 0.45, 0.46 and 0.44, respectively. Therefore, SSR markers can be useful not only in genetic diversity and association analysis of wheat but also in its germplasm for various agronomic traits or mechanisms of tolerance to environmental stresses.

cgSSR marker-based genome-wide association study identified genomic regions for panicle characters and yield in rice (Oryza sativa L.)

BACKGROUND

To improve production efficiency, positive alleles corresponding to yield-related attributes must be accumulated in a single elite background. We designed and used cgSSR markers, which are superior to random SSR markers in genome-wide association study, to identify genomic regions that contribute to panicle characters and grain yield in this study.

RESULTS

As evidenced by the high polymorphic information content value and gene diversity coefficient, the new cgSSR markers were determined to be highly informative. These cgSSR markers were employed to generate genotype data for an association panel evaluated for four panicle characters and grain yield over three seasons. For five traits, 17 significant marker-trait associations on six chromosomes were discovered. The percentage of phenotypic variance that could be explained ranged from 4% to 13%. Unrelated gene-derived markers had a strong association with target traits as well.

CONCLUSION

Trait-associated cgSSR markers derived from corresponding or related genes ensure their utility in direct allele selection, while other linked markers aid in allele selection indirectly by altering the phenotype of interest. Through a marker-assisted breeding approach, these marker-trait associations can be leveraged to accumulate favourable alleles for yield enhancement in rice. © 2022 Society of Chemical Industry.

Parvalbumin: A Major Fish Allergen and a Forensically Relevant Marker

Parvalbumins (PVALBs) are low molecular weight calcium-binding proteins. In addition to their role in many biological processes, PVALBs play an important role in regulating Ca2+ switching in muscles with fast-twitch fibres in addition to their role in many biological processes. The PVALB gene family is divided into two gene types, alpha (α) and beta (β), with the β gene further divided into two gene types, beta1 (β1) and beta2 (β2), carrying traces of whole genome duplication. A large variety of commonly consumed fish species contain PVALB proteins which are known to cause fish allergies. More than 95% of all fish-induced food allergies are caused by PVALB proteins. The authentication of fish species has become increasingly important as the seafood industry continues to grow and the growth brings with it many cases of food fraud. Since the PVALB gene plays an important role in the initiation of allergic reactions, it has been used for decades to develop alternate assays for fish identification. A brief review of the significance of the fish PVALB genes is presented in this article, which covers evolutionary diversity, allergic properties, and potential use as a forensic marker.

A comparison of 25 complete chloroplast genomes between sister mangrove species Kandelia obovata and Kandelia candel geographically separated by the South China Sea

In 2003, Kandelia obovata was identified as a new mangrove species differentiated from Kandelia candel. However, little is known about their chloroplast (cp) genome differences and their possible ecological significance. In this study, 25 whole cp genomes, with seven samples of K. candel from Malaysia, Thailand, and Bangladesh and 18 samples of K. obovata from China, were sequenced for comparison. The cp genomes of both species encoded 128 genes, namely 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes, but the cp genome size of K. obovata was ~2 kb larger than that of K. candle due to the presence of more and longer repeat sequences. Of these, tandem repeats and simple sequence repeats exhibited great differences. Principal component analysis based on indels, and phylogenetic tree analyses constructed with homologous protein genes from the single-copy genes, as well as 38 homologous pair genes among 13 mangrove species, gave strong support to the separation of the two species within the Kandelia genus. Homologous genes ndhD and atpA showed intraspecific consistency and interspecific differences. Molecular dynamics simulations of their corresponding proteins, NAD(P)H dehydrogenase chain 4 (NDH-D) and ATP synthase subunit alpha (ATP-A), predicted them to be significantly different in the functions of photosynthetic electron transport and ATP generation in the two species. These results suggest that the energy requirement was a pivotal factor in their adaptation to differential environments geographically separated by the South China Sea. Our results also provide clues for future research on their physiological and molecular adaptation mechanisms to light and temperature.

The Light Chain Domain and Especially the C-Terminus of Receptor-Binding Domain of the Botulinum Neurotoxin (BoNT) Are the Hotspots for Amino Acid Variability and Toxin Type Diversity

Botulinum neurotoxins (BoNT) are the most potent toxins in the world. They are produced by a few dozens of strains within several clostridial species. The toxin that they produce can cause botulism, a flaccid paralysis in humans and other animals. With seven established serologically different types and over 40 subtypes, BoNTs are among the most diverse known toxins. The toxin, its structure, its function and its physiological effects on the neural cell and animal hosts along with its diversity have been the subjects of numerous studies. However, many gaps remain in our knowledge about the BoNT toxin and the species that produce them. One of these gaps involves the distribution and extent of variability along the full length of the gene and the protein as well as its domains and subdomains. In this study, we performed an extensive analysis of all of the available 143 unique BoNT-encoding genes and their products, and we investigated their diversity and evolution. Our results indicate that while the nucleotide variability is almost uniformly distributed along the entire length of the gene, the amino acid variability is not. We found that most of the differences were concentrated along the protein's light chain (LC) domain and especially, the C-terminus of the receptor-binding domain (H_CC). These two regions of the protein are thus identified as the main source of the toxin type differentiation, and consequently, this toxin's versatility to bind different receptors and their isoforms and act upon different substrates, thus infecting different hosts.

Novel GBS-Based SNP Markers for Finger Millet and Their Use in Genetic Diversity Analyses

Eleusine coracana (L.) Gaertn., commonly known as finger millet, is a multipurpose crop used for food and feed. Genomic tools are required for the characterization of crop gene pools and their genomics-led breeding. High-throughput sequencing-based characterization of finger millet germplasm representing diverse agro-ecologies was considered an effective method for determining its genetic diversity, thereby suggesting potential candidates for breeding. In this study, the genotyping-by-sequencing (GBS) method was used to simultaneously identify novel single nucleotide polymorphism (SNP) markers and genotype 288 finger millet accessions collected from Ethiopia and Zimbabwe. The accessions were characterized at individual and group levels using 5,226 bi-allelic SNPs, with a minimum allele frequency (MAF) of above 0.05, distributed across 2,500 scaffolds of the finger millet reference genome. The polymorphism information content (PIC) of the SNPs was 0.23 on average, and a quarter of them have PIC values over 0.32, making them highly informative. The grouping of the 288 accessions into seven populations based on geographic proximity and the potential for germplasm exchange revealed a narrow range of observed heterozygosity (Ho; 0.09–0.11) and expected heterozygosity (He) that ranged over twofold, from 0.11 to 0.26. Alleles unique to the different groups were also identified, which merit further investigation for their potential association with desirable traits. The analysis of molecular variance (AMOVA) revealed a highly significant genetic differentiation among groups of accessions classified based on the geographic region, country of origin, days to flowering, panicle type, and Al tolerance (p < 0.01). The high genetic differentiation between Ethiopian and Zimbabwean accessions was evident in the AMOVA, cluster, principal coordinate, and population structure analyses. The level of genetic diversity of finger millet accessions varies moderately among locations within Ethiopia, with accessions from the northern region having the lowest level. In the neighbor-joining cluster analysis, most of the improved cultivars included in this study were closely clustered, probably because they were developed using genetically less diverse germplasm and/or selected for similar traits, such as grain yield. The recombination of alleles via crossbreeding genetically distinct accessions from different regions of the two countries can potentially lead to the development of superior cultivars.

Evaluation of the genetic structure of Bromus inermis populations from chemically and radioactively polluted areas using microsatellite markers from closely related species

Hypothesis The ecotoxicological and radiobiological effects can be manifested in a decrease in genetic diversity with an increase in toxic and radiation load, in an increase in the frequencies of rare and/or unique (private) alleles in impact samples, and in a decrease in the differentiation of B. inermis populations within each pollution area.Materials and methods We have selected a collection of primers for Bromus inermis, consisting of 21 microsatellite (SSR) loci from B. sterilis, B. tectorum and Triticum aestivum. The level of toxic load (chemically polluted area) was 4-19 conventional units, and the absorbed dose rate (the Kyshtym accident area) varied from 0.153 to 21.5 μGy h^-1, which is up to two orders higher than the natural background radiation level (≈ 0.1 μGy h^-1).Results Only eight of 21 (38%) of SSR primers showed good transferability and were used for B. inermis population studies from areas of technogenic pollution (heavy metals and radionuclides). We revealed 42 alleles at eight loci, and the number of alleles per locus varied from one to 13 in B. inermis populations. The percentage of polymorphic loci in B. inermis populations was 48.44%, the polymorphism information content (PIC) value was 0.556, and Shannon information index was 0.69 ± 0.3. A total of 22 rare, 14 private and 9 both rare and private alleles were reported for all B. inermis populations. There were no correlations between geographic and genetic distances. Only 6.8% of the genetic variability was distributed among B. inermis populations.Conclusion There was no decrease in genetic diversity ("genetic erosion") found in B. inermis populations growing for a long time under anthropogenic stress. No significant differences in the number of rare and private alleles in the background and impact populations of B. inermis were found. The smooth brome is characterized by low differentiation of the populations. Possible reasons for this phenomenon are discussed.

Mining Indole Alkaloid Synthesis Gene Clusters from Genomes of 53 Claviceps Strains Revealed Redundant Gene Copies and an Approximate Evolutionary Hourglass Model

Ergot fungi (Claviceps spp.) are infamous for producing sclerotia containing a wide spectrum of ergot alkaloids (EA) toxic to humans and animals, making them nefarious villains in the agricultural and food industries, but also treasures for pharmaceuticals. In addition to three classes of EAs, several species also produce paspaline-derived indole diterpenes (IDT) that cause ataxia and staggers in livestock. Furthermore, two other types of alkaloids, i.e., loline (LOL) and peramine (PER), found in Epichloë spp., close relatives of Claviceps, have shown beneficial effects on host plants without evidence of toxicity to mammals. The gene clusters associated with the production of these alkaloids are known. We examined genomes of 53 strains of 19 Claviceps spp. to screen for these genes, aiming to understand the evolutionary patterns of these genes across the genus through phylogenetic and DNA polymorphism analyses. Our results showed (1) varied numbers of eas genes in C. sect. Claviceps and sect. Pusillae, none in sect. Citrinae, six idt/ltm genes in sect. Claviceps (except four in C. cyperi), zero to one partial (idtG) in sect. Pusillae, and four in sect. Citrinae, (2) two to three copies of dmaW, easE, easF, idt/ltmB, itd/ltmQ in sect. Claviceps, (3) frequent gene gains and losses, and (4) an evolutionary hourglass pattern in the intra-specific eas gene diversity and divergence in C. purpurea.

Conservation Priorities Analysis of Chinese Indigenous Pig Breeds in the Taihu Lake Region

Most indigenous pig resources are known to originate from China. Thus, establishing conservation priorities for these local breeds is very essential, especially in the case of limited conservation funds. Therefore, in this study, we analyzed 445 individuals belonging to six indigenous breeds from the Taihu Lake Region, using a total of 131,300 SNPs. In order to determine the long-term guidelines for the management of these breeds, we analyzed the level of diversity in the metapopulation following a partition of diversity within and between breed subpopulations, using both measures of genic and allelic diversity. From the study, we found that the middle Meishan (MMS) pig population contributes the most (22%) to the total gene diversity while the Jiaxing black (JX) pig population contributes the most (27%) to the gene diversity between subpopulations. Most importantly, when we consider one breed is removed from the meta-population, the first two breeds prioritized should be JX pig breed and Fengjing pig breed followed by small Meishan (SMS), Mizhu (MI), and Erhualian (EH) if we pay more attention to the gene diversity between subpopulations. However, if the priority focus is on the total gene diversity, then the first breed to be prioritized would be the Shawutou (SW) pig breed followed by JX, MI, EH, and Fengjing (FJ). Furthermore, we noted that if conservation priority is to be based on the allelic diversity between subpopulations, then the MI breed should be the most prioritized breed followed by SW, Erhuanlian, and MMS. Summarily, our data show that different breeds have different contributions to the gene and allelic diversity within subpopulations as well as between subpopulations. Our study provides a basis for setting conservation priorities for indigenous pig breeds with a focus on different priority criteria.

Molecular Characterization of the Enterohemolysin Gene (ehxA) in Clinical Shiga Toxin-Producing Escherichia coli Isolates

Shiga toxin (Stx)-producing Escherichia coli (STEC) is an important foodborne pathogen with the ability to cause bloody diarrhea (BD) and hemolytic uremic syndrome (HUS). Little is known about enterohemolysin-encoded by ehxA. Here we investigated the prevalence and diversity of ehxA in 239 STEC isolates from human clinical samples. In total, 199 out of 239 isolates (83.26%) were ehxA positive, and ehxA was significantly overrepresented in isolates carrying stx_2a + stx_2c (p < 0.001) and eae (p < 0.001). The presence of ehxA was significantly associated with BD and serotype O157:H7. Five ehxA subtypes were identified, among which, ehxA subtypes B, C, and F were overrepresented in eae-positive isolates. All O157:H7 isolates carried ehxA subtype B, which was related to BD and HUS. Three ehxA groups were observed in the phylogenetic analysis, namely, group Ⅰ (ehxA subtype A), group Ⅱ (ehxA subtype B, C, and F), and group Ⅲ (ehxA subtype D). Most BD- and HUS-associated isolates were clustered into ehxA group Ⅱ, while ehxA group Ⅰ was associated with non-bloody stool and individuals ≥10 years of age. The presence of ehxA + eae and ehxA + eae + stx₂ was significantly associated with HUS and O157:H7 isolates. In summary, this study showed a high prevalence and the considerable genetic diversity of ehxA among clinical STEC isolates. The ehxA genotypes (subtype B and phylogenetic group Ⅱ) could be used as risk predictors, as they were associated with severe clinical symptoms, such as BD and HUS. Furthermore, ehxA, together with stx and eae, can be used as a risk predictor for HUS in STEC infections.

[Research progress and applications of strain analysis based on metagenomic data]

Strain is the fundamental unit in microbial taxonomy. The functional diversity among strains has great influence on host phenotypes. With the development of microbiome research, knowing the composition and functional capacities of complex microbial communities at the strain level has become increasingly valuable in scientific research and clinical applications. This review introduces the principles of bioinformatics algorithms for strain analysis based on metagenomic data, the applications in microbiome research and directions of future development.

Molecular characteristics of eae-positive clinical Shiga toxin-producing Escherichia coli in Sweden

Shiga toxin (Stx)-producing Escherichia coli (STEC) can cause a wide range of symptoms from asymptomatic carriage, mild diarrhea to bloody diarrhea (BD) and hemolytic uremic syndrome (HUS). Intimin, encoded by the eae gene, also plays a critical role in STEC pathogenesis. Herein, we investigated the prevalence and genetic diversity of eae among clinical STEC isolates from patients with diarrhea, BD, HUS as well as from asymptomatic STEC-positive individuals in Sweden with whole-genome sequencing. We found that 173 out of 239 (72.4%) of clinical STEC strains were eae positive. Six eae subtypes (ϵ1, γ1, β3, θ, ζ and ρ) were identified eae and its subtype γ1 were significantly overrepresented in O157:H7 strains isolated from BD and HUS patients. ϵ1 was associated with O121:H19 and O103:H2 strains, and β3 to O26:H11 strains. The combination of eae subtype γ1 and stx subtype (stx 2 or stx 1+stx 2) is more likely to cause severe disease, suggesting the possibility of using eae genotypes in risk assessment of STEC infection. In summary, this study demonstrated a high prevalence of eae in clinical STEC strains and considerable genetic diversity of eae in STEC strains in Sweden from 1994 through 2018, and revealed association between eae subtypes and disease severity.

Evolutionary Pattern of Interferon Alpha Genes in Bovidae and Genetic Diversity of IFNAA in the Bovine Genome

Interferons are secretory proteins induced in response to specific extracellular stimuli which stimulate intra- and intercellular networks for regulating innate and acquired immunity, resistance to viral infections, and normal and tumor cell survival and death. Type 1 interferons plays a major role in the CD8 T-cell response to viral infection. The genomic analysis carried out here for type I interferons within Bovidae family shows that cattle, bison, water buffalo, goat, and sheep (all Bovidae), have different number of genes of the different subtypes, with a large increase in the numbers, compared to human and mouse genomes. A phylogenetic analysis of the interferon alpha (IFNA) proteins in this group shows that the genes do not follow the evolutionary pattern of the species, but rather a cycle of duplications and deletions in the different species. In this study we also studied the genetic diversity of the bovine interferon alpha A (IFNAA), as an example of the IFNA genes in cattle, sequencing a fragment of the coding sequence in 18 breeds of cattle from Pakistan, Nigeria and USA. Similarity analysis allowed the allocation of sequences into 22 haplotypes. Bhagnari, Brangus, Sokoto Gudali, and White Fulani, had the highest number of haplotypes, while Angus, Hereford and Nari Master had the least. However, when analyzed by the average haplotype count, Angus, Bhagnari, Hereford, Holstein, Muturu showed the highest values, while Cholistani, Lohani, and Nari Master showed the lowest values. Haplotype 4 was found in the highest number of individuals (74), and in 15 breeds. Sequences for yak, bison, and water buffalo, were included within the bovine haplotypes. Medium Joining network showed that the sequences could be divided into 4 groups: one with highly similar haplotypes containing mostly Asian and African breeds, one with almost all of the Bos taurus American breeds, one mid-diverse group with mostly Asian and African sequences, and one group with highly divergent haplotypes with five N'Dama sequences and one from each of White Fulani, Dhanni, Tharparkar, and Bhagnari. The large genetic diversity found in IFNAA could be a very good indication of the genetic variation among the different genes of IFNA and could be an adaptation for these species in response to viral challenges they face.

Partitioning heterozygosity in subdivided populations: Some misuses of Nei's decomposition and an alternative probabilistic approach

Nei's decomposition of total expected heterozygosity in subdivided populations into within- and between-subpopulation components, H_S and D_ST , respectively, is a classical tool in the conservation and management of genetic resources. Reviewing why this is not a decomposition into independent terms of within- and between-subpopulation gene diversity, we illustrate how this approach can be misleading because it overemphasizes the within-subpopulation component compared to Jost's nonadditive decomposition based on gene diversity indices. Using probabilistic partitioning of the total expected heterozygosity into independent within- and between-subpopulation contributions, we show that the contribution of the within-subpopulation expected heterozygosity to the total expected heterozygosity is not H_S , as suggested by Nei's decomposition, but H_S /s, with s being the number of subpopulations. Finally, we compare three possible approaches of decomposing total heterozygosity in subdivided populations (i.e., Nei's decomposition, Jost's approach, and probabilistic partitioning) with regard to independence between terms and sensitivity to unequal subpopulation sizes. For the conservation and management of genetic resources, we recommend using probabilistic partitioning and Jost's differentiation parameter rather than Nei's decomposition.

Disentangling the effects of geographic peripherality and habitat suitability on neutral and adaptive genetic variation in Swiss stone pine

It is generally accepted that the spatial distribution of neutral genetic diversity within a species' native range mostly depends on effective population size, demographic history, and geographic position. However, it is unclear how genetic diversity at adaptive loci correlates with geographic peripherality or with habitat suitability within the ecological niche. Using exome-wide genomic data and distribution maps of the Alpine range, we first tested whether geographic peripherality correlates with four measures of population genetic diversity at > 17,000 SNP loci in 24 Alpine populations (480 individuals) of Swiss stone pine (Pinus cembra) from Switzerland. To distinguish between neutral and adaptive SNP sets, we used four approaches (two gene diversity estimates, F_ST outlier test, and environmental association analysis) that search for signatures of selection. Second, we established ecological niche models for P. cembra in the study range and investigated how habitat suitability correlates with genetic diversity at neutral and adaptive loci. All estimates of neutral genetic diversity decreased with geographic peripherality, but were uncorrelated with habitat suitability. However, heterozygosity (H_e ) at adaptive loci based on Tajima's D declined significantly with increasingly suitable conditions. No other diversity estimates at adaptive loci were correlated with habitat suitability. Our findings suggest that populations at the edge of a species' geographic distribution harbour limited neutral genetic diversity due to demographic properties. Moreover, we argue that populations from suitable habitats went through strong selection processes, are thus well adapted to local conditions, and therefore exhibit reduced genetic diversity at adaptive loci compared to populations at niche margins.

The Landscape of Genetic Content in the Gut and Oral Human Microbiome

Despite substantial interest in the species diversity of the human microbiome and its role in disease, the scale of its genetic diversity, which is fundamental to deciphering human-microbe interactions, has not been quantified. Here, we conducted a cross-study meta-analysis of metagenomes from two human body niches, the mouth and gut, covering 3,655 samples from 13 studies. We found staggering genetic heterogeneity in the dataset, identifying a total of 45,666,334 non-redundant genes (23,961,508 oral and 22,254,436 gut) at the 95% identity level. Fifty percent of all genes were "singletons," or unique to a single metagenomic sample. Singletons were enriched for different functions (compared with non-singletons) and arose from sub-population-specific microbial strains. Overall, these results provide potential bases for the unexplained heterogeneity observed in microbiome-derived human phenotypes. One the basis of these data, we built a resource, which can be accessed at https://microbial-genes.bio.

A Natural Variation of Fumonisin Gene Cluster Associated with Fumonisin Production Difference in Fusarium fujikuroi

Fusarium fujikuroi, a member of the Fusarium fujikuroi species complex, stands out as a rice bakanae disease pathogen with a high production of gibberellic acid. Not all, but some F. fujikuroi strains are known to produce a carcinogenic mycotoxin fumonisin. Fumonisin biosynthesis is dependent on the FUM cluster composed of 16 FUM genes. The FUM cluster was detected not only from a fumonisin producing strain, but also from a fumonisin nonproducing strain that does not produce a detectable level of fumonisin. Genetic mapping indicated the causative mutation(s) of fumonisin nonproduction is present in the FUM cluster of the fumonisin nonproducing strain. Comparative analyses of FUM genes between the fumonisin producing and the nonproducing strains and gene complementation indicated that causative mutation of fumonisin nonproduction is not a single occurrence and the mutations are distributed in FUM21 and FUM7. Our research revealed a natural variation in the FUM cluster involving fumonisin production difference in F. fujikuroi.

Monitoring age-related trends in genomic diversity of Australian lungfish

An important challenge for conservation science is to detect declines in intraspecific diversity so that management action can be guided towards populations or species at risk. The lifespan of Australian lungfish (Neoceratodus forsteri) exceeds 80 years, and human impacts on breeding habitat over the last half century may have impeded recruitment, leaving populations dominated by old postreproductive individuals, potentially resulting in a small and declining breeding population. Here, we conduct a "single-sample" evaluation of genetic erosion within contemporary populations of the Australian lungfish. Genetic erosion is a temporal decline in intraspecific diversity due to factors such as reduced population size and inbreeding. We examined whether young individuals showed signs of reduced genetic diversity and/or inbreeding using a novel bomb radiocarbon dating method to age lungfish nonlethally, based on 14 C ratios of scales. A total of 15,201 single nucleotide polymorphic (SNP) loci were genotyped in 92 individuals ranging in age from 2 to 77 years old. Standardized individual heterozygosity and individual inbreeding coefficients varied widely within and between riverine populations, but neither was associated with age, so perceived problems with recruitment have not translated into genetic erosion that could be considered a proximate threat to lungfish populations. Conservation concern has surrounded Australian lungfish for over a century. However, our results suggest that long-lived threatened species can maintain stable levels of intraspecific variability when sufficient reproductive opportunities exist over the course of a long lifespan.

Phylogenetic analysis of human G9P[8] rotavirus strains circulating in Jiangsu, China between 2010 and 2016

Rotavirus A (RVA) is the leading cause of acute viral gastroenteritis in children under 5 years of age worldwide. G9P[8] is a common RVA genotype that has been persistently prevalent in Jiangsu, China. To determine the genetic diversity of G9P[8] RVAs, 7 representative G9P[8] strains collected from Suzhou Children's Hospital between 2010 and 2016 (named JS2010-JS2016) were analyzed through whole-genome sequencing. All evaluated strains showed the Wa-like constellation G9-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1. Furthermore, phylogenetic analysis revealed that the VP7 genes of all strains clustered into lineage G9-III and G9-VI. With the exception of strain JS2012 (P[8]-4), the VP4 sequences of all strains belonged to the P[8]-3 lineage. Sequencing further revealed that amino acid substitutions were present in the antigenic regions of the VP7 and VP4 genes of all strains. Moreover, there were multiple substitutions in antigenic sites I and II of the nonstructural protein 4 (NSP4) genes, whereas the other NSP genes were relatively conserved. In conclusion, our phylogenetic analysis of these 7 G9P[8] strains suggests that RVA varied across regions and time. Therefore, our findings suggest that continued surveillance is necessary to explore the molecular evolutionary characteristics of RVA for better prevention and treatment of acute viral gastroenteritis.

Dynamic Gene-Resource Landscape Management of Norway Spruce: Combining Utilization and Conservation

Traditional gene-resource management programs for forest trees are long-term endeavors requiring sustained organizational commitment covering extensive landscapes. While successful in maintaining adaptation, genetic diversity and capturing traditional growth attributes gains, these programs are dependent on rigid methods requiring elaborate mating schemes, thus making them slow in coping with climate change challenges. Here, we review the significance of Norway spruce in the boreal region and its current management practices. Next, we discuss opportunities offered by novel technologies and, with the use of computer simulations, we propose and evaluate a dynamic landscape gene-resource management in Norway. Our suggested long-term management approach capitalizes on: (1) existing afforestation activities, natural crosses, and DNA-based pedigree assembly to create structured pedigree for evaluation, thus traditional laborious control crosses are avoided and (2) landscape level genetic evaluation, rather than localized traditional progeny trials, allowing for screening of adapted individuals across multiple environmental gradients under changing climate. These advantages lead to greater genetic response to selection in adaptive traits without the traditional breeding and testing scheme, facilitating conservation of genetic resources within the breeding population of the most important forest tree species in Norway. The use of in situ selection from proven material exposed to realistic conditions over vast territories has not been conducted in forestry before. Our proposed approach is in contrast to worldwide current programs, where genetic evaluation is constrained by the range of environments where testing is conducted, which may be insufficient to capture the broad environmental variation necessary to tackle adaptation under changing climate.

Host Specificity of Flagellins from Segmented Filamentous Bacteria Affects Their Patterns of Interaction with Mouse Ileal Mucosal Proteins

Segmented filamentous bacteria (SFB) are known modulators of the mammalian immune system. Currently, the technology for investigating SFB culture in vitro is immature, and as a result, the mechanisms of SFB colonization and immune regulation are not yet fully elucidated. In this study, we investigated the gene diversity and host specificity of SFB flagellin genes. The fliC1 and fliC2 genes are relatively conserved, while the fliC3 and fliC4 genes are more variable, especially at the central and C-terminal regions. Host specificity analysis demonstrated that the fliC1 genes do not cluster together based on the host organism, whereas the fliC3 and fliC4 genes were host specific at the nucleotide and deduced amino acid levels. SFB flagellin protein expression in the ileum mucosa and cecal contents was detected by using fluorescence in situ hybridization (FISH) combined with immunohistochemical (IHC) analysis, immunoblotting, and liquid chromatography-tandem mass spectrometry (LC-MS/MS). Although the purified SFB FliC3 protein originating from both mouse and rat was able to activate Toll-like receptor 5 (TLR5)-linked NF-κB signaling, no host specificity was observed. Interestingly, the patterns of interaction with mouse ileum mucosal proteins were different for mouse FliC3 (mFliC3) and rat FliC3 (rFliC3). Gene Ontology (GO) and KEGG analyses indicated that more adherence-related proteins interacted with mFliC3, while more lysosome- and proteolysis-related proteins interacted with rFliC3. In vitro degradation experiments indicated that the stability of rFliC3 was lower than that of mFliC3 when they were incubated with mouse ileum mucosal proteins. In summary, the gene diversity and host specificity of SFB flagellin genes were investigated, and SFB flagellin expression was detected in gut samples.IMPORTANCE Since SFB genomes contain only one copy of each FliC gene, the diversity of FliC is representative of SFB strain diversity. Currently, little is known regarding the diversity and specificity of members of the group of SFB. The work presented herein demonstrates that select SFB strains, exhibiting unique FliC patterns, are present in a variety of mammalian hosts. SFB fliC genes were found to interact with a number of unique targets, providing further evidence for SFB host selection. Together, this work represents a major advancement in identifying SFB and delineating how members of the group of SFB interact with the host. Future examination of FliC genes will likely enhance our knowledge of intestinal colonization by the gut microbiota.

Genetic diversity of catechol 1,2-dioxygenase in the fecal microbial metagenome

Catechol 1,2-dioxygenase is the key enzyme that catalyzes the cleavage of the aromatic ring of catechol. We explored the genetic diversity of catechol 1,2-dioxygenase in the fecal microbial metagenome by PCR with degenerate primers. A total of 35 gene fragments of C12O were retrieved from microbial DNA in the feces of pygmy loris. Based on phylogenetic analysis, most sequences were closely related to C12O sequences from Acinetobacter. A full-length C12O gene was directly cloned, heterologously expressed in Escherichia coli, and biochemically characterized. Purified catPL12 had optimum pH and temperature pH 8.0 and 25 °C and retained 31 and 50% of its maximum activity when assayed at 0 and 35 °C, respectively. The enzyme was stable at 25 and 37 °C, retaining 100% activity after pre-incubation for 1 h. The kinetic parameters of catPL12 were determined. The enzyme had apparent K_m of 67 µM, V_max of 7.3 U/mg, and k_cat of 4.2 s^-1 for catechol, and the cleavage activities for 3-methylcatechol, 4-methylcatechol, and 4-chlorocatechol were much less than for catechol, and no activity with hydroquinone or protocatechuate was detected. This study is the first to report the molecular and biochemical characterizations of a cold-adapted catechol 1,2-dioxygenase from a fecal microbial metagenome.

Towards Defining Heterotic Gene Pools in Pearl Millet [Pennisetum glaucum (L.) R. Br.]

Pearl millet is a climate resilient crop and one of the most widely grown millets worldwide. Heterotic hybrid development is one of the principal breeding objectives in pearl millet. In a maiden attempt to identify heterotic groups for grain yield, a total of 343 hybrid parental [maintainer (B-) and restorer (R-)] lines were genotyped with 88 polymorphic SSR markers. The SSRs generated a total of 532 alleles with a mean value of 6.05 alleles per locus, mean gene diversity of 0.55, and an average PIC of 0.50. Out of 532 alleles, 443 (83.27%) alleles were contributed by B-lines with a mean of 5.03 alleles per locus. R-lines contributed 476 alleles (89.47%) with a mean of 5.41, while 441 (82.89%) alleles were shared commonly between B- and R-lines. The gene diversity was higher among R-lines (0.55) compared to B-lines (0.49). The unweighted neighbor-joining tree based on simple matching dissimilarity matrix obtained from SSR data clearly differentiated B- lines into 10 sub-clusters (B1 through B10), and R- lines into 11 sub-clusters (R1 through R11). A total of 99 hybrids (generated by crossing representative 9 B- and 11 R- lines) along with checks were evaluated in the hybrid trial. The 20 parents were evaluated in the line trial. Both the trials were evaluated in three environments. Based on per se performance, high sca effects and standard heterosis, F₁s generated from crosses between representatives of groups B10R5, B3R5, B3R6, B4UD, B5R11, B2R4, and B9R9 had high specific combining ability for grain yield compared to rest of the crosses. These groups may represent putative heterotic gene pools in pearl millet.

Genetic polymorphism of 8 Y-STR loci in native population of Isfahan province in central part of Iran

BACKGROUND

Y-chromosome short tandem repeats (Y-STRs) are genetic markers with practical applications in human identification and population studies.

AIM

Here we present the allelic and haplotype frequencies of 8 Y-STR loci most commonly used in forensic medicine in 103 unrelated native males of Isfahan province, central part of Iran.

SUBJECTS AND METHODS

The cases were selected on the basis of strict criteria to assure pure native populations of Isfahan origin. DNA extracted from peripheral blood samples and PCR amplified for each marker. Y-specific STR loci DYS19, DYS385, DYS389I, DYS389II, DYS390, DYS391, DYS392 and DYS393 were included in this study.

RESULTS

The most common alleles for each locus were: DYS19, allele 12; DYS385, allele 12; DYS389I, allele 13; DYS389II, allele 29; DYS390, allele 24; DYS391, allele 10; DYS392, allele 11; and DYS393, allele 13. Gene diversity value was calculated from the allelic frequency for each locus. The average gene diversity was 0.6518. A total of 101 haplotypes were observed in eight Y-specific STR loci, the haplotype diversity was raised to 0.986.

CONCLUSION

The results revealed that a set of eight Y-specific STR loci were able to discriminate most of the male individuals in the population studied. A search through the Y Haplotype Reference Database demonstrated 21 matched haplotypes to 160,693 haplotypes, exclusively with Eurasian-European, Eurasian, and Eurasian-Indo Iranian populations.

Towards plant pangenomics

As an increasing number of genome sequences become available for a wide range of species, there is a growing understanding that the genome of a single individual is insufficient to represent the gene diversity within a whole species. Many studies examine the sequence diversity within genes, and this allelic variation is an important source of phenotypic variation which can be selected for by man or nature. However, the significant gene presence/absence variation that has been observed within species and the impact of this variation on traits is only now being studied in detail. The sum of the genes for a species is termed the pangenome, and the determination and characterization of the pangenome is a requirement to understand variation within a species. In this review, we explore the current progress in pangenomics as well as methods and approaches for the characterization of pangenomes for a wide range of plant species.

Biomarker identification from next-generation sequencing data for pathogen bacteria characterization and surveillance

AIM

The purpose was to develop an analytical pipeline for specific gene analysis and biomarker discovery from next generation sequencing (NGS) data.

MATERIALS & METHODS

As a test case, the fliC gene reference sequences of 24 Salmonella enterica strains of 13 serotypes and NGS reads of 32 serovar Newport, 48 Montevideo and 115 Enteritidis outbreak isolates were retrieved from the National Center for Biotechnology Information database.

RESULTS

Establishment of an analytical pipeline consisting of four steps: reference sequences retrieval and template sequence determination; NGS sequence reads retrieval; multiple sequence alignments and phylogenetic analysis; data mining and biomarker discovery.

CONCLUSION

The pipeline developed provides an effective bioinformatics tool for genetic diversity clarification and marker sequences discovery for pathogen characterization and surveillance.

Allelic diversity for neutral markers retains a higher adaptive potential for quantitative traits than expected heterozygosity

The adaptive potential of a population depends on the amount of additive genetic variance for quantitative traits of evolutionary importance. This variance is a direct function of the expected frequency of heterozygotes for the loci which affect the trait (QTL). It has been argued, but not demonstrated experimentally, that long-term response to selection is more dependent on QTL allelic diversity than on QTL heterozygosity. Conservation programmes, aimed at preserving this variation, usually rely on neutral markers rather than on quantitative traits for making decisions on management. Here, we address, both through simulation analyses and experimental studies with Drosophila melanogaster, the question of whether allelic diversity for neutral markers is a better indicator of a high adaptive potential than expected heterozygosity. In both experimental and simulation studies, we established synthetic populations for which either heterozygosity or allelic diversity was maximized using information from QTL (simulations) or unlinked neutral markers (simulations and experiment). The synthetic populations were selected for the quantitative trait to evaluate the evolutionary potential provided by the two optimization methods. Our results show that maximizing the number of alleles of a low number of markers implies higher responses to selection than maximizing their heterozygosity.

High Genetic Diversity of Microbial Cellulase and Hemicellulase Genes in the Hindgut of Holotrichia parallela Larvae

In this study, we used a culture-independent method based on library construction and sequencing to analyze the genetic diversity of the cellulase and hemicellulase genes of the bacterial community resident in the hindgut of Holotrichia parallela larvae. The results indicate that there is a large, diverse set of bacterial genes encoding lignocellulose hydrolysis enzymes in the hindgut of H. parallela. The total of 101 distinct gene fragments (similarity <95%) of glycosyl hydrolase families including GH2 (24 genes), GH8 (27 genes), GH10 (19 genes), GH11 (14 genes) and GH36 (17 genes) families was retrieved, and certain sequences of GH2 (10.61%), GH8 (3.33%), and GH11 (18.42%) families had <60% identities with known sequences in GenBank, indicating their novelty. Based on phylogenetic analysis, sequences from hemicellulase families were related to enzymes from Bacteroidetes and Firmicutes. Fragments from cellulase family were most associated with the phylum of Proteobacteria. Furthermore, a full-length endo-xylanase gene was obtained, and the enzyme exhibited activity over a broad range of pH levels. Our results indicate that there are large number of cellulolytic and xylanolytic bacteria in the hindgut of H. parallela larvae, and these symbiotic bacteria play an important role in the degradation of roots and other organic matter for the host insect.

Relationships Between Wright's F ST and F IS Statistics in a Context of Wahlund Effect

Waples (2015) has suggested a formula for the Wahlund effect in a case of unequal contribution of samples from genetically different populations that relates Wright's inbreeding coefficient, F IS , and normalized variance in allele frequencies between populations, F ST . I generalize this relationship to a case of multiple alleles and multiple populations not assuming Hardy-Weinberg ratios prior to mixing. This can help to evaluate the impact of a Wahlund effect on heterozygote deficiency relative to other factors such as null alleles, nonrandom mating, or selection. It is suggested that Wahlund effect cannot be an important factor of deviations from Hardy-Weinberg proportions in natural populations in the majority of instances, but it can have a substantial contribution to heterozygote deficiency in a population that has low genetic diversity compared to that among immigrants or in mixed samples that contain comparable fractions of individuals from genetically different populations.

Sequence diversity and molecular evolutionary rates between buffalo and cattle

Identification of genes of importance regarding production traits in buffalo is impaired by a paucity of genomic resources. Choice to fill this gap is to exploit data available for cow. The cross-species application of comparative genomics tools is potential gear to investigate the buffalo genome. However, this is dependent on nucleotide sequences similarity. In this study, gene diversity between buffalo and cattle was determined using 86 gene orthologues. There was approximately 3% difference in all genes in terms of nucleotide diversity and 0.267 ± 0.134 in amino acids, indicating the possibility for successfully using cross-species strategies for genomic studies. There were significantly higher non-synonymous substitutions both in cattle and buffalo; however, there was similar difference in terms of dN- dS (4.414 versus 4.745) in buffalo and cattle, respectively. Higher rate of non-synonymous substitutions at similar level in buffalo and cattle indicated a similar positive selection pressure. Results for relative rate test were assessed with the chi-squared test. There was no significance difference on unique mutations between cattle and buffalo lineages at synonymous sites. However, there was a significance difference on unique mutations for non-synonymous sites, indicating ongoing mutagenic process that generates substitutional mutation at approximately the same rate at silent sites. Moreover, despite of common ancestry, our results indicate a different divergent time among genes of cattle and buffalo. This is the first demonstration that variable rates of molecular evolution may be present within the family Bovidae.

Development of microsatellite markers for Vitex rotundifolia (Verbenaceae), an endangered coastal plant in Lake Biwa, Japan

PREMISE OF THE STUDY

Microsatellite loci were developed for Vitex rotundifolia, an endangered species isolated to Lake Biwa, to investigate its genetic diversity and population structure. •

METHODS AND RESULTS

Ten primer sets were identified in Japanese populations of V. rotundifolia. The number of alleles per locus ranged from one to six and gene diversity per locus ranged from 0.040 to 0.697 between two populations. In addition, all loci could be successfully amplified in V. trifolia. •

CONCLUSIONS

These markers will be useful for studies of genetic diversity and population structure of endangered species isolated to Lake Biwa of V. rotundifolia, to aid in the development of conservation strategies.

Exegeses on maximum genetic differentiation

A canon of population genetics concerns the properties of FST, a descriptor of spatial genetic structure. Interest for FST arose from Wright's early insights linking FST to dispersal parameters as well as to his concept of effective population size (e.g., Wright 1938, 1951). Although there is continued interest in this topic, FST also serves in other applications, such as detecting selected markers in natural populations (Beaumont and Nichols 1996) and more often in routine descriptive works. Remarkably, it is the latter use that seems to attract most discussion. Alternative descriptors have been proposed. Conversely, attempts have been made to draw biological inferences from FST properties that do not depend on biological processes. A reconsideration of its properties under biological scenarios underlines the weaknesses of such approaches.

Allelic diversity and its implications for the rate of adaptation

Genetic variation is usually estimated empirically from statistics based on population gene frequencies, but alternative statistics based on allelic diversity (number of allelic types) can provide complementary information. There is a lack of knowledge, however, on the evolutionary implications attached to allelic-diversity measures, particularly in structured populations. In this article we simulated multiple scenarios of single and structured populations in which a quantitative trait subject to stabilizing selection is adapted to different fitness optima. By forcing a global change in the optima we evaluated which diversity variables are more strongly correlated with both short- and long-term adaptation to the new optima. We found that quantitative genetic variance components for the trait and gene-frequency-diversity measures are generally more strongly correlated with short-term response to selection, whereas allelic-diversity measures are more correlated with long-term and total response to selection. Thus, allelic-diversity variables are better predictors of long-term adaptation than gene-frequency variables. This observation is also extended to unlinked neutral markers as a result of the information they convey on the demographic population history. Diffusion approximations for the allelic-diversity measures in a finite island model under the infinite-allele neutral mutation model are also provided.

Comparative Analysis of the Korean Population of Magnaporthe oryzae by Multilocus Microsatellite Typing

Rice blast fungus, Magnaporthe oryzae, inflicts serious damage to global rice production. Due to high variability of this fungal pathogen, resistance of newly-released rice cultivars is easily broken down. To understand the population structure of M. oryzae, we analyzed the genetic diversity of the Korean population using multilocus microsatellite typing. Eleven microsatellite markers were applied to the population of 190 rice isolates which had been collected in Korea for two decades since the 1980's. Average values of gene diversity and allele frequency were 0.412 and 6.5, respectively. Comparative analysis of the digitized allele information revealed that the Korean population exhibited a similar level of allele diversity to the integrated diversity of the world populations, suggesting a particularly high diversity of the Korean population. Therefore, these microsatellite markers and the comprehensive collection of field isolates will be useful genetic resources to identify the genetic diversity of M. oryzae population.

Allelic diversity and its implications for the rate of adaptation

Genetic variation is usually estimated empirically from statistics based on population gene frequencies, but alternative statistics based on allelic diversity (number of allelic types) can provide complementary information. There is a lack of knowledge, however, on the evolutionary implications attached to allelic-diversity measures, particularly in structured populations. In this article we simulated multiple scenarios of single and structured populations in which a quantitative trait subject to stabilizing selection is adapted to different fitness optima. By forcing a global change in the optima we evaluated which diversity variables are more strongly correlated with both short- and long-term adaptation to the new optima. We found that quantitative genetic variance components for the trait and gene-frequency-diversity measures are generally more strongly correlated with short-term response to selection, whereas allelic-diversity measures are more correlated with long-term and total response to selection. Thus, allelic-diversity variables are better predictors of long-term adaptation than gene-frequency variables. This observation is also extended to unlinked neutral markers as a result of the information they convey on the demographic population history. Diffusion approximations for the allelic-diversity measures in a finite island model under the infinite-allele neutral mutation model are also provided.

Efficiency of conservation management methods for subdivided populations under local adaptation

Computer simulations were used to investigate the efficiency of management methods for the conservation of a structured population when local adaptation exists. A subdivided population, with subpopulations adapted to different optima for a quantitative trait under stabilizing selection, was managed in order to maintain the highest genetic diversity in a 10-generation period. Two procedures were compared. For the first, minimum coancestry contributions were carried out independently for each subpopulation, and random migration of individuals was accomplished thereafter. For the second, minimum coancestry contributions from individuals were globally implemented, including an optimal migration design. This optimal method can be adjusted to control local inbreeding to different extents. Adaptation to local optima implies a reduction in the efficiency of the management methods because of the effective failure in the established migrations. For strong selection, the optimal design can be very inefficient, even more than the random migration scheme because the intended migrants have usually low fitness in the recipient subpopulations. However, for more realistic moderate or weak selection, the optimal method is more efficient than random migration, especially if inbreeding depression on fitness is also taken into account. It is concluded that the optimal management method can be recommended in conservation programs with local adaptation of subpopulations, but this issue should be accounted for when designing the management strategies.

Evolutionary constraints in the β-globin cluster: the signature of purifying selection at the δ-globin (HBD) locus and its role in developmental gene regulation

Human hemoglobins, the oxygen carriers in the blood, are composed by two α-like and two β-like globin monomers. The β-globin gene cluster located at 11p15.5 comprises one pseudogene and five genes whose expression undergoes two critical switches: the embryonic-to-fetal and fetal-to-adult transition. HBD encodes the δ-globin chain of the minor adult hemoglobin (HbA2), which is assumed to be physiologically irrelevant. Paradoxically, reduced diversity levels have been reported for this gene. In this study, we sought a detailed portrait of the genetic variation within the β-globin cluster in a large human population panel from different geographic backgrounds. We resequenced the coding and noncoding regions of the two adult β-globin genes (HBD and HBB) in European and African populations, and analyzed the data from the β-globin cluster (HBE, HBG2, HBG1, HBBP1, HBD, and HBB) in 1,092 individuals representing 14 populations sequenced as part of the 1000 Genomes Project. Additionally, we assessed the diversity levels in nonhuman primates using chimpanzee sequence data provided by the PanMap Project. Comprehensive analyses, based on classic neutrality tests, empirical and haplotype-based studies, revealed that HBD and its neighbor pseudogene HBBP1 have mainly evolved under purifying selection, suggesting that their roles are essential and nonredundant. Moreover, in the light of recent studies on the chromatin conformation of the β-globin cluster, we present evidence sustaining that the strong functional constraints underlying the decreased contemporary diversity at these two regions were not driven by protein function but instead are likely due to a regulatory role in ontogenic switches of gene expression.

Genetic analysis of scattered populations of the Indian eri silkworm, Samia cynthia ricini Donovan: Differentiation of subpopulations

Deforestation and exploitation has led to the fragmentation of habitats and scattering of populations of the economically important eri silkworm, Samia cynthia ricini, in north-east India. Genetic analysis of 15 eri populations, using ISSR markers, showed 98% inter-population, and 23% to 58% intra-population polymorphism. Nei's genetic distance between populations increased significantly with altitude (R(2) = 0.71) and geographic distance (R(2) = 0.78). On the dendrogram, the lower and upper Assam populations were clustered separately, with intermediate grouping of those from Barpathar and Chuchuyimlang, consistent with geographical distribution. The Nei's gene diversity index was 0.350 in total populations and 0.121 in subpopulations. The genetic differentiation estimate (Gst) was 0.276 among scattered populations. Neutrality tests showed deviation of 118 loci from Hardy-Weinberg equilibrium. The number of loci that deviated from neutrality increased with altitude (R(2) = 0.63). Test of linkage disequilibrium showed greater contribution of variance among eri subpopulations to total variance. D('2)IS exceeded D('2)ST, showed significant contribution of random genetic drift to the increase in variance of disequilibrium in subpopulations. In the Lakhimpur population, the peripheral part was separated from the core by a genetic distance of 0.260. Patchy habitats promoted low genetic variability, high linkage disequilibrium and colonization by new subpopulations. Increased gene flow and habitat-area expansion are required to maintain higher genetic variability and conservation of the original S. c. ricini gene pool.

Analysis of expressed sequence tags from the ectomycorrhizal basidiomycetes Laccaria bicolor and Pisolithus microcarpus

•  In an effort to discover genes that are expressed in the ectomycorrhizal basidiomycetes Laccaria bicolor and Pisolithus microcarpus, and in P. microcarpus/Eucalyptus globulus ectomycorrhizas, we have sequenced 1519 and 1681 expressed sequence tags (ESTs) from L. bicolor and P. microcarpus cDNA libraries. •  Contig analysis resulted in 905 and 806 tentative consensus sequences (unique transcripts) in L. bicolor and P. microcarpus, respectively. For 36% of the ESTs, significant similarities to sequences in databases were detected. The most abundant transcripts showed no similarity to previously identified genes. Sequence redundancy analysis between different developmental stages indicated that several genes were differentially expressed in free-living mycelium and symbiotic tissues of P. microcarpus. •  Based on sequence similarity, 11% of L. bicolor unique transcripts were also detected in P. microcarpus. Similarly, L. bicolor and P. microcarpus shared only a low proportion of common transcripts with other basidiomycetous fungi, such as Pleurotus ostreatus and Agaricus bisporus. Such a low proportion of shared transcripts between basidiomycetes suggests, on the one hand, that the variability of expressed transcripts in different fungi and fungal tissues is considerably high. On the other hand, it might reflect the low number of GenBank entries of basidiomycetous origin and stresses the necessity of an additional sequencing effort. •  The present ESTs provide a valuable resource for future research on the development and functioning of ectomycorrhizas.

Molecular diversity of the calcium channel alpha2delta subunit

Sequence database searches with the alpha2delta subunit as probe led to the identification of two new genes encoding proteins with the essential properties of this calcium channel subunit. Primary structure comparisons revealed that the novel alpha2delta-2 and alpha2delta-3 subunits share 55.6 and 30.3% identity with the alpha2delta-1 subunit, respectively. The number of putative glycosylation sites and cysteine residues, hydropathicity profiles, and electrophysiological character of the alpha2delta-3 subunit indicates that these proteins are functional calcium channel subunits. Coexpression of alpha2delta-3 with alpha1C and cardiac beta2a or alpha1E and beta3 subunits shifted the voltage dependence of channel activation and inactivation in a hyperpolarizing direction and accelerated the kinetics of current inactivation. The kinetics of current activation were altered only when alpha2delta-1 or alpha2delta-3 was expressed with alpha1C. The effects of alpha2delta-3 on alpha1C but not alpha1E are indistinguishable from the effects of alpha2delta-1. Using Northern blot analysis, it was shown that alpha2delta-3 is expressed exclusively in brain, whereas alpha2delta-2 is found in several tissues. In situ hybridization of mouse brain sections showed mRNA expression of alpha2delta-1 and alpha2delta-3 in the hippocampus, cerebellum, and cortex, with alpha2delta-1 strongly detected in the olfactory bulb and alpha2delta-3 in the caudate putamen.