KIR2DL4 genetic diversity in a Brazilian population sample: implications for transcription regulation and protein diversity in samples with different ancestry backgrounds.

KIR2DL4 is an important immune modulator expressed in natural killer cells; HLA-G is its main ligand. We have characterized the KIR2DL4 genetic diversity by considering the promoter, all exons, and all introns in a highly admixed Brazilian population sample and by using massively parallel sequencing. We introduce a molecular method to amplify and to sequence the complete KIR2DL4 gene. To avoid the mapping bias and genotype errors commonly observed in gene families, we have developed and validated a bioinformatic pipeline designed to minimize these errors and applied it to survey the variability of 220 individuals from the State of São Paulo, southeastern Brazil. We have also compared the KIR2DL4 genetic diversity in the Brazilian cohort with the diversity previously reported by the 1000Genomes consortium. KIR2DL4 presents high linkage disequilibrium throughout the gene, with coding sequences associated with specific promoters. There are few but divergent promoter haplotypes. We have also detected many new KIR2DL4 sequences, all bearing nucleotide exchanges in introns and encoding previously described proteins. Exons 3 and 4, which encode the external domains, are the most variable. The ancestry background influences the KIR2DL4 allele frequencies and must be considered for association studies regarding KIR2DL4.

Hla-mapper: An application to optimize the mapping of HLA sequences produced by massively parallel sequencing procedures.

A challenging task when more than one HLA gene is evaluated together by second-generation sequencing is to achieve a reliable read mapping. The polymorphic and repetitive nature of HLA genes might bias the read mapping process, usually underestimating variability at very polymorphic segments, or overestimating variability at some segments. To overcome this issue we developed hla-mapper, which takes into account HLA sequences derived from the IPD-IMGT/HLA database and unpublished HLA sequences to apply a scoring system. This comprehends the evaluation of each read pair, addressing them to the most likely HLA gene they were derived from. Hla-mapper provides a reliable map of HLA sequences, allowing accurate downstream analysis such as variant calling, haplotype inference, and allele typing. Moreover, hla-mapper supports whole genome, exome, and targeted sequencing data. To assess the software performance in comparison with traditional mapping algorithms, we used three different simulated datasets to compare the results obtained with hla-mapper, BWA MEM, and Bowtie2. Overall, hla-mapper presented a superior performance, mainly for the classical HLA class I genes, minimizing wrong mapping and cross-mapping that are typically observed when using BWA MEM or Bowtie2 with a single reference genome.

Diversity and Structure of Fungal Communities in Neotropical Rainforest Soils: The Effect of Host Recurrence.

The patterns of the distribution of fungal species and their potential interactions with trees remain understudied in Neotropical rainforests, which harbor more than 16,000 tree species, mostly dominated by endomycorrhizal trees. Our hypothesis was that tree species shape the non-mycorrhizal fungal assemblages in soil and litter and that the diversity of fungal communities in these two compartments is partly dependent on the coverage of trees in the Neotropical rainforest. In French Guiana, a long-term plantation and a natural forest were selected to test this hypothesis. Fungal ITS1 regions were sequenced from soil and litter samples from within the vicinity of tree species. A broad range of fungal taxa was found, with 42 orders and 14 classes. Significant spatial heterogeneity in the fungal communities was found without strong variation in the species richness and evenness among the tree plots. However, tree species shaped the fungal assemblages in the soil and litter, explaining up to 18 % of the variation among the communities in the natural forest. These results demonstrate that vegetation cover has an important effect on the structure of fungal assemblages inhabiting the soil and litter in Amazonian forests, illustrating the relative impact of deterministic processes on fungal community structures in these highly diverse ecosystems.

Mitogenome assembly from genomic multiplex libraries: comparison of strategies and novel mitogenomes for five species of frogs.

Next-generation sequencing continues to revolutionize biodiversity studies by generating unprecedented amounts of DNA sequence data for comparative genomic analysis. However, these data are produced as millions or billions of short reads of variable quality that cannot be directly applied in comparative analyses, creating a demand for methods to facilitate assembly. We optimized an in silico strategy to efficiently reconstruct high-quality mitochondrial genomes directly from genomic reads. We tested this strategy using sequences from five species of frogs: Hylodes meridionalis (Hylodidae), Hyloxalus yasuni (Dendrobatidae), Pristimantis fenestratus (Craugastoridae), and Melanophryniscus simplex and Rhinella sp. (Bufonidae). These are the first mitogenomes published for these species, the genera Hylodes, Hyloxalus, Pristimantis, Melanophryniscus and Rhinella, and the families Craugastoridae and Hylodidae. Sequences were generated using only half of one lane of a standard Illumina HiqSeq 2000 flow cell, resulting in fewer than eight million reads. We analysed the reads of Hylodes meridionalis using three different assembly strategies: (1) reference-based (using bowtie2); (2) de novo (using abyss, soapdenovo2 and velvet); and (3) baiting and iterative mapping (using mira and mitobim). Mitogenomes were assembled exclusively with strategy 3, which we employed to assemble the remaining mitogenomes. Annotations were performed with mitos and confirmed by comparison with published amphibian mitochondria. In most cases, we recovered all 13 coding genes, 22 tRNAs, and two ribosomal subunit genes, with minor gene rearrangements. Our results show that few raw reads can be sufficient to generate high-quality scaffolds, making any Illumina machine run using genomic multiplex libraries a potential source of data for organelle assemblies as by-catch.